Patent Application
20070133696
The general field is a Kalman filter and more specifically an application of Kalman filtering algorithm in an OFDM based communication system.
Orthogonal Frequency Division Multiplexing OFDM has been widely applied in wireless communication systems such as DVB-T/H, 802.11x wireless LAN and 802.16 wireless MAN due to its high bandwidth efficiency and robustness to multipath fading. Due to the fact that the wideband wireless channel is frequency selective and time varying, channel estimation must be performed continuously and the received OFDM subcarriers must be corrected by the estimated CTFs.
In DVB systems, channel estimation is performed by inserting known scattered pilots at predefined subcarrier locations in each OFDM symbol (see ETSI EN 300 744 V.1.4.1 “Digital Video Broadcasting (DVB): Framing Structures, channel coding, and modulation for digital terrestrial television”), normally referred to as “comb-type” pilot channel estimation.
As shown in
For the symbol of index l (ranging from 0 to 67), carriers for which index k belongs to the subset {k=Kmin+3×(l mod 4)+12p|p=int, p≧0, kε[Kmin; Kmax]} are scattered pilots.
Where p is an integer that takes all possible values greater than or equal to zero, provided that the resulting value for k does not exceed the valid range [Kmin; Kmax].
Assume that for current symbol of index l, the Ns inserted scattered pilots according to the above rule are: Xs, s=0, 1, . . . Ns−1, the corresponding received subcarriers at the scattered pilot locations are: Ys, s=0, 1, . . . Ns−1, the channel frequency response at the pilot subcarrier locations can be represented as: Hs, s=0, 1, . . . Ns−1, then the Least-Square estimate of the channel frequency response at the pilot subcarrier locations is given by:
s=0, 1, . . . Ns−1
The above LS estimation is sensitive to noise and ICI, MMSE estimation is known to provide better performance than LS estimation. However, MMSE estimation requires matrix inversion at each iteration, thus not practical for implementation.
In this disclosure, a simplified Kalman filter is provided which reduces the noise effects of the LS estimation. Simulation shows that the simplified Kalman filter is very effective in removing the noise effects, and the overall system performance will be improved by up to 2 dB.
A scalar Kalman filter is applied for a Least-Square estimated value Hs at s. The filter has an input for receiving Hs, a filter equation and an output for the corrected estimated value Hsk for the kth variable. The filter equation is Hsk=KgainSn[k] wherein: correction Sn[k]=S+Kn (Hs−S); prediction of the correction S=KaSn[k]; Kalman filter gain Kn=P/(1+P); minimum predication MSE P=Ka2Pn[k]+Kb; minimum MSE Pn[k]=P(1Kn); and Ka, Kgain and Kb are constants.
A receiver includes an OFDM demodulator, a channel corrector and a channel estimator; and wherein the channel estimator is a Least-Square estimator of a channel frequency response Hs of a subcarrier k. The channel estimator includes the simplified Kalman filter. The constants Kgain, Ka and Kb may be selected as a function of the modulation mode of the subcarrier. The channel estimator processes scattered pilots whose locations s repeats its pattern every r symbols; and the channel estimator includes r*Ns Kalman filters.
The filter gain Kn may also be a constant selected as a function of the modulation mode of the subcarrier. The filter equation is performed in software.
These and other aspects of the present disclosure will become apparent from the following detailed description of the disclosure, when considered in conjunction with accompanying drawings.
In the proposed channel estimation algorithm, for the current OFDM symbol of index l (ranging from 0 to 67), the Least-Square estimate of the channel frequency response at the pilot subcarrier locations
s=0, 1, . . . Ns−1 can be further processed by individual Kalman smoothing filters to reduce the noise and ICI effects before the LS estimate is used in time/frequency domain interpolation.
Since the scattered pilot location repeats its pattern every 4th symbol, there will be a total of (Ns,l mod 4+Ns,(l+1)mod 4+Ns,(l+2)mod 4+Ns,(l+3)mod 4) individual Kalm filters. If a different repeat pattern is used, the number of Kalman filters would match the repeat pattern.
The general form of a scalar Kalman filter is described in the following equations (see “Fundamentals of Statistical Signal Processing Estimation Theory,” Steven M. Kay, PTR Prentice-Hall, Inc., 1993).
Prediction: First order Markov process: ŝ[n|n−1]=aŝ[n−1|n−1]
Minimum prediction MSE: P[n|n−1]=a2P[n−1|n−1]+σu2
Kalman Gain:
Correction: ŝ[n|n]=ŝ[n|n−1]+K[n](x[n]−ŝ[n|n−1]) Where x[n] is the input data at the nth iteration
Minimum MSE: P[n|n]=(1−K[n])P[n|n−1]
The simplifying assumptions of the present design are:
all the scattered pilots are supposed to be uncorrelated, characterized by the first order Markov process;
the measurement noise variance σn2 is supposed to be the same for all carriers;
the signal noise variance σu2 is supposed to be the same for all carriers; and
from the previous assumptions since the Kalman gain converges to a constant after a few iterations, it will be assumed to be a constant.
Based on the simplifying assumptions, every pilot carrier will be filtered independently by a scalar Kalman filter. The scalar Kalman filter equations, for the kth Kalman filter, corresponding to the scattered pilot located at the kth subcarrier, become:
S=KaSn[k]
P=Ka2Pn[k]+Kb
Sn[k]=S+Kn(Ĥs−S)
Pn[k]=P(1−Kn)
Hsk=KgainSn[k]
Where the filtering will take place on symbol bases.
{k=Kmin+3×(l mod 4)+12p|p=int, p≧0,kε[Kmin; Kmax]}.
Ĥs is the LS estimate on the scattered pilot location s, and Hsk is the Kalman filter smoothed output of Ĥs.
In the above equations, Ka, Kb and Kgain are constants and the calculation of Kalman gain factor Kn requires divisions. It is found that Kn will converge to its steady state value over a few OFDM symbols. In order to simplify the implementation, Kn is also set as a constant. The above constants can be set by evaluating the performance in various multipath fading channels and noise conditions. It has also been found that one may set the Kalman constants differently for different modulation modes, such as QPSK, 16 QAM and 64 QAM to achieve better smoothing performance. In a preferred embodiment, the Kalman filter constants are set according to the modulation mode whenever a valid tps frame is decoded.
Although the present disclosure has been described and illustrated in detail, it is to be clearly understood that this is done by way of illustration and example only and is not to be taken by way of limitation. The scope of the present disclosure is to be limited only by the terms of the appended claims.
