1. Introduction
Recent academic and industry interest has shifted from point-to-point multiple-input multiple-output (MIMO) channels to MIMO broadcast channels [1]. In the MIMO broadcast channel, multiple users can be supported simultaneously by sharing spatial resources. The sum capacity region of the MIMO broadcast channel has been extensively studied [2-12]. The optimal transmit strategy that achieves the MIMO broadcast channel capacity is inspired by dirty-paper coding (DPC) [2], but information-theoretic DPC does not give direct realization of practical transmit strategies.
As a practical and simple solution, linear precoding methods such as zero-forcing (ZF) and block diagonalization (BD) have been suggested [13-17]. ZF is based on the channel inversion, but it suffers from excessive transmit power loss due to the weakest stream [13,18]. BD is a generalization of ZF that supports multiple streams per user. BD makes the effective channel block diagonal; hence multiuser interference is cancelled at the transmitter, but the inter stream interference is remained for the decoder processing. Those linear precoding algorithms achieve sum rate throughput close to DPC in many scenarios with far reduced complexity. But the diversity order of those linear techniques is one, resulting in huge error rate performance degradation in higher SNR regime.
To compensate for low diversity order in linear precoding methods, nonlinear precoding or decoding methods can be used to give more diversity order. Nonlinear precoding methods such as [19,20], however, are still too complex compared to linear precoding methods. Complexity issues for the nonlinear decoder have been solved by recent developments on reduced complexity maximum likelihood (ML) decoders [21-26]. As the number of streams per user is smaller than the number of transmit antennas, the usage of nonlinear technique in decoder is advantageous than using nonlinear technique for precoding. A proper joint design of the linear precoder and the nonlinear decoder for the multiuser MIMO systems is, therefore, a promising practical solution to achieve good throughput and diversity order simultaneously.
One most challenging problem of nonlinear decoder in the multiuser MIMO systems is to obtain a post-processing signal to interference plus noise ratio (SINR). Obtaining exact post processing SINR value per stream is important since the obtained SINR plays a key role for the link adaptation. Link adaptation using adaptive modulation and coding (AMC) mechanism for multiuser MIMO uses the effective SINR as an indicator of the link quality. Unlike to the linear receiver where the SINR value can be easily computed, there is no known explicit way to calculate accurate per stream effective SINR after nonlinear decoding, when the SINR level per stream differs from each other. Due to this limitation, one cannot fully take advantage of both linear-based simple transmit strategy and optimal decoding strategy in multiuser MIMO systems so far.
In this paper, we propose a new precoding and postcoding algorithm for multiuser MIMO with ML based nonlinear decoder, when multiple streams are supported per user. We propose a precoding method that makes the effective channel to have equal diagonal element, as in [27,28], to obtain a precise post-processing SINR. The equal channel gain decomposition such as uniform channel decomposition (UCD) and geometric mean decomposition (GMD) in [27,28] were originally suggested for single user MIMO and optimized for the minimum mean square error (MMSE) with successive interference cancellation or ZF-DPC based decoder. We suggest a two-step precoding method using BD and equal channel gain decomposition to eliminate multiuser diversity and to derive SINR easily. Using the proposed two-step precoding at the transmitter, each user independently uses the ML decoder to get full diversity. The proposed algorithm also generates exact per stream SINR to apply for link adaptation.
The main contributions of this paper are summarized as follows.
• We propose a new joint precoding and postcoding method for multiuser MIMO with multiple stream support per user. The proposed algorithm achieves good sum rate and full diversity gain; hence it shows good error rate performance as well as good throughput performance. Both analytical derivation and simulation results of our proposed algorithm are presented.
• Signaling strategies is presented, reducing feedback overhead while minimizing feedback mismatch. We also present the issues when generating link curve using ML decoder and effective SINR for link adaptation. Using the proposed algorithm, the feedback generation mechanism is simplified and the accuracy of the effective SINR for link adaptation is increased.
• We compare our proposed algorithm to several previously proposed single user and multiuser MIMO algorithms. Comparison includes simple singular value decomposition (SVD) based channel decomposition, Orthogonal Beamforming (OBF) [14], BD [16] and UCD [28]. The comparison includes precoder design, decoding method, achievable sum rate, feedback requirement as well as complexity.
The outline of this paper is as follows. In Section 2, we introduce the system model for the MIMO broadcast channel and introduce basic algorithms of BD and UCD. In Section 2, we propose a new precoder and postcoder for the nonlinear decoder. Sum rate analysis, diversity gain analysis and feedback strategy are also presented in the same section. Simulation results and comparison with previously proposed methods are given in Section 4, followed by conclusions and future works in Section 5.
2. System Model and Prior Works
In this section, we introduce the system model and prior works related to this paper. A brief introduction on BD, GMD and UCD is provided in the following few subsections. Throughout this paper, we use upper boldface A and lower boldface a as a matrix and a vector, respectively. We denote the inverse, transpose, Hermitian of matrix A by A-1, AT and A∗, respectively. ‖𝑎‖ denotes vector 2-norm, unless otherwise stated.
2.1 System Model
Consider a MIMO broadcast channel with a single transmitter and U users, as shown in Fig. 1. The transmitter has Nt antennas and the kth user has Nr,k receive antennas, hence the sum of all the receive antennas are 𝑁𝑟=∑𝑘=1 𝑁𝑟,𝑘. Without loss of generality, we let Nr =Nt for simplicity. Let the channel between the transmitter and kth user is 𝑁𝑟,𝑘×𝑁𝑡 matrix Hk, hence the aggregated channel for all U users can be written as 𝑯=[𝑯1𝑇 ⋯ 𝑯𝑈 𝑇]𝑇. Let the transmit symbol vector for kth user is xk and its corresponding transmit precoding matrix is 𝑭= [𝑭1 ⋯ 𝑭𝑈], with size 𝑁𝑡×𝑁𝑡. Let the 𝑁𝑟,𝑘×𝑁𝑡 receive postcoding matrix Wk, then the received vector at the kth user is given by
\(\mathbf{y}_{k}=\mathbf{W}_{k}^{*} \mathbf{H}_{k} \mathbf{F} \mathbf{x}+\mathbf{n}_{k}=\mathbf{W}_{k}^{*} \mathbf{H}_{k} \mathbf{F}_{k} \mathbf{x}_{k}+\sum_{j=1, j \neq k}^{U} \mathbf{W}_{k}^{*} \mathbf{H}_{k} \mathbf{F}_{j} \mathbf{x}_{j}+\mathbf{n}_{k}\) (1)
Fig. 1. System Model
where 𝒙=∑u𝑘=1 𝒙𝑘 is the 𝑁𝑡×1 transmit vector that each of the component is drawn from finite complex alphabet C. The overall received vector over all the users is 𝒚=[𝒚1𝑇 ⋯ 𝒚𝑈𝑇]𝑇. The 𝑁𝑟,𝑘×1 vector n is an additive complex Gaussian noise vector in the kth user, with mean zero and covariance matrix 𝜎𝑛𝑰𝑟,𝑘. Lk is the number of streams sent to the user k, hence ∑𝑈𝑘=1 𝐿𝑘 =𝑁𝑡. We assume that x is normalized such that the average transmit power is PT.
The goal of this paper is to find the precoding matrix F and decoder combining matrix Wk(k = 1, ⋯ , U) to minimize multiuser interference \(\sum_{j=1, j \neq k}^{U} \boldsymbol{W}_{k}^{*} \boldsymbol{H}_{k} \boldsymbol{F}_{j} \boldsymbol{x}_{j}\) in (1) and maximize sum throughput of the multiuser MIMO link.
2.2 Block Diagonalization (BD)
BD is a generalization of ZF linear precoding when the number of streams per user is more than one. Linear preprocessing is performed at the transmitter by multiplying transmit beamformer F to the transmit signal. Receive postcoder Wk for the kth user is used to overcome dimensionality constraint. To eliminate multiuser interference, BD imposes the zero interference constraint, HiFj = 0 for 𝑖 ≠ 𝑗. The zero-interference constraint between users forces Fj to lie in the null space of \(\widetilde{\mathbf{H}}_{j}=\left[\mathbf{H}_{1}^{T} \cdots \mathbf{H}_{j-1}^{T} \mathbf{H}_{j+1}^{T} \cdots \mathbf{H}_{U}^{T}\right]^{T}=\widetilde{\mathbf{U}}_{j} \tilde{\Sigma}_{j}\left[\widetilde{\mathbf{V}}_{j}^{1} \widetilde{\mathbf{V}}_{j}^{0}\right]^{*}\) where \(\widetilde{\mathbf{V}}_{j}^{1}\) and \(\widetilde{\mathbf{V}}_{j}^{0}\) are right singular matrices corresponding nonzero and zero singular values of 𝐇̃𝑗, respectively. Multiplying \(\widetilde{\mathbf{V}}_{j}^{0}\) to the jth user’s channel matrix Hj and applying singular value decomposition (SVD) yields
\(\mathbf{H}_{j} \widetilde{\mathbf{V}}_{j}^{0}=\mathbf{U}_{j}\left[\begin{array}{cc} \Sigma_{j} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{array}\right]\left[\begin{array}{ll} \mathbf{V}_{j}^{1} & \mathbf{V}_{j}^{0} \end{array}\right]^{*}\). (2)
The product of \(\widetilde{\mathbf{V}}_{j}^{0}\) and \({\mathbf{V}}_{j}^{1}\) make the channel Hj free from the multiuser interference since the right singular matrix \({\mathbf{V}}_{j}^{1}\) corresponding nonzero singular vector is obtained by the multiplication of \(\widetilde{\mathbf{V}}_{j}^{0}\) to the channel. Then the BD precoder is obtained as \(F=\left[\widetilde{\mathbf{V}}_{1}^{0} \boldsymbol{V}_{1}^{1} \cdots \widetilde{\mathbf{V}}_{U}^{0} \boldsymbol{V}_{U}^{1}\right] \mathbf{\Lambda}^{\mathbf{1 / 2}}\), where Λ is power scaling matrix. The first Lj column of the left singular matrix Uj is used for decoder combining. The above process can be seen as a decomposition technique of 𝑁𝑟 × Nt physical channel into U parallel (Nr−(U−1)Nt)×Nt channels [29].
2.3 Uniform Channel Decomposition (UCD)
UCD is a matrix decomposition method that makes a MIMO channel into multiple subchannels with the same singular values [28]. It can be considered as a modified QR decomposition since the output is a multiplication of upper triangular matrix and unitary matrices. The precoder for UCD is given as \(\mathbf{V} \mathbf{\Phi} \mathbf{\Omega}^{*}\), where Φ is power loading diagonal matrix as the water-filling solution \(\boldsymbol{\Phi} \triangleq\left(\mu-\alpha \boldsymbol{\Lambda}^{-2}\right)_{+}^{1 / 2} \quad\left([a]_{+} \triangleq \max \{0, a\}\right)\) with sum power constraint \(\operatorname{tr}\left\{{\Phi}^{2}\right\}=P_{T}\) [30]. 𝜴 is 𝑁𝑡 × 𝑁𝑡 unitary matrix to be derived by UCD process and V is the left singular matrix in SVD. Using the extended channel matrix \(\overline{\boldsymbol{H}}=[\boldsymbol{H} \sqrt{\alpha} \boldsymbol{I}]^{T}\) , where 𝛼 = 1/𝑁𝑡 is the inverse of SNR, the channel can be decomposed as
\(\overline{\boldsymbol{H}}_{U C D} \triangleq\left[\begin{array}{c} \boldsymbol{H} \boldsymbol{F} \\ \sqrt{\alpha} \boldsymbol{I}_{N_{t}} \end{array}\right]=\left(\begin{array}{c} \boldsymbol{U} \boldsymbol{\Lambda} \boldsymbol{\Phi} \boldsymbol{\Omega}^{*} \\ \sqrt{\alpha} \boldsymbol{I}_{N_{t}} \end{array}\right)=\boldsymbol{Q} \boldsymbol{R} \boldsymbol{P}^{*}\) (3)
where Q is (𝑁𝑟 + 𝑁𝑡) × 𝑁𝑡 , P is 𝑁𝑡 × 𝑁𝑡 unitary matrices. The 𝑁𝑡 × 𝑁𝑡 upper triangular matrix R with equal singular values \(\bar{r}=\bar{\sigma}=\left(\prod_{k=1}^{N_{t}} \sqrt{\sigma_{k}^{2}+\alpha}\right)^{1 / N_{t}}\) is achieved by symmetric permutation and orthogonalization using the a series of 2 × 2 Givens rotation process. We note that Ω is a multiplication of power loading unitary matrices derived from each Givens rotation process as . Define \(\mathbf{Q}=\left[\widetilde{\mathbf{Q}}_{u} \widetilde{\mathbf{Q}}_{l}\right]^{T}\),where 𝑁𝑟 × 𝑁𝑡 matrix \(\widetilde{\mathbf{Q}}_{u}\) is the first Nr rows of Q and 𝑁𝑡 × 𝑁𝑡 matrix \(\widetilde{\mathbf{Q}}_{l}\) is the last Nt rows of Q. Applying the transmit precoder \(\mathbf{F}=\mathbf{V} \mathbf{\Phi} \mathbf{\Omega}\) and the decoder filter \(\widetilde{\mathbf{Q}}^{*}_{u}\), the output signal can be expressed as
\(\tilde{\mathbf{y}}=\widetilde{\mathbf{Q}}_{u}^{*} \mathbf{H V} \mathbf{\Phi} \mathbf{P} \mathbf{x}+\widetilde{\mathbf{Q}}_{u}^{*} \mathbf{n}=\mathbf{R} \mathbf{x}-\sqrt{\alpha} \widetilde{\mathbf{Q}}_{l}^{*} \mathbf{P} \mathbf{x}+\widetilde{\mathbf{Q}}_{u}^{*} \mathbf{n}\). (4)
The transmitted vector x can be detected by the MMSE with successive interference canceler or Tomlinson-Harashima Precoding (THP). GMD in [27] is a special case of UCD, with \(\boldsymbol{\Phi}=\boldsymbol{I}_{N_{t}}\) and 𝛼 = 0.
2.4 Block Diagonalizable UCD (BD-UCD)
Concatenation of UCD with block diagonalizable matrix are suggested in [31-33]. The main purpose of [31,32] is to make each user have equal rates for its data streams, by using block decomposed GMD and UCD. They allow each user to apply equal rate encoding on its own subchannels. The main purpose of [31,32], however, is to construct a decoder filter that has a block-diagonal form. To nullify multiuser interference in [31,32], a universal postcoder of size 𝑈×𝑁𝑡 should be given to all the users. The method to eliminate multiuser interference at the receive side is nothing but ZF decoder. It differs from the BD in 2.2, which is based on user allocation to the null space of the previously chosen user. In [33], block diagonalized UCD is used only for single user point-to-point channel. It simply pairs two eigenvalues with smaller and larger one by stream ordering and user ML decoding per each pair.
3. A New Two Step Precoding Algorithm
In this section, a new linear transmit-receive joint optimization method is presented. We propose a two-step precoding method based on equal channel gain decomposition with block wise diagonalization to simplify decoding and generating feedback. Sum rate capacity analysis, diversity analysis and comparison with various closed-loop MIMO schemes are also introduced in this section, followed by numerical simulation results in Section 4.
3.1 Transmit Strategies
The goal of the precoder design is to find an appropriate unitary precoding matrix F and decoder combining matrix Wk to make the effective channel \(\overline{\mathbf{H}}\) block diagonal and equal channel gain within each user k, as
\(\mathbf{y}_{k}=\mathbf{W}_{k}^{*} \widetilde{\mathbf{H}}_{k} \mathbf{F} \mathbf{x}+\mathbf{W}_{k}^{*} \mathbf{n}_{k} \triangleq \overline{\mathbf{H}}_{k} \mathbf{x}+\overline{\mathbf{n}}_{k}\) (5)
where 𝐇̃𝑘 is the effective channel with zero multiuser interference. To enforce zero interference criterion between users, a nulling matrix at user k using \(\mathbf{F}_{B D}=\left[\widetilde{\mathbf{V}}_{1}^{0} \boldsymbol{V}_{1}^{1} \cdots \widetilde{\mathbf{V}}_{U}^{0} {V}_{U}^{1}\right]^{T}\) and postcoding matrix Uk can be obtained using block diagonalization described in Section 2.2. Then the effective channel for the kth user is
\(\widetilde{\mathbf{H}}_{k}=\mathbf{U}_{k} \mathbf{H}_{k} \mathbf{F}_{B D}\). (6)
The second precoder 𝑭 = [𝑭1 ⋯ 𝑭𝑈] makes the effective channel per user have equal diagonalization, without hurting the zero-interference criterion. It is straightforward that \(\mathbf{W}_{j}^{*} \widetilde{\mathbf{H}}_{j} \mathbf{F}_{k}=0 \text { yields } \boldsymbol{W}_{j}^{*} \boldsymbol{U}_{j} \boldsymbol{H}_{j} \boldsymbol{F}_{B D} \boldsymbol{F}_{k}=0\) when 𝑗 ≠ 𝑘. By combining (5) and (6), we have the input-output relations for the kth user as
\(\boldsymbol{y}_{k}=\overline{\boldsymbol{H}}_{k} \boldsymbol{x}+\overline{\boldsymbol{n}}_{k}=\boldsymbol{W}_{k}^{*} \boldsymbol{U}_{k} \boldsymbol{H}_{k} \boldsymbol{F}_{B D} \boldsymbol{F} \boldsymbol{x}+\overline{\boldsymbol{n}}_{k}\) (7)
The aggregated effective channel can be written as\(\overline{\boldsymbol{H}}=\left[\begin{array}{lll} \overline{\boldsymbol{H}}_{1}^{T} & \cdots & \overline{\boldsymbol{H}}_{U_{U}}^{T} \end{array}\right]^{T}\), where \(\overline{\boldsymbol{H}}_{k}\) is the kth user’s effective subchannel with equal channel gain with rank Lk.
To compute the second precoder and postcoder for the kth user, we use SVD of \(\widetilde{\mathbf{H}}_{k}=\widetilde{\mathbf{U}}_{k} \widetilde{\Lambda}_{k} \widetilde{\mathbf{V}}_{k}\) and define the precoder F as
\(\mathbf{F}_{k}=\widetilde{\mathbf{V}}_{k} \mathbf{\Phi}_{k} \mathbf{P}_{k}\) (8)
where 𝜱𝑘 is the k-th user’s power loading matrix derived by the waterfilling solution of \(\widetilde{\mathbf{H}}_{k}\) [30] and Pk is the precoding matrix explained in Section 2.3. To prevent error propagation and change of second order channel statistics, we assume Pk be a unitary matrix. Let us define 𝐿𝑘 × 𝐿𝑘 diagonal matrix \(\Sigma_{k} \triangleq \widetilde{\boldsymbol{\Lambda}}_{k} \boldsymbol{\Phi}_{k}=\operatorname{diag}\left\{\sigma_{1}, \sigma_{2}, \cdots, \sigma_{L_{k}}\right\}\). The number of possible streams for the kth user Lk is upper bounded by the dimension of the matrix \(\widetilde{\mathbf{H}}_{k}\), hence 𝐿𝑘 ≤ 𝑁𝑟,𝑘. 𝐿𝑘 is also lower bounded by the minimum numbers of streams for the kth user to have at least one stream, 𝐿𝑘 > 𝑟𝑎𝑛𝑘(𝐇𝑘) + 𝑛𝑢𝑙𝑙𝑖𝑡𝑦(\(\widetilde{\mathbf{H}}_{k}\)) − 𝑁𝑡 . To have more degree of freedom and consequently more diversity gain at the decoder, larger Lk with larger eigenvalues is beneficial. Otherwise, if all of U users have just one stream, block diagonalization simply become an inverse of the aggregated channel matrix, resulting in large transmit power consumption due to the smallest eigenvalue [34]. In such a case, we do not achieve any receive diversity. The sum of all streams for all the users is equal to \(\sum_{k=1}^{U} L_{k}=N_{t}\).
As the diagonal matrix Σk is obtained by \(\widetilde{\boldsymbol{\Lambda}}_{k} \boldsymbol{\Phi}_{k}\), the equal gain decomposition described in Section 2.3 can be applied as
\(\overline{\boldsymbol{H}}_{k}=\left[\begin{array}{c} \sum_{k} \\ \sqrt{\alpha} \boldsymbol{I}_{N} \end{array}\right]=\boldsymbol{Q}_{k} \boldsymbol{R}_{k} \boldsymbol{P}_{k}^{*}, \quad 1 \leq k \leq U\). (9)
Qk is an unitary matrix and the upper triangular matrix Rk has the equal eigenvalues \(\bar{r}=\bar{\sigma}=\left(\Pi_{k=1}^{N_{t}} \sqrt{\sigma_{k}^{2}+\alpha}\right)^{1 / N_{t}}\). The power loading matrix P is derived through UCD method, and the first Lk columns of Qk is used for decoder postcoder of user k. This process is clearly capacity lossless from the following Weyl’s multiplicative majorization conditions.
Theorem 1: Let A be a matrix with rank n. Let the ordered eigenvalues and singular values of A be \(\lambda_{1}, \cdots, \lambda_{n}\) and \(\sigma_{1}, \cdots, \sigma_{n}\), respectively. Then for 1 ≤ k ≤ n,
\(\left|\lambda_{1}\right| \ldots\left|\lambda_{k}\right| \leq \sigma_{1} \ldots \sigma_{k}, \quad\left|\lambda_{1}\right| \ldots\left|\lambda_{n}\right| \leq \sigma_{1} \ldots \sigma_{n}\).
From Theorem 1, we see that the diagonal element of R in (9) is the geometric mean of singular value of H [27]. The permutation matrix for the kth user Pk is an orthonormal matrix, and it can be achieved by QR decomposition. By applying Givens 2 × 2 rotation process up to Lk, we obtain a complete set of the precoding matrix Qk, Rk and Pk in (9). One important point of Qk is that it starts with \(\widetilde{\mathbf{V}}_{k}^{0} \mathbf{V}_{k}^{1}\), not 𝑽𝑘, to preserve zero interference criterion. The first Lk columns of Qk is the decoder postcoder vector of user k, Wk.
After applying block-wise equal channel gain decomposition for every user, we can write Nt × Nt power loading matrix 𝐏 as 𝐏 ≜ 𝑑𝑖𝑎𝑔{𝐏1,𝐏2,⋯,𝐏𝑈}, and consequently have the precoder \(\mathbf{F}=\widetilde{\mathbf{V}} \Phi \mathbf{P}\). Given the precoder F, the kth user’s effective channel \(\overline{\mathbf{H}}_{k}=\mathbf{W}_{k}^{*} \widetilde{\mathbf{H}}_{k} \mathbf{F}\) is obtained. The kth user’s decoder combining matrix \(\mathbf{W}_{k}^{*}\) is given to each decoder as a decoder postcoder. Then the input-output relation of our two-step precoding is finally given by
\(\mathbf{y}_{k}=\mathbf{W}_{k}^{*} \widetilde{\mathbf{H}}_{k} \overline{\mathbf{V}} \mathbf{\Phi} \mathbf{P} \mathbf{x}+\mathbf{W}_{k}^{*} \mathbf{n}=\sum_{k} \mathbf{P}_{k} \mathbf{x}_{L_{x}}+\widetilde{\mathbf{n}}\) (10)
The joint precoder and postcoder design in (10) shows that we may generate an effective channel with equal channel gain without causing interference between users. Our next step is about sum rate capacity and diversity analysis, which are two measurements for the throughput and reliability performance in the MIMO channel.
3.2 Sum Rate Capacity Analysis
In this chapter, the sum rate capacity of our proposed algorithm is analyzed. Recall that the received signal for the kth user in (10) is
\(\mathbf{y}_{k}=\sum_{k} \mathbf{P}_{k} \mathbf{x}_{L_{x}}+\widetilde{\mathbf{n}}_{k}\). (11)
If we let Ux be input covariance matrix and Uk be noise covariance matrix for kth user’s effective noise \(\widetilde{\mathbf{n}}_{k}\), then the mutual information for kth user can be written as
\(l_{k}=\log _{2} \operatorname{det}\left(\boldsymbol{I}_{k}+\sum_{k} \boldsymbol{P}_{k} \boldsymbol{U}_{x}\left(\sum_{k} \boldsymbol{P}_{k}\right)^{*} \boldsymbol{U}_{k}^{-1}\right)\). (12)
Because we maintain zero-interference criterion, the aggregate mutual information for U users is equal to the sum of independent mutual information per user. Also noise covariance matrix Uk is the same to that of AWGN channel because multiplication of unitary postcoder W∗ at the decoder does not change the second order statistics of the noise. Therefore, the total sum rate with U users is given by
\(l_{\{1, \cdots, U\}}=\sum_{k=1}^{U} l_{k}=\log _{2} \operatorname{det}\left(\mathbf{I}_{N_{t} \times N_{t}}+\frac{\Sigma \mathbf{P U}_{x}(\Sigma \mathbf{P})^{*}}{\sigma_{n}^{2}}\right)\). (13)
The maximum sum rate is the mutual information maximizer with given input power constraint, 𝑡𝑟{𝑼𝑥}≤𝑃𝑇. The power allocation matrix P is normalized since it is a unitary matrix. Therefore, the achievable sum rate throughput is
\(C=\max _{\left\{\mathbf{U}_{x}: \operatorname{tr}\left\{\mathbf{U}_{x}\right\} \leq P_{T}\right\}} l_{\{1, \cdots, U\}}=\max _{\left\{\mathbf{U}_{x}: \operatorname{tr}\left\{\mathbf{U}_{x}\right\} \leq P_{T}\right\}} \sum_{k=1}^{U} l_{k}\stackrel{(a)}{=} \log _{2} \operatorname{det}\left(\mathbf{I}+\frac{\Sigma^{2}}{\sigma_{n}^{2}}\right)\). (14)
The equality (a) holds by the determinant identity det(I+AB)=det(I+BA) and capacity lossless property in Theorem 1. Choosing a proper power loading parameter to maximize the channel capacity is equivalent to the standard waterfilling process in [30] and therefore, the sum capacity of our algorithm is equivalent to the block diagonalization with waterfilling solution in [16]. When the number of streams per user goes one, we still achieve the same sum rate throughput, but diversity gain will be decreased as shown in the following section.
3.3 Diversity Gain Analysis
Diversity gain is a tool to measure error rate decay per SNR improvement. Unlike to the open loop MIMO literature where no CSIT is available at the transmitter [36], diversity analysis in closed-loop MIMO is closely related to the transmit power optimization since one can adjust transmit power or transmission rate to prevent error outage event. In case where the data rate is fixed, the diversity analysis is a good measure to see the error rate performance in closed loop literature.
Let us derive an expression for the probability of outage 𝑃𝑒(𝜌) when one is attempting to communicate at rate 𝑟𝑙𝑜𝑔(𝜌) bits/channel use, where r is known as spatial multiplexing gain and ρ is SNR. To make our analysis simple, we assume that the target SNR per user is identical and is simply achieved by proper power allocation. For fixed AWGN channel matrix H, we can formulate the error exponent using sum rate derived in the previous chapter as
𝑃𝑒(𝜌)=𝑃𝑟(𝑙𝑜𝑔 𝑑𝑒𝑡(𝐈+𝜌∑2𝐏2)<𝑟𝑙𝑜𝑔(𝜌)). (15)
From the results of previous section, we know that the sum capacity is achieved when appropriate power is allocated to achieve given SNR requirement. Therefore, we have the following inequality
\(P_{e}(\rho) \leq \operatorname{Pr}\left(\log \left(\prod_{k=1}^{U}\left(1+\rho \lambda_{k} \mathbf{P}_{k}\right)^{L_{k}}\right)<r \log (\rho)\right)\) (16)
where λk denotes the kth user’s eigenvalue. Because of the equal channel gain decomposition, user k has Lk identical eigenvalues λk. Then the outage corresponds to the event \(\log \left(\prod_{k=1}^{U}\left(1+\rho \lambda_{k}\right)^{L_{k}}\right)<r \log (\rho)\) is equivalent to
\(\prod_{k=1}^{U}\left(1+\rho \lambda_{k}\right)^{L_{k}}<\rho^{r}\) (17)
In (17), the U independent eigenvalues \(\lambda_{1}, \cdots, \lambda_{U}\) corresponding to each user can be written as the product of independent Chi-squared distribution χ, each with 2(Nt− Nr,k+1) degrees of freedom, under the high SNR assumption [34,37]. By applying 𝜆𝑖=𝜌−𝛼𝑖 to (17) and high SNR approximation yield the following outrage probability
\(P_{e}(r \log \rho)=\operatorname{Pr}\left(\prod_{i=1}^{U}\left(1+\rho \lambda_{i}\right)^{L_{i}}<\rho^{r}\right)=\operatorname{Pr}\left(\sum_{i}^{U}\left(\left(1-\alpha_{i}\right) L_{i}\right)^{+}<r\right)\). (18)
Our next interest is to formulate of distribution function of the random vector \(\lambda=[\lambda_{1}, \cdots, \lambda_{U}]\). According to [36] Lemma 3] and independent Chi-squared distribution property, the joint pdf of random vector λ is
\(P(\lambda)=A \prod_{i=1}^{U} \lambda_{i}^{L_{i}\left(N_{t}-N_{r, j}\right)} e^{-\sum_{i=1}^{U} \lambda_{i}^{i_{i}}} \prod_{j>i}\left(\lambda_{i}-\lambda_{j}\right)^{2}\). (19)
By applying (18) into (19) with the outage event \(\varepsilon=\left\{\lambda: \sum_{\alpha}\left(\left(1-\lambda_{i}\right) L_{i}\right)^{+}<r\right\}\) yields
\(P_{e}(r \log \rho)=\int_{\varepsilon} p(\alpha) d \alpha=\int_{\varepsilon} \prod_{i=1}^{U} \rho^{-\left(2 i-1-\mid N_{t}-N_{r, i}\right) L_{i} \alpha_{i}} d \alpha\). (20)
From the above equation, it is enough to determine the dominating exponent for the diversity gain. The diversity gain of algorithm can be written as
\(d\left(r^{\prime}\right)=-\lim _{\rho \rightarrow \infty} \frac{\log P_{e}(\rho)}{\log \rho}\) (21)
\(=-\lim _{\rho \rightarrow \infty} \frac{\log \int_{\varepsilon} \prod_{i=1}^{U} \rho^{-\left(2 i-1-\mid N_{t}-N_{R, i}\right) L_{i} \alpha_{i}} d \alpha}{\log \rho}\) (22)
\(=\inf _{\alpha_{i}} \sum_{i=1}^{U} L_{i} \alpha_{i}\left(2 i-1+\left|N_{r, k}-N_{t}\right|\right)\) (23)
Equation (23) shows that when trying to communicate with rate \(r^{\prime}, \text { i.e., } r=r^{\prime}\), the maximum diversity order (inf is achieved when \(\alpha_{1}=\cdots=\alpha_{U-r^{\prime}}=1\) [36]) is
\(d\left(r^{\prime}\right)=\sum_{i=1}^{U-r^{\prime}} L_{i}\left(2 i-1+\left|N_{r, i}-N_{t}\right|\right)\) (24)
\(=\sum_{i=1}^{U-r^{\prime}} L_{i}(2 i-1)+\left(U-r^{\prime}\right)\left(\left|N_{r, i}-N_{t}\right|\right)\). (25)
The results in (25) state that when there is Nt single antenna users, full diversity order \(\left(N_{t}\right.\left.-r^{\prime}\right)\left(N_{r}-r^{\prime}\right)\) is achieved when attempting to communicate at rate \(r^{\prime}\left(r^{\prime} \leq \min \left\{N_{t}, N_{r}\right\}\right)\). If four transmit antennas and two users with two decoder antennas are considered, the maximum achievable rate is \(2\left(2-r^{\prime}\right)\left(4-r^{\prime}\right)\), when attempting to communicate at 𝑟′. This shows that our scheme achieves full receive diversity gain while the maximum multiplexing gain per user cannot exceed R ≤ min (Nt, Nr,k). The numerical results in Section 4 show that our algorithm achieves full diversity gain in given data rate.
3.4 Receiver and Feedback Strategies
In this section, the decoding methods and feedback strategies for the proposed algorithm is analyzed. We assume that the postcoding matrix per user derived in Section 3.1 is given to each decoder by a proper signaling. The result in (10) shows that the composite effective channel matrix is block diagonalizable per user. The kth user has Nr, k × Nt channel with real channel elements, without multiuser interference caused by other users. Since the effective channel is real, we may decode real part and imaginary part of the transmit symbols independently. Hence, the following per symbol ML metric can be applied as
\(\hat{s}_{k}=\arg \min _{s_{k} \in C}\left\|\boldsymbol{y}_{k}-\sum_{k} \boldsymbol{P}_{k} s_{k}\right\|^{2}\) (26)
for real and imaginary part of 𝑠̂𝑘, independently. The decoding complexity is also reduced by the independent decoding, while achieving full receive diversity as analyzed in Section 3.3.
Another advantage of the proposed algorithm is on the feedback signaling for link adaptation at the transmitter. Unlike to the linear decoder such as ZF and MMSE, SINR after ML decoding cannot be explicitly calculated, to our best knowledge. Even though there is an exhaustive search algorithm based on pairwise error probability (PEP) [38], it is too complex and not analogous to multiuser MIMO with link adaptation. As the modulation and coding selection (MCS) per streams is not equal within each user, different modulation order results in different PEP in the multiuser MIMO link. But per stream link adaptation is not required in our algorithm, since the effective channel gain is already equalized for all streams within a user. We may apply link adaptation per user, not per stream, using per user SINR.
The problem to generate effective SINR become simpler since the eigenvalues within each user are equal. The union bound for coded binary transmission and ML decoding can be written as
\(P_{e}(\gamma) \leq \sum_{d=d_{\min }}^{\infty} \alpha_{d} P_{2}(d, \gamma)\) (27)
where γ is effective SINR that we want to derive, αd is the number of codewords with hamming weight d, and P2(d, γ) is PEP with respect to hamming distance d and effective SINR γ. Then the PEP upper bound is obtained using Chernoff bound
\(P_{e}(\gamma) \leq \sum_{d=d_{\min }}^{\infty} \alpha_{d} P_{2}(d, \gamma)<\sum_{d=d_{\min }}^{\infty} \alpha_{d}\left[P_{2}(1, \gamma)\right]^{d} \triangleq P_{e, \text { Chernoff }}(\gamma)\). (28)
From (28), the Chernoff-bounded error probability only depends on the weight distribution of the code and the Chernoff-bounded symbol error probability. If different SINR values are used, (28) may be solved as the Union Chernoff bound for multiple-state channel [39], consequently requires a weighting parameter to adjust the effective SINR to a specific combination of modulation and coding rate. The proposed algorithm makes the work simpler, since multiple identical eigenvalues may reduce the computation from multi-state Chernoff bound to single-state bound. Thus, we have a simple effective SINR for kth user as
\(\gamma_{k, \text { eff }}=-\beta_{M L} \ln \left(\frac{1}{L_{k}} \sum_{i=1}^{L_{k}} \exp \left(-\frac{\gamma_{k, i}}{\beta_{M L}}\right)\right) \approx-\beta_{M L} \ln \left(\exp \left(-\frac{\bar{\gamma}_{k}}{\beta_{M L}}\right)\right)=\bar{\gamma}_{k}\) (29)
where 𝛾̄𝑘 is the received SINR at user k. From (29), we can use effective SINR measured at the decoder before applying ML decoding as an indicator for the link adaptation to the transmitter. If limited feedback is applied as in [14], our strategy is simply feedback the SINR value measured at the decoder, or as a form of MCS indicator by quantizing the channel capacity as
\(Q\left(C_{k}\right)=\left\lfloor C_{k}\right\rfloor \leq L_{k} \log _{2}\left(1+\bar{\gamma}_{k}\right)\) (30)
where ⌊⋅⌋ indicates lower bound in quantization.
Another advantage of the proposed algorithm is on the system level performance evaluation. The effective SINR is used to map the instantaneous channel state information into a specific MCS table satisfying target error rate. It has been a hard problem to obtain per stream effective SINR when ML decoding is used, resulting in performance leakage due to the effective SINR mismatch. Using the proposed methods, each user can generate one representative effective SINR without mismatch. The proposed strategy is the simplest form of the exponential effective signal to interference ratio mapping (EESM) for ML decoder so far, to our best knowledge.
4. Simulation Results
In this section, numerical simulation results for the candidate algorithm are presented. We first compare our algorithm with previously proposed MIMO schemes in Table 1. To illustrate sum rate performance, we present sum rate throughput in various SNR regions in Fig. 2. Error rate performance is presented in Fig. 3 and Fig. 4 to see the diversity gain performance of the proposed algorithm.
Table 1. Comparison of various single-user and multiuser MIMO algorithms
Fig. 2. Sum rate performance comparison between DPC, BD, ZF (channel inversion) and the proposed method for (2, 2, 4) case – 4 transmit antennas, 2 users with 2 receive antennas
Fig. 3. Block error rate performance comparison between ZF (channel inversion), MMSE, SVD without power optimization, BD, BD with ML decoder, BD-GMD with ML decoder and the proposed algorithm in (2, 2, 4) configuration – 4 transmit antennas, 2 users with 2 receive antennas
Fig. 4. Block error rate performance comparison in two antenna configurations (3, 2, 6) and (2, 3, 6) – 6 transmit antennas and 3 users with 2 receive antennas vs 6 transmit antennas and 2 users with 3 receive antennas
Fig. 2 shows the sum rate performance of the proposed algorithm. For four transmit antenna and two users with two receive antenna each, our algorithm shows a good performance compared to conventional ZF precoding methods. Also, the achievable rate region of our algorithm is almost the same to the BD achievable region. This is because our algorithm perfectly cancels the multiuser interference, while maintaining efficient power loading per user. We assume that the number of aggregated receive antennas for all users are equal to the transmit antenna, so no scheduling gain is affected to the performance. If we have more users and thus can be benefited by the multiuser diversity gain (or equivalently scheduling gain), the performance gap between DPC and our proposed algorithm may be shortened.
Fig. 3 and Fig. 4 show the block error rate performance of the proposed algorithm. We assume four transmit antenna and two users with two receive antennas in Fig. 3 and compare the proposed algorithm with ZF, MMSE, SVD, BD, BD with ML and BD-GMD with ML decoder. Conventional BD achieves diversity one as ZF and MMSE, but our proposed algorithm achieves diversity order four. It is better than BD-GMD with ML decoder and BD with ML decoder because UCD achieves full diversity as proven in Section 3.3.
Fig. 4 shows the diversity performance in different antennas configurations. We assume six transmit antennas and two different receive antenna configurations: 3 users with 2 decoder antennas and 2 users with 3 receive antennas. As we analyzed in Section 3.1, configurations with more receive antenna is advantageous in achieving more diversity gain. The simulation results in Fig. 4 reinforce the analysis in Section 3.1 and 3.3.
We compare the proposed algorithm with some previously proposed algorithms in Table 1. The reference algorithms are SVD, OBF [14], BD [16] and UCD [28], which are well known closed-loop MIMO algorithms. We investigate the difference in transmit precoder design, receive filter, effective channel shape, multiuser diversity, decoding method if there is no feed-forward information, achievable sum rate, diversity gain, SINR feedback requirement and finally computational complexity. Table 1 shows that the proposed method has a good tradeoff between complexity and performance. It requires slightly more complexity than BD and UCD but achieves a good sum rate performance and optimal diversity gain. SINR feedback reduction even with ML decoding using single-stage Chernoff bound is another merit of the proposed algorithm.
5. Conclusion
In this paper, a new multiuser MIMO precoder and postcoder design method is proposed. The proposed method uses two step precoding by combining equal rate decomposition under with multiuser interference criterion. We proved that the proposed algorithm can be easily combined with a symbol decodable MLD and can generate simple but accurate effective SINR. We showed that the algorithm achieved good sum rate performance and full diversity gain. Comparison with other algorithms in both SU and MU MIMO was also presented.
We note that the scheduling algorithm of our proposed algorithm is still under investigation. Because diversity order per user is upper bounded by the number of receive antenna per user, the scheduling strategy in terms of maximizing diversity gain is to schedule users with more receive antennas. But when the number of users is high, it is easy to schedule users with single streams, to meet the zero-interference criterion. A qualitative analysis between precoding and beamforming based on the number of users in the scheduling queue will be another interesting research area.
참고문헌
- D. Gesbert, M. Kountouris, R. W. Heath Jr., C. B. Chae, and T. Salzer, "From single user to multiuser communications: Shifting the MIMO paradigm," IEEE Signal Processing Magazine, vol. 24, no. 5, pp. 36-46, 2007. https://doi.org/10.1109/MSP.2007.904815
- G. Caire and S. Shamai, "On the achievable throughput of a multi-antenna Gaussian broadcast channel," IEEE Transactions on Information Theory, vol. 49, no. 7, pp. 1691-1706, July 2003. https://doi.org/10.1109/TIT.2003.813523
- H. Weingarten, Y. Steinberg, and S. Shamai, "The capacity region of the Gaussian multiple-input multiple-output broadcast channel," IEEE Transactions on Information Theory, vol. 52, no. 9, pp. 3936-3964, Sep. 2006. https://doi.org/10.1109/TIT.2006.880064
- A. B. Nogueras, A. Davoodi, P. Kerret, D. Gesbert, N. Gresset, and S. Jafar, "DoF Region of the Decentralized MIMO Broadcast Channel-How many informed antennas do we need?," in Proc. of 2020 IEEE International Symposium on Information Theory, pp. 1468-1473, 2020.
- K. Kim and J. Lee, "Deep CNN based pilot allocation scheme in massive MIMO systems," KSII Transactions on Internet and Information Systems, vol. 14, no. 10, pp. 4214-4230, 2020. https://doi.org/10.3837/tiis.2020.10.016
- T. Yoo and A. Goldsmith, "On the optimality of multiantenna broadcast scheduling using zero-forcing beamforming," IEEE Journal on Selected Areas in Communications, vol. 24, no. 3, pp. 528-541, Mar. 2006. https://doi.org/10.1109/JSAC.2005.862421
- W. Xu and M. Lin, "Downlink Capacity Analysis of Distributed Antenna Systems with Imperfect Channel State Information," KSII Transactions on Internet and Information Systems, vol. 11, no. 1, pp. 253-271, 2017. https://doi.org/10.3837/tiis.2017.01.013
- N. Jindal and A. Goldsmith, "Dirty-paper coding versus tdma for MIMO broadcast channels," IEEE Transactions on Information Theory, vol. 51, no. 5, pp. 1783-1794, May 2005. https://doi.org/10.1109/TIT.2005.846425
- J. Cao, S. Yu, J. Yang, Y. Zhang, and S. Zhao, "Massive MIMO two-hop relay systems over Rician fading channels," KSII Transactions on Internet and Information Systems, vol. 13, no. 11, pp. 5410-5426, 2019. https://doi.org/10.3837/tiis.2019.11.008
- S. Vishwanath, N. Jindal, and A. Goldsmith, "Duality, achievable rates, and sum capacity of Gaussian MIMO broadcast channels," IEEE Transactions on Information Theory, vol. 49, no. 10, pp. 2658-2668, Aug. 2003. https://doi.org/10.1109/TIT.2003.817421
- X. Li, L. Li, and L. Xie, "Achievable sum rate analysis of ZF receivers in 3D MIMO systems," KSII Transactions on Internet and Information Systems, vol. 8, no. 4, pp. 1368-1389, 2014. https://doi.org/10.3837/tiis.2014.04.012
- X. Yunjia, H. Kaizhi, H. Xin, Z. Yi, C. Yajun, and J. Wenyu, "CoMP Transmission for Safeguarding Dense Heterogeneous Networks with Imperfect CSI," KSII Transactions on Internet and Information Systems, vol. 13, no. 1, pp. 110-132, 2019. https://doi.org/10.3837/tiis.2019.01.007
- D. Gore, R.W. Heath, and A. Paulraj, "Transmit selection in spatial multiplexing systems," IEEE Communications Letters, vol. 6, no. 11, pp. 491-493, Nov. 2002. https://doi.org/10.1109/LCOMM.2002.805517
- K. Huang, R. W. Heath, Jr., and J. G. Andrews, "Performance of orthogonal beamforming for SDMA with limited feedback," IEEE Transations on Vehicular Technology, vol. 58, no. 1, pp. 152-164, 2009. https://doi.org/10.1109/TVT.2008.925003
- C. Yuqing, L. Xiaoyan, S. Xixia, and P. Su, "Resource allocation in multi-user MIMO-OFDM systems with double-objective optimization," KSII Transactions on Internet and Information Systems, vol. 12, no. 5, pp. 2063-2081, 2018. https://doi.org/10.3837/tiis.2018.05.009
- Q. Spencer, A. Swindlehurst, and M. Haardt, "Zero-forcing methods for downlink spatial multiplexing in multiuser MIMO channels," IEEE Transactions on Signal Processing., vol. 52, no. 2, pp. 461-471, Feb. 2004. https://doi.org/10.1109/TSP.2003.821107
- L. U. Choi and R. D. Murch, "A transmit preprocessing technique for multiuser MIMO systems using a decomposition approach," IEEE Transactions on Wireless Communications, vol. 3, no. 1, pp. 20-24, Jan. 2004. https://doi.org/10.1109/TWC.2003.821148
- C. B. Peel, B. M. Hochwald, and A. L. Swindlehurst, "A vector-perturbation technique for near capacity multiantenna multiuser communication - part I: channel inversion and regularization," IEEE Transactions on Communications, vol. 53, no. 1, pp. 195-202, Jan. 2005. https://doi.org/10.1109/TCOMM.2004.840638
- C. Windpassinger, R. F. H. Fischer, T. Vencel, and J. B. Huber, "Precoding in multiantenna and multiuser communications," IEEE Transactions on Wireless Communications, vol. 3, no. 4, pp. 1305-1316, July 2004. https://doi.org/10.1109/TWC.2004.830852
- B. M. Hochwald, C. B. Peel, and A. L. Swindlehurst, "A vector-perturbation technique for near capacity multiantenna multiuser communication - part II: perturbation," IEEE Transactions on Communications, vol. 53, no. 3, pp. 537-544, Mar. 2005. https://doi.org/10.1109/TCOMM.2004.841997
- U. Fincke and M. Pohst, "Improved methods for calculating vectors of short length in a lattice, including a complexity analysis," Mathematics of Computation, vol. 44, pp. 463-471, 1985. https://doi.org/10.1090/S0025-5718-1985-0777278-8
- N. Yoo, J. Back, H. Choi, and K. Lee, "ML Symbol Detection for MIMO Systems in the Presence of Channel Estimation Errors," KSII Transactions on Internet and Information Systems, vol. 10, no. 11, pp. 5305-5321, 2016. https://doi.org/10.3837/tiis.2016.11.006
- B. Hassibi and H. Vikalo, "On the sphere-decoding algorithm I: Expected complexity," IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 2806-2818, Aug. 2005. https://doi.org/10.1109/TSP.2005.850352
- K. Higuchi, H. Kawai, N. Maeda, and M. Sawahashi, "Adaptive selection of surviving symbol replica candidates based on maximum reliability in QRM-MLD for OFCDM MIMO multiplexing," in Proc. of IEEE Global Telecommunications conference, vol. 4, pp. 2480-2486, vol. 4, 2004.
- J. Kim, J. Kang, and W. Y. Yang, "Recursive decoding for OFDM systems with multiple transmit antennas," IEICE Transactions on Communications, vol. 88, no. 6, pp. 2664-2667, June 2005.
- J. Jalden, L. Barbero, B. Ottersten, and J. Thompson, "Full diversity detection in MIMO systems with a fixed-complexity sphere decoder," in Proc. of IEEE International Conference Acoustics, Speech and Signal Processing, vol. 3, pp. III-49-III-52, Apr. 2007.
- Y. Jiang, J. Li, and W. Hager, "Joint transceiver design for MIMO communications using geometric mean decomposition," IEEE Transactions on Signal Processing, vol. 53, no. 10, pp. 3791-3803, Oct. 2005. https://doi.org/10.1109/TSP.2005.855398
- Y. Jiang, J. Li, and W. Hager, "Uniform channel decomposition for MIMO communications," IEEE Transactions on Signal Processing, vol. 53, no. 11, pp. 4283-4294, Nov. 2005. https://doi.org/10.1109/TSP.2005.857052
- J. Lee and N. Jindal, "High SNR analysis for MIMO broadcast channels: Dirty paper coding versus linear precoding," IEEE Transactions on Information Theory, vol. 53, no. 12, pp. 4787-4792, Dec. 2007. https://doi.org/10.1109/TIT.2007.909094
- E. Telatar, "Capacity of multi-antenna Gaussian channels," European Transactions Telecommunications, vol. 10, pp. 585-595, Nov. 1999. https://doi.org/10.1002/ett.4460100604
- S. Lin, W. W. L. Ho, and Y.C. Liang, "Block-diagonal geometric mean decomposition (BD-GMD) for multiuser MIMO broadcast channels," IEEE Transactions on Wireless Communications, vol. 7, no, 7 pp. 2778-2789, 2009. https://doi.org/10.1109/TWC.2008.070171
- S. Lin, W. W. L. Ho, and Y. C. Liang, "MIMO broadcast communications using block-diagonal uniform channel decomposition (BD-UCD)," in Proc. of 2006 IEEE 17th International Symposium on Personal, indoor and Mobile Radio Communications, pp. 1-5, Sep. 2006.
- K. J. Lee, B. M. Hochwald, and I. Lee, "Blockwise uniform channel decomposition for MIMO systems," in Proc. of Globecom 2007, IEEE Global Telecommunications Conference, pp. 3154-3158, Nov. 2007.
- A. Edelman, "Eigenvalues and condition number of random matrices," Ph.D. dissertation, Dept. Mathematics, MIT, Cambridge, MA, 1989.
- H. Weyl, "Inequalities between the two kinds of eigenvalues of a linear transformation," National Academy of Sciences of the United States of America, vol. 35, pp. 408-411, 1949. https://doi.org/10.1073/pnas.35.7.408
- L. Zheng and D. N. C. Tse, "Diversity and multiplexing: a fundamental tradeoff in multiple-antenna channels," IEEE Transactions on Information Theory, vol. 49, no. 5, pp. 1073-1096, May 2003. https://doi.org/10.1109/TIT.2003.810646
- T. Anderson, "An Introduction to Multivariate Statistical Analysis," 3rd ed., Wiley, 1958.
- H. Liu, L. Cai, H. Yang, and D. Li, "EESM based link error prediction for adaptive MIMO-OFDM system," in Proc. of 2007 IEEE 66th Vehicular Technology Conference, pp. 559-563, Apr. 2007.
- Ericsson, "System level evaluation of OFDM - further considerations," TSG-RAN WG1 #35, R1-031303, Nov. 2003.