SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Advanced Technologies for LTE Advanced.

Open Access Highly Accessed Open Badges Research

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

Xuewu Dai1, Wuxiong Zhang234, Jing Xu34, John E Mitchell1* and Yang Yang34

Author Affiliations

1 Department of Electronic & Electrical Engineering, University College London, Torrington Place, London, WC1E 7JE, UK

2 Graduate University of Chinese Academy of Sciences, Beijing, 100049, China

3 Shanghai Research Center for Wireless Communications (WiCO), Shanghai, 200050, China

4 Shanghai Institute of Microsystem and Information Technology (SIMIT), Shanghai, 200050, China

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:232  doi:10.1186/1687-1499-2012-232

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2012/1/232

Received:23 August 2011
Accepted:11 June 2012
Published:25 July 2012

© 2012 Dai et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


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.

LTE advanced; Channel estimation; Extended Kalman filter; Pilot-aided-interpolation


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-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 fast-fading 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 time-correlation 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, hk,n) represents a scalar and its bold face lower-case letter represents its corresponding vector (e.g. <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M1">View MathML</a>). 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. <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M2">View MathML</a> (<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M3">View MathML</a>) is the element-wise magnitude of a vector x (matrix A). IN is an N × N identity matrix. Ai,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, Np the number of pilot subcarriers, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M4">View MathML</a> the channel impulse response (CIR) of lth path at kth symbol, referred to as tap, gk the CIR vector at kth symbol time, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M5">View MathML</a>,hk,n the channel frequency response (CFR) at kth symbol time and nth subcarrier, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M6">View MathML</a> the CFR vector at all subcarriers at kth symbol time, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M7">View MathML</a>, hk the CFR vector at Np pilot subcarriers at kth symbol time, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M8">View MathML</a>,ak the time-correlation coefficients of CFR at kth symbol time, xk the vector of transmitted OFDM symbols at pilot subcarriers at kth symbol time and yk the corresponding received OFDM symbol vector of xk.

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.

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

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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M9">View MathML</a>OFDM symbols (the values of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M10">View MathML</a> 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M11">View MathML</a> consecutive subcarriers for one slot (<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M12">View MathML</a> 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.

Table 1. Physical RBs parameters

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 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M18">View MathML</a>, where TCP 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 TCP is longer than the than the channel’s maximum delay, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M19">View MathML</a>.

(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

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M20">View MathML</a>


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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M21">View MathML</a> 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

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M22">View MathML</a>


Hence, the discrete CIR model can be represented by a length-L CIR vector gk

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M23">View MathML</a>


Strictly speaking, gk is only an approximation of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M24">View MathML</a> at kth OFDM symbol (<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M25">View MathML</a>). When the multipath taps do not fall in the discrete sampling grid (i.e., <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M26">View MathML</a>), 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 gk at kth OFDM symbol is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M27">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M28">View MathML</a>, denoting the CFR of nth subcarrier at kth OFDM symbol time is converted from the time-domain CIR via the discrete Fourier transform (DFT)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M29">View MathML</a>


The relationship between the CIR in time-domain and CFR in frequency domain can be described in matrix notation

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M30">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M31">View MathML</a> is the first L columns of the N × N DFT matrices F. And F is denoted by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M32">View MathML</a>


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 second-order 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 first-order AR model for the time-varying channel.

It is easy to verify that the channel coefficients <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M33">View MathML</a> of the time-varying CFR can be modelled by the following dynamic AR process [1,11,12]:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M34">View MathML</a>


where αn represents the time correlation of the channel response between kth and (k + 1)th OFDM symbols at the nth subcarrier. <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M35">View MathML</a> is a mutually independent zero-mean Gaussian complex white noise representing the modelling error.

LTE OFDM reception and channel estimation

In order to estimate the CFR as defined in (4), Np pilot symbols are inserted sparsely among N subcarriers at k-th OFDM symbol duration following the comb pattern shown in Figure 1. Let <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M36">View MathML</a> denote the transmitted pilot vector of Np known pilot symbols at the kth OFDM symbol, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M37">View MathML</a> denotes the vector of the received pilot symbols. After CP removal, the received pilot symbols can be expressed as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M38">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M39">View MathML</a> is an Np × Np diagonal matrix with transmitted pilot symbols xk as its diagonal elements,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M40">View MathML</a>


Here, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M41">View MathML</a> is an additive white complex Gaussian noise with covariance matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M42">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M43">View MathML</a> 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 hk at these Np pilot symbols with as high accuracy as possible. This is an optimization problem described as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M44">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M45">View MathML</a> is an element-wise division with elements yk,n/hk,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, hk, the channel responses at Np pilot subcarriers at kth OFDM symbol are calculated with the common LS estimator or LMMSE estimator; (2) Time-domain interpolation, where the channel responses hk+1 at (k + 1)th OFDM symbol at pilot subcarriers are estimated from hk 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 {hk,hk+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 {hk,hk+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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M46">View MathML</a> in (9), the task for the Kalman interpolator filter can be stated as:

Given the matrix Xk of known transmitted pilot symbols and received signal yk at kth OFDM symbol, to obtain minimum variance estimates of the time-varying multipath CFR hk and interpolate hk to the followed six <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M47">View MathML</a> 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

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M48">View MathML</a>


where hk is the state variable to be estimated, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M49">View MathML</a> is the unknown state transition matrix consisting of the time correlation coefficients αn of channel response. Both vk and wk are mutually independent, zero-mean, Gaussian complex white noises, with covariance <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M50">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M51">View MathML</a>, respectively. It is assumed that vk and wkare independent of the state variable hk. Note that, in this state space model of the CFR, the state transition matrix Ak is unknown and to be estimated together with the state variable hk. Therefore, it is a problem of joint state and parameter estimation. The purpose is to estimate both the channel response hk and channel’s time-correlation matrix Ak from the received pilot symbols yk.

Considering that Ak is a spares matrix in most cases, without loss of generality, we assume Ak has NA<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M52">View MathML</a> unknown entries to be estimated and let a vector ak denote all the NA unknown entries as follows:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M53">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M54">View MathML</a> is an NA × 1 column vector formed by stacking all unknown entries of the matrix Ak in a row-wise order. The time-correlation parameters are now represented by ak which is the vector to be estimated. For the purpose of clarification, Ak is represented by A(ak) explicitly in the following. Assuming a random walk model for the parameter ak, then Equation (12) becomes

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M55">View MathML</a>


where ϵk denotes the process noise of ak and is an independent, zero-mean Gaussian noise with covariance <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M56">View MathML</a>. In order to jointly estimate the state and parameters, a new augmented state zk is defined as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M57">View MathML</a>


and the channel state space model Equation (14) turns into an augmented system

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M58">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M59">View MathML</a> with covariance matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M60">View MathML</a> and f(zk) is the nonlinear state transition function

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M61">View MathML</a>



Since the state transition function f(zn) 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(ak) 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

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M62">View MathML</a>


Applying the standard KF to the model (18) is straightforward. The resulting EKF algorithm for the joint estimation of CFR hk and CFR’s time correlation coefficients ak works in an iterative prediction–correction cycle. The prediction projects forward (in time) the current estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M63">View MathML</a>and error covariance Pk at kth OFDM symbol to obtain the a priori estimates <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M64">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M65">View MathML</a> for the next (k + 1)th OFDM symbol. The correction adjusts the projected estimates <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M66">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M67">View MathML</a> to obtain an improved a posteriori estimate by using an actual measurement of received symbol yk+1 at (k + 1)th OFDM symbol. Here, the subscript <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M68">View MathML</a> corresponds to one-step a priori prediction, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M69">View MathML</a> corresponds to a posteriori correction and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M70">View MathML</a> 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):


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M71">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M72">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M73">View MathML</a>


and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M74">View MathML</a> is the covariance of noises <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M75">View MathML</a>.

2. Correction (once the reception of the OFDM symbol has completed):

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M76">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M77">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M78">View MathML</a>


Here, Kk 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.


In terms of computation complexity, it can be seen that prediction of state error covariance <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M79">View MathML</a> and the update of Kk 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 Ak 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 Ak can be modelled as a diagonal matrix. If both Xk and Ak 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 25Np.

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 Xk is not always available at each iteration. Like most adaptive algorithms, two working modes, namely, training mode and decision-directed interpolation mode, are adopted in the proposed Kalman interpolation filter to address this issue.

thumbnailFigure 2. Channel estimation and interpolation at pilot subcarriers.

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 Xk are known to the estimator, while in the data symbols periods, the transmitted data symbols are estimated as <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M80">View MathML</a>by the decoder and the EKF is fed by the <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M81">View MathML</a>to replace the unknown transmitted symbols Xk. 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M82">View MathML</a>of the transmitted data symbols on the basis of the previous, a priori channel estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M83">View MathML</a>. In the iteration of the OFDM data symbol, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M84">View MathML</a>is also fed to the EKF to calculate a posteriori channel estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M85">View MathML</a> and a priori channel estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M86">View MathML</a>. 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 zk, 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 zk consists of two independent components, ak and hk, their initial values are chosen separately.

For initializing the channel’s time-correlation coefficients ak, we use an identity matrix (<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M87">View MathML</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 hk, 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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M88">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M89">View MathML</a> 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 Xk. As the observation yk is obtained by the DFT at the end of an OFDM symbol duration, the a posteriori CFR <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M90">View MathML</a> is first estimated from yk by using update equations (22)–(24). Then the a priori estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M91">View MathML</a>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 Xk is unknown, Xk is replaced by the decoder’s decision of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M92">View MathML</a>that is supposed to be nearest to Xk. In the decision-directed mode, the prediction and correction processes are the same as the training mode, except Xk is replaced by <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M93">View MathML</a>,.

It is worth noting that, as yk 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M94">View MathML</a>to refine the currently received OFDM symbol, rather than uses the a posteriori CFR <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M95">View MathML</a>.

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 yk is subject to the noise wk, 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M96">View MathML</a>can be calculated by <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M97">View MathML</a>, where Ptx 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M98">View MathML</a>and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M99">View MathML</a>in <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M100">View MathML</a>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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M101">View MathML</a>. At higher SNRs, a better channel estimation is expected and a smaller value is used for <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M102">View MathML</a>.


We now summarize the proposed method for channel estimation in LTE downlink:

Step 1. Initialize <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M103">View MathML</a> when the first pilot symbol is received, make the first a priori prediction <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M104">View MathML</a>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 yk by using DFT

Step 3. Estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M105">View MathML</a> by equalizing yk with previous a posteriori<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M106">View MathML</a>;

Step 4. If yk is pilot symbol, set Xk by the known pilot symbol xk,else set Xk by the estimated data symbol <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M107">View MathML</a>,

Step 5. Correct a posteriori state estimation <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M108">View MathML</a> from yk by (22)–(24).

Step 6. Time-interpolation: Predict a priori state estimate <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M109">View MathML</a>by (19)–(21) for next symbol.

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.

Table 2. The simplified rural area channel model

Table 3. LTE downlink simulation parameters

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.

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

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 time-varying 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 its original QAM constellation position but nearer to another constellation position, the quantization procedure will result in a wrong decision of the data symbol. When the incorrect data symbol decision is fed back to the EKF, it works as an incorrect ‘observation’ resulting in the EKF giving an abrupt change in state estimation. As a result, a sudden jump appears in the CFR estimates and may result in error propagation, making more errors in the following data symbol decision. If these decision errors are infrequent enough, the effects of these errors decay away and the decision-directed equalizer’s performance remains similar to that of the training mode.

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

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.

thumbnailFigure 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.

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, 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.

thumbnailFigure 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.

thumbnailFigure 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.

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

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 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.


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.


Applying the first-order Taylor approximation to the nonlinear state transition function f(zn) around <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M110">View MathML</a>in Equation (17), the state equation (16) becomes


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M111">View MathML</a>

where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M112">View MathML</a> is the linearization error <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M113','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M113">View MathML</a>and

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M114','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M114">View MathML</a>


Here, we assume that A(ak) is independent of hn. Recalling the definitions of ak and A(ak), it is easy to verify that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M115','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M115">View MathML</a>is an Np × NA block-diagonal matrix Hn

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M116','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M116">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M117">View MathML</a> denotes the Kronecker product and the operator <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M118','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M118">View MathML</a> removes these known-columns of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M119','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M119">View MathML</a>. The known column is the <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M120','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M120">View MathML</a>th column when the ith-row-jth-column entry <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M121','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M121">View MathML</a> of A is known. Hence, substituting <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M122">View MathML</a> into Equation (26), we have

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M123','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M123">View MathML</a>


And the linear state space model approximating the AR model (12) is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/232/mathml/M124">View MathML</a>


3GPP: The 3rd Generation Partnership Project (3GPP); AR: autoregressive; CIR: channel impulse response; DFT: discrete Fourier transform; EKF: extended Kalman filter; LMMSE: linear minimum mean square error; LS: least square; LTE: long-term evolution; OFDM: orthogonal frequency-division multiplexing; PSAM: pilot symbol-aided modulation; QPSK: quadrature phase-shift keying; SNR: signal-to-noise ratio.

Competing interests

The authors declare that they have no competing interests.


This study was supported by the EPSRC UK-China Science Bridges: R&D on 4 G Wireless Mobile Communications under grant EP/G042713/1.


  1. B Karakaya, An adaptive channel interpolator based on Kalman filter for LTE uplink in high Doppler spread environments. EURASIP J. Wirel. Commun. Netw. 2009, 1–10 (2009)

  2. JK Cavers, An analysis of pilot symbol assisted modulation for Rayleigh fading channels. IEEE Trans. Veh. Technol. 40(4), 686–693 (1991). Publisher Full Text OpenURL

  3. F Tufvesson, T Maseng, Pilot assisted channel estimation for OFDM in mobile cellular systems. in Proc. of IEEE Vehicular Technology Conference’97 3, 1639–1643 (1997) Phoenix, Az OpenURL

  4. MH Hsieh, CH Wei, Channel estimation for OFDM systems based on comb-type pilot arrangement in frequency selective fading channels. IEEE Trans. Consum. Electron. 44(1), 217–225 (1998). Publisher Full Text OpenURL

  5. J Beek, O Edfors, M Sandell, S Wilson, P Borjesson, On channel estimation in OFDM systems. in Proc. of IEEE Vehicular Technology Conference’95 2, 815–819 (1995) Chicago, Il OpenURL

  6. P Hoecher, S Kaiser, P Robertson, Pilot-symbol-aided channel estimation in time and frequency. in Proc. of IEEE Global Telecommunications Conference’97 Communication Theory Mini-Conference 3, 90–96 (1997) Phoenix, Az OpenURL

  7. WH Chin, DB Ward, AG Constantinides, An algorithm for exploiting channel time selectivity in pilot-aided MIMO systems. IET Commun. 1(6), 1267–1273 (2007). Publisher Full Text OpenURL

  8. O Edfors, J van de Beek, M Sandell, SK Wilson, PO Borjesson, Analysis of DFT-based channel estimators for OFDM. Int. J. Wirel. Personal Commun. 12(1), 55–70 (2000). Publisher Full Text OpenURL

  9. Steepest Ascent Ltd Improving throughput performance in LTE by channel estimation noise averaging, The LTE-Advanced Guide (Online, http://www, 2010), . steepestascent.com/content/mediaassets/pdf/products/LTE_Portal_Article_May_2010.pdf webcite OpenURL

  10. J Hou, J Liu, A novel channel estimation algorithm for 3GPP LTE downlink system using joint time-frequency two-dimensional iterative Wiener filter. in Proc. of 12th IEEE Int. Conf. on Communication Technology (ICCT) 1, 289–292 (2010) Nanjing, China OpenURL

  11. R Iltis, Joint estimation of PN code delay and multipath using the extended Kalman filter. IEEE Trans. Commun. 88(10), 1677–1683 (1990)

  12. K-Y Han, S-W Lee, J-S Lim, K-M Sung, Channel estimation for OFDM with fast fading channels by modified Kalman filter. IEEE Trans. Consum. Electron. 50(2), 443–449 (2004). Publisher Full Text OpenURL

  13. M Huang, X Chen, L Xiao, S Zhou, J Wang, Kalman-filter-based channel estimation for orthogonal frequency-division multiplexing systems in time-varying channels. IET Commun. 1(1), 759–801 (2007)

  14. SK Mitra, Digital Signal Processing: A Computer-Based Approach, 2nd edn. (McGraw-Hill/Irwin, Boston, MA, 2001)

  15. LM Davis, I Collings, R Evans, Coupled estimators for equalization of fast-fading mobile channels. IEEE Trans. Commun. 46(10), 1262–1265 (1998). Publisher Full Text OpenURL

  16. W Li, Estimation and tracking of rapidly time-varying broadband acoustic communication channels, Ph.D. dissertation, Massachusetts Institute of Technology & Woods Hole Oceanographic Institution

  17. 3GPP, Technical specification group radio access network Deployment aspects (release 10), 3GPP TR 25.943V10.0.0,” 3GPP, Technical Report (Online, http://www, April 2011), . 3gpp.org/ftp/specs/html-info/25943.htm webcite OpenURL