Robust Transceiver Design for Reciprocal M*N Interference Channel Based on Statistical Linearization Approximation
11 Robust Transceiver Design for Reciprocal
Interference Channel Based on Statistical Linearization Approximation
Ali D. Mayvan a , Hassan Aghaeinia b , Mohammad Kazemi c a b c Department of Electrical Engineering, Amirkabir University of Technology (Tehran Polytechnic), Tehran, Iran. a Email: [email protected] b Email: [email protected] c Email: [email protected] Corresponding Author: [email protected]
Abstract
This paper focuses on robust transceiver design for throughput enhancement on the interference channel (IC), under imperfect channel state information (CSI).
In this paper, two algorithms are proposed to improve the throughput of the multi-input multi-output (MIMO) IC. Each transmitter and receiver has respectively and antennas and IC operates in a time division duplex mode. In the first proposed algorithm, each transceiver adjusts its filter to maximize the expected value of signal-to-interference-plus-noise ratio (SINR). On the other hand, the second algorithm tries to minimize the variances of the SINRs to hedge against the variability due to CSI error. Taylor expansion is exploited to approximate the effect of CSI imperfection on mean and variance. The proposed robust algorithms utilize the reciprocity of wireless networks to optimize the estimated statistical properties in two different working modes. Monte Carlo simulations are employed to investigate sum rate performance of the proposed algorithms and the advantage of incorporating variation minimization into the transceiver design.
Keyword:
Approximation, Channel State Information, Reciprocity, Robust, Statistical linearization. Introduction
Normally, wireless network scenarios such as interference channel (IC), share the channel among the users, resulting in multi-user interference. A new method, termed interference alignment (IA), leads to the efficient use of communication resources, since it successfully achieves the theoretical bound on the multiplexing gain [1-3]. In this scheme, unwanted signals from other users are fitted into a small part of the signal space observed by each receiver, called interference subspace. The other signal subspace is left free of interference for the desirable signal. The performance of IA scheme is sensitive to channel state information (CSI) inaccuracies. In [4], and [5] the performance of IA under CSI error was quantified. Asymptotic mean loss in sum rate compared to the perfect CSI case was derived [4]. Multiplexing gain is fully achievable when the variance of the CSI measurement error is inversely proportional to the SNR [4]. Performance of IA for SISO and multi-input multi-output (MIMO) IC where the transmitters are provided with the quantized CSI via feedback can be found in [6] and [7]. Other feedback strategies have been presented in [8]. The IA problem has been studied for the cognitive radio networks in [9-11]. The performance of ad-hoc networks using the IA has been evaluated in [12].
Motivation
In order to maximize sum rate of the MIMO interference network, a beamforming strategy based on the interference alignment is used. Such algorithms have been established by progressive minimization of the leakage interference [13, Algorithm 1] and [14]. The Max-SINR algorithm [13, Algorithm 2], minimum mean square error [15] and joint signal and interference alignment [16] are other algorithms. These schemes are established based on availability of perfect CSI. The performance of transceivers is degraded by the CSI error. Different algorithms are proposed to improve the throughput of the IC, under imperfect CSI. In the literature, precoder-decorrelator optimization is proposed for broadcast and point-to-point systems [17–20]. Specifically, in [17], and [20] the authors consider precoding design for multi-input single-output (MISO) broadcast channel, where it is shown that the precoder optimization problem is always convex. In [19] the authors consider precoder-decorrelator design for MIMO broadcast channel using an iterative algorithm based on solving convex sub-problems. On the other hand, in [18] the authors consider a space-time coding scheme for the point-to-point channel with imperfect channel knowledge. However, these existing works cannot be extended to robust transceiver design for the MIMO IC [21]. Researchers have utilized minimum mean square error criterion to design robust transceiver for the MIMO IC [22]. They improved robustness in presence of the channel uncertainty. Authors in [21] proposed a transceiver design that enforces robustness against imperfect CSI as well as providing a fair performance among the users in the interference channel. They adopted worst-case optimization approach to improve robustness and fairness.
Contributions
This paper focused on sum rate improvement of the Max-SINR algorithm under CSI error. In the first proposed scheme, each transceiver adjusts its filter by maximizing the expected value of the signal-to-interference-plus-noise ratio (SINR). Realized SINRs for different realizations of the CSI error matrices, are samples of the SINR random variable. In some cases, this random variable may have large variance. Therefore, the realized SINRs could be very far from the expected value. In the second proposed algorithm, each node adjusts its transmit/receive filter to minimize the variance. This design hedges against variability. Two robust algorithms are designed based on the reciprocity of wireless propagation. There is a high correlation between the original and reciprocal channel’s gains in communication systems working in a time division duplex (TDD) mode. Reciprocity has been exploited by researchers in [23-25]. Algorithms are implemented in a distributed manner with only local channel knowledge required. In other words, each user only needs to estimate the channel between its transmitter and receiver. Monte Carlo simulations demonstrate that the first proposed algorithm, i.e. expectation maximization (EM), achieves higher sum rates compared to the existing algorithms in [21], and [22]. The second proposed algorithm, that is variation minimization (VM), provides a SINR with low variance. Moreover, VM helps to mitigate the effect of the CSI error, but not as satisfactory as EM scheme. Taylor series expansion is exploited to approximate the influence of CSI imperfection on the statistical properties, e. g. mean and variance. Numerical results show that more accurate approximation can be achieved with less error variance / signal power. For practical applications, when estimated mean / variance is used for maximization / minimization, the proposed transceivers will lead to sum rate improvement / SINR with low variance.
Organization
The rest of this paper is organized as follows. System model is studied in section 2. Sections 3 and 4 propose robust transceivers based on expectation maximization and variation minimization, respectively. Simulation results are presented in section 5 and concluding remarks are drawn in section 6. System Model
In the system model under study, each node works in a TDD mode. At the first time slot, nodes on the left hand side send the data to the right hand side nodes, shown in Fig. 1. At the next time slot, nodes on the left hand side receive the data, indicating a change in the roles of the nodes. (This is clear in Fig. 2 following where the reciprocal network is described). The terms, original and reciprocal channels, are used to distinguish between two time slots. Hence, a transmitter in the reciprocal channel plays the role of an original network’s receiver and vice versa.
Fig. 1.
System model of the Original channel.
MIMO IC with K transmitter-receiver pairs is considered in this paper. The transmitter and the * + receiver are equipped with and antennas, respectively, as shown in Fig. 1. The transmitter sends independent data streams [ ] to its intended receiver. Symbol vector has circularly symmetric complex Gaussian distribution with zero mean and covariance matrix , ( ) . True and estimated channel coefficients between transmitter j and receiver k are denoted by and , respectively. In practical scenarios, there is a mismatch between true and estimated channels, as stated below: (1) , ( ) . The received signal at receiver k is: (2) ∑ ( ) , , ( ) . In order to maximize system throughput, a beamforming strategy based on the interference alignment is used. (3) , , - , . Receiver k decodes the transmitted symbol vector using the interference suppression matrix. (4) ̅̅̅̅ , [ ] . For the K -user MIMO IC, reciprocal channel is obtained by switching the role of transmitters and receivers. True and estimated channels indicated by ⃖ and ⃖ , are matrices. Error matrix is denoted by ⃖ with element distribution ( ) . The relation between the original and reciprocal channels’ gains is ⃖ . The ⃖ is precoder and ⃖ indicates receive interference suppression matrices on the reciprocal channel. Since the receivers of the reciprocal channel play the role of the original network’s transmitters and vice versa, therefore ⃖ and ⃖ . Fig. 2 shows the reciprocal network. Fig. 2.
System model. Reciprocal network is obtained by switching the roles of transmitters and receivers in the original channel.
According to the system model, the SINR value for the data stream at receiver is expressed by (where, ‖ ‖ denotes the Euclidian norm.) (5) ‖ ‖ ∑ ∑ ‖ ‖ ‖ ‖ ‖ ‖ . Robust Transceiver Design I
In this section, the first algorithm is formulated. The algorithm starts with arbitrary transmit and receive filters and then iteratively updates these filters to yield solution. The goal is to achieve robust transceiver by optimization problem , - . The iterative algorithm alternates between the original and reciprocal networks ⃖ [ ⃖ ] . Within each network, only the receivers update their filters. In the following, first, approximate expression for the mean value is computed and then the optimization problem is solved. Estimate the Mean of
In (6) , - is computed in terms of the function and the probability density function ( ) ( ) [ ] ∫ ( ) . Unfortunately, a closed form solution cannot be derived for the integration in (6). Hence, the approximate mean should be found. If ( ) is concentrated near its mean, then estimation of the mean value can be expressed by The random variables, and , are used to represent signal-to-interference-plus-noise ratio as a rational function and are defined by: ,
0 ∑ ∑ . / . / . (7) [ ] ( ∑ ) 1 . Where ∑ ∑ and denote, respectively, the estimated covariance matrix of all data streams observed by the receiver k and estimated covariance matrix of the desirable data stream. Since the estimation of mean is common in two algorithms, it is provided in appendix A. Iterative Solution
To obtain columns of , the derivative of (7) with respect to should be obtained and then set equal to zero. Thus, should satisfy the following vector equation (i.e. the derivative of the numerator multiplied by should be equal to the derivative of denominator multiplied by ). (8) [ ] [ ( ∑ ) ] . [ ] [ ] . Eq. 31 in appendix A [ ] [ ( ∑ ) ] . Eq. 32 in appendix A The above vector equation is rearranged as follow by moving the terms involving and to left and to the right hand side: (where, the scalar coefficient is .) (9) ( ) ( ) , (10) ∑ . According to the definition, is the product of scalar value and vector . It is concluded that is in the direction of . Furthermore, only the directions are important. Hence, the scalar factors , ( ) can be removed from (9). Then, the unit vector maximizes (7) is given by (11) ( ) ‖( ) ‖ . Now, consider the reciprocal network. The transmit precoding matrices, ⃖ , are the receive interference suppression matrices from the original network, whose columns are given by (11). The optimal unit column of ⃖ is given by (12) ⃖ ( ⃖ ⃖ ) ⃖ ⃖ ‖( ⃖ ⃖ ) ⃖ ⃖ ‖ . Now, the receive interference suppression matrices in the reciprocal network replace , and then are updated based on them. It is seen from (10) and (11), that and are interdependent. Therefore, prior to repeat steps, should be computed at the end of each iteration. To summarize the iterative procedure the steps are given in Fig. 3. (a)
Maximize SINR Mean Value Optimize Precoding Matrices Optimize Interference Suppression Matrices Transmit Filter Receive Filter (a) (b) Fig. 3 . (a) Schematic view of robust transceiver design I. (b)
Algorithm: Robust Distributed Transceiver Design.
It can be proved that EM filters for the special case , are transmit / receive matrices Pick arbitrary scalar coefficients , ⃖ , and precoding matrices of size to initialize. Switch to the reciprocal channel and set ⃖ . Compute the interference suppression matrices: ( ) ‖( ) ‖ , * + . Compute ⃖ of new receivers: ⃖ ( ⃖ ⃖ ) ⃖ ⃖ ‖( ⃖ ⃖ ) ⃖ ⃖ ‖ , { } . Switch to the original channel and set ⃖ . REPEAT UNTIL ALGORITHM CONVERGES of the Max-SINR algorithm (Proof is provided in appendix B.). In order to implement the algorithm in a distributed manner, receiver k needs to know about and which are locally available. The covariance matrix can be estimated from the autocorrelation of the received signal . Substituting ∑ into (2) yields (13) ( ∑ ) , where, the expectation is computed over error and noise. The receiver j in the reciprocal channel can learn ⃖ in a similar manner. For TDD systems, the transmitters can estimate the channels from the sounding signals received in the reverse link [21, section II-A]. Using MMSE channel prediction, the CSI estimate is obtained, whereas the CSI error is Gaussian distributed and independent from the CSI estimate [21, section II-A]. Proof of Convergence
Now, the convergence of the proposed EM algorithm is demonstrated. Equivalent problem is considered. EM can be written as follow (14) , where, , and ( ∑ ) are matrices and indicates optimization variable. It is shown in [26] that (14) is equivalent to (15) , s. t. . For the equivalent problem, the Lagrangian function is given by ( ) . / . The solution is the eigenvector corresponding to the maximal eigenvalue of , and the Lagrange multiplier is . █ The metric is defined in (16). It is proved here that each step in the algorithm increases the metric. Since it cannot increase unboundedly, this implies that equivalent problem converges and consequently algorithm Fig 3 also converges. It is important to note that the metric is the same for both original and reciprocal networks. (16) ∑ ∑ ( ) . Accordingly: (17) ∑ ∑ ( ) . In other words, given , Step 1 increases the value of (16) over all possible choices of . The filter ⃖ computed in Step 3, based on ⃖ , also maximizes the metric in the reciprocal channel (18). (18) ⃖ ⃖ , ⃖ ∑ ∑ ⃖ . ⃖ / ∑ ∑ ⃖ ⃖ ⃖ ( ⃖ ⃖ ⃖ ) . Since ⃖ and ⃖ , the metric remains unchanged in the original and reciprocal networks, according to following equation: (19) ⃖ ∑ ∑ [ ] ∑ ∑ . [ ( ∑ ) ] / ∑ ∑ ∑ ∑ . Therefore, Step 3 also can increase the value of (16). Since the value of (16) is monotonically increased after every iteration, convergence of the algorithm is guaranteed. Robust Transceiver Design II
Realized SINRs for different realizations of the CSI error matrices, are samples of the SINR random variable. This random variable can has large variance. Hence, the realized SINR, depending on the particular realization of the CSI error matrix, could be very far from the expected value. Therefore, to hedge against such variability, each receiver adjusts its receive interference suppression matrix based on minimizing SINR variance [27]: (20) [ ] , * + . According to ( ) ( ) ( ) , the ( ) is minimized by minimizing ( ) and maximizing ( ) . These two terms may not attain their best values for the same transceivers in some MIMO IC system model, due to the contradiction between them. It will be shown that the VM scheme presents sum data rate (maximize ( ) ) lower than the Max-SINR algorithm but it enables transceivers to hedge against variability for the primary scenario in simulation results. For the second scenario, the VM improves sum data rate superior to the Max-SINR and provides SINR with low variance, simultaneously. In section 4.1, the value of [ ] is approximated by using statistical linearization argument. The iterative solution is similar to Algorithm I. Estimating the Variance of
Lower bound on the SINR is derived in terms of norms of error matrices in [21]. Lower bound is (21) ‖ ‖ ‖ ‖ ∑ ∑ ‖ ‖ ‖ ‖ ∑ ‖ ‖ ‖ ‖ ‖ ‖ . ‖ ‖ is the Euclidian norm of error matrix between transmitter j and receiver k , denotes by in equation (21). ‖ ‖ ⁄ , is the sum of the second power of real independent Gaussians, each having a unit variance. Therefore, it has a Chi-square distribution with degrees of freedom, . Hence, one can conclude [ ] and [ ] [28]. Equation (21), , is a function of random vector , - . It is clear from covariance matrix, ( ) [ ] , that the variance of each element is sufficiently small for practical applications. Besides, the covariance between each two components is zero. It is concluded that the PDF of is concentrated near its mean and it is negligible outside a neighborhood around the mean value. Again by using the statistical linearization argument, first order Taylor series expansion of around the mean value , - , - is employed to yield: (22) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) , Using (22) for approximating the variance, following equation is obtained. (23) [ ] [( ( ) ( ) ( ) [ ]) ] . Since exact computation of [ ] is not feasible, approximate mean ( ) is used (Estimation of mean is provided in appendix A.). Therefore, estimation of variance leads to (24) [ ] ( ) ( ) ( ) . The elements of are given by ( ) ( ) ( )( )( . ∑ / 1 ) , , ( ) ( )( * ( ) ( ∑ ) + )( ( ∑ ) 1 ) , . According to (24) and (25), the estimated variance is [ ] The is a random vector whose component is the derivative of with respect to . In estimating variance both mean vector and covariance matrix are required [29]; hence, lower bound on SINR is chosen because covariance matrix of random vector is known. (26) ( ) ( * ( ) ( ∑ ) + ) (∑ ( ) ) ( ) ( ( ∑ ) 1 ) Iterative Solution
Briefly, should satisfy the following vector equation: To obtain columns of , the equation [ ] should be solved. Thus, should satisfy the following vector equation (The terms involving and are moved to left and the term is moved to the right hand side). (27) . In the above vector equation, , , and denote scalar coefficients and are expanded as follow (28) (∑ ( ) ) ; (∑ ( ) ) ( ∑ ) (∑ ( ) ) ( (∑ ( ) ) ) ( ∑ ) ; ( ( ) (∑ ( ) ) (∑ ( ) ) ) . In (28), the parameters of the scalar coefficients ( , , and ) should be substituted by: (29) * ( ) ( ∑ ) + ; [ ] ; [ ( ∑ ) ] . The unit vector that minimizes (26) is given by (30) ( ) ‖( ) ‖ , . The iterative procedure, Fig. 4, is algorithmically identical to the Algorithm I. Moreover, the distributed implementation explained previously can be applied again. Fig. 4.
Schematic view of Algorithm II. Simulation Results
In this section, simulation results for the sum data rate and variance of SINR are presented. Also, the influence of the mean value approximation, and other influential parameters on the accuracy of the approximation are determined via Monte Carlo simulations. The proposed robust transceiver design algorithms are evaluated through comparison to the following algorithms: 1.
Leakage interference minimization [13]
Minimize Variance Optimize Precoding Matrices Optimize Interference Suppression Matrices Transmit Filter Receive Filter Max-SINR algorithm [13] 3.
Minimum Mean Square Error (MMSE) [15] 4.
Robust MMSE [22] 5.
Worst-case optimization approach [21]. The use of [13], [15] for comparison is not actually fair because these papers design filter with perfect CSI. The reason is that, their sum-rates do not increase linearly with SNR, and they achieve lower sum-rates compared to our proposed algorithm. (Mean loss in sum rate compared to perfect CSI is shown in Fig. 7 e.g. for Max-SINR algorithm.) [21], [22] design robust transceiver for MIMO interference network, under imperfect CSI.
The channel is modeled as Rayleigh fading.
Channel coefficients [ ] are i.i.d. zero mean unit variance circularly symmetric complex Gaussian. MIMO IC has users with antennas at transmitters and receivers and each user transmits data stream denoted by ( ) . Simulation results are also presented for ( ) MIMO IC. 20 error matrices are generated for a true channel and numerical results are averaged over them. Averaging process over the error is repeated for 20 true channels to eliminate the dependence of the numerical results on the true channel, randomly created; In other words the results are obtained after 400 Monte Carlo simulations. The stopping criterion for the convergence of the proposed iterative algorithms is 100 iterations. All numerical results are based on the SINR associated with the imperfect CSI . Throughput Enhancement
Fig. 5 represents the sum rate comparison between the proposed and basic algorithms for ( )
MIMO IC. The filters are designed with the error variance of . It can be observed that the EM scheme achieves higher sum rates compared to all the other schemes over the entire considered SNR range. Proposed EM scheme achieves 7dB SNR gain over the Max-SINR algorithm at providing 14 b/s/Hz sum data rate and etc. Though it does not seem the VM algorithm improves the overall sum data rate as satisfactory as the EM scheme, but it can be roughly concluded that the VM algorithm achieves data rate as much as (slightly lower than) Max-SINR. Sum rate of algorithm 5 is shown in Fig. 5. It presents sum data rate performance as robust MMSE. Average data rate is defined as the average throughput (i.e. the bits/s/Hz successfully delivered to the receiver). Specifically, the throughput of data stream at receiver is given by ( ), where ( ) and ‖ ‖ ∑ ∑ ‖ ‖ ‖ ‖ ‖ ‖ , and ( ) is the actual instantaneous mutual information. The overall sum rate of the system is given by ∑ ∑ . is SNR in the network, since all data streams are of power P and is noise power at all receivers. Fig. 5.
Average sum data rate versus SNR. , , , . Fig. 6 shows sum rate for ( )
MIMO IC. Again, the EM scheme achieves higher sum rates compared to all the other schemes. In comparison with Max-SINR, it cannot achieve data rate higher than 12 b/s/Hz, while the EM scheme improves data rate up to 16 b/s/Hz. Compared to the Max-SINR, the VM helps to mitigate the effect of the CSI error more effectively, as shown in Fig. 6. The algorithm 5 achieves a data rate only superior to the Leakage interference minimization. It should be noted that, this scheme guaranties to provide better worst-case data rate. Fig. 6.
Average sum data rate versus SNR. , , , . The better performances of the proposed schemes compared to other algorithms are accountable to suitably: 1- Both proposed schemes improve the resilience against SINR degradation due to the CSI error. 2- Approximation (When approximate mean or variance is used for maximization or minimization, the proposed transceivers will lead to sum rate improvement).
Hedge against Variability
Average SINR variance in ( ) and ( )
MIMO ICs are reported in Table I and Table II, respectively. It can be observed that proposed VM scheme achieves lower variance compared to other schemes, when the amount of SNR (>10dB) is significant. In case of small SNR (<10dB), the algorithms have variance close to each other. By Considering Fig. 5 and Table I jointly, it is concluded VM cannot achieve data rate higher than and SINR variance lower than Max-SINR simultaneously. This is due to the contradiction between minimizing ( ) and maximizing ( ) in VM algorithm. Table I
Average SINR variance in ( ) interference network with error variance
SNR (dB) Algorithm 2 5 8 12 14 16 18 20 22 24 [ ] Max-SINR ( )
MIMO IC, since variance of the VM approaches other algorithms, as presented in Table II. On the other hand, VM achieves data rate better than Max-SINR as seen in Fig. 6. Therefore, VM presents a balance between minimizing ( ) and maximizing ( ) in this scenario. Table II
Average SINR variance in ( ) interference network with
SNR (dB) Algorithm 2 6 10 12 14 16 18 20 22 24 [ ] Max-SINR
The cost of the proposed robust design methods compared with the contrast schemes is the complexity since the MMSE [15, equation 11], and Robust MMSE [22, equation 4] need the inverse operation of an N-by-N matrix only once to update (or ) at each user per iteration, whereas the proposed algorithms (equations 11 and 30) require d times. In the SINR maximizing algorithm [13, equation 30], contrast scheme 5 [21, equation 13], the transmit and receive filters are column-wise updated, as complex as the proposed algorithms.
Conformity with Performance Evaluation in [4]
Based on equation 21 in [4], it is expected that mean loss in sum rate compared to perfect CSI case increases unboundedly as SNR increases. In addition, based on equation 22 in [4], the achievable multiplexing gain should be equal to zero. These are confirmed in Fig. 7, where for the larger SNR, the wider gap between curves representing perfect CSI and imperfect CSI can be noticed. Slope of the curves or multiplexing gain when is zero, too. Fig. 7.
Average sum rates achieved by the EM and Max-SINR algorithms for ( ) and ( ) interference networks with . Accuracy of Approximation
It is straightforward to say as decreases, the impact of any error in approximating average mutual information diminishes. On the other hand, signal power scales with as it is obvious from (7). Therefore, any approximation error will be attenuated as decreases. But, here the impact of loss due to the mean approximation and other influential parameters on the accuracy of approximation are studied via Monte Carlo simulations. The EM algorithm is used by ( )
MIMO IC to compute precoding and interference suppression matrices. The filters are designed with two CSI error variances and . The numerical and theoretical values of average mutual information of the MIMO IC versus SNR are depicted in Fig. 8. For , the proposed approximation is within 15% and 18% of the true value when and . For similar statement stands when and . The proposed approximation is within 22.16% and 24.81% of the true value, respectively for and , over the entire considered SNR range. Therefore, by decreasing SNR and the theoretical approximation will approach the true value of average mutual information. Theoretical or approximate capacity curve are found by substituting (7) into ∑ ∑ ( ) . [ ] ∑ ∑ ( . ∑ / 1 ) . Fig. 8.
Theoretical prediction of , - and numerical value for ( ) MIMO IC is shown versus SNR. The filters are designed with EM algorithm and two CSI error variances and . Convergence of VM Algorithm
Since it is hard to prove convergence of this problem theoretically (It is mathematically intractable) it is investigated numerically. Convergence of the iterative algorithm is shown numerically by considering fraction of interference leakage to the received signal parameter [22]. Fig. 9 shows parameter for the proposed schemes versus iterations. In ( )
MIMO IC, EM and VM converge after 10 and 50 iterations as Fig. 9 shows. Proposed algorithms converge after 20 iterations in ( )
MIMO IC. Fig. 9.
Sum of fraction of interference leakage to received signal for two interference networks with , and . Convergence behavior of VM is like other algorithms in two considered MIMO IC scenarios.
Parameter decreases and then remains within a neighborhood of final value.
The stopping criterion for the convergence of the iterative algorithms is 100 iterations. Further results regarding the convergence of the VM in larger MIMO ICs are provided for ( ) and ( ) in Fig. 10. Proposed algorithms converge after 10 iterations (left hand figure). EM and VM converge after 10 and 30 iterations in ( )
MIMO IC. Fig. 10. Sum of fraction of interference leakage to received signal in ( ) and ( )
MIMO ICs with , and . Parameter decreases after each iteration and then remains constant around final value. Conclusion
In this paper, two robust algorithms were proposed. In the EM scheme, filters were adjusted based on the problem of expectation maximization of SINR. The other design minimized the variance of SINR to hedge against variability due to the CSI error. Taylor series expansion was exploited to approximate the effect of imperfection in CSI on statistical properties. Proposed robust algorithms utilized the reciprocity of wireless networks to optimize estimated statistical properties in two different working modes. Monte Carlo simulations demonstrated that the EM scheme improves data rate of MIMO IC under imperfect CSI. The VM algorithm provided SINR with low variance. Moreover, it improved sum rate, but not as satisfactory as the EM scheme. Appendix A Estimate the Mean of If ( ) is concentrated near its mean, then can be expressed in terms of and , as the mean values of and , respectively. represents the conditional expected value of (31) [ ] [ ] , and denotes the conditional expected value of (32) [ ] [ ( ∑ ) ] . According to the statistical linearization argument [29], is approximated by a first order Taylor series expansion around mean value ( ) : (33) ( ) ( ) ( ) ( ) ( ) ( ) . In this case, (6) yields (34) [ ] ( ) ∫ ∫( ) ( ) In (31) and (32), the following equality has been used . / . ( ) ∫ ∫( ) ( ) . The value of integrations in (34) is zero. Therefore, estimation of the mean value can be expressed by [ ] . Appendix B Max-SINR: Special Case of Proposed EM Algorithm
This appendix starts with a relation which is needed to proof Max-SINR is a special case of proposed EM algorithm.
Searle identity [30] for matrix and column vector is (35) ( ) . This identity is applied to the column of receive interference suppression matrix of the Max-SINR algorithm (relation 22 in [13]) and it is simplified as follow (36) ( ) ‖( ) ‖ ( ) ‖( ) ‖ , where is covariance matrix of desired data stream at receiver k and is the covariance matrix of all data streams heard. Matrices are with respect to perfect CSI. (37) ∑ ∑ , . Since for we have and , it can be concluded that receive filters of EM design are receive matrices of Max-SINR. Authors Declarations
Competing interests
The authors declare that they have no competing interests.
Funding
There are not any sources of funding for the research to be declared.
Authors' contributions
AD carried out the preliminary studies, participated in the proposed algorithm, and drafted the manuscript. Dr. MK carried out the figure design and adjustment parameters. DR. HA conceived of the study and participated in its design and coordination and helped in drafting the manuscript. All authors read and approved the final manuscript.
Acknowledgements
There is not anyone who contributed towards the article to be acknowledged by authors.
References [1]
Gou, T., Jafar, S. A.: “Degrees of Freedom of the K user M × N MIMO Interference Channel”,
IEEE Trans. Inf. Theory , Dec. 2010, , (12), pp. 6040–6057. [2] Jafar, S., Fakhereddin, M.: “Degrees of Freedom for the MIMO Interference Channel”,
IEEE Trans. Inf. Theory , Jul. 2007, , (7), pp. 2637–2642. [3] Cadambe, V., Jafar, S. A.: “Interference Alignment and the Degrees of Freedom of the K user Interference Channel”,
IEEE Trans. Inf. Theory , Aug. 2008, , (8), pp. 3425–3441. [4] Razavi, S. M., Ratnaraiah, T.: “Asymptotic Performance Analysis of Interference Alignment Under Imperfect CSI”, in
Proc. WCNC , 2014, pp. 532–537. [5]
Guiazon, R. F., Wong, K., Wisely, D.: “Capacity Analysis of Interference Alignment With Bounded CSI Uncertainty”,
IEEE Wireless
Communication Letters , October 2014, , (5), pp. 505–508. [6] Bolcskei, H., Thukral, J.: “Interference Alignment with Limited Feedback”, in
Proc. ISIT , 2009, pp. 1759-1763. [7]
Krishnamachari, R. T., Varanasi, M. K.: “Interference Alignment under Limited Feedback for MIMO Interference Channels”, in
Proc. ISIT , 2010, pp. 619-623. [8]
Ayach, O. E., Heath Jr., R. W.: “Interference Alignment with Analog Channel State Feedback”,
IEEE Transactions on Wireless
Communications , February 2012, , (2), pp. 626–636. [9] Ma, S., Du, H., Ratnarajah, T., Dong, L.: “Robust joint signal and interference alignment in cognitive radio networks with ellipsoidal channel state information uncertainties”,
IET Communications , June 2013, , (13), pp. 1360–1366. [10] Tang, J., Lambotharan, S., Pomeroy, S.: “Interference cancellation and alignment techniques for multiple-input and multiple-output cognitive relay networks”,
IET Signal Processing , February 2013, , (3), pp. 188–200. [11] Rezaei, F., Tadaion, A.: “Interference alignment in cognitive radio network”,
IET Communications , February 2014, , (10), pp. 1769–1777. [12] Luo, Y., Du, H., Xue, J., Ratnarajah, T., Wilcox, D.: “On the performance of asynchronous ad-hoc networks using interference alignment with renewal process”,
IET Communications , May 2013, , (15), pp. 1628–1637. [13] Gomadam, K., Cadambe, V. R., Jafar, S. A.: “A Distributed Numerical Approach to Interference Alignment and Applications to Wireless Interference Networks”,
IEEE Trans. Inf. Theory , Jun. 2011, , (6), pp. 3309–3322. [14] Peters, S. W., Heath Jr., R. W.: “Interference Alignment via Alternating Minimization”, in
Proc. ICASSP , 2009, pp. 2445-2448. [15]
Schmidt, D., Shi, C., Berry, R., Honig, M., Utschick, W.: “Minimum Mean Squared Error Interference Alignment”, in
Proc . Asilomar,
Kumar, K., Xue, F.: “An Iterative Algorithm for Joint Signal and Interference Alignment”, in
Proc. ISIT , 2010, pp. 2293-2297. [17]
M. Botros and T. N. Davidson, “Convex conic formulations of robust downlink precoder designs with quality of service constraints,” IEEE J. Sel. Areas Signal Process., vol. 1, pp. 714–724, Dec. 2007. [18]
A. Pascual-Iserte, D. P. Palomar, A. I. Prez-Neira, and M. A. Lagunas, “A robust maximin approach for MIMO communications with partial channel state information based on convex optimization,” IEEE Trans. Signal Process., vol. 54, pp. 346–360, Jan. 2006. [19]
N. Vucic, H. Boche, and S. Shi, “Robust transceiver optimization in downlink multiuser MIMO systems with channel uncertainty,” in Proc. IEEE ICC’08, 2008. [20]
M. Payaro, A. Pascual-Iserte, and M. A. Lagunas, “Robust power allocation designs for multiuser and multiantenna downlink communication systems through convex optimization,” IEEE J. Sel. Areas Commun., vol. 25, pp. 1392–1401, Sep. 2007. [21] Chiu, E., Lau, V. K. N., Huang, H., Wu, T., Liu, S.: “Robust Transceiver Design for -pairs Quasi-Static MIMO Interference Channels via Semi-Definite Relaxation”, [Online]. Available: http://arxiv.org/abs/1009.4300, 2010. [22]
Shen, H., Li, B., Tao, M., Luo, Y.: “The New Interference Alignment Scheme for the MIMO Interference Channel”, in
Proc. WCNC , 2010. [23]
Rashid-Farrokhi, F., Liu, K. J. R., Tassiulas, L.: “Transmit Beamforming and Power Control for Cellular Wireless Systems”,
IEEE J. Sel.
Areas Commun. , Oct. 1998, , (8), pp. 1437–1450. [24] Song, B., Cruz, R., Rao, B.: “Network Duality and Its Application to multi-user MIMO Wireless Networks with SINR Constraints”, in
Proc.
ICC , 2005, pp. 2684–2689. [25]
Babadi, B., Tarokh, V.: “A Distributed Dynamic Frequency Allocation Algorithm for ad-hoc Networks”, [Online]. Available: http://arxiv.org/abs/0711.3247, 2007. [26]
Chong, E. K. P., Zak, S. H.: “An Introduction to Optimization”, JOHN WILEY & SONS, 2001, pp. 382-383. [27]
Kleywegt, A. J., Shapiro, A.: “Stochastic Optimization” , [Online]. Available:http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.23.4660&rep=rep1&type=pdf, February, 2000. [28] Papoulis, A., Pillai, S. U.: “Probability, Random Variables and Stochastic Processes”, McGraw-Hill, 2002, pp. 162-163. [29]
Kay , S. M.: “Fundamentals of Statistical Signal Processing: Estimation Theory”,
Prentice Hall, pp. 294-299. [30][30]