Implementations of SignalProcessing Algorithms for OFDM Systems
View this Special IssueResearch Article  Open Access
Probabilistic Behavior Analysis of MIMO Fading Channels under Geometric Mean Decomposition
Abstract
Geometric mean decomposition (GMD) has been proposed as a method to realize multiple spatial links with identical gains that are intrinsic to a MIMO channel. In order to simplify system design and implementation based on knowledge regarding probability behavior of MIMOGMD schemes, the main objective of this paper is to statistically characterize the link gains and channel capacities that can be provided via GMD. In particular, closedform univariate and bivariate probability density functions (PDFs) for these metrics under Rayleigh fading are derived using Gamma approximations. By applying these analytical results, the fluctuations of MIMOGMD schemes are examined by modeling both link gains and capacities using finitestate Markov chains (FSMCs).
1. Introduction
Numerous research results have shown that spatial multiplexingbased multipleinput multipleoutput (MIMO) transceivers can boost the data rate of wireless communication systems without using additional power and channel bandwidth. When both transmitter and receiver are equipped with multiple antennas, multiple independent data streams can be concurrently launched into the air by the transmitter and separately decoded at the receiver, with aid of certain signal processing techniques [1]. By combining MIMO with the multicarrier schemes such as orthogonal frequency division multiplexing (the socalled MIMOOFDM), channel capacity can be tremendously increased as the degree of freedom in frequency and spatial domain are jointly utilized.
By assuming that channel state information (CSI) is available at both ends of the link, singular value decomposition (SVD) [2] is commonly considered as the approach to extract spatial links within a MIMO channel. With SVD, the MIMO channel matrix can be transformed into a bank of parallel scalar airpipes (widely termed as eigenmodes in the literature) by applying appropriate beamforming matrices at both transmitter and receiver. Since the eigenmodes realized by SVD are equivalent to singular values of the MIMO channel matrix, their link strengths could be quite different. Therefore, proper bitallocation and power distribution (such as waterfilling) algorithms must be employed to optimize the system capacity. This, however, increases the complexity in adaptive control mechanism, as bit/power allocations among eigenmodes may have to be updated frequently in rapidly timevarying channels. For the sake of convenience, one may simply apply the modulation/coding format or distribute the power uniformly on all eigenmodes. Nevertheless, this may cause degradation on error performance of the system, especially when magnitude of the minimum channel singular value is too small [3].
Geometric mean decomposition (GMD), proposed by Jiang et al. [4], has emerged as an alternative method of designing MIMO transceiver. A set of parallel spatial links with identical gains intrinsic to a MIMO channel can be realized via GMD. Remarkably, these parallel links with the identical gain are equal to the geometric mean of channel channel singular values (eigenmodes). The property of parallel links with identical gain makes GMD useful for adaptive MIMOOFDM systems, as complicated joint bit/power allocation algorithms in frequency and spatial domain can be significantly simplified. In recent years, much attention from both academia and industry has been paid to the GMDbased MIMO transceiver. For instance, GMD has been deemed as a prospective MIMO architecture in some nextgeneration cellular broadband networks, including two of the most important standards in this contextโLTEAdvanced [5] and WiMAX (IEEE 802.16โm) [6]. Note that Jiang et al. have put forward another MIMO transceiver dubbed as uniform channel decomposition (UCD) [7], which also provides spatial links with identical gains. The scope of this paper is focussed on GMD, and the related work on UCD is, however, remained as open problems for future research.
In order to calculate performance metrics such as bit error rate and outage probability of a GMDbased MIMO transceiver operating in mobile fading channels, analytical results on statistical properties of MIMO eigenmode geometric mean are required. To the best of our knowledge, however, only a handful of papers in the existing literature have offered the relevant analysis. In [8], asymptotic behavior of geometric mean of MIMO channel eigenvalues has been examined. Also, the probability density function (PDF) for the link gain of MIMOGMD in special cases with two antennas at both transmitter and receiver can be found in [9], but the results for more general cases are not available. This motivated us to derive closedform expressions of PDFs (either exact or approximated) for systems with arbitrary number of antennas in Rayleigh fading channels. Specifically, apart from the univariate PDFs, we are also interested in the bivariate PDFs for MIMO eigenmode geometric mean and channel capacity, since the statistical behavior relating to channel variation is crucial for systems incorporating with adaptive transmission techniques. For instance, knowledge on how channel quality evolves with time allows engineers to determine the adaptation/feedback rate for MIMOGMD in a more judicious manner. Further, for successive interference cancelation (SIC) detection algorithm in conventional MIMOGMD transceivers [4], the reliability of a data stream is dependent on previously decoded data streams. That is, if the receiver fails to detect a data stream correctly, the data on the remaining streams may not be decoded successfully due to error propagation. Thus, it is possible to gauge how frequently that an error propagation event would occur if timevarying characteristics are known. In correspondence, this paper applies firstorder finitestate Markov chains (FSMC) [10, 11] to model the fluctuations of both link gain and channel capacity under MIMOGMD in baseline cases of Rayleigh fading environments. Note that FSMC has been used to emulate the capacity process of general MIMO wireless channels [12], but channels with MIMOGMD transceiver have not been considered. Also, some researchers have designed the feedback mechanisms for adaptive modulation systems with channels modeled by FSMC [13]. Based on the univariate and bivariate PDFs developed in this paper, the transition probabilities among the discrete states of an FSMC are analytically computed.
In a nutshell, this paper encompasses the following contributions. (1)Derivation of the exact PDFs for link gains and channel capacities of MIMOGMD systems with two antennas at either transmitter or receiver and an arbitrary number of antennas on the other side. (2)Demonstration of Gamma approximations for link gains of MIMOGMD systems with arbitrary numbers of antenna at either transmitter or receiver. (3)Derivation of PDFs for channel capacities of MIMOGMD systems with arbitrary numbers of antenna at either transmitter or receiver, based on Gamma approximation of link gains. (4)Extending the Gamma approximation to provide bivariate PDFs for link gains and channel capacities of MIMOGMD systems with numbers of antenna at either transmitter or receiver. (5)Constructing firstorder FSMCs for both link gains and channel capacities of MIMOGMD systems. The transition probabilities among multiple quantized Markov states are computed based on the derived PDFs.
The rest of the paper is organized as follows. In Section 2, assumptions on channel model and background of the MIMOGMD architecture are briefly reviewed. Then, in Section 3, the closedform expressions of both univariate and bivariate PDFs for eigenmode geometric mean and channel capacity of MIMOGMD schemes are derived. Based on the derived analytical results, FSMCs are constructed in Section 4 to study the time variation of MIMO channel under GMD. Finally, this paper ends with a concise conclusion drawn in Section 5.
2. Background and Assumptions
In this section, we first describe the considered system model. Then, the concept of GMD and the fundamental structure of an FSMC are reviewed subsequently.
2.1. Channel Model
In general, this paper mainly studies a singleuser, pointtopoint MIMO system that has transmit antennas and receive antennas. The MIMO channel, , is therefore an matrix. Here, we denote and . We presume this MIMO system is coupled with OFDM setting, so channel response per subcarrier is basically flat in frequency domain. Therefore, the signal model for each subcarrier can be written as: where is an vector of the received signal, is an vector of the transmitted data symbols, and is an vector of additive complex zero mean Gaussian noise. Since Rayleigh fading is considered, all entries of are modeled as complex Gaussian random variables with zero mean and unit variance, . With SVD, the MIMO channel can be written as where both and are unitary matrices, denotes Hermitian transpose, and is a diagonal matrix with nonzero elements representing the singular values of . The unordered joint density for the singular values of is [14] where . Note that the eigenvalues of the channel correlation matrix, , are equivalent to power gains of the spatial links intrinsic to the MIMO channel, and their joint statistics are governed by the central Wishart distribution as given in [15] It is further assumed that the MIMO channel evolves over time in accordance to where is the zerothorder Bessel function, and are the Doppler frequency and time displacement, respectively, is a random matrix with zeromean, unit variance complex Gaussian elements. For the sake of simplicity, we denote the autocorrelation function as in the remainder of this paper.
2.2. Review of Geometric Mean Decomposition
The authors of [4] suggested designing MIMO transceivers via geometric mean decomposition (GMD), in which the MIMO channel is decomposed as: where and are semiunitary matrices and represents an upper triangular matrix with identical diagonal elements. If the data signal is precoded by at the transmitter and is filtered by at the receiver, the signaling model (1) is tantamount to multiplication between and the upper triangular matrix, . Hence, the transmitted signal could be extracted and decoded using nulling and cancellation procedure as elucidated in [4]. Assuming ideal operations without error propagation, the data signals are sent on parallel spatial links with a common value of amplitude: In other words, the amplitudes of these links are equivalent to the geometric mean of MIMO channel singular values, . Since these spatial links have an identical gain, the โworst subchannelโ problem, in which the overall error performance is dominated by the weakest spatial link, is no longer a concern. Thus, one merit of using GMD is the exemption of complicated bit/power allocations over spatial links. By uniformly allocating the transmission power among all parallel spatial links realized by GMD, the channel capacity can be expressed as: where denotes the signaltonoise ratio (SNR) on each of the parallel links. Also, it is easy to see that the link power gain of these spatial links is the geometric mean of : Conventionally, the channel capacity of a general MIMO system can be written as the sum of multiple random variables; hence, centrallimited theorem (CLT) could be used to approximate the MIMO channel capacity as a Gaussian random variable [16]. For MIMOGMD schemes, however, Gaussian approximation may not be suitable as (8) is simply a scalar integer multiple of a single random variable.
2.3. FiniteState Markov Chain
As one of the objectives in this paper is to model the timevarying behavior of MIMOGMD systems based on FSMC, the fundamental structure of an FSMC is outlined here for pedagogical reasons. To model the fluctuation of a random process using an FSMC, the process is initially quantized into a number of discrete states. Then, the dynamic behavior is described by the transition probabilities among these states. For simplicity, only firstorder FSMC is considered in this paper, which means the transition probability is solely dependent on the most recent observation. In our case, both processes of eigenmode geometric mean and capacity are quantized based on their magnitudes. In particular, the random process of interest is simply partitioned into states by setting threshold levels. In the first state (denoted as ), the process has a value smaller than the lowest threshold level . In the final state (), the process is larger than the highest threshold level . Otherwise, the process is said to be in the state (denoted as ) if its value falls in the range between threshold levels and . For an illustrative example, Figure 1 shows an eigenmode geometric mean process that is partitioned into four states using three threshold levels: , , and . As shown in the later sections of this paper, the transition probabilities can be computed based on both univariate and bivariate PDFs of the process.
3. Derivations of Statistical Distributions
This section targets to derive univariate and bivariate PDFs for both spatial link power gain, , and channel capacity, , for MIMOGMD schemes. For univariate distributions, the exact analytical results are first presented for cases with two spatial links (), and then the cases with more spatial links are dealt with by Gamma approximations.
3.1. Univariate Distributions for Cases with
When either the transmitter or receiver is equipped with only two antennas, only spatial links are realizable as the rank of a MIMO channel is . From (3), the joint PDF for MIMO channel singular values with is Since the geometric mean, , is a function of singular values, Jacobian transform can be applied to obtain the joint PDF for and . We have and hence The PDF for can be acquired by integrating (12) with respect to . After tedious calculation and little algebra, we find the following result: where represents modified Bessel function of the second kind of order 1. The derivation from (12) to (13) could be carried out by certain symbolic manipulation software packages. Note that (13) is a generalized version for cases with ; by setting , (13) can be reduced to the special case result for a (2,2) MIMOGMD system given in [9]. As is the square root of , we may apply simple transformation on (13) to get the PDF for :
Similarly, Jacobian transform can be employed to cope with the channel capacity of MIMOGMD schemes, as the PDF for is now readily available. Based on (8), we have Thus, the PDF for channel capacity with is
The validity of (14) and (16) can be verified through Monte Carlo simulations. In Figure 2, the empirical distributions of eigenvalue geometric mean for (2,2), (2,4), and (2,8) MIMO systems are compared with their respective analytical PDFs that we have derived (14). On the other hand, the comparison between theoretical channel capacity PDFs (16) and simulation results are shown in Figure 3. In all cases, excellent agreements between analytical and simulation results can be observed.
3.2. Univariate Distributions for Cases with
For systems with more than two spatial links, that is, , the similar approach invoked in Section 3.1 could be used. Nonetheless, it is difficult to get closedform results like (14) and (16) because multiple integrals are required in the step of obtaining from as in (13). Fortunately, some previous studies have concluded that MIMO eigenvalues () can be very accurately approximated by Gamma random variables [14]. Moreover, [17] claims that the geometric mean of multiple independent Gamma random variables could be either a Gamma or a mixture of Gamma distributions. Although MIMO eigenvalues are not independent processes, they are weakly correlated. Thus, we may make a hypothesis stating that where the parameters and are shape factor and scale factor of Gamma distribution respectively, and represents Gamma function. Since , the computation for requires the first two moments of . Note that the th moment of can be calculated as where is the joint PDF of all eigenvalues of the channel correlation matrix , as given in (4). To make this paper more selfcontained, we have tabulated the numerical values of , , , and for (4,4), (4,8), and (8,8) cases in Table 1. Hence, one may construct the PDFs for in these MIMO cases accordingly.

The approximations based on Gamma distributions (17) can be further leveraged to find the PDFs for channel capacity in MIMOGMD schemes. Using (15), we have To justify the suitability of Gamma approximations for eigenvalue geometric means, similar comparisons as in Section 3.1 are carried out for (4,4), (4,8), and (8,8) cases. In Figure 4, we can see that Gamma distributions provide excellent approximation for eigenvalue geometric means. Although we did not show it here, it has been found that Gamma distributions also fit the simulation data in cases. Hence, we can claim that, for practical MIMO systems () at least, the distributions of link gains in MIMOGMD can be accurately approximated by Gamma random variables. For channel capacity, we compare the distributions of simulation samples with the capacity PDF (19), which was derived based on Gamma approximation of . The comparisons are shown in Figure 5, and it is clear that the simulation data and computed results are almost aligned with each other. An important observation is that, in contrast to conventional MIMO systems, Gaussian approximation for capacity distribution [16] does not seem to be appropriate in GMDbased MIMO transceiver schemes.
3.3. Bivariate Distributions
Since we have observed that the MIMOGMD link gains can be accurately approximated by a Gamma random variable, it can be conjectured that a bivariate Gamma PDF is well suited for the joint density of their magnitudes at two correlated time instants, and . Thus, by modifying the bivariate Gamma PDF given in [18], we have where and represents the correlation coefficient between and , which is defined as In (21), both and can be calculated by utilizing (4) and (18). Since and are functions of () and (), respectively, the computation for the joint moment, , requires the joint density for MIMO eigenvalues at two time instants, which has been provided in [19]: where with , , , and as before. The joint moment, , can thereby be computed with (22) as Then, can be calculated by substituting (24) into (21). Apparently, the closedform expression for (24) is difficult, if not impossible, to obtain. However, it still can be solved via numerical integrations.
Once the value of is acquired, (20) can be extended to obtain the bivariate PDF for MIMOGMD channel capacity process. Using Jacobian transform, can be written as With PDFs derived in this section, one may construct an FSMC to model the fluctuation of MIMOGMD link gain and capacity, as explained in Section 4.
4. FSMC Construction
In order to investigate the timevarying behavior of MIMOGMD systems, FSMC is adopted to model both link gain and channel capacity. Based on the FSMC structure elaborated in Section 2.3, we delineate how transition probabilities can be analytically calculated using the PDFs derived in Section 3. The transition probabilities from to , denoted as , for the link gain process can be computed as where Note that has been derived as (14) or (17), depending on the value of . Analytical calculations of (27) can be quite difficult, so numerical integrations are used in this work instead. For the sake of simplicity, this paper assumes that the channel variation is slow enough so the process would only transit to one of the adjacent states (from to or ) or stay in the same state (from to ). Therefore, we have For channel capacity in MIMOGMD scheme, the transition probabilities can be computed in a similar fashion.
In order to verify that FSMC is an appropriate tool to model the fluctuation of MIMOGMD channel, several Monte Carlo simulations have been carried out. In particular, simulation results on transition probabilities for both MIMO eigenvalue geometric mean and channel capacity are compared with our calculations by (26) and (29). Note that the threshold levels for state quantization are set arbitrarily in this work. In practical adaptive modulation schemes, for example, the threshold levels for state quantization could be set based on minimum required channel gain for a target error performance. In all simulations, we have โHz, โsec and โdB. In Figure 6, the link gain of a (2,4) MIMOGMD scheme is modeled as an FSMC consisting of four states with , and transition probabilities from both simulation and calculations are plotted. Similarly, we approximate the MIMOGMD capacity process in a (2,2) system using a fourstate FSMC with bps/Hz in Figure 7. In both cases, it is apparent that our calculations can provide very accurate approximations.
5. Conclusion
By applying GMD to a MIMOOFDM system, multiple spatial links with identical gains can be realized within each subcarrier. This assuages the complexity at the transmitter side as the spatial domain bit/power allocation can be simplified. In order to seek for potential opportunities to further decrease the system complexity, this paper has investigated the statistical properties of MIMO systems using GMD. In particular, we have derived PDFs for link gains (MIMO eigenvalue geometric mean) and capacities. Although exact results are only available for cases with two parallel links, Gamma approximations can be used to model eigenvalue geometric means for general MIMO cases with arbitrary antenna configurations. Moreover, these results are extended to derive the bivariate PDFs, which is important for the analysis of timevarying behavior. These results are employed to model the channel variation in MIMOGMD schemes by constructing FSMCs. To be specific, the transition probabilities among states could be computed using the PDFs derived in this paper. This paper proposed a potential approach to predict the fluctuation of MIMO channels under GMD, which may allow the engineers to relent the feedback rate and hence reduce the system burden. The analytical PDF results presented in this paper are for baseline cases of Rayleigh flat fading. For scenarios involve spatial correlation and Ricean fading, statistical properties of MIMOGMD schemes are still open problems and should be addressed in the future.
References
 P. W. Wolniansky, G. J. Foschini, G. D. Golden, and R. A. Valenzuela, โVBLAST: an architecture for realizing very high data rates over the richscattering wireless channel,โ in Proceedings of the IEEE International Symposium on Signals, Systems and Electronics (ISSSE '98), pp. 295โ300, Pisa, Italy, September 1998. View at: Google Scholar
 G. Lebrun, J. Gao, and M. Faulkner, โMIMO transmission over a timevarying channel using SVD,โ IEEE Transactions on Wireless Communications, vol. 4, no. 2, pp. 757โ764, 2005. View at: Publisher Site  Google Scholar
 G. Burel, โStatistial analysis of the smallest singular value in MIMO transmission systems,โ in Proceedings of the International Conference on Signal, Speech, and Image Processing (WSEAS '02), Skiathos, Greece, September 2002. View at: Google Scholar
 Y. Jiang, J. Li, and W. W. Hager, โJoint transceiver design for MIMO communications using geometric mean decomposition,โ IEEE Transactions on Signal Processing, vol. 53, no. 10, pp. 3791โ3803, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 โGMD beamforming for LTEadvanced SU/MUMIMO,โ in R1093799, 3GPP TSG RAN WG1 Meeting 58, Miyazaki, Japan, October 2009. View at: Google Scholar
 โA joint transceiver design for MIMO precoder and antenna selection mechanism,โ in C802.16m08/314r2, IEEE Broadband Wireless Access Working Group TGm Session 55, Macau, China, May 2008. View at: Google Scholar
 Y. Jiang, J. Li, and W. W. Hager, โUniform channel decomposition for MIMO communications,โ IEEE Transactions on Signal Processing, vol. 53, no. 11, pp. 4283โ4294, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 J. Salo, D. Seethaler, and A. Skupch, โOn the asymptotic geometric mean of MIMO channel eigenvalues,โ in Proceedings of the IEEE International Symposium on Information Theory (ISIT '06), pp. 2109โ2113, Seattle, WA, USA, July 2006. View at: Publisher Site  Google Scholar
 D. Tang and M.X. Lv, โAnalysis of multiuser diversity in geometric mean decompositionbased MIMO system,โ in Proceedings of the International Conference on Wireless Communications, Networking and Mobile Computing (WiCOM '07), pp. 61โ64, Shanghai, China, 2007. View at: Publisher Site  Google Scholar
 H. S. Wang and N. Moayeri, โFinitestate Markov channel—a useful model for radio communication channels,โ IEEE Transactions on Vehicular Technology, vol. 44, no. 1, pp. 163โ171, 1995. View at: Publisher Site  Google Scholar
 Q. Zhang and S. A. Kassam, โFinitestate markov model for rayleigh fading channels,โ IEEE Transactions on Communications, vol. 47, no. 11, pp. 1688โ1692, 1999. View at: Google Scholar
 S. Vaihunthan, S. Haykin, and M. Sellathurai, โMIMO channel capacity modeling using Markov models,โ in Proceedings of the 61st IEEE Vehicular Technology Conference (VTC '05Spring), pp. 126โ130, Spring, Stockholm, Sweden, April 2005. View at: Google Scholar
 A. E. Ekpenyong and Y. F. Huang, โFeedbackdetection strategies for adaptive modulation systems,โ IEEE Transactions on Communications, vol. 54, no. 10, pp. 1735โ1740, 2006. View at: Publisher Site  Google Scholar
 P.H. Kuo, Fluctuations of MIMO Wireless Channels: Analysis based on Brownian Motion Theory, Dr. Muller, Saarbrucken, Germany, 2008.
 E. Telatar, โCapacity of multiantenna Gaussian channels,โ European Transactions on Telecommunications, vol. 10, no. 6, pp. 585โ595, 1999. View at: Google Scholar
 P. Smith and M. Shafi, โAn approximate capacity distribution for MIMO systems,โ IEEE Transactions on Communications, vol. 52, no. 6, pp. 887โ890, 2004. View at: Publisher Site  Google Scholar
 R. S. Laurent, โThe distribution of the geometric mean of independent F (or gamma) random variables,โ in Proceedings of the Joint Statistical Meetings (JSM '01), Atlanta, GA, USA, August 2001. View at: Google Scholar
 G. E. Øien, H. Holm, and K. J. Hole, โImpact of channel prediction on adaptive coded modulation performance in Rayleigh fading,โ IEEE Transactions on Vehicular Technology, vol. 53, no. 3, pp. 758โ769, 2004. View at: Publisher Site  Google Scholar
 P.H. Kuo, P. J. Smith, and L. M. Garth, โJoint density for eigenvalues of two correlated complex wishart matrices: characterization of MIMO systems,โ IEEE Transactions on Wireless Communications, vol. 6, no. 11, pp. 3902โ3906, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 PingHeng Kuo and PangAn Ting. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.