Abstract

Variable forgetting factor (VFF) least squares (LS) algorithm for polynomial channel paradigm is presented for improved tracking performance under nonstationary environment. The main focus is on updating VFF when each time-varying fading channel is considered to be a first-order Markov process. In addition to efficient tracking under frequency-selective fading channels, the incorporation of proposed numeric variable forgetting factor (NVFF) in LS algorithm reduces the computational complexity.

1. Introduction

Time-varying frequency-selective fading wireless channels can be modeled by using the tapped-delay-line filter, in which each channel tap coefficient is considered to be an independent autoregressive process [1]. The analytical and simulation results presented in [2] manifest that the first-order Markov channel provides a mathematically tractable model for the time-varying channels. Under such Rayleigh fading environment, the linear least squares algorithm (linear polynomial model-based approach [3]) using variable forgetting factor (LSn-VFF) is developed for channel estimation in [4]. The lag error variance due to time variations and additive white Gaussian noise (AWGN) variance in channel estimation error have a tradeoff relation with each other. The VFF can be determined with the degree of nonstationarity and signal-to-noise ratio by using the well-known least mean square (LMS) algorithm. The LSn-VFF algorithm is reported to perform well only at high signal-to-noise ratios [4]. However, it increases the computational complexity.

Based on an extended estimation error criterion, which accounts for the nonstationarity of signal, a method for determining the numeric variable forgetting factor (NVFF) is presented in [5]. When the signal experiences nonstationarity, NVFF decreases automatically to estimate the global trend quickly using the extended estimation error criterion. On the contrary, NVFF increases under stationary conditions by increasing the memory for accurate estimation. In this correspondence, we propose a channel estimation method using NVFF least squares algorithm combined with polynomial time-varying channel paradigm.

This correspondence is organized as follows. In Section 2, we first describe the time-varying frequency-selective wireless system model. Details about the presented second-order polynomial model-based least squares algorithm using VFF (LSn2-VFF) are given in Section 3, and we also introduce mathematical formulation of NVFF based on the extended estimation error criterion. In Section 4, simulation results are presented to compare the performance of VFF and NVFF under nonstationary environment. Finally, conclusions are given in Section 5.

2. System Model

Let the received signal be𝑦(𝑛)=𝐜𝐻(𝑛)𝐱(𝑛)+𝑣(𝑛),(1) where, 𝐱(𝑛)=[π‘₯(𝑛)π‘₯(π‘›βˆ’1)β‹―π‘₯(π‘›βˆ’πΏ+1)]𝑇 is the transmitted data symbol vector with √π‘₯=Β±1/𝐿, 𝐜(𝑛)=[𝑐0(𝑛)𝑐1(𝑛)β‹―π‘πΏβˆ’1(𝑛)]𝐻 is the frequency-selective time-varying channel coefficient vector, which changes after each symbol period 𝑇𝑠, and 𝑣(𝑛) is the zero-mean AWGN with variance 𝜎2𝑣. The (β‹…)𝐻 is conjugate transposition of the matrix, and 𝐿 is the length of multipath fading channel. Each channel coefficient 𝑐𝑙(𝑛) is an independent stationary ergodic first-order Markov process with correlation coefficient 𝛼=𝐽0(2πœ‹π‘“π‘‘π‘‡π‘ ), where𝑓𝑑 is the maximum Doppler frequency, and 𝐽0(β‹…) is the Bessel function of first-kind and zeroth order [1]. It follows that𝑐𝑙(𝑛)=𝛼𝑐𝑙(π‘›βˆ’1)+π‘žπ‘™(𝑛),(2) where, π‘žπ‘™ is the zero-mean process noise with variance 𝜎2π‘ž. Using estimated channel coefficient vector 𝐰(𝑛)=[𝑀0(𝑛)𝑀1(𝑛)β‹―π‘€πΏβˆ’1(𝑛)]𝐻, the estimated received signal is 𝑦(𝑛)=𝐰𝐻(𝑛)𝐱(𝑛).(3) The estimation error is 𝑒(𝑛)=𝑦(𝑛)βˆ’π‘¦(𝑛) with zero-mean and variance 𝜎2𝑒.

By invoking Taylor’s theorem, the time variations of each channel coefficient is explicitly represented in terms of the polynomial paradigm [3]. It results in π‘π‘™βˆ‘(𝑛)=βˆžπ‘–=0(𝑛𝑖/𝑖!)𝑐𝑙,𝑖(𝑛), where 𝑐𝑙,𝑖 is the time-variation parameter for 𝑙th coefficient such that𝐜(𝑛)=𝐜0𝑛(𝑛)+𝐜1!1𝑛(𝑛)+2𝐜2!2(𝑛)+β‹―,(4) where, πœπ‘–(𝑛)=[𝑐0,𝑖(𝑛)𝑐1,i(𝑛)β‹―π‘πΏβˆ’1,𝑖(𝑛)]𝐻. The least squares algorithm using above second-order channel model is called LSn2 estimation algorithm, which is a modification of least squares algorithm using the first-order channel model (LSn) in [4]. The tracking capability of LSn2 algorithm in the time-varying environment can be further improved by incorporating the variable forgetting factor (LSn2-VFF) without the explicit knowledge of process noise variance.

3. LS Algorithms with Polynomial Channel Model

3.1. LSn2-VFF Algorithm for Time-Varying Channels

The polynomial channel model-based LS-VFF algorithm uses forgetting factor πœ†(𝑛) to update the channel state after each symbol duration. The adaptive weight vector 𝐰(𝑛)=𝐰0(𝑛)+𝑛𝐰1(𝑛)+(𝑛2/2)𝐰2(𝑛) in LSn2-VFF algorithm is implemented by exploiting the least squares algorithm proposed in [4] by Song et al.; it follows that 𝐱𝐗(𝑛)=𝑇(𝑛)𝑛𝐱𝑇𝑛(𝑛)22𝐱𝑇(𝑛)𝑇,𝐖𝐰(π‘›βˆ’1)=𝐻0(π‘›βˆ’1)𝐰𝐻1(π‘›βˆ’1)𝐰𝐻2ξ€»(π‘›βˆ’1)𝐻,(5) where 𝐰𝑖(𝑛)=[𝑀0,𝑖(𝑛)𝑀1,𝑖(𝑛)β‹―π‘€πΏβˆ’1,𝑖(𝑛)]𝐻, πœ€(𝑛)=𝑦(𝑛)βˆ’π–π»(π‘›βˆ’1)𝐗(𝑛),𝐀(𝑛)=𝐏(π‘›βˆ’1)𝐗(𝑛)πœ†(𝑛)𝜎2𝑣+𝐗𝐻,(𝑛)𝐏(π‘›βˆ’1)𝐗(𝑛)𝐖(𝑛)=𝐖(π‘›βˆ’1)+𝐀(𝑛)πœ€(𝑛)βˆ—,1𝐏(𝑛)=ξ€Ίπœ†(𝑛)𝐏(π‘›βˆ’1)βˆ’π€(𝑛)𝐗𝐻.(𝑛)𝐏(π‘›βˆ’1)(6)

The minimum mean square error (MMSE) in channel estimation is 𝐽=𝐸[|𝐜(𝑛)βˆ’π°(𝑛)|2]. For unknown process noise variance, VFF is updated asξ‚»πœ‡πœ†(𝑛)=πœ†(π‘›βˆ’1)+𝜎2𝑣𝐃Re𝐻(π‘›βˆ’1)𝐗(𝑛)πœ€(𝑛)βˆ—ξ€»ξ‚Όπœ†maxπœ†minξ€Ί,(7)𝐃(𝑛)=πˆβˆ’π€(𝑛)𝐗𝐻(𝑛)𝐃(π‘›βˆ’1)+𝐌(𝑛)𝐗(𝑛)πœ€(𝑛)βˆ—πœŽ2𝑣,𝐌(𝑛)=πœ†(𝑛)βˆ’1ξ€Ίπˆβˆ’π€(𝑛)𝐗𝐻(𝑛)𝐌(π‘›βˆ’1)πˆβˆ’π—(𝑛)𝐀𝐻+ξ€·πœ†(𝑛)(𝑛)πœŽπ‘£βˆ’2ξ€Έβˆ’1𝐀(𝑛)𝐀𝐻(𝑛)βˆ’πœ†(𝑛)βˆ’1𝐏(𝑛),(8) where ΞΌ is the step size [6], which controls convergence and stability of the LMS algorithm in (7). The variation range for VFF is [πœ†min,πœ†max] to ensure the bounded nonnegative value of VFF, which increases the computational complexity of LSn2-VFF algorithm. However, it requires the knowledge of 𝜎2𝑣 at receiver.

3.2. LSn2-NVFF Algorithm for Time-Varying Channels

Equations (5)–(6) are used in combination with NVFF to develop LSn2-NVFF algorithm. The speed of adaptation is proportional to the asymptotic memory length 𝑁=(1βˆ’πœ†)βˆ’1 [5]. The memories corresponding to πœ†max and πœ†min are denoted by 𝑁max and 𝑁min, respectively. If process noise variance is small in comparison to variance of AWGN that is 𝜎2π‘žβ‰ͺ𝜎2𝑣 (see the appendix), then 𝜎2π‘’β‰ˆπœŽ2𝑣. The extended estimation error is determined by1𝑍(𝑛)=π‘€π‘€βˆ’1ξ“π‘š=0||||𝑒(π‘›βˆ’π‘š)2.(9)

To ensure that the averaging in above equation is not obscuring the nonstationarity introduced by time-varying channel, 𝑀 is kept smaller than minimum asymptotic memory length, that is, 𝑀β‰ͺ𝑁min. The NVFF is determined by using the extended estimation error in (9). It follows thatπœŽπ‘(𝑛)=2𝑒𝑁maxπ‘β‰ˆπœŽ(𝑛)2𝑣𝑁max𝑍,(𝑛)πœ†(𝑛)=1βˆ’(𝑁(𝑛))βˆ’1.(10)

It takes relatively long time for accurate parameter estimation for a value of NVFF close to unity, when the signal experiences stationarity [5]. Therefore, 𝑁max controls the speed of adaptation. However, a small value of NVFF appears beneficial under nonstationary environment, which is bounded (lower) by πœ†(𝑛)β‰₯πœ†min to guarantee positive nonzero values of NVFF.

4. Simulation Results

For simulations, the BPSK independent and identically distributed data is considered to be input. The presented results are based on the ensemble average of 250 independent simulation runs. Note that we have kept 𝜎2π‘ž=0.001, 𝜎2𝑣=0.01 and 𝑀=3. The channel tracking performances of LSn2-VFF and LSn-VFF algorithms are compared in Figure 1 for 𝐿=3 and 𝑓𝑑𝑇𝑠=0.001, where πœ†max, πœ†min and step-size ΞΌ are empirically chosen as 0.99, 0.75 and 0.005, respectively, for all cases (see [4]). The actual channel coefficient is denoted as β€œTRUE” in Figures 1 and 2. The tracking performance results presented in Figure 1 depict that LSn2 algorithm (with second-order channel model) combats lag noise more efficiently than LSn algorithm, but at the cost of increased computational complexity. It is apparent from the simulation results shown in Figures 1 and 2 that LSn-VFF and LSn-NVFF algorithms supersede LSn2-VFF and LSn2-NVFF algorithms.

Both variable forgetting factors overwhelm the loss in tracking capability caused due to first-order channel model. The simulation results are in good agreement with previous studies [4]. For smoothly fading channels, LSn-NVFF algorithm reduces tracking weight error relatively more as compared to LSn-VFF algorithm (as shown in Figure 3). At 𝑓𝑑𝑇𝑠=0.01, the tracking performances of both adaptive algorithms are observed to be approximately equal. However for fast fading channels, LSn-VFF outperforms LSn-NVFF algorithm due to the inaccurate extended estimation error used in determining NVFF.

5. Concluding Remarks

The LSn-NVFF algorithm not only performs better than LSn-VFF algorithm, but also precludes the need of LMS algorithm in variable forgetting factor updating at each iteration, which in turn reduces computational burden. It is inferred from simulation results that the higher-order polynomial model-based LS algorithms (e.g., LSn2) in conjunction with NVFF are not providing any additional advantage. However, the linear least squares algorithm using NVFF is found to be efficient under slow and smoothly time-varying fading channels.

Appendix

Using (1) and (3), the estimation error is𝑒(𝑛)=𝑦(𝑛)βˆ’=𝑦(𝑛)(𝐜(𝑛)βˆ’π°(𝑛))𝐻𝐱(𝑛)+𝑣(𝑛).(A.1) The above equation can be simplified using (2) as𝑒(𝑛)=(π›Όπœ(π‘›βˆ’1)+πͺ(𝑛)βˆ’π°(𝑛))𝐻𝐱(𝑛)+𝑣(𝑛).(A.2) Under optimum conditions, it is assumed that 𝐰(𝑛)β‰…π›Όπœ(π‘›βˆ’1). It leads to𝑒(𝑛)β‰…πͺ𝐻(𝑛)𝐱(𝑛)+𝑣(𝑛).(A.3) Therefore, the estimation error variance is 𝜎2𝑒=𝜎2π‘£ξƒ©πœŽ1+2π‘žπœŽ2𝑣ξƒͺ.(A.4) However, 𝜎2π‘’β‰ˆπœŽ2𝑣 for 𝜎2π‘ž/𝜎2𝑣β‰ͺ1.