Channel Parameters Identification Based on IMM Algorithm for Variant Correlation Channel

In wireless communication systems, correct knowledge of the correlation of a fading channel is essential for channel estimation. Both the reliability of the estimated channel impulse response (CIR) and the adjustment of an adaptive communication system need the accurate correlation information, which is difficult to identify especially when changing. By modeling the fading channel as a hybrid dynamic system, a channel estimation algorithm based on Interacting Multiple Model (IMM) is presented with the consideration of time-variant channel correlation. Applying the IMM algorithm, the proposed channel estimator can identify the channel correlation.With the accurate information of channel correlation, the proposed algorithm is capable of performing accurate estimation on the fading wireless channel with time-variant or time-invariant correlation. Our simulations demonstrate that the IMMbased channel estimation algorithm has good performance in estimating CIR as well as in identifying the channel correlation.


Introduction
Information transmission with high data rates and reliable performance is required in wireless communication systems.However, the performance of the communication systems suffers from the signal distortion caused by wireless channel.As a fundamental technology to ensure communication performance, a channel estimation algorithm is required to measure the channel parameters and reduce the influence on the communication system.
In the application of wireless control and wireless sensor systems, high performance communication between maneuvering objects is required.Obviously, the impulse response of the channel is time-variant because of the Doppler effect.Furthermore, the channel correlation which governs the way that the channel varies is also time-variant because of the frequent change in the speed of the moving object.This scenario can be treated as a channel that varies in different modes.In this case, the performance of channel estimation is affected not only by the variance of channel impulse response but also by the changing of channel modes.
In the adaptive communication systems, system parameters, signal modes, and transmission modes can be adjusted according to the channel quality.Therefore, the statistical information of the channel, such as channel correlation, reflecting how the channel changes, is required at each moment to adjust communication parameters.Without the channel statistical information directly provided by the channel estimator, the system parameters can only be adjusted according to the channel information calculated indirectly, for example, bit error rate (BER).Therefore, tracking the channel correlation is important for the channel estimation algorithm to achieve good performance.Channel estimation algorithms that provide CIR as well as the statistical information should be developed.This has not been systematically studied.
Among channel estimation algorithms, Linear Minimum Mean Square error (LMMSE) is widely used, since it is optimum in minimizing the Mean Square Error (MSE) of the estimated channel parameters in the presence of Additive White Gaussian Noise (AWGN).It is shown that LMMSE is very attractive for the channel estimation in [1][2][3].However, 2 Mathematical Problems in Engineering the computational complexity of LMMSE is very high due to the large amount of information operated comprehensively in the estimation algorithm.
The Kalman filter (KF), as an approximation of the Optimal Bayesian filter, is used in a wide range of engineering and econometric applications because of its high accuracy and efficiency on parameter estimating in dynamic systems [4].In wireless channel estimation, the Kalman filter has been applied with the wireless channel described by autoregressive (AR) model [5][6][7].These methods are functional with low computational complexity and the ability of fast tracking the varying channel in the situation that the speed of the channel fading is constant.
In the case of a maneuvering receiver, channel correlation varies into different modes because of the changing Doppler shift corresponding to the relative speed between the transmitter and the receiver.Without information on channel correlation, traditional channel estimation algorithms are limited in performance in such situations, because they cannot effectively respond to the changes in channel mode [8,9].
The Multiple Model (MM) filter was developed to solve the problem of system mode changes [4].The Interacting Multiple Model (IMM) algorithm is one of the most efficient approaches among MM algorithms and is widely used in object tracking [10] and other hybrid dynamic systems [11,12].In communication systems the IMM algorithm can provide accurate channel parameters by applying two Kalman filters to the static and moving states of the receiver [13].
The contribution of this paper is to propose a new channel estimation algorithm based on the IMM algorithm.This algorithm characterizes the fading channel using state-space models and describes the dynamic channel correlation with multiple models.The proposed method is capable of tracking the CIR and identifying the channel correlation which changes according to the maneuverability of the receiver.
This paper is organized as follows.The signal transmission of a communication system as well as the Rayleigh fading multipath wireless channel is modeled mathematically in Section 2. In Section 3, the processes of channel correlation identification and channel estimation based on IMM theory are presented with theoretical proof after the introduction of the KF.The importance of channel correlation identification and the applicability of the proposed algorithm are also discussed in Section 3.After presenting and analyzing a number of simulation results in Section 4, the conclusions are summarized in Section 5.

Transmission System Model.
The transform and transmission of data symbols are illustrated in Figure 1.In the transmitter, a sequence composed of binary data [] is modulated in baseband and transforms into a sequence of complex signals  [𝑘].After being modulated by carrier with frequency   , [] transforms into waveform () and transmitted over a Rayleigh fading channel ℎ(, ), where  is the delay and  denotes current time.In the receiver, after carrier demodulation and sampling, [] is obtained from the receiving signal and can be expressed as follows: where where  is the number of the paths and  channel paths have the significant multipath energy.The time delay of each path is   .ℎ[, ] is the complex gain of the th path in time   .And the superscript " T " in (2) denotes transpose.
In mobile communication systems, the time-variance of the CIR taps is caused by the motion of the transmitter and/or receiver and is quantified by the maximum Doppler shift.When the number of multiple reflective paths is large and there is no line-of-sight signal component, the envelop of every tap of CIR can be statistically described by a Rayleigh probability density function, and the phases of each tap are uniformly distributed in [0, 2].According to the property of the Rayleigh distribution, the real and imaginary components of each tap of CIR are independently and identically distributed (i.i.d.) zero mean Gaussian.As a result, the CIR at time  − 1 and time  can be represented as a first-order AR model [5] in the following form: where k[] is the complex noise vector.Making  = 1, (3) can be written as The real and imaginary components of each element of k[] are zero mean i.i.d.Gaussian.The covariance matrix of k[] is determined by the autocorrelation of h[].For the wide sense stationary uncorrelated scattering model [14,15], the autocorrelation function of the multipath Rayleigh fading channel is [15] E where the superscript " * " denotes the conjugate transposition operation;  = 1/(∑   −/ ) is a normalization constant. 0 () is the zeroth-order Bessel function of the first kind.The maximum Doppler frequency shift in hertz  d is shown as follows: where V is the relative speed between the transmitter and the receiver,  0 is the emitted frequency, and  denotes the speed of wave.

Interacting Multiple Model (IMM) Based Channel Estimation
Assuming that the coherence time of the channel is  coh , to estimate the CIR, the same training sequences with the length of   are sent every   , where   ≤  coh .The received signal of the th training sequence can be written in the form of matrix as follows: where   × 1 vector, y[] denotes the th received training signal after being sampled and can be defined as w[] in ( 7) is the complex noise which distributes independently with zero means in both real and imaginary parts and covariance matrix R[].D in ( 7) is the transmitted training signal matrix.Assuming that the transmitter sends the same training sequence every time, D is a constant matrix which can be written as where the element [] denotes the training sequence.Observing ( 4) and ( 7), it can be found that they are in the typical form of state-space model of a dynamic system.Equation ( 4) is the stochastic difference equation.h[] is the vector to be estimated of the dynamic system.Vector k[] is the process noise with zero mean and covariance Q[] which is denoted as Equation ( 7

Kalman Filter (KF) Theory. Separating the real and image parts h[𝑘] and y[𝑘], they can be written as
where Re{⋅} and Im{⋅} are the real part and the image part of the complex value.Then ( 4) and ( 7) can be written as where D ex is denoted as The recursion of KF is given by the following equations.
(1) The predicted mean and covariance matrix: (2) The predicted measurement, innovation covariance matrix, and Kalman gain: (3) The posterior mean (estimated value) and covariance matrix: where "̂" denotes the estimated value; the superscript " − " of the time  denotes the predicted value obtained from the estimated value in time  − 1; P[] is the covariance matrix of From the equations of KF, it can be noticed that, to ensure the usability of KF, the covariance matrixes of process noise and measurement noise should be known.However, in real cases, they can only be calculated statistically rather than measured.In that case, noise covariance, say process noise, can be denoted as With (17), KF can provide an accurate noise covariance when the noise covariance is constant.

IMM Based Channel Estimation.
In hybrid systems, the dynamic parameter to be estimated varies in different modes.
In this case, the IMM estimator is one of the best compromises available between complexity and performance, because of its low computational requirements and the accuracy which is almost the same as that of many other algorithms with much higher complexity [10].As a powerful approach to adaptive estimation in dynamic hybrid systems, IMM estimator is functional of solving problems of structural and parametric changes [4,11,16].
In wireless communications, CIR varies randomly.In addition, the channel features, such as channel correlation, change from time to time.As a result, channel estimation can be achieved by applying IMM estimation method.
(1) Channel Correlation Identification.In mobile communications, the CIR varies with the speed reflected by the channel correlation, which is described in (5).On the one hand, the information of channel correlation is needed by the channel estimator to estimate the CIR accurately.On the other hand, channel correlation is important information for the upper layer of the adaptive communication networks to adapt the signal transmission scheme according to the channel feature.From the channel stochastic difference equation denoted by (4), the process noise vector k[] reflects the change of the CIR.And Q[], the covariance matrix of k[], reflects the channel correlation, which can be derived as ) where R ℎ [,   ] and E ℎ [] are the correlation matrix and the covariance matrix of CIR and are defined as follows: The derivation of ( 18) is shown in Appendix A. Seen from ( 5), the gains of each tap are independent.In ( 18), E ℎ [] is in the following form: In practice, E ℎ [] can be obtained by using the estimated CIR.
The elements of R ℎ [,   ] can be calculated by ( 5) and be written as where  d [] is the maximum Doppler frequency shift at time .As E ℎ [] and R ℎ [,   ] are both diagonal matrixes, Q[] is also diagonal and can be denoted as Seen from (23), ( 22), and (6),  d is caused by the relative speed between the transmitter and the receiver and governs the channel correlation, which reflects the varying speed of the CIR.As a result, both Q[] and R ℎ [,   ] can be identified by establishing multiple models on Doppler shift  d and conducting IMM estimation based on these models.Models are established by dividing  d into a series of values which are indexed by  ()  d .When mode  () d is in effect at time , the system mode is denoted as where M is the set of all modal states.  is the amount of models in M. Making Y  the set of all the measurements y[] in (7) from the initial time to time , Y  is denoted as The optimal estimation of R ℎ [,   ] in the MMSE sense is the expectation based on the posterior distribution.The estimated value of R ℎ [,   ] is given by where R() ℎ [,   ] is the estimated value of R ℎ [,   ] based on model  and is defined as () [] is the posterior probability that model  is affected at time  and is given by In (26), R() ℎ [,   ] can be calculated by ( 22) using the value of  d in each model.Consequently, the problem of R ℎ [,   ] identification is reduced to the problem of calculating  () [], which is solved in the following sections.
(2) Channel Estimation.Seen from the discussion above, multiple models are established on Doppler shift  d .The mode changes are governed by a Markov chain with the transition probability matrix Γ, whose element Γ () is denoted as Seen from (29), Γ () is the transition probability from mode  in time  − 1 to mode  in time  and can be obtained statistically.The optimal estimation of CIR can be denoted as ĥ where ĥ() [] is the estimated value of h[] based on model  and is defined as ĥ Y  and  () [] are defined by ( 25) and (28).In optimal Bayes' theorem, where Pr{y } is the prediction model probability, which can be obtained by As a result,  () [] can be described in the following form: Seen from ( 30) and (34), in optimal Bayes filters, every time CIR is estimated, the history of measurements Y  is needed.The longer time that the estimator works for, the larger the amount of calculation is required.In view of this, by suboptimally approximating Bayes filter, fixed memory algorithms have been proposed [4,17].The generalized pseudo-Bayesian filter of order  (GPB) approaches the optimal algorithm taking a memory history of  − 1 steps into account.Hence, in channel estimation    filters are needed in each time .IMM estimator performs as well as GPB2, with the lower computational requirements than GPB1.Thus the IMM algorithm has enjoyed remarkable success as one of the most cost-effective schemes for the estimation in hybrid systems.
The IMM channel estimator is recursive.  Kalman filters, corresponding to each mode, run in parallel in each cycle, which consist of four steps: interacting/mixing, filtering, and combination.The structure of IMM channel estimator is also shown in Figure 2, where ĥ() [], P () [], and R ()  ℎ [,   ] are the estimated CIR by KF  at time , its covariance, and its correlation; ĥ(0) [] and P (0) [] are the mixed condition for KF  at time  and its covariance; ĥ[], P[], and R ℎ [,   ] are the combined estimated CIR at time , its covariance, and its correlation;  () [ − ] and  () [] are the predicted and posterior probabilities for KF ;  () [] is the likelihood function of KF .In each cycle, the estimated CIR and model probability  () [] are firstly interacted with each other and mixed for each mode under the assumption that this mode is in effect at the current time.Therefore the initial condition for each filter based on certain model is obtained.Secondly, the interacted/mixed initial conditions are sent to their particular model-based filters.In channel estimation, KFs based on different models are applied in the filtering step.The estimated CIRs and the updated model probabilities based on each model are obtained.After the filtering step, the outputs of each filter are combined as the final result at the current time.Seen from the upper description and Figure 2, as a weighting factor, the model probability  () [] is important to the interacting/mixing and combination steps.In IMM algorithm, the model probabilities are given by Mathematical Problems in Engineering [k] ĥ(N  ) [k] P (1) [k] where [] is the normalization constant and N[u () []; 0, S [𝑘]] denotes the probability of u () [] in the normal distribution with zero mean and the covariance of S ()   [𝑘].u () [] and S ()  [] are the residual and its covariance at time  for the model  based KF and illustrated in the summary of IMM based channel estimator in Appendix B.

Discussions.
The data format suitable for IMM based channel estimation algorithm on the basis of single carrier communication system is shown in Figure 3(a).Training sequences are sent frequently.The time interval between the two training sequences is   , with the assumption that the channel correlation R ℎ [,   ] is big enough and the influence of CIR variance during   can be omitted.In an adaptive communication system, R ℎ [,   ], which reflects the speed of channel variance, can be used to adapt the transmitter.And   can be adjusted according to the channel correlation to achieve better reliability.Moreover, with the help of software radio, adaptive communication can also change its method of communication, such as modulation and coding, with respect to the channel correlation.
From the previous analysis, although presented on the basis of single carrier communication systems, the IMM based channel estimation algorithm can also be applied in multiple carrier communication systems, such as OFDM systems in two ways, training mode and pilot mode, whose training symbols are arranged as in Figures 3(b) and 3(c).In training mode, training symbols are sent in all subcarriers at the same time; the IMM channel estimation is functional by setting the elements of matrix S in (9) as the sampled signal at the receiver.In pilot mode training symbols are sent in some of the subcarriers with data symbols sent in other subcarriers at the same time.A multiple carrier communication is equivalent to the signal transmitted in parallel.Therefore, IMM based channel estimation can be applied separately in subcarriers with training symbols.In each subcarrier, the wireless channel can be seen as a Rayleigh flat fading channel.The CIR at each subcarrier can be estimated by IMM based channel estimation algorithm, and then these estimates are interpolated via different methods to obtain the general CIR.

Simulation Results and Performance Analysis
In this section, the performance of both KF and IMM based channel estimation algorithms is compared and analyzed  The estimated value is more accurate than KF when the model is set precisely.A similar result can be found in BER which is shown in Figure 6. Figure 6 is the relation between BER and   / 0 with different channel estimation algorithms applied in the receiver.It can be concluded that the channel estimation algorithm with better MSE performance can lead to a better BER performance.
In the second case, the maximum Doppler shift  d switches between 100 Hz and 200 Hz every 150 ms.The IMM channel estimator is carried using two models ( (1)  d = 100 Hz, The channel correlation R ℎ [,   ] is analyzed first.Figure 7 4 and 5 in the first case, the KF performs much worse when R ℎ [,   ] is time-variant, because a reliable time-variant R ℎ [,   ] cannot be obtained by KF.Consequently, the BER performances of the algorithms are affected by the accuracy of the estimated CIRs. Figure 11 shows the BERs in different measurement noise levels.It is the same as the pervious analysis that the IMM channel estimation algorithm performs better than the KF method in BER when R ℎ [,   ] changes in time.

Summary and Conclusion
Correlation of a fading channel is important information for channel estimation.The performance of traditional channel estimation with the consideration of time-variance channel correlation is limited.In this paper, by modeling wireless fading channel as a hybrid dynamic system, the channel estimation algorithm based on IMM is presented.The proposed algorithm can identify the channel correlation with the help  of the IMM algorithm.With accurate information of channel correlation, the proposed algorithm is capable of performing accurate estimation on the fading wireless channel with timevariant or time-invariant correlation.The results of a number of simulation experiments show that the proposed algorithm is efficient with good performance in estimating CIR as well as the channel correlation.

A. Derivation of the Process Noise Covariance
From (4), the following equation can be obtained: Calculating the expectation and taking (10) and (20) into the upper equation, 2) The following equation can be obtained from (A.

Figure 1 :
Figure 1: Diagram of a communication system.
) is the measurement equation, where y[] is the measurement of h[], measured by matrix D and polluted by the measurement noise vector w[] with covariance R[].

Figure 3 :
Figure 3: Typical arrangement of training symbols.

Figure 4 :
Figure 4: MSEs of CIR amplitude with constant maximum Doppler frequency.

Figure 8 :
Figure 8: MSEs of channel correlation with time-variant maximum Doppler frequency.

Figure 9 :
Figure 9: MSEs of CIR amplitude with time-variant maximum Doppler frequency.