About this Journal Submit a Manuscript Table of Contents
International Journal of Antennas and Propagation
Volume 2012 (2012), Article ID 246348, 11 pages
http://dx.doi.org/10.1155/2012/246348
Research Article

Effect of Polarization on the Correlation and Capacity of Indoor MIMO Channels

1Information and Communications Research Laboratories, Industrial Technology Research Institute, Hsinchu 31040, Taiwan
2Department of Electronic Engineering, National Taipei University of Technology, Taipei 10608, Taiwan
3Graduate Institute of Communication Engineering, National Taiwan University, Taipei 10617, Taiwan

Received 17 June 2011; Revised 30 August 2011; Accepted 31 August 2011

Academic Editor: Yan Zhang

Copyright © 2012 Hsun-Chang Lo et al. 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.

Abstract

We describe a simple multiple-input/multiple-output (MIMO) channel measurement system for acquiring indoor MIMO channel responses. Four configurations of the polarization diversity antenna, referred to as , , and , were studied in terms of the capacity of indoor MIMO systems. Measurements were taken for a MIMO system in the 2.4 GHz band. In addition, the channel capacity, singular-value decomposition, and correlation coefficient were used to explain the effects of various polarization schemes on MIMO fading channels. We also propose an analysis method for polarization channel capacity; this method includes the normalization of the received power and polarization effect for different polarization schemes. The validation of the model is based upon data collected in both light-of-sight (LOS) and non-light-of-sight (NLOS) environments. From the numerical simulation results, the proposed analysis method was close to measurements made in an indoor environment.

1. Introduction

The multiple-input/multiple-output (MIMO) system proposed by Foschini [1] is a remarkable structure that achieves large capacity through parallel channeling. However, the capacity of an MIMO system is highly dependent on the correlation properties of the channels. It is well known that because of the high isolation of orthogonal polarizations, using a cross-polarization antenna scheme symmetrically at both transmitting and receiving ends can provide a higher channel capacity than a conventional antenna polarization scheme. Moreover, part of the existing infrastructure uses cross-polarized antenna elements, which could be used to support both polarization diversity and polarization multiplexing [2]. The channel capacity for MIMO systems has been investigated theoretically [35]. In addition, polarization diversity that employs both vertical and horizontal polarizations can make the fading correlation between copolarized and cross-polarized channels sufficiently low, regardless of the antenna spacing [68].

Recently, many researchers have examined multiple polarizations for an MIMO antenna system. Andrews et al. [9] discussed a wireless MIMO link that provides six independent signals with three electric dipoles and three magnetic dipoles to take full rank of the channel. They assumed an antenna model with ideal polarizations and a rich scattering environment and confirmed their analysis using three orthogonal electric dipoles. Li and Yu [10] compared the MIMO correlation coefficient properties and channel capacities, and displayed different polarization combinations are an efficient way for enhancing channel capacity. Svantesson [11] studied the effect of multipath angular spread and antenna radiation patterns on the channel capacity and showed that the capacity increase is due to a combination of polarization and pattern diversity. Weichselberger and Özcelik [12]. introduced a formulation for narrowband MIMO channel model characterized by eigenvectors of covariance matrices and coupling matrices to fulfill the deficiencies. The structure of formulation resembles virtual channel representation [13], in which steering vectors and virtual channel representation matrix are involved.

On the other hand, one way of reducing the level of correlation between antenna elements is to provide sufficient interelement spacing. However, for the practical application of MIMO systems to a wireless local area network, the spacing between adjacent antenna elements cannot be too large. Given this constraint, the most probable solution is a dense MIMO array with nonnegligible correlation. In addition, the use of polarization diversity may be a solution for obtaining a more compact antenna array layout since another diversity dimension can be provided to the MIMO radio channel. Therefore, modeling the polarization diversity technique is an interesting topic of study in MIMO radio channels, and it necessitates the construction of realistic MIMO radio channel models featuring both space and polarization diversity [1416].

Dissimilar channel environments and different incoming wave angular distributions result in a distinct spatial correlation among antenna array elements, which are the main parameters that affect MIMO channel characteristics [17]. To investigate the system performance and correlation properties, we use various antenna polarization combinations in different environments for the same antenna element spacing.

The validation of the model is supported by measurement results. Using two measurement setups having several transmitting and receiving elements, four polarization schemes have been investigated in several different LOS and NLOS rooms as indoor environments. The parameters of the MIMO model are extracted from the measurement data and fed to the model to compare simulation results with the measurement results.

In this paper, we propose an analysis method, which utilizes a correlation coefficient for both transmission and reception and the average receiving power, to calculate the MIMO channel capacity in polarization systems. In addition, we derive the optimal received-power matrix for both additive white Gaussian noise (AWGN) and multipath fading channels. Furthermore, measurement results show that the proposed analysis method outperforms traditional method in polarization systems.

The remainder of this paper is organized as follows. In Section 2, the experiments on the channel measurement setup and antenna configuration are described in detail. In Section 3, we discuss the concepts of the MIMO system, including the singular-value decomposition (SVD) method for finding channel gain, subchannel correlations, and channel capacity. The normalization method for comparing the effect of various polarization schemes is mentioned in this section. In Section 4, the channel capacity and correlation coefficient of the measured channel are calculated, and some meaningful results are revealed. Finally, we present conclusions in Section 5.

2. Measurement Setup

MIMO measurements with an setup, where and are the numbers of elements at the transmitter () and receiver (), respectively, are conducted with . A simplified sketch of the MIMO measurement setup is shown in Figure 1.

246348.fig.001
Figure 1: MIMO measurement system.

To measure the response of the single-input/single-output (SISO) channel, a two-port vector network analyzer is used. Both and are attached to a 1-to-3 switch with switching times of 15 s and 5 s for and , respectively. Furthermore, a signal generator is required to control the status of the switches. A personal computer is used to acquire the measured data and to remotely control the signal generator so that the correct clock signal is sent through the general-purpose interface bus. During the measurement, one of the transmitting antennas is switched on, and three receiving antennas are activated one after the other. Next, the second transmitting antenna is switched on, and the control circuit at the receiver activates the receiving antennas successively. The process is repeated for all the transmitting antennas. We assume that the indoor channel is quasistatic so that we can measure each radio link by controlling the switching circuits. The measuring time for each link is 5 s, so it takes about 45 s to complete a MIMO channel matrix measurement. The parameters of our measurement are listed in Table 1.

tab1
Table 1: Parameters of measurement.
2.1. Measurement Environment

The measurement sites are located in Electrical Engineering Building II at the National Taiwan University campus. Two dissimilar environments are selected for their different characteristics. Figure 2 shows the equipment and setting for indoor channel measurements. The size of the room is about . The transmitting antenna array is fixed at the center of the classroom, and the receiver is placed at eight locations. These locations lie on the circumference of a circle with a radius of 3 m that is centered at the center of the room; this means that all eight locations are equidistant from the transmitter. Each receiving antenna array is moved around a 1-m2 area to obtain channel responses for sixteen positions. As a result, we obtain a total of 128 MIMO channel matrices at the end of the experiment. Figure 3 shows the placement of equipment in two rooms on the same floor for the NLOS environment. There is a concrete wall with a thickness of 10 cm between the two rooms. In Figure 3, the locations of the receiving and transmitting antenna arrays are denoted as “X” and “●”, respectively. The height of the antenna arrays at and is 1.2 m. The measurements are made in the absence of bodily movement to investigate time-stationary picocell environments.

246348.fig.002
Figure 2: The experimental layout for indoor LOS environment.
246348.fig.003
Figure 3: The experimental layout for indoor NLOS environment.

For both environments, a large number of receiving locations are used to gather more statistical information for the environment. In addition, to remove the effect of attenuation or loss resulting from equipment, we calibrate the data by dividing the measured frequency response by that measured for 1-m - separation.

2.2. Antenna Configuration

In this study, we investigate different polarization configurations: vertical (V), horizontal (H), and slanting (Y) polarizations (angle of slant = ). The relative orientations of the transmitting antenna array and receiving antenna array are shown in Figure 4. Whenever the receiving antenna array is moved, the same relative orientation between the transmitting antenna array and the receiving antenna is maintained for the same polarization type. Figure 5 shows photos of the four different polarization combinations, in which each antenna has an omnidirectional pattern and a 5-dBi antenna gain. At both transmitting and receiving ends, the adjacent antennas are spaced half wavelengths apart. Transmission is carried out at 2.4 GHz center frequency with a frequency span of 500 MHz over 801 individual frequencies (points).

fig4
Figure 4: Four antenna configurations which are denoted as (a) , (b) , (c) , and (d) .
fig5
Figure 5: Photos of the four different polarization combinations denoted as (a) , (b) , (c) , and (d) in practical measurement environment.

3. Narrowband MIMO Channel Model Analysis

In this section, we analyze the measurement results obtained for the MIMO radio channels used in narrowband systems. Although our measurement system is wideband, we only use data in the frequency range 2.35–2.45 GHz, corresponding to 161 frequency points in the central region, because the antenna possesses a narrowband operational bandwidth. Thus, our discussion deals with narrowband analysis. Moreover, to analyze the capacity and eigenvalue of the measured data, we first need to normalize the data. For each location, let denote the channel matrix measured at the th position ( = 1–16) and th frequency point.

3.1. Normalization for Different Polarizations

We discuss the relationship between polarization and channel capacity. The method for normalizing the configuration is defined as follows: where is an channel matrix at one position and one frequency point. is the square of the Frobenius norm of . Moreover, denotes the th position (where the total number of positions in one location is denoted as ) in each location and denotes the th frequency point (where the total number of frequency points is denoted as ).

The normalization method is different for the case. We use the case for the same location as the baseline. We first calculate the mean power of all antenna elements in the case over all the frequency points and positions in the same location. We then normalize all the channel matrices by the square root of the mean power for the case: where and denote the channel matrices at the th position and th frequency point for the case and case, respectively. The normalization methods for the other cases are the same as the method used for normalization.

3.2. Eigenvalue and Capacity for Equal Power

In analyzing the capacity of an MIMO system, transmitting elements and receiving elements are assumed, and the received signal is expressed as where is an channel matrix, is the transmitted signal, and is the AWGN.

For a fixed channel realization, the channel capacity has the following constraints: the transmitter has no channel state information and the transmitted power is equally allocated to each transmitting element. Therefore, the capacity can be expressed as [1] where det denotes the determinant of a matrix, and denotes a Hermitian transpose. is an identity matrix and is the average received signal-to-noise ratio.

To investigate the characteristics of , we can perform an SVD of to diagonalize and determine the eigenvalues. The SVD expansion of any matrix is written as [18] where and are unitary matrices, which means that is nonnegative and diagonal with entries specified by where diag is a vector consisting of the diagonal elements of are the nonzero eigenvalues of ; .

The columns of and are the eigenvectors of and [19], respectively. The SVD (8) shows that the channel matrix can be diagonalized to a number of independent orthogonal subchannels, where the power gain of the th channel is [20], Thus, (7) can be rewritten as [20] where are nonzero eigenvalues of in (10).

Therefore, the channel capacity is affected not only by the maximal value of but also by the minimal value of. We can define the condition number as

If a channel has a low condition number, then its correlation is low, its diversity is high, and thus its capacity is high. The channel is then said to be “well conditioned.” Otherwise, the channel is referred to as being “ill conditioned” [21].

3.3. Subchannel Correlation

We modify the method of analysis of channel capacity for different polarization schemes. A different polarization configuration has a different received power owing to cross-polarization. Therefore, we consider the effect of the cross-polarization ratio (XPR) in different polarization schemes. Moreover, the use of the XPR for polarization antenna systems has been previously studied.

The XPR is defined as the ratio of the copolarized average received power to the cross-polarized average received power. The XPR has been used in the capacity analysis of different polarization combinations [22].

To obtain XPRs for different polarization cases, we normalize the scheme so that the average normalized power matrix is equal to one. The average received power matrix is defined as From (13), we obtain the matrix for different polarization schemes. Next, we obtain as the normalized power matrix and use the case in the same location as the baseline. Other polarization schemes refer to the scheme. Finally, the proposed MIMO channel matrix is written as where is an matrix containing independent and identically distributed (0,1) elements, and and are transmitting and receiving correlation matrices respectively. Then, (or ) is the correlation coefficient of the th and th transmitting (or receiving) antennas.

For a measured MIMO channel matrix as denoted in (3), we can thus define the transmitter and receiver correlation coefficients as Subsequently, the correlation coefficients are calculated as follows:

4. Measurement Results

Our experiment focuses on indoor measurements. We intend to gain a deeper understanding of the effects of different polarizations from our experimental results.

The first and the second MIMO scenario performances are illustrated in terms of the cumulative distribution functions (CDFs) of their eigenvalues, as shown in Figures 6 and 7, in which three subchannel gains are plotted as straight lines. We find that regardless of whether the LOS or NLOS environment is used, the largest eigenvalue () in the schemes is smaller than of other polarization schemes owing to the significantly lower channel gain. Moreover, is higher for the scheme than for the , , and schemes; thus, the channel capacity of the combination in indoor environments is the smallest among the capacities of all the schemes. However, all the eigenvalues are closer to each other in the case than in the other polarization cases. Hence, the channel is the most well conditioned.

fig6
Figure 6: The plots of CDF versus eigenvalues when the polarization system employing (a) , (b) , (c) , and (d) cases in LOS indoor environment.
fig7
Figure 7: The plots of CDF versus eigenvalues when the polarization system employing (a) (b) (c) (d) cases in NLOS indoor environment.

Let us consider the channel correlation matrices listed in Tables 2 and 3. We observe that the copolarized cases have greater correlation coefficients than the other polarization cases. In addition, if the polarization between the transmitting antenna and receiving antenna is orthogonal, then the correlation coefficient is the lowest among all the polarized pairs under the same measurement environment. From this viewpoint, we study the polarization scheme. The correlation coefficient () is smaller than () and () owing to orthogonal polarization transmission. The same relationship holds for the and polarization schemes. Consequently, the correlation coefficients of the nondiagonal matrix elements are lower; the channel of the MIMO systems is more uncorrelated; the channel has a higher capacity.

tab2
Table 2: Correlation matrix for different polarization schemes in LOS environment.
tab3
Table 3: Correlation matrix for different polarization schemes in NLOS environment.

We also examine the correlation for direct path power which defines the first path in time domain among the subchannels in the four polarization schemes. As shown in Tables 4 and 5, if the polarization between transmitting antenna and receiving antenna is orthogonal, then the correlation coefficient of the direct path power corresponds to the lowest subchannel gain among all the polarized pairs. All matrix components of are 1 because of the direct path powers are close to each other. In explaining the effects of polarization combinations, we use the notation P-Q, where P and Q denote the transmitting antenna and receiving antenna, respectively. Considering the normalized power matrix () in Tables 6 and 7, we note that the V-V copolarized case and the Y-Y (“V”-“H”) cross-polarized case are very different. The cross-polarized combination shows three subchannels with considerably high gain and remaining subchannels with a lower gain; in contrast, the copolarized combination shows high gains for all the subchannels, with the gains being close to each other, although the scheme provides good isolation as the antenna elements are orthogonal for the transmission and reception pairs. Nevertheless, as we observe normalized received power matrices, it is almost the loss of all the subchannel gains. This phenomenon shows the benefit of adopting orthogonal polarizations. Table 8 lists the gains of different polarization schemes; these gains help to clarify the above phenomenon. As a result, in using orthogonal polarization, we lose part of the channel gain but achieve greater isolation. The same results are found for the scheme.

tab4
Table 4: Correlation matrix of direct path power analysis for different polarization schemes in LOS environment.
tab5
Table 5: Correlation matrix of direct path power analysis for different polarization schemes in NLOS environment.
tab6
Table 6: Normalized power for different polarization schemes in indoor LOS environment.
tab7
Table 7: Normalized power for different polarization schemes in indoor NLOS environment.
tab8
Table 8: Average normalized power for various polarization schemes.

We now present results from the analysis of the measurement data in CDFs. The channel capacities of the four polarization schemes obtained from measurement and simulation analyses are, respectively, shown as straight and dotted lines in Figure 8 and Figure 9. The four polarization schemes—, , , and —are indicated by the symbols, , , ○, and □, respectively. The signal-to-noise ratio is set at 10 dB to plot the CDF of the capacity statistical properties. From the preliminary analysis result shown in Figure 8, we find that the structure has the best capacity among all the investigated antenna array configurations. Moreover, when the normalized power matrix is considered with the analysis of channel capacity, the scheme has the worst channel capacity among the polarization combinations.

246348.fig.008
Figure 8: The plots of CDF versus capacities obtained by measurement analysis (solid line), proposed simulation analysis (dotted line) for dissimilar polarization systems in indoor LOS environment.
246348.fig.009
Figure 9: The plots of CDF versus capacities obtained by measurement analysis (solid line), proposed simulation analysis (dotted line) for dissimilar polarization systems in indoor NLOS environment.

By observing the CDF of capacity analysis under different conditions, it is found that the simulation analysis method, which makes use of (15), is close to a realistic measurement environment. This method is useful for application to polarization schemes. As seen in Figure 8, using the scheme in an indoor environment provides higher channel capacity owing to antenna isolation and greater total received power. As a result, we need to consider both isolation and correlation properties for describing an MIMO system that utilizes orthogonal polarized antennas. Using only one property is insufficient in determining the system performance. Therefore, the correlation property is not a sufficient measure when considering parallel transmission with orthogonal polarized schemes.

For the NLOS condition, as shown in Figure 9, we find that the scheme has a higher capacity than the scheme. As for the analysis result in Figure 7, has the same power level for both the and schemes, but the other eigenvalues in the scheme are closer to each other. Therefore, we have verified our proposed simulation analysis method and found it suitable for approximating polarization transmission and reception systems. Consequently, in attaining higher channel capacity in MIMO systems, there is a tradeoff between antenna isolation and total received power in each transceiver structure.

5. Conclusion

We introduced the concept of the MIMO system and factors that determined the channel capacity, including eigenvalues and correlation coefficients. We studied the effect of space and polarization diversity on the MIMO system in detail. Furthermore, we designed a series of experiments to determine the correlation between antenna polarization and the channel characteristics. From measurement data, the channel capacity and correlation coefficient were used for explaining the effects of various polarization schemes in the MIMO channels. Moreover, the concepts of normalizing the received power and the polarization effect were described in modifying the numerical analysis of the polarized channel capacity. In addition, we found that the performance of an MIMO system exploiting a copolarized antenna combination can be described simply using spatial correlation properties, but when adopting a cross-polarized antenna combination, both the isolation and correlation properties were needed to fully describe the system performance. Consequently, we found and verified the proposed algorithm of (15) which was close to realistic environment in polarization systems.

Acknowledgment

This work was supported by the National Science Council, Taiwan, under the Grant of NSC95-2219-E-002-003 and NSC97-2218-E027-007.

References

  1. G. J. Foschini, “Layered space-time architecture for wireless communication in a fading environment when using multi-element antennas,” Bell Labs Technical Journal, vol. 1, no. 2, pp. 41–59, 1996. View at Scopus
  2. K. Yu, M. Bengtsson, B. Ottersten, D. McNamara, P. Karlsson, and M. Beach, “Modeling of wide-band MIMO radio channels based on NLoS indoor measurements,” IEEE Transactions on Vehicular Technology, vol. 53, no. 3, pp. 655–665, 2004. View at Publisher · View at Google Scholar · View at Scopus
  3. L. Dong, H. Choo, R. W. Heath, and H. Ling, “Simulation of MIMO channel capacity with antenna polarization diversity,” IEEE Transactions on Wireless Communications, vol. 4, no. 4, pp. 1869–1872, 2005. View at Publisher · View at Google Scholar · View at Scopus
  4. 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 · View at Google Scholar · View at Scopus
  5. D.-S. Shiu, G. J. Foschini, M. J. Gans, and J. M. Kahn, “Fading corre- lation and its effect on the capacity of multielement antenna system,” IEEE Transactions on Communications, vol. 48, pp. 502–513, 2000.
  6. J. P. Kermoal , et al., “Experimental investigation of the joint spatial and polarisation diversity for MIMO radio channel,” in Proceedings of the 4th International Symposium on Wireless Personal Multimedia Communication (WPMC '01), pp. 147–152, September 2001.
  7. P. Kyritsi, D. C. Cox, R. A. Valenzuela, and P. W. Wolniansky, “Effect of antenna polarization on the capacity of a multiple element system in an indoor environment,” IEEE Journal on Selected Areas in Communications, vol. 20, no. 6, pp. 1227–1239, 2002. View at Publisher · View at Google Scholar · View at Scopus
  8. R. U. Nabar, H. Bölcskei, V. Erceg, D. Gesbert, and A. J. Paulraj, “Performance of multiantenna signaling techniques in the presence of polarization diversity,” IEEE Transactions on Signal Processing, vol. 50, no. 10, pp. 2553–2562, 2002. View at Publisher · View at Google Scholar · View at Scopus
  9. M. R. Andrews, P. P. Mitra, and R. Decarvalho, “Tripling the capacity of wireless communications using electromagnetic polarization,” Nature, vol. 409, no. 6818, pp. 316–318, 2001. View at Publisher · View at Google Scholar · View at Scopus
  10. H. J. Li and C. H. Yu, “Correlation properties and capacity of antenna polarization combinations for MIMO radio channel,” in Proceedings of the IEEE International Antennas and Propagation Symposium, pp. 503–506, San Antonio, Tex, USA, June 2003. View at Scopus
  11. T. Svantesson, “On capacity and correlation of multi-antenna systems employing multiple polarizations,” in Proceedings of the IEEE Antennas and Propagation Society International Symposium, pp. 202–205, San Antonio, Tex, USA, June 2002. View at Scopus
  12. W. Weichselberger and H. Özcelik, “A novel stochastic MIMO channel model and its physical interpretation,” in Proceedings of the 6th International Symposium on Wireless Personal Multimedia Communications (WPMC '03), Yokosuka, Japan, October 2003.
  13. A. M. Sayeed, “Deconstructing multiantenna fading channels,” IEEE Transactions on Signal Processing, vol. 50, no. 10, pp. 2563–2579, 2002. View at Publisher · View at Google Scholar · View at Scopus
  14. S. R. Saunders, Antennas and Propagation for Wireless Communication systems, John Wiley & Sons, New York, NY, USA, 1999.
  15. V. Jungnukel, , V. Pohl, and C. von Helmolt, “Capacity of MIMO Systems with closely spaced antennas,” in Proceedings of the IEEE Semiannual Vehicular Technology Conference (VTC '02), chapter 3.2., Birmingham, UK, 2002.
  16. V. R. Anreddy and M. A. Ingram, “Capacity of measured ricean and rayleigh indoor MIMO channels at 2.4 GHz with polarization and spatial diversity,” in Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '06), vol. 2, pp. 946–951, April 2006. View at Scopus
  17. J. P. Kermoal, L. Schumacher, P. E. Mogensen, and K. I. Pedersen, “Experimental investigation of correlation properties of MIMO radio channels for indoor picocell scenarios,” in Proceedings of the 52nd Vehicular Technology Conference (VTC '00), vol. 1, pp. 14–21, September 2000. View at Scopus
  18. V. Erceg, H. Sampath, and S. Catreux-Erceg, “Dual-polarization versus single-polarization MIMO channel measurement results and modeling,” IEEE Transactions on Wireless Communications, vol. 5, no. 1, pp. 28–33, 2006. View at Publisher · View at Google Scholar · View at Scopus
  19. I. E. Telatar, “Capacity of multi-antenna gaussian channels,” Tech. Rep. BL0112170-950615-07TM, Office Official Pubilications of the European Communities, AT&T Bell Labs., 1995.
  20. D. Gesbert, M. Shafi, D. S. Shiu, P. J. Smith, and A. Naguib, “From theory to practice: an overview of MIMO space-time coded wireless systems,” IEEE Journal on Selected Areas in Communications, vol. 21, no. 3, pp. 281–302, 2003. View at Publisher · View at Google Scholar · View at Scopus
  21. T. S. Rappaport, Wireless Communications: Principles and Practices, Prentice Hall, Upper Saddle River, NJ, USA, 2nd edition, 2002.
  22. P. Soma, D. S. Baum, V. Erceg, R. Krishnamoorthy, and A. J. Paulraj, “Analysis and modeling of Multiple-Input Multiple-Output (MIMO) radio channel based on outdoor measurements conducted at 2.5 GHz for fixed BWA applications,” in Proceeding of IEEE International Conference on Communications, pp. 272–276, May 2002. View at Scopus