Abstract

We address the problem of estimating time-varying fading channels in filter bank multicarrier (FBMC/OQAM) wireless systems based on pilot symbols. The standard solution to this problem is the least square (LS) estimator or the minimum mean square error (MMSE) estimator with possible adaptive implementation using recursive least square (RLS) algorithm or least mean square (LMS) algorithm. However, these adaptive filters cannot well-exploit fading channel statistics. To take advantage of fading channel statistics, the time evolution of the fading channel is modeled by an autoregressive process and tracked by Kalman filter. Nevertheless, this requires the autoregressive parameters which are usually unknown. Thus, we propose to jointly estimate the FBMC/OQAM fading channels and their autoregressive parameters based on dual optimal Kalman filters. Once the fading channel coefficients at pilot symbol positions are estimated by the proposed method, the fading channel coefficients at data symbol positions are then estimated by using some interpolation methods such as linear, spline, or low-pass interpolation. The comparative simulation study we carried out with existing techniques confirms the effectiveness of the proposed method.

1. Introduction

Wireless multicarrier (MC) communication systems are parallel data transmission techniques in which high data rates can be achieved by the simultaneous transmission over many orthogonal subcarriers [1]. Using MC communications, a wide-band frequency-selective fading channel is divided into many narrow-band frequency nonselective flat fading subchannels, facilitating channel estimation and equalization as well as time synchronization and narrow-band interference mitigation. In addition, the division of the whole bandwidth into many subchannels provides scalability and flexibility when configuring the communication link.

The most widely used multicarrier modulation technique is orthogonal frequency division multiplexing (OFDM) with cyclic prefix (CP) [2]. Due to the various advantages of OFDM, it has become the key physical layer transmission technology adopted in current broadband communication systems such as wireless local area networks (WLAN), worldwide interoperability for microwave access (WiMAX), and long term evolution (LTE) as well as in digital video and audio broadcasting (DVB and DAB) [2, 3]. The OFDM system offers simple implementation using the inverse fast Fourier transforms (IFFT) and the fast Fourier transform (FFT) pairs in the modulator and demodulator, respectively. However, the large side lobes of the frequency response of the filters that characterize the subcarrier channels result in significant interference among subcarriers. As an alternative to OFDM, filter bank multicarrier (FBMC) modulation, which is also implemented based on IFFT/FFT pairs, has several advantages over OFDM [4, 5]. Firstly, it does not require CP resulting in better use of the allocated spectrum. Secondly, instead of using a rectangular window, a prototype filter designed with the Nyquist pulse shaping principle is adopted, which can reduce greatly the spectral leakage problem of OFDM. This results in negligible intersymbol interference (ISI) and intercarrier interference (ICI). Thirdly, the combination of filter banks with offset quadrature amplitude modulation (FBMC/OQAM), known also as OFDM/OQAM, leads to a maximum data transmission rate [6].

Fading channel estimation and equalization techniques are crucial for coherent symbol detection at the receiver. In the framework of multicarrier systems, many channel estimation techniques were proposed in the literature particularly for OFDM systems [711]. They can be classified as training sequence/pilot symbols based techniques and blind methods. In this paper, as blind techniques require a large amount of data and have higher complexity than training/pilot based techniques, we will focus our attention on the latter techniques.

Compared to OFDM, few FBMC/OQAM studies have addressed the problem of channel estimation and equalization [1218]. Indeed, when classical channel estimation methods used for OFDM systems are applied directly to FBMC/OQAM system, an intrinsic intersymbol-interference is observed. This results in system performance degradation. To reduce this intrinsic interference, several methods are recently proposed, including both preamble-based [1214] and scattered pilots-based [15, 16] channel estimators. In [17], a two-stage minimum mean square error (MMSE) equalizer is proposed to cope more efficiently with intercarrier interference (ICI). A fractionally spaced adaptive decision feedback equalizer (DFE) based on the least mean square (LMS) algorithm was proposed in [18] where the input of each equalizer comprises only the output of each subcarrier.

When dealing with Kalman filter based channel estimators, several approaches were developed for multicarrier systems such as multicarrier code division multiple access (MC-CDMA) and OFDM [10, 11, 19]. Indeed, the Kalman filter combined with an autoregressive (AR) model to describe the time evolution of the fading channel is shown to track channel variations and to provide superior bit error rate (BER) performance over the standard LMS and recursive least square (RLS) estimators [19]. Nevertheless, the AR parameters are unknown and should be estimated. To jointly estimate the autoregressive process and its parameters from noisy observations, one of the authors of this paper has recently proposed dual Kalman filters based structure [11, 20]. This structure consists of two cross-coupled Kalman filters where the first filter uses the latest estimated AR parameters to estimate the fading process, while the second filter uses the estimated fading process to update the AR parameters. In this paper, we propose to investigate the relevance of the dual Kalman filters for the estimation of FBMC/OQAM time-varying fading channels based on pilot symbols. Once the fading channel coefficients at pilot symbol positions are estimated by the proposed method, the fading channel coefficients at data symbol positions are then to be estimated by using some interpolation methods such as linear, spline, or low-pass interpolation.

The remainder of the paper is organized as follows. In Section 2, we explain the FBMC/OQAM system model. The estimation of the fading channels based on dual Kalman filters is introduced in Section 3. Simulation results are reported in Section 4. Conclusion remarks are drawn in Section 5.

2. FBMC/OQAM System Model

There are different filter bank multicarrier structures that are proposed in the literature. Among them are discrete wavelet multitone [21], filtered multitone (FMT) [22], cosine modulated multitone (CMT) [23], modified discrete Fourier transform (MDFT) [24], exponentially modulated filter bank (EMFB) [25], and FBMC/OQAM [6]. The last is the most popular filter bank scheme that results in high bandwidth efficiency.

The fundamental parts of the FBMC/OQAM system are the synthesis filter bank (SFB) at the transmitter and the analysis filter bank (AFB) at the receiver, arranged in the so-called transmultiplexer (TMUX) configuration as shown in Figures 1 and 2. It should be noted that the SFB and AFB pairs can be efficiently implemented using FFT and IFFT of size aided by polyphase filtering structures. When using uniformly modulated filter banks [26], a prototype filter of length is shifted in frequency to generate subchannels which cover the whole system bandwidth. Thus, the discrete-time baseband signal at the output of the FBMC transmitter with OQAM modulation can be expressed as where is the sample index at SFB output/AFB input, is the sample index at OQAM preprocessing output/postprocessing input, is the number of subcarriers in the filter bank, and

In (3), are real-valued data symbols at subcarrier , transmitted at a rate of where the signaling period is defined as with being the subcarrier spacing. The pair of symbols and represent, respectively, the in-phase and quadrature parts of the complex-valued input symbol . This complex-to-real conversion operation is considered to introduce an up-sampling by two. The inverse operation of real-to-complex conversion is performed at the receiver and results in down-sampling the signal by two. Thus, the complex-to-real operation performs the following mapping: where is the sample index at OQAM preprocessing input/postprocessing output. Note that the complex-valued input symbols belong to QAM alphabet with the real and imaginary parts being interleaved by a time offset of , resulting in OQAM modulation. To maintain orthogonality, is multiplied by .

In order to achieve high spectral efficiency, complex modulated filter banks are usually used, which means that all subchannel filters are frequency shifted versions of the prototype filter . So, the synthesis filter can be expressed as where . The length of the prototype filter depends on the size of the filter bank and the number of FBMC/OQAM symbol waveforms that overlap in time as .

The analysis filter is simply a time-reversed and complex-conjugated version of the corresponding synthesis filter:

It should be noted that the design of the prototype filter must satisfy perfect reconstruction (PR) conditions or at least provide nearly perfect reconstruction (NPR) characteristics. However, the PR property is only achieved under the condition of ideal transmission channel. As interferences in the wireless channel are unavoidable, there is no way to have PR conditions. Thus, prototype filters are designed to satisfy NPR characteristics.

In this paper, we used frequency sampling technique [27] or windowing based approach [28] to design NPR prototype filter. In these methods, the prototype filter coefficients can be obtained using a closed-form representation that includes only a few adjustable design parameters. The impulse response coefficients of the filter are obtained according to the desired frequency response, which is sampled on a uniformly spaced frequency point . So, the finite impulse response (FIR) of the low-pass prototype filter can be expressed as [29] where , , and for , for, .

When polyphase filtering structures are used with the FFT and IFFT pairs to implement the FBMC/OQAM system, the synthesis polyphase filter at the transmitter is defined as At the receiver, the analysis filter polyphase is given by

The transmitted FBMC/OQAM signal in (1) is assumed to go through a time-varying frequency-selective Rayleigh fading channel with background additive white Gaussian noise (AWGN). Assuming that the allocated bandwidth for each subcarrier is less than the coherence bandwidth of the wireless channel, each subcarrier undergoes frequency nonselective flat-fading. Thus, the received signal sample over the th subcarrier for the FBMC/OQAM symbol can be expressed as where is a complex valued fading process over the th subcarrier for the th FBMC/OQAM symbol and is an additive white Gaussian noise (AWGN) process. The noise processes are assumed to be mutually independent and identically distributed zero-mean complex Gaussian processes, with equal variances . The fading process is assumed to be a zero-mean complex Gaussian process with uniformly distributed phase on and with Rayleigh distributed envelope . The variances of the fading processes are all assumed to be equal to .

After processing the received signal with the analysis filter bank block, the resulting signal at the input of channel estimation and equalization is given by where with denoting the convolution operator and being a Gaussian noise process with variance . Under the assumption of nearly perfect reconstruction (NPR) conditions of the prototype filter, it follows that . Thus, (12) is reduced to

The statistical properties of the fading process are given by its power spectrum density (PSD) and autocorrelation function (ACF). Indeed, the PSD of is defined by the well-known U-shaped band-limited Jakes spectrum with maximum Doppler frequency [30]: where with is the mobile speed and is the wave length of the carrier wave. The corresponding normalized discrete-time autocorrelation function (ACF) hence satisfies where is the zero-order Bessel function of the first kind, is the symbol period, and denotes the Doppler rate.

To exploit the statistical properties of the fading channel given by its PSD and ACF, the fading process can be modeled by a th order AR process, denoted by and defined as follows [31]: where are the AR model parameters and denotes the zero-mean complex white Gaussian driving process with equal variance over all subcarriers. Increasing the AR mode order will provide better fitting between the resulting model and the theoretical PSD and ACF [11, 31]. However, the computational complexity of the resulting channel estimation algorithm will also increase. Therefore, a compromise between the accuracy of the model and the computational complexity of the estimation algorithm has to be found. When the Doppler rate is available at the receiver, the AR parameters can be computed by solving the so-called Yule-Walker (YW) equations. However, as the Doppler rate is usually unknown, we propose to complete the joint estimation of the fading process and its AR parameters based on dual Kalman filters.

3. Dual Kalman Filters Based Channel Estimation

The estimation of the fading process along the th FBMC/OQAM symbol will be performed in two steps. Firstly, the fading process at the pilot symbol position is estimated using dual Kalman filters. Secondly, the fading process at data symbol position will then be estimated by using some interpolation methods such as linear, spline, or low-pass interpolation.

When using comb-type pilot arrangement [8] as shown in Figure 3, pilots are uniformly inserted into as follows: where is the pilot symbol with , and is the pilot interval spacing with being the number of subcarriers.

Thus, using known comb-type pilot symbols , we propose to jointly estimate the fading process and its AR parameters based on dual Kalman filters as shown in Figure 4. Indeed, the first Kalman filter in Figure 4 uses the pilot symbol , the output of the analysis filter bank , and the latest estimated AR parameters to estimate the fading process , while the second Kalman filter uses the estimated fading process to update the AR parameters.

3.1. Fading Process Estimation at Pilot Symbol Positions

In this subsection, our purpose is to estimate the fading process over the FBMC/OQAM symbol based on Kalman filtering with known pilot symbols. To this end, let us define the state vector as follows:

Note that, for the sake of simplicity and clarity of presentation, the time index subscript is dropped. Then, (16) can be written in the following state-space form: where In addition, given (13) and (18), it follows that where .

Hence, (19) and (21) represent the state-space model dedicated to the FBMC/OQAM symbol fading channel system (13) and (16). A standard Kalman filtering algorithm can then be carried out to provide the estimation of the state vector given the set of observations . To this end, let us introduce the so-called innovation process which can be obtained as follows: The variance of the innovation process can be expressed as where the so-called a priori error covariance matrix can be recursively obtained as follows: The Kalman gain is calculated in the following manner: The a posteriori estimate of the state vector and the fading process are, respectively, given by The error covariance matrix is updated as follows:

It should be noted that the state vector and the error covariance matrix are initially assigned to zero vector and identity matrix, respectively; that is, and .

Equations (22)–(27) can be evaluated providing the availability of the AR parameters that are involved in the transition matrix and the driving process variance . They will be estimated in the next subsections.

3.2. Estimation of the AR Parameters

To estimate the AR parameters from the estimated fading process , (26) are firstly combined such as the estimated fading process is a function of the AR parameters: where and . In addition, the variance of the process is given by

When the channel is assumed stationary, the AR parameters are invariant and satisfy the following relationship:

As (28) and (30) define a state-space representation for the estimation of the AR parameters, a second Kalman filter can be used to recursively estimate as follows: where the Kalman gain and the update of the error covariance matrix are, respectively, given by with initial conditions and .

3.3. Estimation of the Driving Process Variance

To estimate the driving process variance , the Riccati equation is first obtained by inserting (24) in (27) as follows: Taking into account that is a symmetric Hermitian matrix, one can rewrite (25) in the following manner: By combining (33) and (34), can be expressed as follows: where is the pseudoinverse of .

Thus, we propose to estimate recursively as follows: where the variance of the innovation process is replaced by its instantaneous value and is the forgetting factor. It should be noted that can be either constant or time-varying (e.g., ).

3.4. Fading Process Estimation at Data Symbol Positions

Once the fading process at pilot symbol positions is estimated by the dual Kalman filters, the fading process at data symbol positions will then be estimated by using some interpolation methods such as linear, spline, or low-pass interpolation [8].

(1) Linear Interpolation. Using linear interpolation method, the channel estimates at data positions, , are given by where is the estimated fading process at pilot symbol position.

(2) Low-Pass Interpolation. The low-pass interpolation is carried out by inserting zeros into the data symbol positions of the original sequence and then applying a special low-pass filter. In this paper, we use the Matlab function interp.

(3) Spline Interpolation. Spline interpolation produces a smooth and continuous polynomial fitted to the estimated fading process at pilot symbol positions. In this paper, we use the Matlab function spline.

3.5. Fading Channel Equalization

Once the fading process at pilot and data symbols is estimated using the proposed approach, channel equalization can be performed by multiplying (13) with a normalized version of the complex conjugate of the channel estimate as follows:

4. Simulation Results

In this section, we carry out a comparative simulation study on the estimation of FBMC/OQAM fading channels between the proposed dual Kalman filters based channel estimator and the conventional LMS and RLS channel estimators. In addition, we test the performance of three interpolation techniques, namely, low-pass, spline, and linear interpolation. Furthermore, we investigate the effect of the number of pilot symbols and Doppler rate on the BER performance.

In all of our simulations, the fading channels are generated according to the autoregressive based method presented in [31]. The autoregressive model order in the proposed dual Kalman filters based estimator is set to . The step size for LMS algorithm is set to and the forgetting factor for RLS algorithm is set to .

Figure 5 shows the envelope of estimated fading process using the various channel estimators with low-pass interpolation, SNR = 15 dB, , , and. One can notice that the dual Kalman filters based estimator provides much better estimation than the LMS and RLS based ones.

Figure 6 shows the BER performance versus SNR for the FBMC/OQAM system when using the various channel estimators with , pilots, , and using low-pass interpolation. According to Figure 6, the proposed dual Kalman filters based estimator outperforms the conventional LMS and RLS estimators with the performance difference increasing as the SNR increases. Figure 7 displays the BER performance of the FBMC/OQAM system versus Doppler rate when using the various channel estimators with low-pass interpolation, , and pilots. One can notice that the dual Kalman filters based estimator yields lower BER performance than the LMS and RLS estimators with the BER performance difference increasing as the Doppler rate increases. This confirms that the dual Kalman filters based estimator can exploit the fading channel statistics and can track fast variations of rapidly varying fading channels.

The effect of changing the number of pilot symbols on the BER performance of the FBMC/OQAM system when using dual Kalman filters based channel estimator with low-pass interpolation is shown in Figure 8. Indeed, increasing the number of pilot symbols will improve the BER performance of the system.

Figure 9 shows the BER performance of the FBMC/OQAM system when using dual Kalman filters with the various interpolation methods, , pilots, and . It is confirmed that the low-pass interpolation yields the best BER performance.

5. Conclusions

This paper addresses the estimation and equalization of FBMC/OQAM fading channels based on pilot symbols. A dual Kalman filters based structure is proposed for the joint estimation of the fading process and its AR parameters over each subcarrier at pilot symbol positions. The fading process at data symbol positions is then obtained by using some interpolation methods such as linear, spline, or low-pass interpolation. The simulation results showed that the proposed dual Kalman filters based estimator outperforms the LMS and RLS based ones. In addition, the low-pass interpolation is confirmed to outperform both spline and linear interpolation.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.