Abstract

A novel interference alignment (IA) scheme based on the errors-in-variables (EIV) mathematic model has been proposed to overcome the channel state information (CSI) estimation error for the MIMO interference channels. By solving an equivalently unconstrained optimization problem, the proposed IA scheme employing a weighted total least squares (WTLS) algorithm can obtain the solution to a constrained optimization problem for transmit precoding (TPC) matrices and minimizes the distortion caused by imperfect CSI according to the EIV model. It is shown that the design of TPC matrices can be realized through an efficient iterative algorithm. The convergence of the proposed scheme is presented as well. Simulation results show that the proposed IA scheme is robust over MIMO interference channels with imperfect CSI, which yields significantly better sum rate performance than the existing IA schemes such as distributed iterative IA, maximum signal-to-interference-plus-noise ratio (Max SINR), and minimum mean square error (MMSE) schemes.

1. Introduction

The importance of the interference management has been emphasized in [1, 2]. It is of paramount importance to conceive efficient interference management schemes for multiuser wireless networks. Interference alignment (IA) has been proposed as a powerful and promising technique to mitigate the interference in multiuser interference channels [13]. In [1], Jafar and Shamai first characterized the degrees of freedom (DoF) for two-user multiple input multiple output (MIMO) interference channels. The earlier studies [13] of his group were focused on the DoF for various distributed systems. Since then, IA technique has been used to structure interfering signals to occupy a reduced-dimensional interference subspace at the receivers and thus maximize the multiplexing gain (or the sum rate of system).

Recently, various algorithms for IA and IA-inspired schemes in different scenarios [48] have been proposed. Among them, maximum signal-to-interference-plus-noise ratio (Max SINR) scheme [5] so far has been regarded as one of the most effective schemes, without being proven to be convergent. Some other kinds of iterative schemes, such as iterative CJ08 scheme [5] and MMSE scheme [9], are also widely used in IA. And most of these former works rely upon perfect channel state information (CSI). However, realistic practical channel estimation (CE), feedback, quantization, and so forth have errors. Therefore, CSI at the transmitter and receiver is far from being perfect. It is revealed in [10] that the performance of IA scheme is sensitively degraded due to the uncertainties of CSI. Lately, several works on IA schemes with imperfect channel knowledge have focused on this aspect. In [11], the authors investigate an iterative beamforming design algorithm so as to maximize the sum outage rate for MIMO interference channels. Specifically, in [12, 13], the authors propose a minimum mean squared error (MMSE) scheme that takes CSI error into account to align the interference. This MMSE approach is further extended to the robust MSE-based iterative transceiver designs in [14], which effectively improve the BER performance of MIMO interference system. A robust lattice alignment method designed by using the existing conventional Iterative IA algorithm is presented in [15] for quasistatic MIMO interference channels with imperfect CSI as well. In addition, [16, 17] investigate the effect of CSI mismatch and develop several adaptive interference alignment schemes with CSI error.

In contrast to previous works above, we, in this paper, take the errors-in-variables (EIV) model into consideration from the essence of IA with imperfect CSI and employ an unbiased parameter estimation technique in such a linear measurement error model instead of the conventional method of minimizing projector distances of interference subspaces. Hence, the estimation problem can be easily converted into the problem of solving a linear system of equations. Our goal is to overcome the distortion caused by the imperfect CSI at transmitters and optimize the sum rate performance under a given and feasible DoF. To this end, an EIV-based system model of interference management is established. Meanwhile, a weighted total least squares- (WTLS-) based estimation algorithm is proposed to update the transmit precoding (TPC) matrices for the implementation of accurate IA. Unlike the Max SINR scheme, our proposed WTLS-based (or EIV-based) IA scheme is proven to be convergent. Simulation results show that the proposed WTLS-based IA scheme can improve IA performance under the scenarios with different variances of error for direct link estimation and interference links estimations. Against this background, the major contributions of this paper are summarized as follows:(1)We firstly introduce the EIV model into IA. More explicitly, we set out to establish an EIV-based system model of an IA scheme with imperfect CSI for the sake of exploiting the unbiased WTLS parameter estimation algorithm to estimate TPC matrix in such a linear measurement error model.(2)We propose a WTLS-based IA scheme to convert the IA-based constrained optimization problem under the condition of imperfect CSI into an equivalently unconstrained optimization problem in the transmitted signal estimation. More explicitly, the unconstrained optimization problem is derived by minimizing analytically over the correction of the measurement data matrix, and an iterative algorithm for the solution of the unconstrained optimization problem is presented. What is more, our WTLS-based IA scheme is proven to be convergent by employing the large data matrix.The rest of this paper is organized as follows. In the next section, the EIV-based system model is first introduced. In Section 3, we present the proposed WTLS-based IA scheme with CSI uncertainties and discuss the convergence issues. In Section 4, we demonstrate the proposed algorithm with numerical simulations. At last, we conclude with Section 5.

Notation. The superscript , superscript , and superscript stand, respectively, for conjugate, transpose, and Hermitian transpose. Matrices and vectors are set in bold-face letters. or denotes the th row of , and or denotes the th column of . is the identity matrix and is the all-zero matrix. is the th element of . denotes the covariance matrices of a random vector . , , and denote -norm, the expectation, and the determinant of , respectively. is the space of complex matrices. is complex Gaussian distribution with mean and covariance matrix .

2. EIV-Based System Model

In this paper, we consider -user IA scheme over MIMO interference channels. As shown in Figure 1, each transmitter and receiver is equipped with and antennas, respectively. The channel output at receiver is defined as follows:where represents the received signal vector at receiver and is the transmitted signal vector at transmitter . represents the channel fade coefficients from transmitter to receiver for , where each element is assumed to obey an independent and identically distributed (i.i.d.) complex Gaussian random process with zero mean and unit variance . is noise vector at receiver , which is complex-valued additive white gaussian noise (AWGN) with elements normal distributed as . is the TPC matrix with unit-norm and linearly independent columns at transmitter , where is the DoF to meet feasibility of IA scheme at transmitter . represents the data streams from transmitter , where the transmitted streams are i.i.d. such that .

Furthermore, the data streams received at receiver can be described aswhere is receiving filter matrix with unit-norm and linearly independent columns at receiver .

In most IA schemes, perfect CSI is assumed. However, CSI is far from being perfect in realistic practical system. We assume that only imperfect estimation of global CSI is available at each terminal. The feedback channels are assumed to be error-free. The imperfect CSI model can be quantized according towhere is the estimated channel matrix from transmitter to receiver and is the CSI error matrix between the true and available information. This error matrix is assumed to be i.i.d. zero-mean complex Gaussian, andwhere the parameters and indicate the variance of error for direct link estimation and interference links estimations [10], respectively. From (3) to (6), we know and . If , , then the imperfect CSI model reduces to the case that perfect CSI is available at terminals, and when , , it represents that no CSI is available terminals. The parameters and can potentially have different values corresponding to different accuracy of channel estimation of the direct links and that of the interference links.

Due to imperfect CSI, then (2) can be written as

From (7), we hardly find that the IA scheme with CSI error is relative with EIV model. However, when the last two terms on the right-hand part of (7) are moved to the left-hand part, it will be clear that

Obviously, from above (8), the terms of interference and noise are considered as the measurement errors of the interested parameter , while the CSI error term is considered as the measurement error of the interested parameter . Definitely, this IA system model with imperfect CSI in (8) satisfies the EIV model [18, 19].

Briefly, the above EIV-based system model can be described by a corrected system of equations [18, 19] , where and is the parameter to be estimated, while and are the perturbations of and , respectively.

As a result, the transmitted signal estimation problem for the IA with imperfect CSI can be defined as a constrained optimization problem: an appropriate cost function depending on the data is minimized over the estimated parameters.

The classical method, so-called generalized total least squares (TLS) [18, 2022], is proposed as a parameter estimation technique for the EIV model when all elements in channel matrices are perturbed by i.i.d. Particularly, the WTLS [19, 2325] became popular in the mathematics field during this decade since its estimator has a better statistical accuracy under more general noise assumptions. So in the next section we will focus our attention on the WTLS in a realistic CSI error-based interference channel.

3. WTLS-Based Interference Alignment Scheme with Imperfect CSI

In this section, the proposed WTLS-based IA scheme will be investigated to solve the IA-based constrained optimization under the condition of imperfect CSI. While minimizing analytically over the correction matrix, the constrained optimization problem can be changed into an equivalently unconstrained optimization problem. And an effective way to set the weights will also be considered. At the end of this section, an iterative algorithm will be presented for the sake of solving the WTLS-based optimization.

3.1. IA-Based Constrained Optimization with CSI Uncertainties

Looking back on (8), we assume that perfect IA is available and is the true value of its approximation estimation such thatwhere we start with arbitrary TPC matrix and receiving filter matrix. Here we define a measurement data matrix with CSI errors at th receiverwhere and . According to (3), the measurement data matrix can be given bywhere , , , and . is the true value of data matrix , and is measurement error. And , is th independent column of ; it is with covariance matrices (covariance information of measurement error)

We define an extended matrix of the approximate estimation

An alternative expression for the EIV-based system model (8) iswithprovided that the initial data matrix is available and the error weights information is , corresponding to each row of data matrix , which totally depends on the quality of CSI. When the data matrix is confirmed, we set out to seek for a minimal correction to compensate for the measurement errors , so that the WTLS optimization [23] in the presence of a mixture of perfect CSI and CSI error can be defined aswhere . The optimization variables are and . Let be an optimal point of the WTLS optimization problem. Thereupon, is an optimal WTLS estimation of the true value and is an optimal WTLS estimation of the true data .

The first step of WTLS optimization is to minimize the cost function. The minimization with respect to the correction can be carried out as follows:

For a certain , is fixed and the constraint of (17) can be treated as linear equations to optimize variables . The set of linear equations in can be expressed as

Let us define residual matrix . The th row of the is denoted by ; that is, . As a result, the constraint of (17) can be treated the same way asand thus the constrained optimization problem can be separated according to different . In this way, the constrained minimization problem (17) is transformed to independent optimization problemsand the optimal correction isEmploying (21) and the solution of can be written as

Hence, the solution of (17) becomes the well-known function [18, 26] as follows:where .

As a result, the WTLS constrained optimization problem (17) is changed into an equivalently unconstrained optimization problem as follows:

3.2. Unconstrained Optimization

The second step of the WTLS optimization is to solve the unconstrained optimization problem (24). However, generally, the analytic solution is unavailable [18, 23]. Thus, we propose a numerical method to solve the unconstrained optimization problem.

From the unconstrained optimization problem (24), we have the necessary condition for (23)where is conjugate gradient [27] for the unconstrained optimization problem (24).

According to the matrix theory [27, 28], the conjugate gradient which is presented in the appendix iswhere , , , and are, respectively, the th column of , , , and . Note that, in the EIV-based IA system, and are vectors; thus the solution of is univariate case. In this case, is a single component matrix, while , , and are all single component vectors. The data weights are

Since there is no closed-form solution for (25), we use an iterative algorithm to solve the problem. Let represent the approximate value of the th iterative result and let represent the value of the th iterative result. Thus, we define the approximation function of (24) as follows:

Substituting the expression of (the expression of is detailed in the appendix) into (28), thus (28) can be rewritten as

Hence, the approximation is defined as the solution of the following equation:

Substituting (29) into (30), we obtain a system of linear equations with respect to ,which is equivalent to a standard system of linear equations where

On the basis of (29)–(32), we repeat the iterative process and obtain a series of values and to approximate to accurate value. The iteration continues until , and is a given tolerance.

After the iterative process above, we have the final , approximate to accurate value. Assume that

Finally, the estimated TPC matrix can be readily obtained. Hence, the series of estimated values are optimized to the true value of the TPC matrix by using WTLS-based IA optimization method.

The procedures of the above iteration are outlined in the following algorithm which is presented to solve (24).

Algorithm 1 (WTLS-cased IA optimization algorithm). (1)Start with arbitrary TPC matrix and receiving filter matrix .(2)Set the data matrix with CSI error , the measurement error , and the extended matrix of the approximate estimation . Caculate the error weights information . According to (4) and (5), give the variances of CSI error, while a convergence tolerance of WTLS estimation algorithm should also be given.(3)According to (13), obtain an initial value of WTLS estimation.(4)Begin iteration, and set the iteration counter .(5)Let , for .(6)Compute the residual matrix .(7)Define the iterative approximation function as and work out .(8)Increase the iteration counter .(9)Repeat steps 5 through 8 until .(10)The solution of optimal WTLS estimation is .(11)Finally, according to (13) and constraint of TPC matrix [5], obtain the TPC matrix . Let .(12)Calculate from or MMSE method [5].

In this WTLS-based IA optimization scheme, the sum rate result for stream of user is given by (34) where the in (34) is the desired signal, and the first, second, third, and fourth terms in (35) represent the interference purely by channel uncertainty, interstream interference, other user interferences, and AWGN, respectively:where is the TPC for stream of user and is the power of noise at user . The numerator of (34) is the desired signal power, while the denominator of (34) is interferences plus noise [11]. And the first, second, third, and fourth terms in of the right-hand side of (35) represent the internal perturbation purely by CSI error, interstream interference, other user interferences, and AWGN, respectively. Hence, all sum rate of -user IA scheme can be derived easily as follows:

3.3. Convergence Issue

In this part, we present the convergence issue of the above iterative algorithm for the WTLS optimization in the previous part. The convergence issue equals WTLS estimation which converges to or which converges to in probability [23, 29, 30].

Define the variance of each element of measurement error matrix We allow some of to be equal to 0. The following assumptions to are satisfied:(I)For a fixed , every and every satisfy .(II)There exists such that, for every , .(III)The elements of are nonzero vector.(IV)There exists with such that .(V), and is infinite.(VI), and is infinite.Then the WTLS optimization estimation converges to almost surely (a.s.), as is large enough. The details of the proofs of convergence and consistency can be found in [23]. And the convergence, availability, and consistency of WTLS optimization estimation are well investigated [18, 23, 29, 30].

From above statements, we find that if is large enough, estimation converges to with probability tending to 1. That means the number of rows of the data matrix must be so large that estimation converges to with 100% probability. Mathematically, it seems feasible. However, realistically, it is unfeasible in interference management of multiuser MIMO communication system. Since the DoF cannot be infinite. A compromise is reached by the extending the row spacef of the data matrix where and and are the arbitrary initialized TPC matrix and receiving filter matrix, respectively. The relevant evaluation of will be discussed in the next section.

4. Numerical Results

In this section, simulations are conducted to evaluate the performance of our WTLS-based IA scheme presented in Section 3 for multiuser MIMO interference system with imperfect CSI.

Without loss of generality, the CSI errors are assumed to be i.i.d. zero-mean complex Gaussian, and Gaussian input is considered. size of data matrix is fixed in Figures 2 to 5; that is, the value of in (38) is fixed (here we set ), while the WTLS estimation convergence tolerance . The parameters and indicate the variance of CSI errors for direct link and interference links, respectively. For given and , we first generate and randomly according to a zero-mean complex Gaussian distribution from (4) and (5). After and are generated as such, we generate and according to (3) to (6), and the true channel is determined by (3). In each plot presented in this section, the sum rate performance is computed via Monte Carlo simulations using 1000 independent channel realizations.

Figure 2 shows the sum rate performance of different IA schemes with imperfect CSI for the cases and , respectively, where we have users and antennas at each node. The theoretical maximum DoF which is 4 can be achieved, which means one of the transmitters can obtain 2 DoF and the other two are with 1 DoF. That is to say, this system is unsymmetrical MIMO interference network. It can be observed that the sum rates of all schemes increase as the variance of CSI error decreases. However, the sum rates of all schemes have the error-floor effect because of the internal perturbation purely by CSI error and the interstream interference. From Figure 2, it is seen that the proposed WTLS-based IA scheme outperforms the Max SINR [5], MMSE [13, 14], and Iterative IA [5] (i.e., CJ08 algorithm) schemes for all SNR under different variances of CSI error scenarios. The Max SINR and MMSE schemes show good performance almost comparable to our proposed WTLS-based IA scheme at low SNR. However, as the SNR increases, the sum rates of Max SINR and MMSE schemes increase at slower rate than that of the WTLS-based IA scheme. Particularly, the sum rate of Max SINR scheme degrades to that of MMSE scheme as the SNR increases (these two schemes themselves converge as SNR increases), while both of them still outperform Iterative IA scheme. Generally, Figure 2 shows that there is a considerable gain by using our WTLS-based IA scheme especially under a serious CSI error scenario. Meanwhile, the number of iterations in WTLS-based IA scheme is not more than 10 [19]. It is much less than 100 iterations in both Max SINR and Iterative IA schemes and still less than 16 iterations [13, 14] in MMSE scheme (this MMSE scheme is slightly more complicated as Lagrange multipliers need to be computed to design the transmit beamformers [31]).

Figure 3 portrays the sum rate performance of different IA schemes with imperfect CSI for the cases and , respectively, where we have users and antennas at each node. The theoretical maximum DoF is 6 which can be achieved, which means each transmitter is with 2 DoF. That is to say, this system is symmetrical MIMO interference network. Comparing to Figure 2, although the sum rates of all schemes increase as the variance of CSI error decreases, the increasement of the sum rates of all schemes is limited due to the more inter-stream interference with imperfect CSI. It can be observed that the proposed WTLS-based IA scheme still outperforms the other schemes for all SNR under different variances of CSI error scenarios. Our proposed WTLS-based IA scheme can achieve a better performance than the other schemes even at low SNR. Generally, the proposed WTLS-based IA scheme is proven to have good performance under the scenario of symmetrical multiple interstream interference.

Figure 4 illustrates the achievable sum rate of the proposed WTLS-based IA scheme for different variances of CSI error, where we have users and antennas at each node. It is seen that the achievable sum rate decreases as the variance of CSI error increases. When the variance of CSI error becomes sufficiently small, the achievable sum rate surely approaches the one with perfect CSI. The starlike and dotted lines represent the achievable sum rate of the proposed WTLS-based IA scheme with different variances of CSI error for direct link and interference links , respectively. It can be observed that the sum rate decays as the variance of CSI error only for direct link or interference links increases. The increasement of the variance of CSI error for direct link has a much worse effect on the achievable sum rate than that of the variance of CSI error for interference links . This leads to a conclusion that transmitter can achieve a desired transmission rate by allocating minimum required resources according to CSI.

Figure 5 describes achievable per-stream rate of the proposed WTLS-based IA scheme versus SNR in -user interference channels, where we have . As shown in Figure 5, the achievable per-stream rate monotonically increases, and it saturates in the high SNR regime. However, because of the negative effect of interuser interference, the achievable per-stream rate monotonically decreases as the number of the users increases. It can be also observed that the per-stream rate performance of single stream case for each user is much better than the case of the multiple streams case for each user. The reason is that interstream interference has a negative effect on the rate performance under the CSI error scenario. Generally, under a CSI error scenario (especially a serious CSI error scenario), the more users or antennas the communication system has, the worse the rate performance will be.

Figures 6(a) and 6(b) show the sum rate performance of the proposed WTLS-based IA scheme with different sizes of data matrix for and , respectively, where we have users and antennas at each node. It can be observed that the sum rate increases as the size of data matrix magnifies. In other words, when the value of in (38) increases, the sum rate performance of the proposed WTLS-based IA scheme increases. From both Figures 6(a) and 6(b), it can be found that the sum rate performance increases more slowly when the value of is larger than 64. The sum rate begins to converge as . These observations confirm that if the number of rows of the data matrix is large enough, estimation converges to with probability tending to 1.

5. Conclusions

In this paper, we have considered the sum rate performance of different IA schemes under the condition of CSI uncertainties. We have established an EIV-based system model of IA with imperfect CSI and presented an IA scheme based on WTLS error analysis algorithm for turning the imperfect CSI-based IA optimization (i.e., constrained optimization problem) into an equivalently unconstrained optimization problem in the transmitted signal estimation. Meanwhile, the convergence of the proposed WTLS-based IA scheme is presented as well. Numerical results have shown that the proposed WTLS-based IA scheme can effectively improve the sum rate performance of the MIMO interference channel system compared with several existing IA schemes under different scenarios of CSI error. Numerical results have also shown that the larger the number of rows of the data matrix is, the higher the accuracy of estimation is.

Appendix

Conjugate Gradient of X

, , and are all the functions of argument . For the sake of convenience, they can be expressed by as follows:

Obviously, is a scalar function. We can treat as function with respect to and . According to the conjugate gradient method for unconstraint optimization and the chain rule of matrix derivative [27, 28], we havewhere is a vector, while is matrix. Note that, in the EIV-based IA system, and are vectors; thus the solution of is univariate case. In this case, is a single component matrix, while , , and are all single component vectors. According to the derivative rules of function with vector argument [27, 28], the in the first term of (A.4) can be rewritten asAccording to the derivative rules of function with matrix argument relevant to matrix inverse [28], in the second term of (A.4) can be written as

Then (A.4) can be rewritten as follows:

We can derive (A.7) by deducing the terms of and . Firstly, the derivative term of can be derived from (A.2)

Secondly, derivative term of can be derived from (A.3)Finally, substituting (A.8) and (A.9) into (A.7), we obtain the conjugate gradient of

Disclosure

Zhengmin Kong and Shixin Peng should be considered co-first authors.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Zhengmin Kong and Shixin Peng contributed equally to this work.

Acknowledgment

This work was supported in part by the National Natural Science Foundations of China under Grant 61201168.