Wireless communication systems have evolved substantially over the last two decades. The explosive growth of the wireless communication market is expected to continue in the future, as the demand for all types of wireless services is increasing. Due to their ability to provide high data rates for multimedia applications, Orthogonal Frequency Division Multiplexing (OFDM) is gaining a strong interest for wide-area, local area, and personal area networks. In OFDM, carrier frequencies are chosen in such a way that there is no influence of other carriers in the detection of the information in the carrier of interest when the orthogonality of the carrier is maintained. Maintaining this orthogonality requires some special care for the frequency and symbol timing synchronizations. Cyclic extension of the symbols is usually used to relax the requirements for symbol synchronization.
In wireless OFDM systems, channel estimation is an integral part of the coherent receiver design as the performance of the receiver is greatly affected by the quality of the channel estimation. Extensive studies on the topic can be found in the literature. In OFDM systems, channel estimation is frequently employed in the frequency domain after taking the Discrete Fourier Transform (DFT) of the time synchronized digital samples. There are numerous approaches for estimating the Channel Frequency Response (CFR) over the OFDM subcarriers. The direct Least-Squares (LS) estimation assumes the channel over each subcarrier to be independent. However, in practice, the CFR is often oversampled over these subcarriers, and the estimated coefficients are correlated. On the other hand, the noise in these subcarriers can be independent. By exploiting the correlation of CFR over OFDM subcarriers, the noise can be reduced significantly, and hense the channel estimation accuracy can be improved. Assuming that the channel frequency correlation and noise variance are known, Minimum Mean-Square Error (MMSE) filtering of the LS estimates has been shown to provide optimal performance under Additive White Gaussian (AWGN).
Many of the channel estimation approaches proposed in the literature assume perfect symbol timing. However, in practice, the symbol timing used in OFDM systems is not perfect. As such, the symbol timing is often intentionally shifted towards the Cyclic Prefix (CP) so that any possible error in symbol timing that might create the loss of orthogonality can be avoided. Even though this intentional bias in synchronization avoids the loss of orthogonality of the carriers and intercarrier-interference, it results in the effective CFR to be less correlated due to the additional carrier-dependent phase shift. Synchronization errors in the receiver cause a linear phase rotation at the output of the DFT block. The correlation between the channel coefficients at different subcarriers is weakened due to this phase rotation. As a result, the performance of MMSE channel estimation degrades significantly since the noise averaging effect will be reduced.
Accordingly, what is needed in the art is an improved system and method for performing Minimum Mean-Square Error (MMSE) channel estimation in an Orthogonal Frequency Division Multiplexing (OFDM) channel under synchronization errors.
In accordance with a particular embodiment of the present invention, one approach to solving the problem of synchronization errors is by estimating the timing offset and removing the linear phase rotation caused by it. Fortunately, the linear phase of the estimated channel in the frequency domain is mainly due to the timing offset with less effect from the noise and the actual multi-path channel. So by estimating and removing the linear phase of the estimated channel the system performance can be improved, especially for medium and high signal-to-noise ration (SNR) values.
In accordance with an embodiment of the present invention, a method for improving the Minimum Mean-Square Error (MMSE) channel estimation in an Orthogonal Frequency Division Multiplexing (OFDM) channel under synchronization errors includes the steps of receiving data over a multipath OFDM channel, estimating a timing offset for the channel, wherein the timing offset results from synchronization errors in the channel, estimating a linear phase rotation resulting from the synchronization errors in the channel, wherein the linear phase rotation is dependent upon the estimated timing offset, estimating a channel frequency response for the channel using a direct least-squares estimation, removing the estimated linear phase rotation from the estimated channel frequency response estimate, filtering the channel frequency response estimate for the channel using the MMSE channel estimation and then adding the estimated linear phase rotation back into the filtered channel frequency response estimate.
In a specific embodiment, the step of estimating a timing offset for the channel further comprises approximating the channel linear phase to the nearest value in C=[C0C1 . . . CN-CP-l], to obtain C{tilde over (d)} where {tilde over (d)} is the timing offset estimate.
In an additional embodiment, another approach to remedy the effect of timing offsets is by increasing the channel correlation. Since the timing offset weakens the correlation between the channel's frequency coefficients, by removing the linear phase corresponding to different timing offset values and choosing the one that results in the channel with maximum correlation we can get a very good estimate of the timing offset. Moreover, even with a wrong timing offset assumed, we would still be left with a highly correlated channel thus improving the channel estimation performance. Another advantage to this approach is its independence of the noise level, since the noise is uncorrelated between different frequency subcarriers as long as their orthogonality is maintained.
As such, in accordance with an additional embodiment of the present invention, the step of estimating a timing offset resulting from the synchronization errors of the channel further includes estimating a timing offset by identifying the timing offset that results in a maximum correlation of the channel in the frequency domain. In a particular embodiment, the timing offset estimating further includes calculating the timing offset {tilde over (d)} where:
where R{tilde over (H)}(Δ)=E{{tilde over (H)}k†{tilde over (H)}k−Δ}, is the frequency-domain channel correlation function with a frequency separation Δ, {tilde over (d)} is the timing offset estimate, {tilde over (H)}=Θd−1
In accordance with the present invention, filtering the channel frequency response for the channel using the MMSE channel estimation further includes estimating the channel frequency response using:
{tilde over (H)}MMSE=Θ{tilde over (d)}F†Θ{tilde over (d)}−1ĤLS
where,
F†={tilde over (R)}HH ({tilde over (R)}HH+σ2I)−1
and where, {tilde over (H)}MMSE is the MMSE channel estimation, ĤLS is the least-squares estimate of H, F is the MMSE filter, I is N×N identity matrix, RHH=E{{tilde over (H)}{tilde over (H)}†}, {tilde over (H)}=Θ{tilde over (d)}−1
According, the embodiments of the present invention are effective in improving the system performance, thus allowing the system to operate at higher noise levels with less errors and better channel estimation.
As such, the present invention provides a novel system and method for MMSE channel estimation that takes synchronization errors (intentional or not) into account. The proposed channel estimation in accordance with the present invention improves the noise averaging capability and takes advantage of channel correlation fully by removing the effect of synchronization errors during the estimation process.
For a fuller understanding of the invention, reference should be made to the following detailed description, taken in connection with the accompanying drawings, in which:
According to a particular embodiment of the invention, we consider an N subcarrier OFDM system with Xk as the frequency domain transmitted symbol at the kth subcarrier. The time domain complex baseband data sequence x=[x0x1 . . . xN−1]T is obtained at the output of the Inverse Discrete Fourier transform (IDFT) block with X=[X0X1 . . . XN−1]T as the input. Before transmission, a cyclic prefix of length NCP is added at the beginning of the data sequence to eliminate the Inter-Symbol Interference (ISI) and preserve the orthogonality of the frequency subcarriers. The data is then transmitted over a multipath channel. The time domain channel impulse response is described as a pulse train
Where l is the total number of paths, αi(t) is complex-valued Gaussian random variable with zero mean, and τi is the delay of the ith path. The CP length NCP is chosen such that 0<τl−1<NCPTS where TS is the sampling interval so that ISI will is eliminated. The received signal is sampled and fed into an N-point DFT block. The output of the DFT block is
Y=XH+W (2)
Where Y is the received vector, X is a diagonal matrix containing the transmitted signal, H is the channel vector, and W is a complex zero-mean Gaussian noise with variance σ2.
When there are no synchronization errors (zero timing offset) the MMSE channel estimation is as follows:
ĤMMSE=F†ĤLS (3)
and
F†=RRHH(RHH+σ2I)−1 (4)
Where F is the MMSE filter, I is N×N identity matrix, RHH=E{HH†} is the channel autocorrelation matrix (the superscript (.)†denotes the Hermitian transpose), and ĤLS is the LS estimate of H,
ĤLS=X−1Y=H+X−1W (5)
However, a timing synchronization error will cause an offset of d samples at the input of the receiver DFT block. For no ISI, the timing offset d should be 0≦d≦NCP−l−1. Since d can have negative values, an intentional timing offset is usually added at the receiver to guarantee that d24 0. When the timing offset is within this range, the equivalent CFR at subcarrier k is
The received frequency domain signal including the synchronization error can then be presented as
Y=X
where Θd is a diagonal matrix containing the phase rotation exp(j2πdk/N), k=0,1, . . . N−1. Due to this timing offset, the effective channel will appear to be changing at a faster rate than the actual channel and the correlation between the channel coefficients at different subcarriers will be weaker. As a result, the performance of the MMSE estimation will degrade.
As previously described in the background of the invention, synchronization errors in the OFDM receiver will result in a linear phase rotation that is dependent on the timing offset d. The present invention provides a system and method to estimate this phase rotation and reverse its effect.
The channel phase at subcarrier k
∠(
where Cd=j2πd/N, and d=0,1, . . . NCp−l−1. So, the linear phase component of the effective channel is equal to (Cd+φ) where φ is the linear phase introduced by the actual channel Hk. However, the phase of the actual channel is random with a uniform distribution between 0 and π, which means that Cd is the dominant component in the linear phase of the effective channel. By approximating the channel linear phase to the nearest value in C=[C0C1 . . . CN-CP-l], we obtain C{tilde over (d)} where {tilde over (d)} is the timing offset estimate. The autocorrelation matrix is then given by
{tilde over (R)}HH+E{{tilde over (H)}{tilde over (H)}†} (9)
where {tilde over (H)}=Θ{tilde over (d)}−1
{tilde over (H)}MMSE=Θ{tilde over (d)}F†Θ{tilde over (d)}−1ĤLS (10)
where
F†={tilde over (R)}HH({tilde over (R)}HH+σ2I)−1 (11)
As shown in (10) first the estimated phase rotation is removed from the channel LS estimate Θ{tilde over (d)}−1ĤLS and then the MMSE estimation filter F is applied. Again, if the timing offset estimation is perfect, F will be the same filter used when there are no synchronization errors. Finally, Θ{tilde over (d)} is added. In other words, the timing offset d estimation and then the actual channel MMSE estimation steps are performed separately.
Since the synchronization errors weaken the correlation between the channel coefficients at different subcarriers, the present invention will choose the timing offset resulting in a channel with maximum correlation. Usually, the maximum correlation is obtained when {tilde over (d)}=d. However, even if an error has been made is estimating d, the MMSE estimation will then be performed on a highly correlated version of the channel with less variation, thus improving the MMSE overall performance. First, the timing offset d is found that results in the channel with maximum correlation in the frequency domain
where
R{tilde over (H)}(Δ)=E{{tilde over (H)}k†{tilde over (H)}k+Δ} (13)
is the frequency-domain channel correlation function with a frequency separation Δ, and {tilde over (H)}=Θd−1
The system performance is expressed in terms of Mean Square Error (MSE). First, the case of no synchronization error is considered. The mean square error at a subcarrier k is,
εk=E{|Ĥk−Hk|2} (14)
The average MSE is found to be,
where Tr(.) denotes the trace of a matrix and SH is a diagonal matrix containing E{Hk|2}(which is equal to the diagonal elements of RHH). By substituting (4) into (15) and using the fact that RHH†=RHH,
The channel autocorrelation matrix is obtained as follows,
where p=[p0p1 . . . pl−1] is the channel Power Delay Profile (PDP) vector. The average MSE with synchronization errors can be derived using (16) and (17) with a small modification. The channel PDP will be exchanged with the effective power delay profile pE, where pE is the convolution of the channel PDP and the timing offset Probability Density Function (PDF).
To evaluate the performance of the method of the present invention, a new error term γd is defined, where
γd=d−{tilde over (d)} (18)
The significance of this term is that after removing the estimated phase, the estimated CFR at subcarrier k is
Note the resemblance between (6) and (19). Hence, to evaluate the average MSE for one of the proposed methods, first pγ is obtained, where pγ is the convolution of the channel PDP and the PDF of γd. Next, equations (16) and (17) with pγ instead of p are used to obtain the average MSE. The PDF of γd can be obtained through computer simulations.
In an exemplary embodiment, an OFDM system with N=64 frequency subcarriers and a CP of length NCp=16 samples is considered. Accordingly, a 5-tap channel with an exponential PDP such that the power of the ith path is given by
The timing offset d caused by the synchronization block is approximated as a truncated Gaussian random variable with zero mean and a variance of 4 samples2, where −5≦d≦5. The Gaussian approximation is based on the timing offset's statistical distribution of some of the existing synchronization techniques. In the receiver, an intentional timing offset of 5 samples is added to guarantee non-negative overall timing offset. The average MSE is obtained and plotted against the ratio of energy per bit to the spectral noise density (Eb/N0) as illustrated in
The timing offset PDF for γd at different Eb/N0 values is obtained to evaluate the semi-analytical performance of the methods in accordance with the present invention.
Simulation results are used to confirm the semi-analytical results shown in
Note that at MSE=2×10−3a loss of about 5 dB is caused by the synchronization errors. As shown, a gain of 3.5 dB is achieved using the maximum correlation method in accordance with the present invention and a gain of 2 dB is achieved using the phase estimation algorithm in accordance with the method of the present invention. Furthermore, at low Eb/N0 values, the performance of the phase estimation method in accordance with the present invention degrades. This is expected since for lower Eb/N0, the contribution of the noise to the phase of the estimated channel is higher. However, this is not the case for the maximum correlation method in accordance with the present invention, since the noise (regardless of its level) is uncorrelated across the OFDM subcarriers as long as the orthogonality is maintained. In
The present invention considers MMSE channel estimation for OFDM systems under synchronization errors. The present invention provides a system and method to combat the effect of synchronization errors and restore the performance of the MMSE estimation. In accordance with one embodiment of the present invention, the linear phase rotation is estimated—caused by the timing offset in the receiver—and is removed before applying the MMSE estimation. In an additional embodiment, the timing offset is removed that results in the maximum correlated channel. Both embodiments exhibit a significant improvement in the MMSE estimation performance.
Additionally, a derivation is presented for the performance of the two embodiments. Simulation results are used to check the semi-analytical results. Note that the phase estimation algorithm gain improves for lower noise levels because of the noise effect on the linear phase estimation. On the other hand, the maximum correlation algorithm gain is constant even for high noise levels since the noise is uncorrelated between subcarriers. This, however, comes at the cost of higher computational complexity.
It will be seen that the advantages set forth above, and those made apparent from the foregoing description, are efficiently attained and since certain changes may be made in the above construction without departing from the scope of the invention, it is intended that all matters contained in the foregoing description or shown in the accompanying drawings shall be interpreted as illustrative and not in a limiting sense.
It is also to be understood that the following claims are intended to cover all of the generic and specific features of the invention herein described, and all statements of the scope of the invention which, as a matter of language, might be said to fall therebetween. Now that the invention has been described,
This application claims priority from U.S. Provisional Patent Application No. 60/754,817, filed on Dec. 29, 2005, and titled, “System and Method for Improved Channel Estimation for Wireless OFDM Systems”.
Number | Date | Country | |
---|---|---|---|
60754817 | Dec 2005 | US |