Abstract
The Cholesky decompositionblock diagonalization (CDBD) interference alignment (IA) for a multiuser multiple input multiple output (MUMIMO) relay system is proposed, which designs precoders for the multiple access channel (MAC) by employing the singular value decomposition (SVD) as well as the mean square error (MSE) detector for the broadcast Hermitian channel (BHC) taken advantage of in our design. Also, in our proposed CDBD IA algorithm, the relaying function is made use to restructure the quasieigenvalue decomposition (quasiEVD) equivalent channel. This approach used for the design of BD precoding matrix can significantly reduce the computational complexity and proposed algorithm can address several optimization criteria, which is achieved by designing the precoding matrices in two steps. In the first step, we use Cholesky decomposition to maximize the sumofrate (SR) with the minimum mean square error (MMSE) detection. In the next step, we optimize the system BER performance with the overlap of the row spaces spanned by the effective channel matrices of different users. By iterating the closed form of the solution, we are able not only to maximize the achievable sumofrate (ASR), but also to minimize the BER performance at a high signaltonoise ratio (SNR) region.
1. Introduction
Recently, wireless relay networks which are capable of improving the power efficiency, as well as the network coverage, have been studied with a lot of interest because relaying transmission is a promising technique which can be applied to extend the coverage or increase the system capacity. The capacity achieved by a pointtopoint MIMO network has been shown to increase linearly with the minimum number of transceiver’s antennas [1, 2]. Therefore, by employing multiple antennas at the transmitter or the receiver, the system can significantly improve the transmission reliability.
If multiple antennas are applied at both the transmitter and receiver sides, the channel capacity can be enhanced linearly with the minimum number of transmit and receive antennas [3].
Relay precoder designs for such a system have been reported in [4–6]. The problem of designing optimal beamforming vectors for multicasting is hard in general, mainly due to its nonconvex nature. In [4], the authors propose a transceive precoding scheme at the relay node by using zeroforcing (ZF) and MMSE criteria with certain antenna configurations. The information theoretic capacity of the multiantenna multicasting channel is studied in [5] with a particular focus on the scaling of the capacity and achievable rates as the number of antennas and users approaches infinity. In [6], the authors develop one algorithm to compute the globally optimal beamforming matrix at the relay node and characterize the system capacity region.
Most of the works mentioned above assume the availability of perfect channel state information (CSI) at the relay node [7, 8]. In practice, the CSI available at the relay node is usually imperfect due to different factors such as estimation error, quantization, and feedback delay. Interference alignment (IA) is proposed to achieve the maximum degree of freedom (DOF) for the user interference channels [9]. It designs the signals transmitted by all users with perfect CSI in such a way that the interfering signals at each receiver fall into a reduceddimensional subspace. In order to implement IA scheme in the slow fading environment, multiple channels can be used for multiple carriers or multiple antennas [10]. Since these resources are limited, IA scheme with time extension is still efficient to support multiple users. In the past decades, researches on information theory have been exploring the capacity regions of Gaussian interference channels [11, 12]. In the user interference channel, it is proved that the IA scheme can provide the following capacity for each user: Thus, in highSNR regime, the capacity scales linearly with the number of users.
In this paper, we consider the problem of jointly designing the precoders and the relay transformation matrix for a oneway relay MIMO relay system, where all nodes have multiples antennas. Our goal is to use BHC and BD precoding design to decouple MUMIMO channel into a set of parallel independent SUMIMO channels and CDBD algorithm to reduce the computational complexity. In particular, the leakage interference is minimized in order to achieve interference alignment. By iterating the closedform solution and precoding design, we reach the maximum sumofrate capacity and better performance in BER as shown in simulations.
The organization of the paper is as follows: Section 2 describes a general system model for the pairs oneway relay system, the definition of quasiEVD, and global CSI. In Section 3, we propose an iterative CDBD algorithm and optimal precoder design. In Sections 4 and 5, we discuss the ASR, DOF, and computational complexity for efficient channel model. The simulation results are presented to show the good performance of the proposed algorithm for the pairs relayaided system in Section 6, and Section 7 concludes the paper.
Notation. For matrix , , rank , , , , and denote the trace, rank, determinate, transpose, conjugate transpose, and inverse of , respectively. and denote the space of matrices with complex and real entries. stands for the expectation and are the diagonal matrix whose elements on the diagonal are .
2. System Model
In this section, we propose the oneway relay system, whose key idea to structure the quasiEVD channel is using the relay function to cancel the unitary matrices of multiple access channel (MAC) and broadcast hermitian channel (BHC).
2.1. Protocol Description
Consider pairs interference single relayaided system that proceeds in two phases, which are multiple access channels (MAC) and broadcast hermitian channel (BHC) as shown in Figure 1, where transmitter and receiver are equipped with antennas, and the relay node has antennas. The channel coefficients and define links from the source to relay and relay to the destination , where and (decodable condition). The received signal at relay in the MAC phase is given by where represents the additive white Gaussian noise (AWGN) vector with zero mean and variance . The transmitted signal form to relay is obtained by the precoding matrix ; that is, for , where is the transmitted signals form user and is date stream. The proposed precoder can be obtained in two steps as follows: , which will be further discussed in Section 3. The term is subject to a power constraint, with , where is the transmit power at .
In the BHC phase, relay sends which is combined with the linear precoding matrix , to as follows: where the relay precoding matrix is subset of relay filter . We assume that the maximum transmission power at relay node is ; that is, where we have used the assumption that the source signals and the relay noise are independent with each other. Then, the relay broadcasts to the destination nodes and the received signals at can be written as where denotes the additive noise vector at with . Due to the received signal given by (5), the destination can detect the message by the MMSE criterion or where is an linear decode matrix at .
2.2. QuasiEVD and Global CSIT
We assume that the global channel state information (CSI) and the designed precoding matrices are perfectly known at all the nodes; thus, the channel coefficient can be denoted as SVD decomposition or Hermitian of SVD. In our proposed system, the channel matrices may be defined as follows:(a)MAC phase: ,(b)BHC phase: ,
where and are unitary matrices. and are eigen value matrices, where is the element of eigenvalues.
In addition, we propose the channel gain matrix which has its singular value matrix in its middle as well as its eigen matrix and unitary matrix in its right or left side appropriately, which results in the new diagonal matrix. This kind of structure is called quasiEVD. Firstly, we show a result which is helpful to define the quasiEVD equivalent channel as follows: where .
First, we proceed by reviewing the feasibility conditions of interference alignment and cancellation. Next, we turn to structure of the quasiEVD diagonal channel and the problem of the optimization of the precoders and MSE detectors.
3. Optimal Filters Design and CDBD Algorithm
3.1. Interference Alignment and Cancellation
As shown in [15], the IA scheme is a linear precoding technique to align interference in reduced dimensional signal subspace at each receiver. The feasibility conditions for MIMO interference channel (IC) consist of the one interferencefree constraint and a signal space rank constraint. The perfect IA requirements for all areAn efficient distributed algorithm to find matrices and are derived in [16] by using the channel reciprocity. The condition (8a) guarantees that all the interfering signals at destination are aligned in a subspace of dimensions and can be zeroforced by . Condition (8b) guarantees that destination is able to decode all intended data streams successfully. If conditions (8a) and (8b) are satisfied, then the effective channel is free from interference; the structure is feasible for the given DOF .
3.2. Effective Equivalent Diagonal Channel
Due to the SVD of channel in Section 2, the equivalent channel for the total system can be described as where is the relay precoding matrix. To eliminate the quasiEVD channel, we adopt the relay precoding matrix defined as If has full rank, are also with full rank. It implies that both pseudoinverses of and exist. In order to get the optimal leakage interference, the relay filter should satisfy the constraint By substituting (10) into (11), the abovementioned equation can be written as Obviously, the relay function results in optimal leakage interference condition. In order to achieve the optimal leakage interference, it should satisfy the constraint as follows: where , is the relay power constraint shown in (4). Therefore, when interference alignment is feasible, the objective function in (13) can be minimized. By using relay function and (7), we may structure a quasiEVD channel as Subsequently, this efficient channel for the pair of user in total system can be shown in Figure 2.
Therefore, constitutes the useful signal space in which it is expected to observe all symbols transmitted by user , while is the space where all interference is observed. In addition, to make the leakage interference zero, the relaying function can be inserted at the relay.
The total interference leakage at the destination is given by [17] where is the power constraint shown in (4). Based on equivalent channel, (15) can be rewritten as For the perfect interference alignment, the leakage interference should be zero, which means that . This equation is equivalent to the zeroforcing at which is elegantly employed to achieve a good performance in the proposed scheme. The channel state information is perfectly known at every node; the optimization problem in (15) can be written as where is the transmit power at relay. It shows that the optimization problem contains only and ; we will further discuss details in next section.
3.3. Global Optimal Precoder and Detector Design
The proposed optimal precoder design involves two steps, that is, MMSE detector design at destination and optimal precoding design at transmitter. It contains two phases as follows.
3.3.1. MMSE Detector Design
For the abovementioned parameters, the sum of leakage interference can be reshaped as and it may be given as follows by denoting that : which is the optimal MMSE decoder design proved in Appendix A. Therefore, the minimum is equivalent to sum of least dominant eigenvalues of .
3.3.2. Optimal Precoding Design and Iterative Algorithm
Based on MMSE detector , precoding matrices at source nodes should be collaboratively designed. To simply discuss the optimization problem, we assume that the noises are with same variance; that is, . By using optimal MSE detector design shown in (19), the MSE matrix of the signal waveform estimation at receiver can be denoted as or where .
Lemma 1. The optimal precoding matrices design is a convex optimization in highSNR region. For proof see Appendix B.
By applying the MMSE inversion to the combined channel matrix, we have where is the combined equivalent channel matrix; that is, and is the regularization factor. Considering a highSNR case, it can be shown that approaches zero and we have . This means the off diagonal block matrices of converge to zero with high SNR. In addition, we exclude the th pair user’s channel matrices and define and as Thus, the equivalent excluded channel may be denoted as Obviously, the matrix is approximately in the null space of which can be expressed as Considering the SVD of , we have where and are unitary matrices and is eigen value matrix. Since and are invertible, we have Thus, satisfies the BD constraint to balance the interference and the noise term. Therefore, the first step precoding design is completed with result . On the other hand, the interference generated to the other users is determined by . Thus, the final precoder for user may be obtained as After the precoding process, the MUMIMO channel is decoupled into a set of parallel independent SUMIMO channels by the BD precoding. In order to decode the desired signals at the corresponding receivers, the following constraints should be satisfied [9]: where the precoder is subject to the signal space. We can optimize the precoder matrix tailored to individual rate. Consequently, the total leakage interference is As the variance of noises and is small enough in the wireless systems, the convexity can be ensured by substituting (10) and (27) into (29). While it is hard to derive a closedform solution for (29), it can be efficiently solved using the optimal package provided in [18]. Therefore, the minimum is equal to the sum of the least dominant eigenvalues of ; therefore, the optimal precoder and decoder design are completed.
The proposed relayaided interference alignment algorithm is given in Algorithm 1. By employing the minimization technique, it can iteratively update the coding vectors at transmitters, the zeroforcing vectors at receivers, and relaying function at relay to minimize the total leakage interference.

4. Performance Analysis
In this section, we carry out an analysis of the performance of proposed system. We consider a performance analysis in terms of BER, achievable sum of rate (ASR).
For the RBD precoding [13], the residual interference is not zero between the users which is the solution in highSNR region shown as follows: By comparing (26) and (30), we can see that the impact of our proposed precoding would be smaller than that of the conventional RBD precoding algorithm.
Assuming that there exist intersections between desired signal channel and interference signal channel, the following equation will be satisfied: where is the transmitted signals from user . After spanned interference signals into one dimension, we can full cancel them [19]. Therefore, the observations at the relay in (2) can yield where denote column vector of total effective MAC channel matrix with size . Consequently, after the relay filter , the effective propagation of total system is structured and the observations of user for MMSE precoding under the highSNR scenario can be obtained as Consequently, the factor that with rank and with rank , it is simple that the normalization factors and can be written as where the quantity , , , and are the th singular value of , th singular valve of , energy of th, and th stream of , respectively. From (34), the received SNR for th date of user is obtained as Then, the SR upper bound for th user can be calculated as It shows that contains only normalization factors and . The maximum value of is achieved only and only if ; thus, the ASR for total system at highSNR region can be expressed as Therefore, the total achievable DOF for this network can be represented as the sum of DOF for each link [20]. Consider where denotes the DoF for the transmission from user to user .
5. Computational Complexity Analysis
In this section, we will compare the computational complexity of proposed scheme and prior works. We use the total number of floating point operations (PLOPs) to measure the computational complexity. According to [21], the required FLOPs of each matrix operation are described as follows:(i)multiplication of and complex matrices: ;(ii)multiplication of and complex matrices: ;(iii)SVD of and complex matrix where only is obtained: ;(iv)SVD of and complex matrix where only and are obtained: ;(v)SVD of and complex matrix,where only , , and are obtained: ;(vi)inversion of an real matrix using GaussJordan elimination: ;(vii)Cholesky factorization of an complex matrix: /3.
For the conventional RBD method [13], the authors consider a multiuser MIMO downlink precoding system with a base station communicating with users simultaneously. For the nonregenerative MIMO relay systems [14], the authors consider a 3node MIMO relay, where multiple antennas are equipped at the source , the relay , and the destination . We compare the required FOLPs of each precoding algorithm for proposed method, conventional RBD, and nonregenerative MIMO relay system in Tables 1, 2, and 3, respectively, where we assume that and .
For instance, the case denote a system with user , each user with antennas, and total transmit antennas is . The required FLOPs of the proposed method, conventional RBD, and the nonregenerative MIMO relay system are counted as 34638, 40824, and 45306, respectively. From the results, we can see that the reduction in the number of FLOPs and the proposed method precoding are 15.15% and 23.55% as compared to the conventional RBD and the nonregenerative MIMO relay systems. Thus, the proposed algorithms exhibit lower complexity than the conventional RBD and the nonregenerative MIMO relay system approaches, and the complexity advantage grows as , , and increase.
6. Simulation Results
In this section, we show the performance of the proposed scheme in terms of the computation complexity, achievable sumofrate (ASR), and BER performance with some simulation results.
Using Tables 1, 2, and 3, we give the calculated results of FLOPs of the alternative methods in Figures 3 and 4. In the first comparison shown in Figure 3, we consider the case that . We set and express the computation cost as a function of .
In Figure 4, we fix user and while the computation cost as a function of . For conventional RBD method, the orthogonal complementary vector with dimension is obtained; it requires times SVD operations and if we only want to compute , the computational is not efficient. In Step 5, after we got efficiency channel , the second SVD operation should be carried out with dimension , where is the rank of .
For nonregenerative MIMO relay system method, to simply discuss computational complexity, only the indirect link part algorithm is shown. In Steps 1 and 2, two SVD operations are required for the channels from the source to relay and relay to the destination Two variances and are needed to structure as shown in Step 5. Finally, SVD and diagonalize .
For the proposed algorithm, the second precoding matrix is structured by using Cholesky decomposition instead of SVD operation and the first precoding matrix is calculated by SVD of , but only eigenvalue matrices are obtained. Obviously, the proposed method shows a clear advantage in comparisons.
In Figures 5 and 6, we compare the sumofrate of various MUMIMO schemes under full CSI known at each node. The total capacity is obtained by using [22] and the ASR of proposed method is computed using (35), (36), and (37). Figures 5 and 6 illustrate the sumofrate as a function of SNR for and cases, respectively.
In Figures 5 and 6, the nonregenerative MIMO relay systems show a better sumofrate than others at high SNRs, because its capacity includes direct links form source to the destinations and indirect links via relay. The RBD precoding with SVD provides higher ASR than BD at whole SNRs. It is clear that the ASR of our proposed precoding algorithm is lower than the BR at low SNRs, but at highSNR regime, it is higher than SVDRBD and almost same as RBD.
In Figure 7, we compare the BER performance of BDwater filling, RBD, SVDRBD, and proposed method, where QPSK modulation is applied. The proposed algorithm achieved better performance than existing precoding algorithms. As shown in Figure 7, the global optimal scheme in Section 3.3 is evaluated, the reason is that the precoding matrix restricts the interference between the users close to zero while the other precoding algorithm is . The performances significantly improve with increase of SNR.
7. Conclusion
In this paper, motivated by the structure of the quasiEVD based channel in the relayaided system, we have demonstrated a novel iterative algorithm. Our goal is to achieve the maximum sumofrate and the minimum leakage interference. To minimize leakage interference, we use interference alignment to minimize the overlap of the row spaces spanned by the effective channels of different users. The design of the precoding matrix presented in this paper is general, which also can target minimum BER and reduce the computational complexity. In the first step, we use the Cholesky and the singular value decomposition to design the second part of precoder and solve the optimization problem for the total system with MMSE detector. In the next step, we apply the MMSE inversion to the equivalent channel to minimize the BER, which completes the first part of the precoder design. According to the precoding processes, the MUMIMO channel is decoupled into a set of parallel independent SUMIMO channels. Simulation results show that the proposed algorithm outperforms the existing techniques.
Appendices
A. The proof of Optimal MSE
The MSE at receiver can be further expressed as where we have assumed that the signals and noise are independent with each other. Based on (10), the derivation of optimal MSE detection matrix is equivalent to solving the following equation: where . To evaluate the efforts of the result, can be further developed as follows by applying singular value decomposition (SVD) shown in Section 2.2 on BCH channel: where is the eigenvalues of . Then, the closedform expression of can be obtained, which can be expressed as this completes the proof.
B. The proof of Lemma 1
Proof. In the highSNR region, the objective function can be expressed approximately as Since the matrix in the above function is Hermitian and positive definite, we can decompose this matrix using Cholesky factorization as where is an upper triangular matrix. Thus, the MSE can be rewritten as Using equivalent channel , can be denoted as , replace into (B.1), we can rewrite (B.1) as When MSE of the signal waveform estimation is adopted as the optimal problem in (20) which is solved in [23], the precoding matrices at source can be designed as where is a diagonal matrix as power constraint, . Replacing the precoding matrix into , the optimization problem is obtained as where is structured shown in (7); that is, and is the diagonal elements of matrices . Similar to Lemma 2 in [24], is convex if and only if is convex and nonincreasing with and is a concave function of . The Hessian matrices of is which is seminegative definite; it holds that is a concave function of . Thus, Lemma 1 has been proven.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work presented in this paper was supported in part by MEST 2012002521 NRF, BK 21+ Korea, and Double Innovation Plan Eastern Liaoning University.