Assisted Channel Estimation with Dynamic Allocation Pilot Subcarriers for OFDM Communications

In this paper, a transmission system based on OFDM (Orthogonal Frequency Division Multiplexing) technique and channel estimation is proposed. The idea is to enhance the assisted channel estimation without data rate loss. A solution based on frequency diversity is proposed. With the same number of pilots, the global system performances are better for higher mobility speed. The main results will be presented in the case of high mobility context.


Introduction
The channel estimation is a fundamental step in a transmission system whose performances are directly related to the severity of this step in the reception process.The estimation quality is highly dependent on variations of the channel propagation.In the literature, one can find three main types of channel estimation techniques: supervised or assisted techniques, blind techniques and semiblind techniques with decision feedback [1].The blind estimation technique does not require insertion of pilot subcarriers in the frame.It uses some statistical properties of the received signal (i.e.redundancy due to the insertion of cyclic prefix, higher order statistics) and directly addresses the problem of loss of spectral efficiency [2][3][4][5].However, this technique often requires significant computation time and isn't appropriate for this case of high mobility.The estimation technique with a semiblind decision feedback using both known data (i.e.pilot subcarriers) and estimated data sent with a forward link.This estimation is very sensitive to the detection of errors on the estimated symbols [1,4], and in particular in the context of the high mobility with rapid channel variations.The assisted channel estimation is based on inserting, in the OFDM frame, symbols which are known by the receiver and called pilot symbols.The associated subcarriers are called the subcarrier pilots.Training sequence corresponds to the set of pilot subcarriers dedicated to the channel estimation.It consists in a pattern robust to the frequency selectivity of the channel.It is suitable for quasi-stationary channels for which the channel coefficients vary little from one OFDM symbol to another.The insertion of pilot "comb type" consists in reserving some pilot subcarriers only to channel estimation for all transmission duration.This pattern is robust to the channel temporal variations from an OFDM symbol to another [6].The main disadvantage of the channel assisted estimating technique is the spectral efficiency loss due to the insertion of pilots between the transmitted data.However, we chose this technique despite the loss of spectral efficiency.In the broadband challenge for high mobility, the hostile characteristics nature of mobile propagation channel requires the use of pilots.The rapid variations of the channel require a minimum of channel knowledge by pilots in order to get a correct detection datum.
In this paper, we will focus on two aspects which influence strongly the channel estimation: pilot's geometry used and the algorithms applied to estimate the transfer function gains associated to data subcarriers.In the following section, we rely on the work of [7,8], using a model of BEM (Basis Expansion Modeling) channel development.Then, we propose a pilot subcarriers configuration to improve the estimation for rapid changes in time channel.Our objective is to enhance the spectral efficiency with the same number of pilots in order to increase the data rate.
The paper is organized as follows.Section II formulates the data model for OFDM system.Section III pre-sents the chosen propagation channel model and its Basis Expansion Modeling (BEM) approach.In Section IV, this approach is applied to channel estimation with a dynamic allocation of pilot subcarriers.Simulation results are shown in Section V. Finally, we conclude the presented study in section VI.
Notation: x is a real scalar, is a vector, is a matrice.

System Model
The multipath effect introduced by the propagation channel causes interferences between symbols when using a transmission type series.A good solution to overcome the effects of the channel is to use a parallel or multicarrier transmission.The OFDM technique consists to convert the serial coded data stream on parallel blocks and to transmit simultaneously these blocks over several orthogonal carriers using the Inverse Fast Fourier Transform (IFFT).In the transmitter side, the data stream is grouped into blocks of N data symbols, called OFDM symbol, and are represented by the vectors   i X k ( : are frequency and OFDM symbol index).Then, an IFFT is performed on each data symbols block: where   i X k represents the complex data (M-PSK or M-QAM) which modulates the frequency k f , and . b is the symbol OFDM duration.By sampling the continuous signal at time sampling e , we obtain the following equation: The received signal is generally the sum of a convolution of the channel impulse response It is implicitly assumed that the transmitter and the receiver are perfectly synchronized: where is the symbol of the circular convolution operator and    i V n are samples of additive white gaussian noise (AWGN).At the receiver side, the cyclic prefix is removed and the data symbols are obtained by performing the FFT operation: and H represent the channel matrices respectively in time and frequency domain for each OFDM symbol.The channel matrix H is not diagonal because of the channel temporal variations during an OFDM symbol (i.e.

  k
, k z H , transmitted and received carriers).It is expressed using the following equation [9,10]: where represents the Fourier transform of the impulse response with is the total number of paths in the multipath fading channel.

Propagation Channel Model
Since many years, different approaches for mobile radio channels modeling and simulation have been developed in the literature.But, no specific channel model is proposed for the mobile communications taking into account very high velocities such those encountered for very high speed train.The model of Clarke [11], which is the mathematical model of reference, and the simplified version of Jakes [12] are the most widely used to model variations in the signal channels called Rayleigh fading.
Simulation of a fading channel with Doppler calls for the generation of a Gaussian random process in order to mimic the sampled version of the channel waveform.These simulations are widely used in the process of designing efficient and robust next generation wireless communication and broadcasting systems.One important aspect to be considered is the accuracy of such simulation methods.Statistically, the provided discrete channel response samples should meet certain requirements related to the spectrum and time-domain variations.Hence, some quality measures have been introduced such as the probability density function (PDF), the cumulative distribution function (CDF), the autocovariance function (ACF), the level crossing rate (LCR) and the average duration of fades (ADF).Generated samples should meet the requirements in terms of quality measures during the predefined simulation period.The task of designing a channel simulator should tackle these issues while keeping a reasonable computation complexity and memory requirements.
State-of-the-art Rayleigh fading simulators can be divided into three main categories: the sum of sinusoids (SOS) method, filtering white Gaussian variables method and the inverse fast Fourier transform (IFFT) method.In [13], Pop and Beaulieu showed that the Jakes simulator is not wide-sense stationary in some cases and proposed other modifications.However, although the number of sinusoids is considered close to infinity, the higher order statistics of the simulator do not follow the theoretical statistical model of Clarke.This was also shown in [14,15].In [16], Xiao and Zheng have proposed a new model for a sum of sinusoids Rayleigh fading channel that generates a non-deterministic signal with statistical properties converging to those of Jakes model even with a small number of sinusoids s N .The following equations give the fading process of this simulator:       where  and  are random variables statistically independent and uniformly distributed for any number of sinusoids s N .n  , n  represent respectively the arrival angle in the direction of movement of the receiver and the initial phase.This model is selected for the proposed system simulations presented in this paper.In the next section, we describe the approximation of the complex channel gains by a development in basis functions.The goal is to model the temporal channel variations during an OFDM symbol and to determine their evolutions during to the transmission duration.

Basis Expansion Modeling
The idea of Basis Expansion Modeling (BEM) is described in [8,17,[18][19][20][21][22].The use of such a development is particularly useful to reduce the number of channel parameters to be estimated on a symbol time.For a block transmission during symbols, the channel is characterized by samples, is the total number of broadband channel paths according to the representation given by the relation (5) in which each path has a different realization at each time.The BEM model represents the total samples of the channel complex gains where is the matrix approximation of the basis functions and is the order the BEM development.
contains the Q coefficients corresponding to the basic samples N   l h n and l  is the modeling error that's often neglected and dependent on the number .For the basis function matrix , there are several approaches: the Discrete Spheroidal Sequences (DSS-BEM), the Karhunen-Loeve (KL-BEM), the Complex Exponential (CE-BEM) or the Polynomial (P-BEM) [8,18,[23][24][25].Subsequently, we will focus on the P-BEM.

Q B
Considering the modeling error l  equal to zero (i.e.

   
), the correlation matrix of the basis coefficients is defined by: .
Regarding the channel estimation that is modeled by a BEM development, we only need to estimate the basis coefficients.We chose to model the channel path by a P-BEM [26].In this case, has the following expression where n

Proposed Channel Estimation Technique Channel Estimation with Dynamic Allocation of Pilots
In this section we consider the estimation of a timevarying channel with pilot subcarriers whose positions change from one OFDM symbol to another.We perform time channel estimation instead of the frequency estimation.Indeed, the last estimation based on conventional methods such as ML (Maximum Likelihood, MMSE (Minimum Mean Square Error) or LS (Least Square), gives poor performances because the estimated channel matrix is no longer diagonal in this context [8,9,26].
We propose to vary the positions of the pilot carriers in the OFDM symbol, in order to scan the whole spectrum of the OFDM symbol (see Figure 1).The aim is, firstly, to avoid the estimation of the channel on the same pilots positions, which could be lost by fading, and secondly, to obtain an adaptive estimation of the channel.This is the proposed technique called Dynamic Allocation of Pilots (DAP).
Figure 2 illustrates the operating mechanism of the DAP.We note here that before applying the IFFT signal   i X k , it consists of complex payload (QPSK or 16QAM) and pilots  BEM basic functions for the channel model in order to reduce system complexity while maintaining an acceptable level of performance.Equation ( 14) can be rewritten as [27,28].So the symbol can be written as:

 
or and Assuming that is the number of pilot subcarriers is the space between two successive pilot subcarriers.It is then possible to reformulate as the sum of two matrices such as: where corresponds to the pilots and in Equation ( 14), we obtain: Now, we will use the pilot positions to estimate the channel impulse response .Let h corresponding to the positions of the pilots.We can build linear equations using the following equation: To find the impulse response as a least squares solution of the system, we'll replace it with its basis coefficients .Omitting the indices, the Equation ( 18) becomes [27,28]: hence the expression of the estimated channel h becomes:

Simulation Results
In the following, we consider an OFDM system.The total number of subcarriers is 128 in which subcarriers are deployed for data transmission and subcarriers are used like pilots.We set to zero the middle DC carrier and the lateral subcarriers placed on left and on right of the OFDM symbol spectrum.A Cyclic Prefix (CP) with samples is used and the OFDM symbol duration is 72 μs.
For the channel model simulation, we have chosen the Zheng model with a Jakes spectrum for mobile radio channel with respect to the motion between the transmitter and the receiver.Assuming that the synchronization between the transmitter and the receiver is perfect at receiver, the Bit Error Rate (BER) is evaluated for three normalized Doppler frequencies F 1 = 0.0003, F 2 = 0.0006 and F 3 = 0.0008 corresponding to a mobile speed v of 140, 280 and 350 kmph at a carrier frequency f 0 = 2.5 GHz.   which does not exceed 0.25 dB.This can be due to the size of the constellation of 16 QAM.

Conclusion
In this paper we proposed a solution to minimize the Doppler Effect damage and consequence on the performances of OFDM communication systems.OFDM technology provides efficient estimation of time-variant channels.In order to estimate the channel, pilot subcarriers are inserted into the OFDM symbol.The channel is estimated in the time domain and the temporal variations of each path are modeled by BEM.A pilot structure based on varying progressively the positions of pilot subcarriers from a symbol to another, is proposed.The pilot subcarrier positions are allocated dynamically during the transmission of the OFDM symbols with the DAP.The estimated time-varying channel is better and therefore the system performance will be improved in presence of the Doppler Effect.Using the DAP, the results show that the time-varying channel estimation is improved and the severe consequences of the Doppler Effect are reduced.

Figure 1 .
Figure 1.OFDM frame structure with pilot carrier offset.

Figure 2 .
Figure 2. Block diagram of the DAP.

Figures 3 and 4
Figures 3 and 4 illustrate the evolution of BER with increasing the value of the normalized Doppler frequency.Concerning the QPSK modulation, in a channel where the temporal variations are quantified by F 1 and F 2 and energy per bit 25 b o E dB N  , The DAP improvement is about 1 or 2 dB.The DAP proposed method tracks channel variations better than the conventional case with fixed pilots.On the other hand, for the value F 3 (350 kmpl), the DAP offers a lower gain of 0.5 dB.This is explained by the rapid variations in the channel and changing the pilot positions doesn't achieve the same performance as the values for F 1 and F 2 .Concerning the 16 QAM and for three values of the normalized Doppler frequency, the effect of mobility on performance is significant.The DAP is limited by a gain

Figure 3 .
Figure 3.Comparison of classical and DAP channel estimations for QPSK modulation.

Figure 4 .
Figure 4. Comparison of classical and DAP channel estimations for 16 QAM modulation.