Abstract

Channel estimation problem is one of the key technical issues in sparse frequency-selective fading multiple-input multiple-output (MIMO) communication systems using orthogonal frequency division multiplexing (OFDM) scheme. To estimate sparse MIMO channels, sparse invariable step-size normalized least mean square (ISS-NLMS) algorithms were applied to adaptive sparse channel estimation (ACSE). It is well known that step-size is a critical parameter which controls three aspects: algorithm stability, estimation performance, and computational cost. However, traditional methods are vulnerable to cause estimation performance loss because ISS cannot balance the three aspects simultaneously. In this paper, we propose two stable sparse variable step-size NLMS (VSS-NLMS) algorithms to improve the accuracy of MIMO channel estimators. First, ASCE is formulated in MIMO-OFDM systems. Second, different sparse penalties are introduced to VSS-NLMS algorithm for ASCE. In addition, difference between sparse ISS-NLMS algorithms and sparse VSS-NLMS ones is explained and their lower bounds are also derived. At last, to verify the effectiveness of the proposed algorithms for ASCE, several selected simulation results are shown to prove that the proposed sparse VSS-NLMS algorithms can achieve better estimation performance than the conventional methods via mean square error (MSE) and bit error rate (BER) metrics.

1. Introduction

High-rate data broadband transmission over multiple-input multiple-output (MIMO) channel has become one of the mainstream techniques for the next generation communication systems [1]. The major motivation is due to the fact that MIMO technology, as shown in Figure 1, is a way of using multiple antennas to simultaneously transmit multiple streams of data in wireless communications [2] and hence it can bring considerable improvements such as data rate, reliability, and energy efficiency. In fact, coherent receivers require accurate channel state information (CSI) since the received signals are distorted by multipath fading transmission. The accurate estimation of channel impulse response (CIR) is a crucial aspect and challenging issue in coherent modulation and its accuracy has a significant impact on the overall performance of the communication system.

During last decades, there exist many channel estimation methods which were proposed for MIMO systems [311]. All these methods are categorized into two groups. The first group contains the linear channel estimation methods, for example, least squares (LS) algorithm, based on the assumption of dense CIRs. By applying these approaches, the performance of linear methods depends only on the size of MIMO channel. Note that narrowband MIMO channel may be modeled as dense channel model because of its very short time delay spread. Accurately, broadband MIMO channel is often modeled as sparse channel model [1214]. A typical example of sparse channel is shown in Figure 2. It is well known that linear channel estimation methods are relatively simple to implement due to their low computational complexity [38]. Unfortunately, their main drawback is the failure to exploit the inherent channel sparsity. The second group is the sparse channel estimation methods which use compressive sensing (CS) theory [15, 16]. Interested authors are recommended to refer to [17]. Basically, optimal sparse channel estimation often requires that its training signal satisfies restrictive isometry property (RIP) [18] in high probability. However, designing the RIP-satisfied training signal is a nonpolynomial (NP) hard problem [19]. Also, there exist some proposed methods which are stable with the cost of extra computational burden, especially in time-variant MIMO systems. For example, sparse channel estimation method using Dantzig selector was proposed for double-selective fading MIMO systems [10]. Indeed, the proposed method needs to be solved by linear programming which incurs high computational complexity. To reduce the computational cost, sparse channel estimation methods using greedy iterative algorithms were also proposed in [9, 11]. But their complexity still depends on the number of nonzero taps of MIMO channel.

Unfortunately, the mentioned proposed methods do not have adaptive estimation capability. Adaptive sparse channel estimation (ASCE) methods using sparse invariable step-size (ISS) least mean square algorithms (ISS-LMS) were proposed in [20] for single-input single-output (SISO) channels. However, conventional ISS-LMS methods have two main drawbacks: (1) sensitive to random scale of training signal and (2) unstable in low signal-to-noise ratio (SNR) regime.

To overcome the two harmful factors on channel estimation and extend their applications to estimate MIMO channels, sparse ISS normalized least mean square (ISS-NLMS) algorithms, for example, zero-attracting ISS-NLMS (ZA-ISS-NLMS) and reweight ZA-ISS-NLMS (RZA-ISS-NLMS), were proposed in [21]. It is well known that step-size is a critical parameter which controls the estimation performance, convergence rate, and computational cost. Different from conventional sparse ISS-NLMS algorithms [21], zero-attracting variable step-size NLMS (ZA-VSS-NLMS) algorithm was proposed for ASCE to improve estimation performance in sparse multipath single-input single-output (SISO) systems [22]. Unlike the previous works, this paper proposes two sparse VSS-NLMS algorithms for estimating sparse MIMO channels. The main contribution of this paper is summarized as follows. First, we derive the lower bound of proposed MIMO channel estimator for introducing the research motivation. Second, we extend the proposed VSS-ZA-NLMS for estimating SISO channels in [22] to MIMO channels. Third, a reweighted ZA-VSS-NLMS (RZA-VSS-NLMS) is proposed to further improve the estimation performance of MIMO channels. In addition, we explain the reason why sparse VSS-NLMS algorithms can achieve better performance than conventional sparse ISS-NLMS ones. Finally, Monte Carlo based computer simulations are conducted to confirm the effectiveness of our proposed algorithms via two metrics: bit error rate (BER) and mean square error (MSE).

The remainder of this paper is organized as follows. A baseband MIMO system model is described and problem formulation is presented in Section 2. In Section 3, sparse ISS-NLMS algorithms are overviewed. In Section 4, sparse VSS-NLMS algorithms are proposed and a figure example is also given to explain the difference between ISS and VSS based algorithms. Simulation results are presented in Section 5 in order to assess the proposed methods. Finally, we conclude the paper in Section 6.

Notations. Capital bold letters and small bold letters denote matrices and row/column vectors, respectively. The discrete FOURIER transform (DFT) matrix is denoted by with entries , ; , , and denote the transpose, conjugate transpose, matrix inversion, and absolute operations, respectively; denotes the expectation operator; assume any vector ; and denote -norm, that is, , and -norm, that is, ; is a component-wise function which is defined as for , for , and for , where denotes any component in vector ; represents the channel estimator of .

2. System Model

A frequency-selective fading MIMO communication system using OFDM modulation scheme is considered in Figure 3. Initially, frequency domain signal vector , , is fed to inverse discrete Fourier transform (IDFT) at the th antenna, where is the number of subcarriers and is the number of transmit antennas. Assume that the transmit power is normalized as . The resultant vector is padded with cyclic prefix (CP) of length to avoid interblock interference (IBI). After CP removal, the received signal vector at the th antenna for time is written as , where . As shown in Figure 4, the received signal vector and input signal vector are related by where collects all of the input signal vectors from different antennas at the transmitter, is an additive white Gaussian noise (AWGN) variable with distribution , and th received multiple-input single-output (MISO) channel vector is written as and the matrix-vector form of system model (1) is also written as where received signal vector , noise vector , and channel matrix can be represented, respectively, as follows: where , , is assumed to be equal -length sparse channel vector from receiver to th antenna. In addition, we also assume that each channel vector is only supported by dominant channel taps.

3. Overview of Sparse ISS-NLMS Algorithms

According to the system model in (1), the th updating estimation error can be written as where denotes an MISO channel estimator of the ; denotes receive error vector at the th adaptive update; and is the received signal at the th receive antenna.

3.1. ISS-ZA-NLMS

According to (5), the cost function of ISS-ZA-LMS [23] at the th antenna of the receiver can be constructed as where is a regularization parameter to balance the square estimation error and sparse penalty of . Hence, the corresponding update equation of ISS-ZA-LMS [23] for MIMO channel estimation is derived as for , where and is the ISS. To mitigate random scaling of input signal , based on the ISS-ZA-LMS algorithm in (7), the update equation of improved ISS-ZA-NLMS [20, 23] was proposed as

3.2. ISS-RZA-NLMS

It is well known that ISS-ZA-LMS cannot distinguish between zero taps and nonzero taps as it gives the same penalty to all the taps which are often forced to be zero with the same probability; therefore, its performance will degrade in less sparse systems. Motivated by the reweighted -norm minimization recovery algorithm [24], Chen et al. have proposed a heuristic approach to reinforce the zero attractor which was termed as the ISS-RZA-LMS [25]. The cost function of ISS-RZA-LMS is given by where is the regularization parameter and reweighted factor is the positive threshold. In computer simulation, the threshold is set as which is also suggested in previous papers [26, 27]. The th coefficient of ISS-RZA-LMS channel estimator is then updated as for , where . According to (10), hence, ISS-RZA-NLMS [20, 23] was proposed as Note that the sparse penalty term in (11) replaces these channel coefficients under the threshold as zero.

3.3. Drawback of the Sparse ISS-LMS Algorithms

Comparing the standard ISS-NLMS algorithm [28], sparse ISS-NLMS algorithms have a common ability of exploiting channel sparsity. Without the loss of generality, we derive the steady-state mean square error (MSE) performance of the ISS-ZA-NLMS [23] as for the typical example to illustrate the drawbacks of the sparse ISS-NLMS algorithms. Assuming that denotes the sparse MIMO channel estimator, under the independence assumption, in [25], the steady-state MSE of ISS-ZA-NLMS estimator was derived as where and . To exploit the channel sparsity, should be selected in the range so that . According to (12), the lower bound of depends on the three factors: . However, and are determined by the input signal and additive noise , respectively. Only selecting the smaller step-size can further achieve better MSE performance. However, if small step-size is adopted, it will incur slow convergence speed (i.e., high computation complexity) on overall adaptive channel estimation. Hence, it is expected that large step-size is used in the case of large MSE to accelerate the convergence speed, while small step-size is used in the case of smaller MSE to improve the steady-state MSE performance. Assume denotes the th update MIMO channel estimator using sparse VSS-NLMS algorithms. As , the lower bound of steady-state MSE of sparse VSS-NLMS algorithms is derived as To simultaneously achieve higher convergence speed and lower steady-state MSE performance, we propose sparse VSS-NLMS algorithms for estimating MIMO channels in the next section.

4. Proposed Sparse VSS-NLMS Algorithms for Estimating MIMO Channels

Recall that the ISS-ZA-NLMS algorithm in (8) does not make use of the VSS rather than ISS. Inspirited from the VSS-NLMS algorithm which has been proposed in [29], to improve estimation performance of MIMO channels, sparse VSS-NLMS algorithms are proposed. Unlike the sparse ISS-NLMS algorithm, sparse VSS-NLMS algorithms are time-variant with respect to the accuracy of updating estimators.

4.1. VSS-ZA-NLMS

At time , based on the previous research on the ISS-ZA-NLMS and VSS-NLMS algorithms, VSS-ZA-NLMS algorithm is proposed as follows: where is the VSS which is given by where is a positive threshold parameter which is related to received signal-to-noise ratio (SNR), . According to (15), the range of VSS is given as , where is the maximal step-size of gradient descent. Theoretically, the maximal step-size is less than 2 to ensure the adaptive algorithm stability [28]. Please note that in (15) is given by where is a smoothing factor to trade off VSS and estimation error.

4.2. VSS-RZA-NLMS

The th channel coefficient of is then updated by Then, the matrix-vector form of (15) can also be expressed as Please note that the second term in (16) attracts the channel coefficients , whose magnitudes are comparable to to zeros. For estimating MIMO channels, two proposed sparse VSS-NLMS algorithms are summarized in Algorithm 1.

Input                (1) and ;
       (2) and ;
       (3) for VSS-ZA-NLMS;
       (4) and for VSS-RZA-NLMS;
Output       Channel estimator .
Initialization       ;
        ;
        .
While         or   Do
        ;
        ;
        ;
        ;
        ;
        ;
        ;
        for VSS-ZA-NLMS in (14) or
       
       for VSS-RZA-NLMS in (18);
End        

Remark 1. To better understand the difference between ISS and VSS, based on (8), (11), and (15), it is worth mentioning that step-size for sparse ISS-NLMS algorithm is invariable but the step-size for sparse VSS-NLMS algorithm is variable as depicted in Figure 5, where the maximal step-size and ISS are set as and , respectively. From the figure, one can easily find that ISS is kept invariant. Unlike the ISS, VSS decreases adaptively as the estimation performance increases and vice versa. In other words, sparse VSS-NLMS algorithms adopting VSS for adaptive gradient descend; large step-size is adopted to speed up convergence rate for reducing computational complexity; small step-size is adopted to ensure algorithm stability in the case of high-accuracy estimator for further improving estimation performance.

5. Computer Simulations

To confirm the effectiveness of the proposed methods, two metrics, that is, MSE and BER, are adopted for performance evaluation. Channel estimators are evaluated by average MSE which is defined by and system performance is evaluated by the BER metric which adopts different data modulation schemes, such as phase shift keying (PSK) and quadrature amplitude modulation (QAM). The results are averaged over 1000 independent Monte-Carlo runs. The length of each channel vector is set as equal length with and corresponding number of dominant taps is set to . Each dominant channel tap follows random Gaussian distribution as and their positions are randomly decided within the length of . In addition, MISO channel vector is subject to . The received SNR is defined as , where is the power of received signal. Computer simulation parameters are listed in Table 1. Based on the research work in [30], it is worth mentioning that threshold parameters of sparse VSS-NLMS algorithms are adopted for 5 dB and for 10 dB and 20 dB, respectively.

In the first example, average MSE performance of proposed methods is evaluated in the case of and in Figures 6, 7, 8, 9, 10, and 11 under three SNR regimes, that is, 5 dB, 10 dB, and 20 dB. To confirm the effectiveness of the proposed method, we compare it with previous methods, that is, ISS-NLMS [28], VSS-NLMS [29], and sparse ISS-NLMS [23, 25]. In addition, to achieve a better steady-state estimation performance, regularization parameters for sparse VSS-NLMS algorithms, that is, VSS-ZA-NLMS and VSS-RZA-NLMS, are adopted from [27], which depend on the number of nonzero taps of a channel. In the case of different SNR regimes, for example, 5 dB, 10 dB, and 20 dB, as shown in Figures 6, 7, 8, 9, 10, and 11, two proposed methods achieved better estimation performance than sparse ISS-NLMS ones.

As it can be observed from Figures 5 and 6, since sparse VSS-NLMS algorithms take advantage of the channel sparsity as for prior information, hence they achieve better estimation performance than standard VSS-NLMS algorithm, especially in a very sparse channel case, for example, . The sparse VSS-NLMS algorithms can exploit much more sparse information for sparser channel. It is obviously observed that the performance gaps between proposed methods and previous methods in Figure 6 () are bigger than the gaps in Figure 7 ().

In the second example, system performance using proposed channel estimators is also evaluated with respect to BER performance. Two kinds of signal modulation schemes, that is, multiple PSK and multiple QAM, are considered. Received SNR is defined by , where is the average received power of symbol and is the noise power. In Figure 12, multiple PSK schemes, that is, QPSK, 8PSK, and 16PSK, are considered for data modulation and system performance was evaluated. One can find that there is no big performance difference using QPSK and 8PSK due to high transmission. If the higher-order data modulation is adopted, much better BER performance will be achieved when compared with previous methods, that is, ISS-NLMS, VSS-NLMS, and sparse ISS-NLMS algorithms. In Figure 13, multiple QAM schemes, that is, 16QAM, 64QAM, and 128QAM, are considered for data modulation. One can easily find that the proposed method can achieve a better estimation than previous methods. In addition, we also compare the system performance with respect to different modulation schemes, PSK and QAM. In Figure 14, 16PSK and 16QAM are adopted as for a typical example of performance evaluation. We can find that 16QAM based system performance is better than 16PSK based system using the same channel estimators.

6. Conclusion

Traditional adaptive MIMO channel estimation methods utilize sparse ISS-NLMS algorithms using ISS. One of the main disadvantages of the traditional methods is the inability to balance the convergence speed and the estimation accuracy. In this paper, two sparse VSS-NLMS algorithms were proposed for estimating MIMO channels. Unlike the traditional sparse ISS-NLMS algorithms, the proposed algorithms utilized VSS which can change adaptively as the estimation error. Simulation results were provided to confirm the effectiveness of the proposed methods in three aspects: convergence speed, estimation performance, and system performance. First, convergence speed of the proposed methods using VSS is faster than ISS-NLMS based methods due to the fact that VSS for adaptive gradient descent is more efficient than ISS. In other words, VSS can well balance that fast convergence speed is dominant in the case of large estimation error while high accuracy is dominant in the case of small estimation error. Second, the proposed adaptive estimators can achieve better MSE gain than previous methods in different SNR regimes especially for sparser channels. At last, system performance using the proposed channel estimators can also achieve better BER performance than previous methods especially in high-order modulation signal based systems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper. Authors of the paper do not have a direct financial relation that might lead to a conflict of interests for any of the authors.

Acknowledgments

The authors would like to extend their appreciation to the anonymous reviewers and editor for their constructive comments. This work was supported by the National Natural Science Foundation of China under grants (No. 61201273 and No. 61201275) as well as fundamental research funds for the central universities (ZYGX2013J026).