Robust transceiver design for reciprocal M × N interference channel based on statistical linearization approximation

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, M and N 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.


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][2][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][10][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][18][19][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][24][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.respectively.In practical scenarios, there is a mismatch between true and estimated channels, as stated below: ( , ) .
The received signal at receiver k is: ( In order to maximize system throughput, a beamforming strategy based on the interference alignment is used. ( Receiver k decodes the transmitted symbol vector using the interference suppression matrix.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 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 1 The random variables, and , are used to represent signal-to-interference-plusnoise ratio as a rational function and are defined by: .
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 ). ( 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 .) .
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.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 ⃖ .
EAT UNTIL ALGORITHM CONVERGES of the Max-SINR algorithm (Proof is provided in appendix B.).
In Accordingly: 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). ( Since ⃖ and ⃖ , the metric remains unchanged in the original and reciprocal networks, according to following equation: 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]: 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 The elements of are given by ( 25) , .
According to ( 24) and ( 25), the estimated variance is The is a random vector whose component is the derivative of with respect to . 3 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.( 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.

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
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 4 .

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 5 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. 4 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 ∑ ∑ .
5 is SNR in the network, since all data streams are of power P and is noise power at all receivers.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.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.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.

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.) .

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.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

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.

Fig. 1 .
Fig. 1.System model of the Original channel.
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 .
Fig. 2. System model.Reciprocal network is obtained by switching the roles of transmitters and receivers in the original channel.
. 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).

Fig. 7 .
Fig. 7. Average sum rates achieved by the EM and Max-SINR algorithms for ( ) and ( ) interference networks with .

Fig. 8 .
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 .

Fig. 9 .
Fig. 9. Sum of fraction of interference leakage to received signal for two interference networks with, and

Fig. 10 .
Fig. 10.Sum of fraction of interference leakage to received signal in ( ) and ( ) MIMO ICs with (16)r to implement the algorithm in a distributed manner, receiver k needs to know 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 channelThe 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.
MIMO IC, since variance of the VM approaches other algorithms, as presented in TableII.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.
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 Theoretical or approximate capacity curve are found by substituting (7) into 6