Kalman interpolation filter for channel estimation of LTE downlink in high-mobility environments

The estimation of fast-fading LTE downlink channels in high-speed applications of LTE advanced is investigated in this article. In order to adequately track the fast time-varying channel response, an adaptive channel estimation and interpolation algorithm is essential. In this article, the multi-path fast-fading channel is modelled as a tapped-delay, discrete, finite impulse response filter, and the time-correlation of the channel taps is modelled as an autoregressive (AR) process. Using this AR time-correlation, we develop an extended Kalman filter to jointly estimate the complex-valued channel frequency response and the AR parameters from the transmission of known pilot symbols. Furthermore, the channel estimates at the known pilot symbols are interpolated to the unknown data symbols by using the estimated time-correlation. This article integrates both channel estimation at pilot symbols and interpolation at data symbol into the proposed Kalman interpolation filter. The bit error rate performance of our new channel estimation scheme is demonstrated via simulation examples for LTE and fast-fading channels in high-speed applications.


Introduction
Channel estimation plays an important role in communication systems and, particularly, in the 3GPP Long-Term Evolution (LTE) which aims at continuing the competitiveness of the 3G Universal Mobile Telecommunications System technology.Orthogonal frequency-division Multiple Access (OFDM) is considered as one of the key technologies for the 3GPP LTE to improve the communication quality and capacity of mobile communication system.As the support of high mobility is required in 3GPP LTE systems, the signals at the OFDM receivers are likely to encounter a multi-path, fast time-varying channel environment [1].Thus, good channel estimation and equalization at the receiver is demanded before the coherent demodulation of the OFDM symbols.In mobile communication, since the radio channel is modelled by some dominant spare paths and is represented by path taps, the channel estimation is to estimate and track the channel taps adaptively and efficiently.
In wideband mobile communications, the pilot-based signal correction scheme has been proven a feasible method for OFDM systems.The 3GPP LTE standard employs a Pilot Symbol-Aided Modulation (PSAM) scheme but does not specify the methods for estimating the channel from the received pilot and data signals.
In the 3GPP LTE downlink, pilot symbols, known by both the sender and receiver, are sparsely inserted into the streams of data symbols at pre-specified locations.Hence, the receiver is able to estimate the whole channel response for each OFDM symbol given the observations at the pilot locations.Pilot-symbol-aided channel estimation has been studied [2][3][4] and the common channel estimation techniques are based on least squares (LS) or linear minimum mean square error (LMMSE) estimation [5].Note that most pilot-symbol-aided channel estimators, including those mentioned above, work in the frequency domain.LS estimation is the simpler algorithm of the two as it does not use channel correlation information.The LMMSE algorithm makes use of the correlation between subcarriers and channel statistic information to find an optimal estimate in the sense of the minimum mean square error.
In the literature, based on these two basic estimators, various methods are proposed to improve the performance of the channel estimation.As the LS and LMMSE estimators only give the channel estimate at the pilot symbol, most current work on pilot-aided channel estimation considers interpolation filters where channel estimates at known pilot symbols are interpolated to give channel estimates at the unknown data symbols.Since the 3GPP LTE downlink pilot symbols are inserted in a comb pattern in both the time and the frequency domain, the interpolation is a 2D operation.Although some 2D interpolation filters have been proposed [6], presently, interpolation with two cascaded orthogonal 1D filters is preferred in 3GPP LTE.This is because the separation of filtering in time and frequency domains by using two 1D interpolation filters is a good trade-off between complexity and performance.Various 1D interpolation filters have been investigated.Examples are linear interpolation, polynomial interpolation [7], DFT-based interpolation [8], moving window [9] and iterative Wiener filter [10].
From a system point of view, the channel estimation is a state estimation problem, in which the channel is regarded as a dynamic system and the path taps to be estimated are the state of the channel.It is known that Kalman filter (KF) provides the minimum mean square error estimate of the state variables of a linear dynamic system subject to additive Gaussian observation noise [11].By considering the radio channel as a dynamic process with the path taps as its states, the KF has shown its suitability for channel estimation in the time domain [1].In the frequency domain, Kalman-based channel estimator in OFDM communication has also been studied [1,12,13].For example, in [1,12], a modified KF is proposed for OFDM channel estimation where the time-varying channel is modelled as an autoregressive(AR) process and the parameters of the AR process are assumed real and within the range [0.98, 1] for slow-fading channels.However, in the high-mobility environment, these parameters are relative large (e.g. in the 200 km/h environment, they are complex values with magnitudes varying in [0, 1.5]) representing a fastfading channels.
The difference between the KF in [12] and the one proposed in this article is that the former estimated the parameters of AR by a gradient-based recursive method separately, rather by the linear KF.Whereas, we derive an extended Kalman filter (EKF) for jointly estimating the channel response and the parameters of the AR model simultaneously.In addition, the parameters of the AR model are assumed time-invariant and known in priori by solving Yule-Walker equation in [1].The authors of [13] only considered the comb-type pilot patterns in which some subcarriers are full of pilot symbols without unknown data.As a result, the KF in [13] requires continuous stream of pilot symbols and is not suitable for 3GPP LTE, as the 3GPP LTE employs a scattered pattern where the pilot symbols are distributed sparsely among the data streams.
Although the KF-based channel estimation for LTE uplink has been reported recently [1], there has been no KF-based joint estimation of both time-varying channel taps and the time-correlation coefficients of 3GPP LTE downlink in frequency-time domain.This article focuses on the major challenge of scattered pilot-aided channel estimation and interpolation for a time-varying multipath fast-fading channel in 3GPP LTE downlink.An AR process is used to model the time-varying channel.Both the taps of the multipath and the timecorrelation coefficients are jointly estimated by treating the channel as a nonlinear system.Then, a combined estimation and interpolation scheme is present under the EKF framework.
The main contribution of the proposed method is (1) both the time-correlation coefficients and channel taps are estimated simultaneously in the framework of EKF; (2) no assumption on the upper/lower boundaries of the time-correlation coefficients to achieve a good tracking of fast-fading channel in high-mobility scenario; (3) applicable to preamble pilot patterns, comb-type pilot patterns and scattered pilot patterns.
This article is organized as follows: Section "System model" gives an overview of the LTE 3GPP downlink system and formulates its channel estimation problem.In Section "EKF for channel estimation", an EKF is derived by using a first-order Taylor approximation for the joint estimation of channel taps and time-correlation coefficients at pilot symbols.Section "EKF for channel interpolation" describes the combined estimation and interpolation scheme and summarizes the proposed algorithm.Simulation results of the proposed Kalman interpolation filter are presented and its performance is demonstrated in Section "Simulation results and performance analysis".

Notation and terms
Unless specified otherwise, an italic letter (e.g.T, h k,n ) represents a scalar and its bold face lower-case letter represents its corresponding vector (e.g. ).A bold face upper-case letter (e.g.A) represents a matrix.The subscriber k denotes the time index of an OFDM symbol, n denotes the index of subcarriers in the frequency domain, l denotes the lth path of the radio channel.x j j ( A j j) is the elementwise magnitude of a vector x (matrix A).I N is an N × N identity matrix.A i,j denotes the entry at the ith row and the kth column of A.
L denotes the total number of possible paths in a radio channel, referred to as channel length, N denotes the total number of subcarriers, N p the number of pilot subcarriers, g k;l the channel impulse response (CIR) of lth path at kth symbol, referred to as tap, g k the CIR vector

System model
Figure 1 describes the LTE downlink baseband system used in this article.Here, we only consider baseband processing and omit all analogue components, higher layer protocols and application processing.The baseband processor receives the digitized signal as complex samples from the analogue-to-digital convertors and posts the decoded data stream to the higher layer protocol and the application processor.

Pilot symbols in LTE downlink
As depicted in Figure 1, a radio frame of the LTE downlink has duration of 10 ms and consists of ten subframes each of 1 ms.Each subframes has two 0.5-ms time slots with each slot consisting of N DL symb OFDM symbols (the values of N DL symb for various configurations are given in Table 1).The transmitted downlink signal is represented as a time-frequency resource grid.Each small box within the grid represents a single subcarrier for one symbol period and is referred to as a resource element.Note that in MIMO applications, there is a resource element mapping graph for each transmitting antenna.A resource block (RB) is defined as consisting of N RB sc consecutive subcarriers for one slot (N DL symb OFDM symbols).An RB is the smallest unit of bandwidth-time resource allocation assigned by the base station scheduler, and the specification for the parameters of one RB is shown in Table 1.
In order to successfully receive a data transmission, the receiver must estimate the CIR to mitigate the multi-path interference.In packet-oriented networks (like IEEE 802.11), a physical preamble is used to facilitate this purpose.In contrast to 802.11, LTE makes use of PSAM, where known reference symbols, referred to as pilot symbols, are inserted into the stream of data symbols, as shown in Figure 1.Generally, there are three kinds of time-frequency allocation pattern of pilot symbols, namely, entirely known OFDM symbols, pilot subcarriers and scattered pilots. 3GPP LTE adopts a scattered pattern involving the sparse insertion of known pilot symbols in a data symbol stream.For example, in the scenario of a single transmitting and a single OFDM Symbol k 1 2 3 4 5 6 7 8 9 ... k k  receiving antenna, pilot symbols are transmitted at the first and the fifth OFDM symbols of each slot at the pilot subcarriers.In the frequency domain, reference signals are spread over every six subcarriers.The effect of the channel response on the known pilot symbols can be computed directly by calculating the attenuation of each pilot symbol [5].For the remaining unknown data symbols, interpolation has to be used to estimate the channel response among adjacent pilot symbols.A simple way of performing this interpolation is the linear approximation in both time and frequency.The concept of PSAM in OFDM systems allows the use of both the time and frequency correlation properties of the channel to improve the channel estimation.Therefore, an efficient channel estimation procedure may apply a complicated 2D time-frequency interpolation or a combination of two simple 1D interpolations [6] to provide an accurate estimation of the channel states for each OFDM symbol.

Channel model
In this article, we consider an LTE downlink system with N subcarriers over a Rayleigh-fading channel.For the purpose of analysis, the following notation and assumptions are taken in this article.
(1)The system bandwidth is B = 1/T, where T is the duration of one time-chip.The duration of one OFDM symbol is T s ¼ NT þ T CP , where T CP is the duration of cyclic prefix (CP) for every OFDM symbol.
(2)The number of possible path is L and the maximum delay due to multi-path is (L -1)T.(3)The length of CP is carefully designed to eliminate inter-symbol interference between consecutive OFDM symbols.That is T CP is longer than the than the channel's maximum delay, ðL À 1ÞT < T CP .(4)The Rayleigh-fading channel varies in consecutive OFDM symbols, but is assumed constant within one OFDM symbol.
The time-varying multi-path channel can be represented in the continuous time-domain function by a collection of paths gðt; τÞ ¼ where the lth path is represented by a tap with complex amplitude α l (t) at time instant t and a delay τ l .The impulse response of the physical channel consists of independent Rayleigh-fading impulses, uniformly distributed over the length of the CP.
In the OFDM implementation of the 3GPP LTE, the transmitted and received signals are sampled for D/A and A/D conversion with an interval of chip duration T, the CIR (1) in the continuous time-domain is converted into an equivalent discrete channel model with sampling interval T. We define g k;l ¼ α l ðkT Þ as representing the complex magnitude of the lth path with delay lT during the kth OFDM symbol.The equivalent discrete model of the radio channel ( 1) is therefore Hence, the discrete CIR model can be represented by a length-L CIR vector g k Strictly speaking, g k is only an approximation of g t; τ ð Þ at kth OFDM symbol ( kTs < t < ðk þ 1ÞTs ).When the multipath taps do not fall in the discrete sampling grid (i.e., τ l 6 ¼ lT ), the discrete-time CIR vector will be infinite in length.However, the pulse's energy decays quickly outside the neighbourhood of the original pulse location [5,14], it is still feasible to capture the impulses with a length-L vector.In this study, we assume that the tails of the impulse response function are negligible beyond L samples, which is also the assumption made in OFDM to justify that no ISI occurs.
In the frequency-domain, the frequency response of the CIR g k at kth OFDM symbol is where h k ½n f g; n 2 1; 2; . . .N f g , denoting the CFR of nth subcarrier at kth OFDM symbol time is converted from the time-domain CIR via the discrete Fourier transform (DFT) The relationship between the CIR in time-domain and CFR in frequency domain can be described in matrix notation where F L 2 C NÂL is the first L columns of the N × N DFT matrices F. And F is denoted by It has been shown that time-varying path taps in a fading channel can be modelled by an AR process [11,15], which is applicable to general fading channels, and in particular to mobile communication.Examples include the first-order AR model in [1,11,16] and the secondorder AR model [15].Although the first-order AR model is just an approximation to the actual statistics of the random radio propagation process, it is more realistic than those models assuming constant channel parameters (identity matrix) or using linear interpolation.Furthermore, the use of higher-order models will lead to higher computational costs, which may not be justified by the performance improvement.Compared to the higher-order model, a lower-order model may reduce the overall computational complexity significantly with only a relatively small performance sacrifice.Here, we are concerned with the basic derivation of the proposed Kalman interpolator filter in LTE downlink.As shown in our following derivation, higher-order models can also be incorporated into the proposed scheme with only minor modifications.For the purpose of analysis, we restrict ourselves to a firstorder AR model for the time-varying channel.
It is easy to verify that the channel coefficients h k of the time-varying CFR can be modelled by the following dynamic AR process [1,11,12]: where α n represents the time correlation of the channel response between kth and (k + 1)th OFDM symbols at the nth subcarrier.v k;n is a mutually independent zeromean Gaussian complex white noise representing the modelling error.

LTE OFDM reception and channel estimation
In order to estimate the CFR as defined in (4), N p pilot symbols are inserted sparsely among N subcarriers at kth OFDM symbol duration following the comb pattern shown in Figure 1.Let x k ¼ x k;1 ; x k;2 ; . . .; denote the transmitted pilot vector of N p known pilot symbols at the kth OFDM symbol, y k ¼ y k;1 ; y k;2 ; . . .; y k;N p Â Ã T 2 C N p denotes the vector of the received pilot symbols.After CP removal, the received pilot symbols can be expressed as where X k ¼ diagðx k Þ is an N p × N p diagonal matrix with transmitted pilot symbols x k as its diagonal elements, Here, w k 2 C N p is an additive white complex Gaussian noise with covariance matrix σ 2 w I N p and h k 2 C N p is the CFR at pilot subcarriers at k-th OFDM symbol.
The goal of channel estimation is to estimate the whole CFR for all data carriers from h k at these N p pilot symbols with as high accuracy as possible.This is an optimization problem described as arg min where y k =h k is an element-wise division with elements y k,n /h k,n .
It is worth noting that, as the pilot symbols in LTE downlink are inserted into the data symbols sparsely in a frequency-time scatter pattern, the channel response at data symbols are typically interpolated from the channel estimates at pilot symbols.As shown in literature, if the OFDM symbol is short compared with the coherence time of the channel, the time correlation between the channel attenuation of consecutive OFDM symbols is high.There is also a substantial frequency correlation between the channel attenuation of adjacent subcarriers.For a better channel estimation at data symbols, both of these time and frequency correlation properties of the fading channel can be exploited by the channel estimator.
Generally, as illustrated in Figure 1, the whole process of such a pilot-aided channel estimation includes three steps: (1) Estimation at pilot symbols, where, h k , the channel responses at N p pilot subcarriers at kth OFDM symbol are calculated with the common LS estimator or LMMSE estimator; (2) Time-domain interpolation, where the channel responses h k+1 at (k + 1)th OFDM symbol at pilot subcarriers are estimated from h k by tracking the parameters of each path.(3) Frequency-domain interpolation, where the channel responses at all N subcarriers are estimated by interpolating or smoothing these estimates {h k, h k+1 ,. ..} at pilots subcarriers.This article integrates the first two steps into one framework called the Kalman interpolator filter.

EKF for channel estimation
In this section, we are interested in deriving a minimum variance estimator/interpolator for the channel response {h k, h k+1 ,. ..} at pilot subcarriers from the observation of sparse pilot symbols.We present a combined estimation and interpolation scheme, where the time correlation among consecutive OFDM symbols is taken into account to estimate the CFR at the known pilot symbols and then to interpolate to estimate the CFR at the unknown data symbols at the pilot subcarriers.The proposed scheme is based on the idea of Kalman filtering to improve the accuracy of the estimation and interpolation.More specifically, recalling the LTE reception model y k ¼ X k h k þ w k in (9), the task for the Kalman interpolator filter can be stated as: Given the matrix X k of known transmitted pilot symbols and received signal y k at kth OFDM symbol, to obtain minimum variance estimates of the time-varying multipath CFR h k and interpolate h k to the followed six g data symbols at the pilot subcarriers until the next pilot symbol (k + 7-th OFDM symbol) is received.

Augmented state space model
Considering a time-varying channel described in Equation (8), the CFR at pilot subcarriers can be described as a state space model where h k is the state variable to be estimated, A k 2 C N p ÂN p is the unknown state transition matrix consisting of the time correlation coefficients α n of channel response.Both v k and w k are mutually independent, zero-mean, Gaussian complex white noises, with covariance Q v ≜σ 2 w I N p and Q w ≜σ 2 w I N p , respectively.It is assumed that v k and w k are independent of the state variable h k .Note that, in this state space model of the CFR, the state transition matrix A k is unknown and to be estimated together with the state variable h k .Therefore, it is a problem of joint state and parameter estimation.The purpose is to estimate both the channel response h k and channel's time-correlation matrix A k from the received pilot symbols y k .
Considering that A k is a spares matrix in most cases, without loss of generality, we assume A k has N A N A < N p Â N p À Á unknown entries to be estimated and let a vector a k denote all the N A unknown entries as follows:

column vector formed by stacking
all unknown entries of the matrix A k in a row-wise order.The time-correlation parameters are now represented by a k which is the vector to be estimated.For the purpose of clarification, A k is represented by A(a k ) explicitly in the following.Assuming a random walk model for the parameter a k , then Equation ( 12) becomes where E k denotes the process noise of a k and is an independent, zero-mean Gaussian noise with covariance In order to jointly estimate the state and parameters, a new augmented state z k is defined as and the channel state space model Equation ( 14) turns into an augmented system where EKF Since the state transition function f(z n ) in the augmented state model ( 16) is a nonlinear function and an EKF has to be used to estimate the augmented states.The development of the EKF basically consists of two procedures: linearizing the augmented model ( 16) and applying the standard KF to the linearized model.The linearization procedure is included in the Appendix where the derivation of the EKF algorithm for a general matrix A(a k ) is demonstrated.The basic concept is to form the Taylor approximation of the nonlinear transition function.The resulting linear state space model approximating the AR model ( 12) is Applying the standard KF to the model ( 18) is straightforward.The resulting EKF algorithm for the joint estimation of CFR h k and CFR's time correlation coefficients a k works in an iterative prediction-correction cycle.The prediction projects forward (in time) the current estimate ẑk and error covariance P k at kth OFDM symbol to obtain the a priori estimates ẑkþ1 k j and P kþ1 k j for the next (k + 1)th OFDM symbol.The correction adjusts the projected estimates ẑkþ1 k j and P kþ1 k j to obtain an improved a posteriori estimate by using an actual measurement of received symbol y k+1 at (k + 1)th OFDM symbol.Here, the subscript k þ 1 k j corresponds to one-step a priori prediction, k k j corresponds to a posteriori correction and k k j is denoted by k for the purpose of short notation.More specifically, the filtering algorithm is presented as follows: 1. Prediction (before receiving a OFDM symbol): where and Correction (once the reception of the OFDM symbol has completed): Here, K k is the Kalman gain of the EKF.The EKF makes use of a first-order Taylor approximation of the state transition and thus does not approach the true minimum variance estimate when the linearization error is non-negligible.Nevertheless, the resulting EKF is a practical approximation to the minimum variance estimator when the state equation is nonlinear, and will be shown to provide a good performance in time-varying channel estimation.Furthermore, the EKF has successfully been applied to the problem of joint channel state and parameter estimation in [11,16], and thus it seems reasonable to apply EKF to the time-varying channel estimation.

Remark
In terms of computation complexity, it can be seen that prediction of state error covariance P kþ1 k j and the update of K k consumes the major amount of computation.Fortunately, in general, cross-path coupling is confined within a small neighbourhood, and thus the off-diagonal elements of A k representing the coupling between multiple paths are small and may be neglected.As shown in the AR model ( 8) of time-varying channel, the channel's time-correlation matrix A k can be modelled as a diagonal matrix.If both X k and A k are diagonal matrices, the number of complex multiplications and additions is be reduced to a great extent.More specifically, the number of multiplication and division operations in Equations ( 19)-( 24) is 25N p .

EKF for channel interpolation
In this section, the proposed EKF is further extended to interpolate the CFR estimate to unknown data symbols and the whole estimation and interpolation process of the proposed EKF is summarized.
Figure 2 illustrates the block diagram of the baseband channel model and the proposed Kalman interpolation filter for LTE downlink channel equalizer.The EKF works in an iterative prediction-correction manner and, in the application of LTE downlink channel estimation, each iteration corresponds to the duration of an OFDM symbol.However, due to the fact that the known pilot symbols are inserted sparsely into the unknown data symbols, the coefficient matrix X k is not always available at each iteration.Like most adaptive algorithms, two working modes, namely, training mode and decisiondirected interpolation mode, are adopted in the proposed Kalman interpolation filter to address this issue.
The estimator is trained during these periods when a pilot symbol is received.Then it switches to an interpolation mode, in which a decision-directed method is applied to estimate the channel response until the next pilot symbol is received.During the training period, the transmitted symbols X k are known to the estimator, while in the data symbols periods, the transmitted data symbols are estimated as Xk by the decoder and the EKF is fed by the Xk to replace the unknown transmitted symbols X k.Indeed, the channel estimator is fed with one pilot symbol and six estimates of the data symbols in one LTE slot.The proposed Kalman interpolator filter method yields an adaptive algorithm and can be implemented recursively.
At each iteration, the equalizer and the decoder compute an estimate Xk of the transmitted data symbols on the basis of the previous, a priori channel estimate ĥk kÀ1 j .In the iteration of the OFDM data symbol, Xk is also fed to the EKF to calculate a posteriori channel estimate ĥk k j and a priori channel estimate ĥkþ1 k j .By exchanging their estimates, both EKF and equalizer are able to improve their performance iteratively.This is particularly useful at these iterations of unknown data symbols.

Initialization by LS estimation
Although a KF is able to convergence under any reasonable initial value of the state variable z k , a good initial condition will reduce the duration of convergence.Generally, if the initial value of the state variable is set to the neighbourhood of the true value, a faster convergence can be obtained.Since the state variable z k consists of two independent components, a k and h k , their initial values are chosen separately.
For initializing the channel's time-correlation coefficients a k , we use an identity matrix (A a 0 ð Þ ¼ I N A ) assuming the channel response at next OFDM symbol is the same as the current OFDM symbol.Although an identity matrix represents a time-invariant channel, an identity matrix would be the best choice of the channel's initial condition, given we have no a priori knowledge about the channel.
For CFR h k , we shall use the conventional version of an LS estimation to get the initial value.When the first group of pilot symbols is received, the LS method is performed as follows: where ĥ0;LS 2 C N p is the initial CFR estimate.A more complicated LMMSE estimator using the channel's frequency correlation may be applied to obtain a more accurate initial estimate of the CFR.It should be pointed out that the EKF is initialized until the first group of pilot symbol is received.

Trained estimation
After the state variable is initialized, the EKF works iteratively either in the training mode or in the interpolation mode.During the pilot symbols, the EKF switches to the training mode, where the known pilot symbol forms the matrix X k .As the observation y k is obtained by the DFT at the end of an OFDM symbol duration, the a posteriori CFR ĥk k j is first estimated from y k by using update equations ( 22)-( 24).Then the a priori estimate ĥkþ1 k j is calculated by Equations ( 19)-(21) for next OFDM symbol.

Decision-directed interpolation
During periods where the pilot symbol is not available, the EKF switches to decision-directed interpolation mode to continue adaptation.For these data symbols, as the transmitted symbol X k is unknown, X k is replaced by the decoder's decision of X k that is supposed to be nearest to X k .In the decision-directed mode, the prediction and correction processes are the same as the training mode, except X k is replaced by X k ,.
It is worth noting that, as y k is only available at the end of the current symbol duration, the correction process has to be carried out at the end of the symbol duration.Thus, the equalizer uses the a priori channel estimate ĥk kþ1 j to refine the currently received OFDM symbol, rather than uses the a posteriori CFR ĥk k j .

Selection of the covariance matrices
In most applications of Kalman filtering, it is difficult to measure the variance of noises.In practice, the covariance matrices are a priori approximated by applying the best available knowledge and tuned empirically in the application.As shown in the state space model ( 12), the channel measurement y k is subject to the noise w k , the additive white complex Gaussian noise in the wireless channel.Since the transmission power and signal-to-noise ratio (SNR) are usually available in a communication system, the elements of the variance matrix Q w can be calculated by σ 2 w ¼ P tx 10 S N R=10 , where P tx is the transmission power measured in Watts and SNR is in dB.Presuming a small process variance and linearization errors in (18), the values of σ 2 v and σ 2 E in Q v are empirically selected from {0.1, 0.01, 0.001} according to the SNRs.At low SNRs, the channel estimate is less accurate due to large observation noise and thus a larger value is used for Q u .At higher SNRs, a better channel estimation is expected and a smaller value is used for Q u .

Summary
We now summarize the proposed method for channel estimation in LTE downlink: Step 1. Initialize ½a T 0 h T 0 T when the first pilot symbol is received, make the first a priori prediction ½a T 1j0 h T 1j0 T for next OFDM symbol and set k = 1; When a new OFDM symbol (kth symbol) has been received, repeat the following steps 2-6.
Step 2. Calculate y k by using DFT Step 3. Estimate xk by equalizing y k with previous a posteriori h k kÀ1 j ; Step 4. If y k is pilot symbol, set X k by the known pilot symbol x k ,else set X k by the estimated data symbol xk , Step 5. Correct a posteriori state estimation ½a T kjk h T kjk T from y k by ( 22)-(24).
Step 7. Frequency-interpolation: The CFR at data subcarriers for next symbol is interpolated using a DFT-based interpolation [8].
Step 8. k = k + 1, wait for next symbol and goes back to step 2.
It can be seen that, the proposed KF-based channel estimation scheme is a combination of the estimator (for pilot symbols) and the interpolator (for data symbols).When the pilot symbol is available at kth iteration, a direct observation of the channel state is obtained and the EKF works at the training mode giving the optimal estimate of CFR in the sense of minimum variance.In the followed six {k + 1, k + 2,. ..,k + 6} data symbols, the EKF interpolates the CFR in decision-directed model until the next pilot symbol (k + 7-th OFDM symbol) is received.

Simulation results and performance analysis
In this section, simulation is performed to validate the performance of the proposed Kalman interpolation filter for LTE downlink systems.A simplified rural area model defined by 3GPP [17] is adopted to configure the Rayleigh channel with additive white Gaussian noise and the parameters are listed in Table 2.The LTE downlink simulation parameters are listed in Table 3.The total number of subcarriers is 512 with 300 of them used for data/pilot transmission and a quadrature phase-shift keying (QPSK) modulation employed.For simplicity, the raw bits randomly generated are not coded with turbo coding schemes.In the single input single output scenario, 100 of the 300 subcarriers are used for carrying pilot symbols during the pilot OFDM symbol time period.Three speeds of user equipment are simulated, namely, 50 200 and 300 km/h.For each speed, the simulation is repeated 10 times (10 runs) in order to obtain reliable statistics and each simulation run simulates the transmitting/receiving four LTE downlink subframes containing 56 OFDM symbols.
Figure 3 illustrates a plot of the channel surface for the urban channel model at a moving speed of 200 km/h, where the Doppler frequency is around 480 Hz.This plot shows the time-varying and frequency-selective nature of the channel gain and provides an image of the true values of the CFR.Studying the channel surface indicates that fluctuations in the frequency are clearly visible, but relatively smoothly varying in time, implying an AR process would be able to represent the channel's time correlation.The channel surface also suggests that a linear interpolation may not be good for such a nonlinear CFR.
The simulations are carried out at different noise levels with the SNR varying from 0 to 40 dB at a step size of 5 dB. Figure 4 shows an example of the CFR estimation errors at the 100 pilot subcarriers in one simulation run (with 4 subframes containing 56 OFDM symbols) at SNR = 20 dB and moving speed 20 km/h. Figure 4a depicts the CFR surface estimation error given by the proposed EKF scheme and Figure 4b depicts the  estimation errors of the LS scheme, where the improvement of CFR estimation in the proposed scheme can be seen clearly and the mean square error (MSE) of the EKF is 0.066 and that of the LS is 0.09.The smaller CFR estimation error demonstrates the proposed EKF's ability to filter the noises in observation and to track the timevarying channel parameters.Particularly, towards the edge of the LTE downlink subframe (i.e. at the 14th, 28th, 42th, 56th OFDM symbols), a larger estimation error occurs in the LS estimation which can be seen clearly in Figure 4b.This is caused by the extrapolations errors in the LS scheme as no pilot symbols are inserted at the edge of each subframe.However, it is worth noting that the proposed EKF's estimation errors may have peaks at some data symbols, due to the fact that the incorrect data symbol decision is fed back to the EKF in the decision-directed mode.If a large deviation occurs and thus makes the received OFDM data symbol far from  Figure 5 shows the average CFR estimation MSEs of the LS and EKF schemes at different SNRs.It can be seen that the EKF achieves a smaller estimation errors and gives a better CFR estimation.
The BER performances are plotted in Figures 6, 7 and 8. Figure 6 is for low-speed environment (50 km/h), Figures 7 and 8 are for high-speed environment (200 and 300 km/h), respectively.In the BER comparisons, Figure 6 BER performance comparison at 50 km/h (A0 denotes a perfect channel estimation algorithm where the actual CFR is known to the receiver).Blue solid line with triangle: the BER performance of LS method.Red solid line with star: the BER performance of the proposed EKF method.Blue solid line with circle: the BER performance of LMMSE method.Blue solid line with square: the BER performance of the perfect channel estimation method.The perfect channel estimation means that the actual CFR is known to the receiver.
the popular (LMMSE) algorithm [5] is also employed.Note that since the proposed EKF method is based on the LS estimation, we explicitly denote it by 'EKF with LS' in the legends of these figures.As a performance benchmark, the BER performances of a perfect channel estimation algorithm (denoted by A0 in Figures 6, 7 and 8) are also depicted, where the perfect channel estimation refers to the actual CFR being known by the receiver in advance.
As expected, A0 gives the best performance among all of the three methods, since it has the perfect CFR.The BER performance of A0 can be regarded as the BER's lower bound.Obviously, the LS method has the poorest BER performance in all these three scenarios and the LMMSE is able to improve the BER performance.It can be seen that BERs of the proposed Kalman interpolation filter fall between the LMMSE's performances and the performances of perfect channel, although the EKF shows a slightly higher BER than LMMSE in low SNRs (i.e.0 and 5 dB).It is worth noting that the EKF is always better than the LS method.This is to be expected since the concept behind the observation equation in the proposed EKF method is the Figure 7 BER performance comparison at 200 km/h (A0 denotes a perfect channel estimation algorithm where the actual CFR is known to the receiver).Blue solid line with triangle: the BER performance of LS method.Red solid line with star: the BER performance of the proposed EKF method.Blue solid line with circle: the BER performance of the LMMSE method.Blue solid line with square: the BER performance of the perfect channel estimation method.The perfect channel estimation means that the actual CFR is known to the receiver.same as the LS method, where it assumes the CFRs at adjacent pilot subcarriers are independent.Nevertheless, compared to the LS estimation, the proposed Kalman interpolation filter shows a significant improvement.This is particularly obvious at high SNRs and high-speed environment.As seen in Figure 7, when using the proposed EKF instead of the LS estimator, a gain in SNR up to 8 dB can be obtained for certain BERs (e.g.0.002) at high-speed application.The average SNR gain is about 3-5 dB.

Conclusions
This article focuses on channel estimation and interpolation for a time-varying multipath fading channel in 3GPP LTE downlink.The time-varying radio channel is modelled as an AR process represented in state space form and an EKF is developed for the purpose of both channel estimation at pilot symbols and interpolation at data symbols.The time-varying channel estimation is a joint state and parameter estimation problem, where both the channel taps and AR parameters need to be estimated simultaneously to achieve an accurate channel estimate.We convert the state model into an augmented system and a corresponding EKF is proposed.Furthermore, the interpolation channel estimate at data symbols are also integrated into the EKF and the proposed Kalman interpolation filter shows a good performance of estimating a time-varying channel in the 3GPP LTE downlink.

Appendix
Applying the first-order Taylor approximation to the nonlinear state transition function f(z n ) around ẑn n j in Equation ( 17), the state equation ( 16) becomes where μ n is the linearization error μ n ¼ f ẑnjn À Á À F n ẑnjn and Here, we assume that A(a k ) is independent of h n .Recalling the definitions of a k and A(a k ), it is easy to verify that @Aða T n Þ @a n h n is an N p × N A block-diagonal matrix H n where denotes the Kronecker product and the operator SðHÞ removes these known-columns of I N p h T n .The known column is the ði À 1ÞN p þ j th column when the ith-row-jth-column entry A i;j of A is known.Hence, substituting z n ¼ ẑn n j into Equation (26), we have And the linear state space model approximating the AR model ( 12) is

Figure 1
Figure 1 LTE downlink frame structure and the time-frequency allocation of pilot symbols (one transmitting antenna scenario).

2
Channel estimation and interpolation at pilot subcarriers.

Figure 3
Figure 3 Channel surface of an LTE radio channel at a moving speed of 200 km/h.

Figure 4
Figure 4 CFR estimation errors (a) for the proposed EKF scheme and (b) for the LS scheme.

Figure 5
Figure 5 Average mean square estimation errors at various SNRs (200 km/h).Blue solid line with triangle: the average MSE of channel estimates by LS method.Red solid line with star: the average MSE of channel estimates by the proposed EKF method.

Figure 8
Figure8BER performance comparison at 300 km/h (A0 denotes a perfect channel estimation algorithm where the actual CFR is known to the receiver).

Table 1
Physical RBs parameters

Table 2
The simplified rural area channel model

Table 3
LTE downlink simulation parameters