Dynamic rate-adaptive MIMO mode switching between spatial multiplexing and diversity

In this article, we propose a dynamic multiple-input multiple-output (MIMO) mode switching scheme between spatial multiplexing and diversity modes, which also includes adaptive modulation. At each transmission, we select the modulation level and the MIMO mode that maximize the spectral efficiency while satisfying a given target bit error rate. The dynamic MIMO mode scheme considers instantaneous spectral efficiency whereas the conventional static scheme considers only the average SNR. As for adaptive modulation, a new method is proposed to compute the SNR thresholds for adaptive modulation in each MIMO mode, and it can avoid the computational difficulty of the conventional Lagrangian (optimal) method at high average SNR. To deal with the case where the rates of the two MIMO modes are the same, we also propose a new measure based on the BER exponent, which has lower computational complexity than a conventional measure. Numerical results show that the proposed dynamic mode switching improves over the conventional static mode switching in terms of average spectral efficiency.


Introduction
Today's wireless communication systems demand high data rate and spectral efficiency with increased reliability. Multiple-input multiple-output (MIMO) systems have been popular techniques to achieve these goals because increased data rate is possible through spatial multiplexing scheme [1] or improved diversity order is possible through transmit diversity scheme (e.g., space-time block code, STBC) [2]. Other ways are link adaptation techniques, where transmission parameters such as modulation and coding are dynamically adapted to the varying channel condition [3]. A typical link adaptation technique is adaptive modulation in which an adequate modulation level is selected by means of the current signal-to-noise ratio (SNR).
Recently, adaptive modulation schemes in conjunction with MIMO techniques have been investigated [4][5][6][7][8][9][10][11]. The prior study in the literature mainly tried to maximize the average spectral efficiency (ASE) for only one MIMO mode, either spatial multiplexing [8] or transmit diversity [5][6][7]11]. In [12], the mode switching between diversity and multiplexing was first proposed. But the authors focused on the situation where both MIMO *Correspondence: junglee@snu.ac.kr School of Electrical Engineering and Computer Science, Seoul National University, Seoul, Korea modes have equal spectral efficiency without considering adaptive modulation, so that they showed the result that spatial multiplexing is preferred in low SNR region. The mode switching scheme combined with adaptive modulation was proposed in [9,10], but the analysis was focused on the static mode switching which depends only on the average SNR.
In this article, we propose a dynamic MIMO mode switching scheme which considers instantaneous channel condition in conjunction with rate adaptation. Although the adaptive modulation part is based on the existing methods [5,[7][8][9]11], we compare the performance of the existing techniques, and also propose a sub-optimal method to obtain the SNR thresholds for the average BER constraint. Its complexity is lower than that of the optimal method using a Lagrange multiplier, but the performance degradation is negligible. The proposed mode switching scheme is based on the instantaneous spectral efficiency (ISE). In case the ISE's of the two modes are equal, an additional rule is necessary for mode selection. Although the Demmel condition number proposed in [12] can be a choice, we propose a new method which has lower complexity than the Demmel condition number without performance loss.
This article is organized as follows. In the section of System overview, we outline the system and the channel model as well as the structure of the considered http://jwcn.eurasipjournals.com/content/2012/1/238 MIMO mode switching system. The dynamic MIMO mode switching scheme combined with adaptive modulation is then proposed in the section of rate-adaptive MIMO mode switching. In the section of Simulation results, we compare the performance of the proposed algorithm with that of the existing methods. Since numerical methods are necessary in order to get the SNR thresholds for the average BER constraint, we shows detailed results in this section. Finally, conclusions are drawn in the last section.

System overview
We consider a MIMO system with M transmit antennas and N receive antennas. The block diagram of the proposed system is shown in Figure 1. The system consists of a transmitter with a switch between a multiplexing and a diversity modulator, a receiver unit with the corresponding pair of receivers, a modulation level and mode selector, and a low rate feedback path. At the receiver side, the modulation level and the MIMO mode are selected according to the current channel condition. The information about the selected modulation level and the MIMO mode is sent to the transmitter through the feedback path. The transmitter then switches the MIMO mode with the modulation level based on the feedback information.
Suppose that the N × M flat fading channel matrix H has i.i.d. complex Gaussian random entries. The (i, j)th entry [ H] i,j = h ij is distributed as CN (0, 1) a . The channel is assumed to be quasi-static (channel coefficients do not change during one time interval, and change independently in the next interval). The input-output relation for the MIMO channel is given by

Rate-adaptive MIMO mode switching
In this article, we propose a new rate-adaptive MIMO mode switching algorithm. The goal is to maximize the ASE while satisfying a given bit error ratio constraint. The proposed algorithm can be summarized by the following three steps.
1. Calculate the post-processing SNR in each MIMO mode. 2. Decide the modulation order in each MIMO mode. 3. Decide one MIMO mode based on a given selection rule.
For analysis, we consider a linear receiver for the spatial multiplexing mode, and orthogonal space-time block codes (OSTBC) for the diversity mode. In Step 2, we analyze several adaptive modulation techniques subject to an instantaneous BER constraint as well as an average one. In Step 3, we propose the mode selection rule based on the ISE as well as the rule which can be applied to the case when both of the two MIMO modes have the same data rate.

Post-processing SNR calculation (Step 1)
The post-processing SNR at the receiver is calculated separately in each MIMO mode with a given detection algorithm. At first, in the spatial multiplexing mode, the post-processing SNR of the mth (m = 1, 2, . . . , M) output data stream of the zero forcing (ZF) receiver, denoted as γ m,ZF , is given by ( [13], Eq. (7.43)) and the SNR of the minimum mean-square error (MMSE) receiver, denoted as γ m,MMSE , is given by ( [13], Eq. (7.49)) The post-processing SNR of the OSTBC system, denoted as γ OSTBC , is given by where ζ is the code rate of the OSTBC.

Decision of the modulation order (Step 2)
Using the post-processing SNR obtained in Step 1, we can choose an appropriate modulation order for each MIMO mode which enhances spectral efficiency without exceeding a given target BER at the receiver. Since adaptive modulation for one MIMO mode with a given BER constraint has been studied in [4][5][6][7][8]11], we take a similar approach of the literature. For analysis, we consider a discrete The receiver chooses the constellation M l whenever θ l ≤ γ < θ l+1 . If γ < θ 1 , data transmission is suspended for the corresponding channel since the respective BER constraint cannot be satisfied. Moreover, the maximum SNR threshold is set to infinity, i.e., θ L+1 = ∞.

SNR thresholds for instantaneous BER constraint
An easy way to set the switching thresholds θ l 's is to use the instantaneous BER (I-BER). In this approach, the BER of every reception has to be less than or equal to the target BER δ 0 . In order to meet the constraint, the BER for a QAM in AWGN channels can be used. Although the exact BER expressions for M-QAM are shown in [14], they are not easily inverted with respect to the SNR, so that a numerical method is necessary. Instead, in the adaptive modulation literature [5][6][7][8], an exponential function form is used, which is given by where a l = 0.2 and c l is a constellation specific constant defined as [4] c l = 6 5·2 l −4 for rectangular QAM (odd l) 3 2(2 l −1) for square QAM (even l) .
If we want a more accurate form than the above approximation, we can find the modulation specific constants a l and c l numerically using curve-fitting methods [11]. Table 1 shows those values of M-QAM's which are used in [11]. Inverting (5) with respect to γ , the switching threshold is determined as Although it is simple, I-BER approach keeps the instantaneous BER at all time instants below the target BER δ 0 . This is so conservative that the average BER (A-BER) is lower than δ 0 . In order to make the A-BER be equal to δ 0 , SNR thresholds should be lowered. Therefore, there is potential for improving the ASE by adjusting the switching threshold of each modulation.

SNR thresholds for average BER constraint
Generally, the ASE η for one channel use is given by where b l = log 2 M l is the number of bits corresponding to the lth modulation and p l is the probability that the postprocessing SNR falls into the lth bin, given by where f (γ ) is the probability density function (pdf ) of the post-processing SNR. The A-BER can be denoted as the average number of error bits N e,avg divided by the average number of transmitted bits N b,avg . It is observed that N b,avg = η by the definition in (8) and N e,avg is given by where P e (l) is the A-BER when the SNR falls into the lth bin, given by Since it has already been known that the pdf 's of the two MIMO modes have Gamma distributions, using the above formulas, the ASE and the A-BER can be obtained as closed forms as (26) and (30) in case of ZF spatial multiplexing system and (34) and (38) in case of OSTBC system, respectively. See Appendix for details. The pdf of γ m,MMSE can be well approximated to a Gamma distribution and to http://jwcn.eurasipjournals.com/content/2012/1/238 a generalized Gamma distribution [15]. Thus, though it is approximation, the analysis of the MMSE receiver can be done with the same procedure as that of the ZF receiver.
Optimal method In A-BER approach, the goal is to maximize the ASE under the constraint that the A-BER should be lower than or equal to δ 0 . Defining the set of adjustable switching thresholds as = θ l | l = 1, 2, . . . , L c , the optimization problem can be formulated as This problem can be solved with a Lagrange multiplier. Since P e,avg is denoted as changing the constraint P e,avg ≤ δ 0 into N e,avg ≤ δ 0 η for convenience, the Lagrangian of (12) is defined as Differentiating (14) with respect to θ l and equating to zero, the following relationship for l = 2, 3, . . . , L is obtained as According to the relationship, once θ 1 is chosen, all the other θ l 's are uniquely determined. Thus, the optimal SNR thresholds can be found numerically by adjusting θ 1 only.
Descending search method (sub-optimal) Intuitively, as the SNR range assigned to a high order modulation increases, the ASE increases while the BER also increases at the same time. Thus, in order to maximize the ASE, the SNR threshold θ l in descending order (l = L, L − 1, . . . , 1) has to be lowered as much as possible. Since the A-BER has to be kept below δ 0 , we take the constraint that P e (l) ≤ δ 0 , ∀l. The detailed algorithm is as follows: Find θ l such that P e (l) = δ 0 4: if θ l then 5: if P e (l) < δ 0 then 7: Exit for loop. 9: else 10: Switch off {M l }.

11:
θ l ← θ l+1 12: end if 13: end if 14: end for In the proposed algorithm, a certain modulation can be completely switched off because the solution that satisfies P e (l) = δ 0 may not exist at all or only higher order modulations may satisfy the BER constraint on the whole SNR region.
Although the A-BER approach has higher computational complexity than the I-BER approach, the thresholds can be calculated off-line, and the A-BER approach is still practical.

Dynamic MIMO mode switching (Step 3)
After the modulation order in each MIMO mode is chosen, the ISE can be calculated. The ISE of the spatial multiplexing system with ZF receiver, denoted as R ZF , can be written as where b m,ZF is the number of bits corresponding to the selected modulation on the mth subchannel. Likewise, the ISE of the OSTBC system, denoted as R OSTBC , can be expressed as where b OSTBC denotes the number of bits corresponding to the selected modulation.
If R ZF and R OSTBC are different from each other, the mode selection rule is as follows: if R ZF > R OSTBC , spatial multiplexing is chosen for the next transmission mode, and vice versa. In case of R ZF = R OSTBC , a general rule is to select a mode that gives lower BER, and the MIMO mode can be chosen based on the following two methods. One method is the Demmel condition number approach which was proposed in [12]. The Demmel condition number κ D is defined as where λ min (H) denotes the minimum singular value of H, and spatial multiplexing is preferred if where d min,ZF is the minimum Euclidean distance of the transmit constellation of the spatial multiplexing system with a ZF receiver, and d min,OSTBC is the minimum Euclidean distance of the OSTBC system. Another method is to use the I-BER of each MIMO mode, which can be measured by the exponent of the BER Equation (5). In other words, a MIMO mode which has lower I-BER than the other can be chosen. Assuming that the approximation of the I-BER in (5) and the constellation specific constants in [4] are used, the following rule can be derived. Spatial multiplexing is preferred if where γ min,ZF = min m γ m,ZF , which means that the worst stream of the spatial multiplexing mode is used to calculate the I-BER. Although the Demmel condition number can be used for mode switching, it has higher complexity than (20) because singular value decomposition is necessary to get κ D in (18), whereas (20) only uses the parameters which are already obtained from the previous steps. Therefore, the proposed measure is desirable in terms of computational complexity.

Simulation results
In simulation, a simple MIMO system is considered whose

Adaptive modulation in a MIMO mode SNR thresholds for the I-BER constraint
The thresholds are obtained in the following three ways.
In the first method, the thresholds are obtained from numerical search using the exact BER expressions of M-QAM in AWGN channels. In the second method, the thresholds are from (7) with constellation specific constants defined in . In the third method, the thresholds are also from (7) with the constants of a l = 0.2 and c l defined in (6). As shown in Table 2, the second thresholds are quite close to the first, so the first and the third thresholds are used for the I-BER simulation.

SNR thresholds for the A-BER constraint
It is observed that, in (30) and (38) of Appendix, the A-BER of each MIMO mode depends on the three parameters, M, N, and ρ. Thus, if M and N are fixed, the optimal SNR thresholds change only according to ρ. We obtained the optimal thresholds using a Lagrange multiplier, and obtained the sub-optimal thresholds using the descending search method described in the previous section by changing ρ from 0 to 40 dB with the interval of 1 dB. Tables 3 and 4 show the optimal and the sub-optimal thresholds of the spatial multiplexing system with ZF receiver, and Tables 5 and 6 show those   Tables 3 and 5, searching for the optimal SNR thresholds fails at high average SNR. This is because lower modulation orders cannot be switched off even when using only one higher http://jwcn.eurasipjournals.com/content/2012/1/238  order modulation for the whole SNR region produces the A-BER below δ 0 , and also because the curve fitting approximation of the BER is used [11]. Note that it is assumed that all the constellations are used in the optimal method. In [11], the last valid thresholds with a lower SNR is used to tackle this situation, but it causes some performance loss. Instead, in the suboptimal method, since it is possible to switch off some constellations, the thresholds can always be obtained. Note that only 64-QAM is used for the entire instantaneous SNR region with high average SNR as shown in Tables 4 and 6. Figures 2 and 3 shows the ASE and the A-BER of the spatial multiplexing system with ZF receiver, and Figures 4  and 5 show those of the OSTBC system. Note that 'Z' , 'H' , 'LM' , and 'DS' stand for the Zhou method [4], the Huang method [11], the Lagrange multiplier method, and the descending search method, respectively. The Zhou method means that the constellation specific constants http://jwcn.eurasipjournals.com/content/2012/1/238 of (6) are used to get the switching thresholds for adaptive modulation, and the Huang method means that the constellation specific constants of Table 1 are used. The A-BER approach has SNR gain of 2-3 dB compared to the I-BER method at the same ASE. In the I-BER approach, even though the exact BER function in AWGN channels is used, the A-BER is still below δ 0 so that it is too conservative. In the A-BER approach, the optimal method ('H+LM' in the figure), in which constellation specific constants defined in , and Lagrange multiplier are used, gives the best ASE performance while satisfying the A-BER constraint. In case that the constants of (6) are used ('Z+LM'), even though Lagrange multiplier is used, gives lower performance. This means that using an accurate BER function in (13) is important. The sub-optimal method ('H+DS') also gives as good performance as the optimal method. Since the optimal thresholds are not available in high average SNR range (ZF: ≥33 dB, OSTBC: ≥23 dB), the last valid thresholds are used (ZF: 32 dB, OSTBC: 22 dB) as mentioned in [11]. At low average SNR, the optimal method performs better in terms of ASE than the sub-optimal method. At high average SNR, the suboptimal method is slightly better in terms of ASE than the optimal method because the search of the optimal SNR thresholds fails at high average SNR, and the last valid threshold with a lower SNR was used. Thus, we can use a hybrid scheme which uses the optimal thresholds at low average SNR, and the sub-optimal thresholds at high average SNR.

ASE and BER performance
Performance of the dynamic rate-adaptive MIMO mode switching Figures 6 and 7 shows the ASE and A-BER of the proposed scheme. The SNR thresholds of the proposed scheme are obtained from the A-BER approach using the hybrid approach described in the previous section. At low average SNR, the optimal thresholds are used, whereas the sub-optimal thresholds are used at high average SNR where the optimal thresholds cannot be computed (i.e., ZF: ≥ 33 dB, OSTBC: ≥23 dB). As shown in Figure 6, the two ASE curves that only one MIMO mode is used cross at about 17 dB. Thus, the mode switching rule of the static method is as follows: if ρ < 17, the diversity mode (OSTBC) is selected. Otherwise, the spatial multiplexing mode (ZF) is selected. It is observed that the proposed dynamic switching scheme has 1-2 dB SNR gain over the static switching scheme in the 10-26 dB SNR range while keeping the A-BER close to the target value as shown in Figure 7. It is also observed in Figure 6 that the proposed mode switching criterion has the same performance as the Demmel condition number (DCN) method. Therefore, the proposed criterion is preferable in terms of computational complexity when the ISE of each MIMO mode is equal to each other.

Conclusion
In this article, we proposed a dynamic rate-adaptive MIMO mode switching scheme between spatial multiplexing and diversity modes. The proposed dynamic scheme shows better ASE performance than the static MIMO mode switching scheme with adaptive modulation. In the conventional static mode, the MIMO mode changes depending on the average SNR instead of the instantaneous channel condition. We also proposed a suboptimal algorithm to find the SNR thresholds with the A-BER constraint in the high average SNR range. When the two MIMO modes have the same ISE, we can use a new mode switching criterion based on the BER exponent. The BER exponent method shows the same performance as the Demmel condition number method with lower computational complexity. Simulations show that the performance of the proposed algorithm is close to that of the optimal case. The proposed algorithms appear to be promising as a practical MIMO mode switching technique. http://jwcn.eurasipjournals.com/content/2012/1/238 Endnote a In this article, we use boldface lowercase letters to denote vectors, boldface uppercase to denote matrices. s k is the kth element of the vector s while [ H] i,j is the element in the ith row and jth column of the matrix H. H F is the Frobenius norm on H. b We use (·) H for conjugate transpose, and E to denote expectation. c θ L+1 is fixed as infinity.

Spatial multiplexing system with zero-forcing receiver
It has been shown in [16] that the pdf of γ m,ZF in (2) the ASE η ZF is given by Since it is assumed that channel is i.i.d., the pdf 's of all the γ m,ZF 's are identical.
Thus, dropping the subscript m, the ASE can be rewritten as η ZF can be expressed as a closed form as follows: where P e,ZF (l) is given by (30)

Orthogonal space-time block coding system
It is known that the pdf of γ OSTBC in (4) is distributed as [7] f Since the pdf is also Gamma distributed as that of ZF case, the derivation procedure is similar. Defining the set of adjustable thresholds OSTBC as OSTBC = {θ 1 , θ 2 , . . . , θ L }, the ASE η OSTBC is given by Substituting (31) into (33), η OSTBC can be expressed as a closed form as follows: where D OSTBC = NM. Likewise, the A-BER P e,avg,OSTBC can be denoted as where P e,OSTBC (l) is given by where D O means D OSTBC .