Advanced
An Efficient Identification Algorithm in a Low SNR Channel
An Efficient Identification Algorithm in a Low SNR Channel
Journal of the Korea Institute of Information and Communication Engineering. 2014. Apr, 18(4): 790-796
Copyright © 2014, The Korea Institute of Information and Commucation Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License(http://creativecommons.org/li-censes/by-nc/3.0/) which permits unrestricted non-commercial use, dist ribution, and reproduction in any medium, provided the origin al work is properly cited.
  • Received : February 11, 2014
  • Accepted : March 20, 2014
  • Published : April 30, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
지원 황
Department of Information Technology, Chonb uk National University, Jeonju 561-756, Korea
주필 조
Department of Radiocommunication Engineering, K unsan National University, Kunsan 573-701, Korea
stefano@kusan.ac.kr

Abstract
Identification of communication channels is a problem of important current theoretical and practical concerns. Recently proposed solutions for this problem exploit the diversity induced by antenna array or time oversampling. The method resorts to an adaptive filter with a linear constraint. In this paper, an approach is proposed that is based on decomposition. Indeed, the eigenvector corresponding to the minimum eigenvalue of the covariance matrix of the received signals contains the channel impulse response. And we present an adaptive algorithm to solve this problem. Proposed technique shows the better performance than one of existing algorithms.
Keywords
I. INTRODUCTION
In HOS-based methods, because the performance index as the optimization criterion is nonlinear with respect to estimation parameters and these methods require a large amount of data samples. These methods have the disadvantage that their computational complexity may be large. See, for example, [1] and references therein. Since the seminal work by Tong et al. the problem of estimating the channel response of multiple FIR channel driven by an unknown input symbol has interested many researchers in signal processing and communication fields. This is achieved by exploiting assumed cyclostationary properties, induced by oversampling or antenna array at the receiver part [1 , 2] . The basic blind channel identification problem involves a channel model where only the observation signal is available for processing in the identification channel. Earlier blind channel identification approaches mostly depend on higher order statistics (HOS), because the second order statistics (SOS) does not contain phase information for stationary signal [3 - 4] . Most communication channels are time-varying in practice. Therefore, the algorithms should be able to track the change of the channel impulse response. Moreover, in a fast fading channel, the multipath channels in wireless communications vary rapidly, and we only have a few data samples corresponding to the same channel characteristics. Blind channel identification technique has been developed in adaptive algorithm based on vector-correlation method [8 , 9 , 11] . But most algorithms neglected the effect of channel noise.
Most notations are standard: vectors and matrices are boldface small and capital letters, respectively; the matrix transpose, the complex conjugate, the Hermitian, and convolution are denoted by (·) T , (·) * , (·) H and ⊗, respectively; I P is the P × P identity matrix; E (·) is the statistical expectation.
This paper is organized as follows. In section II, we review the basic assumption and identification issues. And the existing adaptive algorithms of the block LS methods are described also. A novel blind channel identification technique based on eigenvlaue decomposition and adaptive implementation are proposed in section III. Simulation results with real measured channel are performed in section IV. Section V concludes our results.
Ⅱ. Basic Assumption and Issues
In this paper, consider a special case, when the channel output is two times oversampled or there are two antennas at the receiver, this is equivalent to two channel representation ( M =2). From the Fig. 1 , in the absence of noise, it is apparent that the output of each subchannel is
PPT Slide
Lager Image
Then
PPT Slide
Lager Image
Let 𝒳( t ) be the signal at the output of a noisy channel
PPT Slide
Lager Image
PPT Slide
Lager Image
두 서브채널간 상관관계 Fig. 1 The cross relation between two subchannels
PPT Slide
Lager Image
M 서브채널을 갖는 등가 SIMO 모델 Fig. 2 Equivalent SIMO model with M subchannels
where s ( k ) denotes the transmitted symbol at time kT , h ( t ) denotes the continuous-time channel impulse response, and v ( t ) is additive noise. As shown in [3] , the single channel system can be considered as the multichannel system by the sampling the received signal at a rate faster than the input symbol rate. The source signal s ( n ) then passes through M equivalent symbol rate linear filters. And as shown in Fig. 2 , x i ( n ) denotes the output from the ith channel with the noisy FIR channel impulse response { h i ( n )}, which is driven by the same input s ( n ). Clearly, for linearly modulated communication signals, x i ( n ) , a i ( n ) , s ( n ) , v i ( n ) , and h i ( n ) are related as follows
PPT Slide
Lager Image
where L is the maximum order of the M channels.
The blind identification problem can be stated as follows: Given the observation of channel output { x i ( n ) , i =1,⋯, M ; n = L ,⋯, N }, determine the channels and further recover the input signals { s ( n ) }. As in classical system identification problems, certain conditions about the channel and the source must be satisfied to ensured identifiability. We assume the following throughout in this paper about the channel and source conditions.
  • A1) Subchannels do not share common zeros, or in other words, they are coprime.
The assumption that L is known may be practical. To address this problem, there are three approaches [5] . First, channel order detection and parameter estimation can be performed separately. Second, some statistical subspace methods require only upper bound of L . Third, channel order detection and parameter estimation can be performed jointly.
Ⅲ. Proposed Scheme
As described in [5] , to avoid the trivial solution to minimization problem a proper condition must be selected. In this section, a new approach is proposed that is based on eigenvalue decomposition. Indeed, the eigenvector corresponding to the minimum eigenvalue of the covariance matrix of the received signals contains the channel impulse response. This approach is based on the unit norm constraint that is apart from the linear constraint introduced in the previous section [6] .
- A. Concept of the Proposed Scheme
Number equations consecutively with equation numbers in parentheses flush with the right margin, as in (1). We assume that the channel is linear and time invariant within small time interval; therefore, we have the following relation as described in (4)
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and the channel impulse response vector of length L are defined as
PPT Slide
Lager Image
The covariance matrix of the two received signals is given by
PPT Slide
Lager Image
Consider the 2L´1 vector as follows:
PPT Slide
Lager Image
From (5) and (8), it can be seen that R x h = 0 , which means that the vector h is the eigenvector of the covariance matrix R x corresponding to the eigenvalue 0.
Moreover, if the two channel impulse response h 1 and h 2 have no common zeros and the autocorrelation matrix of the source signal s ( n ) is full rank, which is assumed in the rest of this paper, the covariance matrix R x has one and only one eigenvalue equal to zero. Consider the noisy channel case as described in (2) and let M =2. It follows from (1) that
PPT Slide
Lager Image
where x ( n )=[ x 1 T ( n )) x 2 T ( n )] T and v ( n )=[ v 1 T ( n )) v 2 T ( n )] T .
If the correlation matrix of the vector x ( n ) is denoted by R x , a direct of conclusion of (10) will be
PPT Slide
Lager Image
We note from (11) that h is the eigenvector of the correlation matrix R x and
PPT Slide
Lager Image
is the corresponding eigenvector of R x . The knowledge of
PPT Slide
Lager Image
can be obtained as a by product if wanted.
PPT Slide
Lager Image
- B. Adaptive Algorithm
In practice, it is simple to estimate iteratively the eigenvector corresponding to the minimum eigenvalue of R x , by using an algorithm similar to the Frost algorithm that is a simple constrained LMS algorithm [7] .
Minimizing the quantity h H R x h with respect to h and subject to ∥ h 2 = h H h = 1 will give us the optimum weight h opt .
Let us define the error signal
PPT Slide
Lager Image
where x ( n )=[ x 1 T ( n ) x 2 T ( n )] T . Note that minimizing the mean square value of e ( n ) is equivalent to solving the above eigenvalue problem. Taking the gradient of e ( n ) with respect to h ( n ) gives
PPT Slide
Lager Image
and we obtain the gradient-descent constrained LMS algorithm:
PPT Slide
Lager Image
where μ , the adaptation step-size, is a positive constant. Substituting (13) and (14) into (15) gives
PPT Slide
Lager Image
and taking statistical expectation after convergence, we get
PPT Slide
Lager Image
which is what is desired: the eigenvector h (∞) corresponding to the smallest eigenvalue E [| e ( n )| 2 ] of the covariance matrix R x .
In practice, it is advantageous to use the following adaptation scheme
PPT Slide
Lager Image
The algorithm (18) presented above is very general to find the eigenvector corresponding to the smallest eigenvalue of any matrix R x . If the smallest eigenvalue is equal to zero, which is the case here, the algorithm can be simplified as follows:
PPT Slide
Lager Image
and
PPT Slide
Lager Image
IV. Simulation Results
Computer simulations were conducted to evaluate the performance of the proposed algorithm in comparison with existing algorithms. In all the simulations, two channel SIMO model is assumed. This means two times oversampling or two sensors at the receiver in real situation. The input signal is 4-QAM. For simplicity of comparison, we assumed that the channel order L is known. The performance index is achieved by examination the root mean square error (RMSE) that is defined as [4] .
PPT Slide
Lager Image
where N t is number of Monte Carlo trials, and
PPT Slide
Lager Image
is the estimate of the channels from the i th trial.
We used real-measured microwave channel. The shortened length-16 version of an empirically measured T /2-spaced digital microwave radio channel ( M =2) with 230 taps, which we truncated to obtain a channel with L =7. The Microwave channel chan1.mat is founded at http://spib.rice.edu/spib/ microwave.html.
채널 계수Table. 1Channel Coefficients
PPT Slide
Lager Image
채널 계수 Table. 1 Channel Coefficients
The shortened version is derived by linear decimation of the FFT of the full-length T /2-spaced impulse response and taking the IFFT of the decimated version (see [10] for more details on this channel). The channel coefficients for both sets of channels are listed in Table 1 . A total number of 50 independent trials were performed. All algorithms were initiated at h (0)=[1, 0, ..., 0, 1, 0, ..., 0] T with the step size μ =0.01.
Fig. 3 shows the RMSE of the channel estimates from existing algorithms and the proposed algorithm. From these figures, we can see that the proposed algorithm always performs better than others. By inspection, we can observe that RMSE values of the proposed method are decreased more or less 6-10 dB, and 1-2 dB under 20 dB, and 10 dB, respectively. Clearly, we can observe the significant improvement of the proposed algorithm over existing algorithms.
PPT Slide
Lager Image
제안 알고리즘과 기존 알고리즘의 RMSE 비교 (a) SNR=20dB and (b) SNR=10dB Fig. 3 RMSE comparison of the proposed and existing algorithms (a) SNR=20dB and (b) SNR=10dB
V . Conclusion
In this paper, an approach to channel identification has been presented. The method is based on eigenvalue decomposition. The eigenvector corresponding to the minimum eigenvalue of the covariance matrix of the received signals contains the channel impulse response. And we use a simple constrained LMS algorithm to estimate iteratively the eigenvector corresponding to the minimum eigenvalue. In comparison with algorithms, the proposed one seems to be more efficient in a low SNR channel and much more accurate.
BIO
황지원(Jeewon Hwang)
1995년: 전북대학교 전자공학과 공학박사
2003년: 뉴질랜드 CPIT 방문교수
1992년 ~ 현재 전북대학교 IT정보공학과 교수
※관심분야 : 컴퓨터구조, Cognitive Radio, LTE 요소기술
조주필(Juphil Cho)
2001년 : 전북대학교 전자공학과 공학박사
2000년 ~ 2005년 : ETRI 이동통신 연구단 선임연구원
2006년~2007년 : ETRI 초빙연구원
2011년 : 미국 USF, 교환교수
2005년~ 현재 : 군산대학교 전파공학과 부교수
※관심분야 : Cognitive-Radio, 주파수 융합기술, LTE
References
Haykin S. 1996 Adaptive Filter Theory Prentice-Hall Englewood Cliffs, NJ
Baccala L. A. , Roy S. 1994 “A new blind time-domain channel identification method based on cyclostationarity,” IEEE Signal Processing Letters http://dx.doi.org/10.1109/97.295342 1 (6) 89 - 92    DOI : 10.1109/97.295342
Tong L. , Xu G. , Kailath T. 1994 “Blind identification and equalization based on second order statistics: A time domain approach,” IEEE Trans. Inform. Theory http://dx.doi.org/10.1109/18.312157 40 (2) 340 - 349    DOI : 10.1109/18.312157
Xu G. , Liu H. , Tong L. , Kailath T. 1995 “A least-squares approach to blind channel identification,” IEEE Trans. Signal Processing http://dx.doi.org/10.1109/78.476442 43 (12) 2982 - 2983    DOI : 10.1109/78.476442
Tong L. , Perreau S. 1998 “Multichannel blind identification: from subspace to maximum likelihood methods,” Proc. IEEE 86 (10) 1951 - 1968
Gesbert D. , Duhamel P. 2000 “Unbiased blind adaptive channel identification and equalization,” IEEE Trans. Signal Processing 48 (1)
Frost III O. L. 1972 “A algorithm for linearly constrained adaptive arrays,” Proc. IEEE 60 (8) 926 - 935
Heath Jr. R. W. , Halford S. D. , Giannakis G. B. 1997 “Adaptive blind channel identification of FIR channels for viterbi decoding,” Proc. 31th Asilomar Conf. Signals, Syst., and Comput
Higa Y. , Ochi H. , Kinjo S. , Yamaguchi H. 1999 “A gradient type algorithm for blind system identification and equalizer based on second order statistics,” IEICE Trans. Fundamentals E32-A (8) 1544 - 1551
Endres T. J. , Halford S. D. , Johnson C. R. , Giannakis G. B. 1998 “Simulated comparisions of blind equalization algorithms for cold startup applications,” Int J. Adaptive Contr. Signal Process http://dx.doi.org/10.1002/(SICI)1099-1115(199805)12:3<283::AID-ACS479>3.0.CO;2-E 12 (3) 283 - 301    DOI : 10.1002/(SICI)1099-1115(199805)12:3<283::AID-ACS479>3.0.CO;2-E
Goekel D. L. , Hero A. O. , Stark W. E. 1998 “Data-recursive algorithms for blind channel identification in oversampled communication systems,” IEEE Trans. Signal Processing 2217 - 2220