Research Article  Open Access
Xiaoming Chen, "Spatial Correlation and Ergodic Capacity of MIMO Channel in Reverberation Chamber", International Journal of Antennas and Propagation, vol. 2012, Article ID 939104, 7 pages, 2012. https://doi.org/10.1155/2012/939104
Spatial Correlation and Ergodic Capacity of MIMO Channel in Reverberation Chamber
Abstract
It has previously been shown that a reverberation chamber can conveniently be used to measure ergodic multipleinput multipleoutput (MIMO) capacity for overtheair (OTA) tests. However, the MIMO channel in the reverberation chamber has not been fully studied before. In this paper, the spatial correlation of the MIMO channel in the chamber is studied by comparing the measured channel with two popular MIMO channel models. It is shown that the models can accurately predict the ergodic MIMO capacity of the channel in the reverberation chamber, but not the outage capacity (especially at high signaltonoise regime). It is verified that the capacity estimation error is due to the fact that the measured MIMO channel in the chamber does not satisfy multivariate normality (MVN), which causes the capacity error increases additively with MIMO size and multiplicatively with signaltonoise (SNR).
1. Introduction
Multipleinput multipleoutput (MIMO) systems have drawn considerable popularity, over the past decade, due to their performanceenhancement capability in multipath environments [1]. Lots of work has been carried out for measuring MIMO capacity in reallife (outdoor and indoor) multipath environments [2–6]. As opposite to the reallife measurements, the reverberation chamber is being considered as a strong candidate for standardization of overtheair (OTA) measurements for characterization of MIMO terminals due to its fast, repeatable, and economical measurements. The ergodic MIMO capacity can be measured in reverberation chambers [7–11]. While most of the reverberation chamber works have been focused on characterizing multiport antennas, the MIMO channel in the chamber has not been fully studied yet.
The aim of this paper is to study the measured channel in the chamber by comparing it with two wellknown channel models, (i.e., Kronecker model and fullcorrelation model [1]) and physically explain the models’ discrepancies with measurements. Some different properties of reverberation chamber and reallife multipath environments are also discussed. Most previous reallife measurements only used Kronecker model to compare with measurements [2–6]. Although there exist other sophisticated models, it is well known that the fullcorrelation model offers the best accuracy at the expanse of the most complexity. Therefore, we include both Kronecker model and the fullcorrelation model in this paper. It is found that both models have the same performance in terms of capacity estimation, which is a bit surprising. The reasons for this are discussed in Section 4. It is also found that, although both models can well predict the ergodic capacity (with only slight overestimation at high SNR regime), they fail to estimate the outage capacity, or the cumulative distribution function (CDF) of the capacity. Moreover, instead of capacity underestimations as indicated by reallife measurements [4], the present paper shows that Kronecker model tends to overestimate capacities based on measurements in a reverberation chamber. The reasons for this are discussed in Section 5.
This work is of particular interest for OTA characterization of MIMO terminals in reverberation chambers, because it helps to understand the channel conditions under which the passive and/or active MIMO measurements have been conducted in the chamber.
2. Measurement Description
It has been shown that the ergodic MIMO capacity of a multiantenna system can be easily determined based on the reverberation chamber measurement [7, 9–11]. The reverberation chamber is basically a metal cavity with many excited modes that are stirred to create a Rayleighfading environment [12]. The chamber used in the present paper is Bluetest HP reverberation chamber with a size of (see Figure 1). It has two plate mode stirrers, a turntable platform, and three antennas mounted on three different walls (referred to as wall antennas thereafter). The wall antennas are actually wideband halfbowtie (or triangular sheet) antennas. In the measurements, the platform (with a radius of 0.3 m), on which the MIMO terminal under test was mounted, was moved to 20 positions equally spaced by 18°, and, for each platform position, the two plates simultaneously moved to 10 positions (equally spanned on the total distances that they can move along the walls). All the mechanical (stepwise) moves were controlled by a computer. At each stirrer position and for each wall antenna, a full frequency sweep was performed by the vector network analyzer (VNA), during which the channel transfer functions at different frequencies were sampled. The frequency step was set to 1 MHz always. Therefore, for each transmit and receive antenna element pair, there are 200 channel transfer function samples per frequency point.
In order to calibrate out the longterm fading, or attenuation, in the chamber (so that only shortterm fading came into play) [1], a reference measurement needed to be performed first, where the average power transfer function was measured using a reference antenna with known radiation efficiency. The reference level, , was obtained by dividing the average power function by the radiation efficiency of the reference antenna. Then, the multiport antenna under test, in this case a sixmonopole array (see also Figure 1), was measured (actually this procedure was repeated for each of the monopole while the others are terminated with 50 ohm). During this measurement, the three wall antennas were assumed to be the transmit antenna elements; the monopole array was assumed to be the receive antenna. The monopoles have physical length of 8.3 cm (resonating at around 900 MHz). The ground plane has a radius of 14 cm. The monopoles are uniformly mounted on the ground plane in a circle, where adjacent monopoles have a separation of 4.6 cm. The small separation is necessary to have reasonably large correlations, in order to effectively compare different channel models [4].
The resulting channel matrix is a function of frequency and stirrer positions. For convenience, we introduce the following notation for the normalized measured channel matrix where the reference level, , is described above. Note that the total radiation efficiency of the wall antenna is also calibrated out by (1).
Assume that the receiver has perfect channel state information and that transmitted power is equally allocated among transmitting antenna elements, the ergodic MIMO capacity can be computed from the measured channel matrices by [1] where is signaltonoise ratio (SNR) and the expectation can be approximated by averaging over all channel samples.
3. MIMO Channel
3.1. Channel Characterization
Wireless channel that can be assumed as widesense stationary uncorrelated scattering (WSSUS) [13] is usually characterized by its coherence bandwidth, coherence time, and coherence distance, or equivalently by their reciprocal counterparts: delay spread, Doppler spread, and angular spread, respectively [1]. Since the channel was sampled (by the VNA) stepwisely, that is, in static condition at each stirrer position, the Doppler spread of the channel is negligible. The angular distribution in the chamber is shown to be threedimensional (3D) uniform, that is, isotropic [14]. The coherence bandwidth of the channel at the frequency of interest is around 12 MHz [15]. In practice, wireless channels are seldom WSSUS. Fortunately, most of them can be assumed as quasiWSSUS, that is, the channel statistics do not change within certain time and frequency intervals. These intervals are defined in [16] as stationarity time and stationarity bandwidth. Since channel was sampled under static condition, its stationarity time is infinite. Under the assumption of correlation underspread, the stationarity bandwidth is larger than 20 MHz (i.e., at least 10 times larger than the coherence bandwidth). In the reverberation chamber, every subchannel (for each transmit and receive antenna element pair) is assumed to have the same timefrequency statistics.
Note that different fadingtype environments can be emulated involving reverberation chamber. Holloway [17] showed that loading the chamber with electromagneticabsorbing objects can generate Rician fading. It was shown in [10] that two cascaded reverberation chambers can be used to emulate doubleRayleigh (keyhole) fading. More outdoorlike fading can be generated by connecting a reverberation chamber to an anechoic chamber [18]. As pointed out in [19], for MIMO channels, the most important aspect is the spatial correlation of the channel. Therefore, the main focus of this paper will be on the spatial correlation of the channel and its effects on capacity in an unloaded (with little electromagneticabsorbing objects) reverberation chamber.
3.2. Channel Models
Assume a singleuser narrow band MIMO system, consisting of transmit antennas and receive antennas, in a frequency flat Rayleighfading environment. The output vector of the MIMO system can be expressed as where is the channel matrix, is vector of transmit signals, and is zero mean white Gaussian noise vector. This paper will focus on narrowband models only, knowing that they can be readily extended to wideband model by assuming channel taps with different delays (with each tap modeled as a narrowband channel) [13].
A general, socalled, fullcorrelation channel model is given by [1] where is independent and identically distributed (i.i.d.) zero mean spatially white × channel matrix, is the operator stacking the matrix into a vector columnwise, , and is the covariance matrix that can be expressed as where superscript is Hermitian operator. is normalized so that its Frobenius norm satisfies . This normalization has the same physical meaning as (1), that is, path loss calibration. The fullcorrelation model is acknowledged as the most accurate model. However, it suffers from huge covariance matrix size and analytical intractability.
The Kronecker model assumes separable , that is, where denotes Kronecker product and and are covariance matrices at the transmit and receive sides, respectively. The superscript is transpose operator. Under this assumption, can be represented using the socalled Kronecker model [1], This model has been verified in [2, 3, 5] by reallife measurements. This model is relatively simple and analytically tractable. Furthermore, it allows independent optimizations of transmit and receive MIMO terminals. Therefore, the Kronecker model becomes the most popular MIMO channel model. However, it was shown in [4, 6] that Kronecker model leads to inaccurate capacity estimation. It was pointed out in [4] that the Kronecker model only rendered correct capacity when the antenna number at either side is no larger than three and that it tends to underestimate capacity otherwise.
4. Measurement Results
From Section 1, it is known that there are only 200 MIMO channel samples at one single frequency, which is too few to support reliable estimation. One simple way of increasing channel samples is to treat the channel samples at different frequencies as different channel realizations (or samples). This methodology has been used in [2] for reallife measurements. In a reverberation chamber, it is usually referred as frequency stirring or electronic stirring [20] (to increase independent samples). However, the frequency stirring bandwidth has to be carefully chosen so that more independent (or less correlated) samples can be included without changing the channel statistics. From Section 3, the coherence bandwidth of the channel is around 12 MHz, while the stationarity bandwidth is larger than 20 MHz. In practice, the antenna bandwidth will also affect the channel characteristics, since the (effective) channel also includes the antennas. Hence, the frequency stirring bandwidth should be larger than coherence bandwidth but smaller than stationarity bandwidth and antenna bandwidth. As a result, an empirical frequency stirring bandwidth of 8 MHz (with 1MHz frequency step) is chosen, which is limited by the monopole bandwidth. As a result, there are 1600 MIMO channel samples for capacity evaluation.
The full channel covariance matrix is estimated from as where is the number of samples. Similarly, estimations of the covariance matrices at the transmit and receive sides are, respectively, Note that covariance matrices calculated using (8) and (9) include both antenna efficiencies and correlation coefficients. The corresponding ergodic capacity of the fullcorrelation model is Similarly, the corresponding ergodic capacity of the Kronecker model is Apart from comparing capacity, the discrepancy between the fullcorrelation and Kronecker models can be examined using the Kronecker correlation error defined in [2, 6], To examine the Kroneker model error as a function of the antenna element number, we fix the three wall antennas and increase the number of monopoles from two to six (the monopoles are always chosen from adjacent ones). The Kronecker correlation errors are then calculated and plotted in Figure 2. It can be seen that the Kronecker correlation model error increases with increasing receive monopole number. Similar result was shown in [2] based on reallife measurements.
From Figure 2, it is tempting to anticipate that the Kronecker model has worse capacityestimation performance than the fullcorrelation model for a large antenna element number in a reverberation chamber, just as the cases for reallife multipath environments [4]. However, by comparing capacities based on the Kronecker model, fullcorrelation model, and measurements in the reverberation chamber (see Figures 3 and 4), it is shown that the Kronecker correlation error obtained in the chamber is a bit misleading in the sense that even when it is as high as 10%, the capacities obtained using both channel models are almost the same, in spite of the deviation from measured ones at high SNR regime.
(a)
(b)
(c)
(a)
(b)
(c)
Figure 3 shows the comparison of ergodic capacities of the Kronecker and fullcorrelation models against measurements for the (two adjacent monopoles and three wall antennas), (four adjacent monopoles and three wall antennas), and (all six monopoles and three wall antennas) MIMO systems. As can be seen, the Kronecker model gives almost identical ergodic capacity as the fullcorrelation model; both models well predict ergodic capacity from measured channel in the chamber. Both Kronecker and fullcorrelation models slightly overestimate the measured ergodic capacity as the monopole number exceeds three. The ergodic capacity of the corresponding i.i.d. channels are also plotted for each case. The capacity degradation of the measurement from the corresponding i.i.d. channel is due to antenna correlations and efficiencies.
Figure 4 shows the CDFs of the capacities from the Kronecker model, fullcorrelation model, and measurements for the same , , and MIMO systems with SNRs of 10 dB and 25 dB. For all cases, it is seen that the Kronecker and fullcorrelation models give almost the same result. However, both models overestimate the measured capacity CDFs for more than three antenna elements, especially at high SNR regime, opposite to the capacity underestimation of the Kronecker model observed in [4] based on reallife measurements. The capacity overestimation of the Kronecker model was also observed in [2] for a MIMO system based on reallife measurements. Note that it is shown in [21] that Kronecker model could also overestimate capacity.
It is shown that the Kronecker model has the same performance as the fullcorrelation model because that multibounce rich scattering property of the chamber makes the correlations at transmit and receive sides separable [22]. And, due to this property, the reverberation chamber can be used to characterize the MIMO performance of a MIMO terminal independent of the other MIMO end, which is highly desired in OTA MIMO terminal tests.
5. Multivariate Normality Test
The good agreement between the Kronecker model and the other two advanced models means that, unlike in reallife multipath environments, the discrepancy of capacity of the Kronecker model with the measured one in the reverberation chamber is not due to the Kronecker structure. Instead, it is probably because the entries of the measured MIMO channel matrix are not jointly Gaussian (or normal). It has been found that the HenzeZirkeler’s test [23] has a good overall power against other alternatives to MVN [6, 24]. Therefore, it is applied to the measured channel matrices in the chamber.
Let be the vector of the th channel sample, then the test statistic is where The mean and variance of are, respectively, Let the null hypothesis be that is MVN, that is, the test statistic is approximately lognormally distributed. The probability that the null hypothesis is rejected although it is true is called significance level [25]. If the complementary CDF (CCDF) at for a lognormal distribution ( value), is no smaller than the significance level, satisfies MVN. Otherwise, does not satisfy MVN. A good significance level for HenzeZirkeler’s MVN test is found to be 0.05 [25], which is also used in this paper.
Applying the HenzeZirkeler’s test to the measured channel matrices, it is found that every subchannel coefficient of is still in (complex) normal distribution, but they are together not jointly normal, that is, does not satisfy MVN. This finding is a bit surprising, since it is well known the singleinput singleoutput (SISO) channel in a reverberation chamber is in (complex) normal distribution [12]. And it is usually assumed that MIMO channels in the chamber are jointly normal. To find out what is the reason of the nonMVN, the submatrices of are tested separately. The corresponding values are listed in Table 1. By comparing the values with the significant level (0.05), it is found that the subchannel vectors of the monopoles (for one wall antenna) are jointly normal up to five monopole elements, while the subchannel vectors of the wall antennas (for one monopole) are only jointly normal for two wall antennas, that is, the subchannels of the three wall antennas not jointly normal.

It is found that the channels with large MIMO sizes show strong nonMVN and that channels with MIMO small size (e.g., ) approximately satisfy MVN. The different MVN properties for the monopole array and wall antennas are due to the fact that the monopole array is uniformly and circularly distributed (with the same polarization), while the wall antennas are located arbitrarily on three different walls with orthogonal polarizations (which is necessary to keep a desired polarization balance). Thus, the monopole array satisfies antenna stationarity (i.e., antenna correlation depends only on the antenna element separation, not the element positions) [19], while wall antennas do not. The antenna stationarity of an antenna array enables its correlation matrix to converge, which guarantees MVN of the channels seen by the array (cf. weakly converging Gaussian vectors theorem [26]). Note that a wellstirred reverberation chamber is polarization balanced, that is, there is no polarization preference in the chamber. Thus, there is no need to introduce a polarization matrix to the Kronecker model as for the polarizationunbalanced cases [1].
Due to the nonMVN of , neither of the abovementioned channel models can estimate the MIMO capacity accurately, because all of them involve (i.i.d. zero mean complex Gaussian matrix). It can be seen from the MIMO capacity formula (2) that the channel model error is additively increased by MIMO size and multiplicatively increased by SNR. This explains why the capacity estimation error increases with either increasing MIMO size or increasing SNR. And, therefore, the nonMVN is the main contribution to capacity estimation errors of the channel models in the reverberation chamber. It seems that the only way to circumvent this problem is to give up the MVN assumption in channel modeling. However, the ultimate goal of channel modeling is to find a model that can predict the actual channel correctly and, at the same time, can be used with reasonably low complexity from communication and informationtheoretic viewpoints. Since MVN is virtually the only mathematically tractable multivariate distribution [26], any model giving up the MVN assumption might not be very much useful from communication and informationtheoretic point of views, even though it might offer better estimation accuracy. Therefore, no effort is exerted in finding a model with better accuracy in this paper.
Note that the nonMVN may not be the main contribution for channel model errors for reallife multipath environments, where the fullcorrelation model can outperform the Kronecker model. The equal performance of all the models for the reverberation chamber measurement is because of the multibounce rich scattering property of the chamber, as discussed in Section 4.
6. Conclusion
In this paper, the MIMO channel in a reverberation chamber is studied by comparing the measurements with different channel models. It is found that both models have the same performance in terms of capacity estimation for the reverberation chamber measurements and that all of them can well predict the ergodic capacity up to six antenna elements with only slight overestimation at highSNR regime. However, the models fail to predict the CDFs of the capacities for more than three antenna elements, especially at high SNR regime. The reason for this is because of the nonMVN of the MIMO channel in the chamber. Since all the models involves i.i.d. complex Gaussian channel, there will be modeling errors due to the nonMVN of the measured channel. And the channel modeling errors will additively increase with MIMO size and multiplicatively increase with SNR, for MIMO capacity estimations. The equal performance of the Kronecker model and the fullcorrelation model (for reverberation chamber measurements) implies that the correlations at transmit and receive sides can be treated separately. This is actually very desirable, since it allows characterizing the performance of a MIMO terminal independently (without the effect of the other MIMO side) by doing an OTA test in the chamber, which in turn allows fair comparisons of different MIMO terminals.
Acknowledgments
This work has been supported by The Swedish Governmental Agency for Innovation Systems (VINNOVA) within the VINN Excellence Center Chase. The author would like to thank Dr. Li Yang and the anonymous reviewers for their helpful comments to increase the quality of this paper.
References
 A. Paulraj, R. Nabar, and D. Gore, Introduction to SpaceTime Wireless Communication, Cambridge University Press, Cambridge, Mass, USA, 2003.
 K. Yu, M. Bengtsson, B. Ottersten, D. McNamara, P. Karlsson, and M. Beach, “Modeling of wideband MIMO radio channels based on NLoS indoor measurements,” IEEE Transactions on Vehicular Technology, vol. 53, no. 3, pp. 655–665, 2004. View at: Publisher Site  Google Scholar
 J. P. Kermoal, L. Schumacher, K. I. Pedersen, P. E. Mogensen, and F. Frederiksen, “A stochastic MIMO radio channel model with experimental validation,” IEEE Journal on Selected Areas in Communications, vol. 20, no. 6, pp. 1211–1226, 2002. View at: Publisher Site  Google Scholar
 H. Özcelik, M. Herdin, W. Weichselberger, J. Wallace, and E. Bonek, “Deficiencies of “Kronecker” MIMO radio channel model,” Electronics Letters, vol. 39, no. 16, pp. 1209–1210, 2003. View at: Publisher Site  Google Scholar
 D. Chizhik, J. Ling, P. W. Wolniansky, R. A. Valenzuela, N. Costa, and K. Huber, “Multipleinputmultipleoutput measurements and modeling in Manhattan,” IEEE Journal on Selected Areas in Communications, vol. 21, no. 3, pp. 321–330, 2003. View at: Publisher Site  Google Scholar
 T. Svantesson and J. W. Wallace, “Tests for assessing multivariate normality and the covariance structure of MIMO data,” in Proceedings of the IEEE International Conference on Accoustics, Speech, and Signal Processing, Hong Kong, April 2003. View at: Google Scholar
 K. Rosengren and P.S. Kildal, “Radiation efficiency, correlation, diversity gain and capacity of a sixmonopole antenna array for a MIMO system: theory, simulation and measurement in reverberation chamber,” in Proceedings of the IEEE Microwaves, Optics and Antenna, vol. 152, pp. 7–16, February 2005. View at: Publisher Site  Google Scholar
 J. F. Valdés, M. A. Fernandez, A. M. Gonzalez, and D. A. Hernandez, “The influence of efficiency on receive diversity and MIMO capacity for Rayleighfading channels,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 5, pp. 1444–1450, 2008. View at: Publisher Site  Google Scholar
 L. GarciaGarcia, B. Lindmark, N. Jaldén, and C. Orlenius, “MIMO capacity of antenna arrays evaluated using radio channel measurements, reverberation chamber and radiation patterns,” IET Microwaves, Antennas and Propagation, vol. 1, no. 6, pp. 1160–1169, 2007. View at: Publisher Site  Google Scholar
 O. Delangre, P. de Doncker, M. Lienard, and P. Degauque, “Coupled reverberation chambers for emulating MIMO channels,” Comptes Rendus Physique, vol. 11, no. 1, pp. 30–36, 2010. View at: Publisher Site  Google Scholar
 X. Chen, P.S. Kildal, J. Carlsson, and J. Yang, “Comparison of ergodic capacities from wideband MIMO antenna measurements in reverberation chamber and anechoic chamber,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 446–449, 2011. View at: Publisher Site  Google Scholar
 J. G. Kostas and B. Boverie, “Statistical model for a modestirred chamber,” IEEE Transactions on Electromagnetic Compatibility, vol. 33, no. 4, pp. 366–370, 1991. View at: Publisher Site  Google Scholar
 R. Vaughan and J. B. Andersen, Channels, Propagation and Antennas for Mobile Communications, Intelligent Energy Europe, London, UK, 2003.
 D. A. Hill, “Plane wave integral representation for fields in reverberation chambers,” IEEE Transactions on Electromagnetic Compatibility, vol. 40, no. 3, pp. 209–217, 1998. View at: Google Scholar
 X. Chen, P.S. Kildal, C. Orlenius, and J. Carlsson, “Channel sounding of loaded reverberation chamber for overtheair testing of wireless devices: coherence bandwidth versus average mode bandwidth and delay spread,” IEEE Antennas and Wireless Propagation Letters, vol. 8, pp. 678–681, 2009. View at: Publisher Site  Google Scholar
 G. Matz, “Characterization of nonWSSUS fading dispersive channels,” in Proceedings of the International Conference on Communications (ICC '03), pp. 2480–2484, Anchorage, Alaska, USA, May 2003. View at: Google Scholar
 C. L. Holloway, D. A. Hill, J. M. Ladbury, P. F. Wilson, G. Koepke, and J. Coder, “On the use of reverberation chamber to simulate a rician radio environment for the testing of wireless devices,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 11, pp. 3167–3177, 2006. View at: Publisher Site  Google Scholar
 H. Fielitz, K. A. Remley, C. L. Holloway, Q. Zhang, Q. Wu, and D. W. Matolak, “Reverberationchamber test environment for outdoor urban wireless propagation studies,” IEEE Antennas and Wireless Propagation Letters, vol. 9, Article ID 5406096, pp. 52–56, 2010. View at: Publisher Site  Google Scholar
 M. Herdin, Nonstationary indoor MIMO radio channels, Ph.D. thesis, Vienna University of Technology, Vienna, Austria, 2004.
 D. A. Hill, “Electronic mode stirring for reverberation chambers,” IEEE Transactions on Electromagnetic Compatibility, vol. 36, no. 4, pp. 294–299, 1994. View at: Publisher Site  Google Scholar
 V. Raghavan, J. H. Kotecha, and A. M. Sayeed, “Why does the Kronecker model result in misleading capacity estimates?” IEEE Transactions on Information Theory, vol. 56, no. 10, Article ID 5571876, pp. 4843–4864, 2010. View at: Publisher Site  Google Scholar
 A. S. Y. Poon, D. N. C. Tse, and R. W. Brodersen, “Impact of scattering on the capacity, diversity, and propagation range of multipleantenna channels,” IEEE Transactions on Information Theory, vol. 52, no. 3, pp. 1087–1100, 2006. View at: Publisher Site  Google Scholar
 N. Henze and B. Zirkler, “A class of invariant consistent tests for multivariate normality,” Communications in Statistics—Theory and Methods, vol. 19, no. 10, pp. 3595–3618, 1990. View at: Publisher Site  Google Scholar
 A. TrujilloOrtiz, R. HernandezWalls, K. BarbaRojo, and L. CupulMagana, HZmvntest: HenzeZirkler's Multivariate Normality Test, 2007.
 R. B. Agostino and M. A. Stephen, GoodnessofFit Techniques, Marcel Dekker, New York, NY, USA, 1986.
 I. I. Gikhman and A. V. Skorokhod, Introduction to the Theory of Random Processes, Dover publication, New York, NY, USA, 1996.
Copyright
Copyright © 2012 Xiaoming Chen. 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.