- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Antennas and Propagation
Volume 2012 (2012), Article ID 625165, 8 pages
Joint Gain/Phase and Mutual Coupling Array Calibration Technique with Single Calibrating Source
1Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
2Temasek Laboratories, NTU, Singapore
3School of Electrical and Electronic Engineering, NTU, Singapore
Received 9 July 2012; Revised 19 September 2012; Accepted 20 September 2012
Academic Editor: Hon Tat Hui
Copyright © 2012 Wei Chen 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.
An iterative-based method for joint gain/phase and mutual coupling array calibration is proposed in this paper. It estimates the array gain/phase and mutual coupling coefficients with a set of simultaneous equations formed by using the beam pattern property of the array. Only one calibrating source with known direction is requiblue to obtain the unique estimate. The effectiveness of this approach is illustrated by simulation results and by experimental data collected with an antenna array operating in high-frequency radio band.
Smart antenna technology is widely used in mobile communication. In a Smart Antenna System (SAS), an antenna array is used in the base station to transmit and receive the mobile’s signal directionally so that interferences are rejected and the system’s channel capacity is enhanced. But in many practical situations, the performance of the antenna array is degraded due to gain/phase imbalances between the array sensors, unknown mutual coupling coefficients, and unknown sensor locations. This implies that an estimate of the true array manifold is necessary before the processing of the array data can be carried out. In fact, most of the array processing algorithms available in the literature, such as the MUSIC direction finding algorithm  and the minimum-variance distortionless response (MVDR) adaptive beamforming algorithm , are devised under the assumption that the array manifold of the array processing system is known exactly. Therefore, antenna array calibration is indispensable for estimating the array manifold accurately. Antenna array calibration mainly compensates for the following factors, namely, sensor location uncertainty, gain/phase mismatch between array sensors, and mutual coupling effects between array sensors.
There are several methods for array calibration. The first type of methods uses signals from known sources and known directions for calibration. Friedlander and Weiss  presented a method for estimating the gain and phase of each sensor and the mutual coupling in the receiving array, but it relies on the calibrating signals from several sources (although it does not require to know their actual impinging directions). The second type of methods injects an equal phase signal to all channels of the antenna array. In , Ertel et al. compensated the gain and phase imbalance in the uplink system of an antenna array but without considering the effects of mutual coupling. The third type of methods uses network analyzer for calibration. Dandekar et al.  gave a mathematical model of smart antenna and calibrated the smart antenna array by network analyzer measurements and computational electromagnetic (CEM) simulations. But as an off-line method, it has limitations in compensating mutual coupling effects which change due to thermal effects, the environment around the sensor array (e.g., the effect of metal objects near an antenna array on its beam pattern), and other factors. Therefore, this paper is motivated to attempt to suggest an alternative calibration technique to estimate the array manifold parameters. Moreover, recent research on array calibration can be found in [6, 7]. It is worth noting that the calibration methods discussed here refers to the calibration method for the antenna array when it is in receive mode . As widely reported in references [9–11], there exists other alternative means for dealing with the mutual coupling that belongs to mutual coupling compensation techniques.
The basic idea of the technique presented in this paper is to form a set of linear equations using the null and peak characteristics of the beam pattern, from which the errors or the true steering vector that contains the manifold information of the array sensors can be estimated. The proposed technique only uses one calibrating source with known direction to obtain the unique estimate. This paper is organized as follows. Section 2 contains the signal model. Section 3 describes the gain/phase calibration technique. Section 4 describes the joint gain/phase and mutual coupling calibration technique. Section 5 gives some simulation results to show the performance of the proposed technique, as well as its application to experimental data collected with an antenna array operating in high-frequency band and Section 6 concludes the paper.
2. Signal Model
Consider an antenna array having sensors and the azimuth direction of arrival of (measublue with respect to the array normal axis for the case of linear array), and let . In the ideal situation, the received signal of the array can be described as where and denotes the received signal at th element of the array, is the transmitted signal, is the steering vector for the direction , and is the noise vector.
However, under the practical conditions, the received signal of the array may be described as where is the mutual coupling complex matrix, is the real diagonal matrix whose elements represent the sensor gains, and is an complex diagonal matrix whose elements represent the sensor phases .
Let , where the complex diagonal matrix represents the sensor gain/phase matrix. From (3), matrix is given by For a linear uniform array, matrix is a Toeplitz matrix. For a uniform circular array, matrix is a circulant matrix. This paper uses the properties of Toeplitz matrix and circulant matrix (see [3, Lemmas 1, 2, and 3]) to calculate the mutual coupling matrix as well as the sensor gain/phase matrix . It is worth noting that the assumption that the coupling matrix exhibits the Toeplitz and circulant matrix structure also implies that the coupling between two adjacent array elements is independent of the presence of other nonadjacent array elements. This also means the coupling matrix for an edge element is identical to that of an element located in the center of the array. This condition is only true under the minimum scattering array principle .
3. Estimating the Gain/Phase Vector
3.1. The Characteristics of Beam Patterns
For the ideal condition expressed by (1), the covariance matrix of received signals is where represents the expectation operation and the superscript “” represents the conjugate transpose operation.
The characteristics of the beam patterns are discussed as follows for the uniform linear array and the uniform circular array, respectively.
3.1.1. Uniform Linear Array
Decomposing the covariance matrix , we can obtain the standard beam pattern that is expressed as where is the normalized eigenvector of the signal subspace, which is expressed as ; is the normalized steering vector, which has the expression ; is the intersensor spacing and is the wavelength.
Moreover, in the -space the beam pattern can be expressed as
For the purpose of illustration, Figure 1 shows the beam pattern of a 10-element uniform linear array with half wavelength intersensor spacing and the direction of arrival is . We can see that the beam pattern has one peak and nulls. The peak of the pattern occurs when while the nulls (more detailed derivation and explanation on the location of the nulls in the ULA beam pattern is available ) occur when
3.1.2. Uniform Circular Array
Decomposing the covariance matrix , we can get the uniform excitation beam pattern that is expressed by where is the uniformly eigenvector of the signal subspace, which is expressed as ; is the normalized steering vector which is expressed as ; and is the wavelength. denotes the circular radius of the array geometry.
Since the polar angle is not consideblue here, we can simplify the expression in . In the -space, the beam pattern is expressed as where is the Bessel function of the first kind of order . The variables and have the following expressions, respectively:
And in (11), the term is the principle term while the other terms are the residuals. The main lobe behavior is adequately described by the term . And if the interelement spacing on the arc satisfies the inequality then the residuals will be negligible.
Hence, we can use to describe the characteristics of the standard beam pattern for uniform circular array.
Figure 2 shows the beam pattern of a 10-element uniform circular array with one-third wavelength intersensor spacing on the arc and the looking direction is . We can get the beam pattern value corresponding to any angle by the Bessel function table .
The proposed technique for calibration is based on the knowledge of one calibrating source with known direction. We can use some values of the beam pattern to get a set of simultaneous equations for estimating the gain/phase and mutual coupling matrix. In the following part of this paper, we will use the peak and the nulls of beam pattern for uniform linear array’s calibration. The calibration for uniform circular array can be performed in a similar way as the uniform linear array, by using other values of the beam pattern.
3.2. Gain/Phase Calibration Technique
The true steering vector at any look direction for a gain/phase perturbed array can be expressed as where is the nominal steering vector; whose elements are the coefficients of the gain/phase mismatch; the superscript “” represents the transpose operation.
Moreover, from (8) and (9) we have the following equations Rewriting the above equations in matrix form, we obtain where is an matrix given by and is a vector given by As the nominal steering vectors according to different DOA are independent with each other, matrix is full rank and has an inverse matrix. The gain/phase mismatch vector is given by From (19), under asymptotic condition, the first element of should yield a value of one as the first array element is chosen as a reference element and constrained as one. However, in actual implementation, the data covariance matrix available is not really perfect, especially when the number of snapshots taken to form the covariance matrix is small. In this case, the following constrained minimization formulation would be more appropriate: where denotes the 2-norm operation. Equation (20) yields the following solution: where
4. Jointly Estimating the Gain/Phase Vector and Mutual Coupling Matrix
Once the DOA of the source is known, which means is known, reasonable estimates of and may be obtained by minimizing the cost function The proposed minimization algorithm is based on a two-step procedure. First, given , we minimize over the gain/phase parameters. Then given , we minimize over the mutual coupling matrix components. These minimization steps can be repeated until converges.
(i)Set the iteration counter to zero: . (ii)Select initial values for the gain/phase matrix and initial value for the mutual coupling matrix . Usually the initial values are based on some previous knowledge (e.g., last measublue values or pblueictions based on the idealized model). (iii)Compute the covariance matrix from the discrete time-domain samples of the received signal at the antenna array. (iv)Perform eigen-analysis and construct signal subspace .
4.2. Two-Step Procedure
Step 1 (estimating the gain/phase vector). Fixing the matrix , we now minimize with respect to the gain and phase of each of the sensors. From (23), we have where We can minimize (24) using the constrained minimization formulation expressed by (20). The solution is similar to the formulation (21): where Then, we can compute the gain-phase matrix from the vector given by (26):
Step 2 (estimating the mutual coupling matrix). In this step, we hold the sensor gain-phase fixed and find the mutual coupling matrix that minimizes the cost function . Using the lemmas  depicted in the appendix, we obtain where and is the mutual coupling coefficients vector whose elements are related to the coupling matrix , as expressed by The formulation of is given by Lemma A.3  in the appendix.
We can minimize (29) using the constrained minimization formulation expressed by (20). The solution is similar to the formulation in (21), which is, where Then, we can reconstruct the mutual coupling matrix form the vector given by (32).
4.3. Convergence Check
Here, we compute using the estimated sensor gain-phase and the mutual coupling matrix. If ( is a preset threshold), then we will update the iteration counter and continue with the two-step procedure. Otherwise, we will stop the process. The algorithm performs the iterations until converges. Note that at each step, the cost function should blueuce so that
4.4. Extension to Azimuth-Elevation Case
Although the proposed calibration method is derived for linear geometry array, the extension from azimuth-only (1D) to azimuth-elevation (2D) case is possible by reformulating the azimuth-elevation steering vector into its decoupled form [17, 18]. Using this decoupled formulation, the proposed calibration can be applied individually on the coupling matrix associated with the azimuth as well as elevation components of the steering vector.
5. Simulation and Experimental Results
To illustrate the behavior of the algorithm, consider a narrowband half-wavelength spaced uniform linear array of 10 omnidirectional sensors. Noise components are complex white Gaussian with zero mean. The following parameters are used in the simulations: SNR = 20 dB; Snapshots = 100; the known source’s location: ; coupling coefficient: ; coupling coefficients are independently distributed according to . gains: ; phases: .
In the presence of gain/phase imbalances and unknown mutual coupling coefficients, the ideal beam pattern characteristics of a uniform linear can no longer be observed. This is depicted clearly in Figure 3 for this particular realization. In comparison, we have also plotted the ideal beam pattern computed given the gain/phase imbalances and mutual coupling coefficients, as well as the after-calibration beam pattern given by where and are the estimated coupling matrix and gain/phase imbalances. These results demonstrate that the proposed two-step procedure in Section 4.2 is effective in estimating jointly the gain/phase and coupling matrix that satisfies the beam pattern’s nulls and peak characteristics.
In addition, we have also tested the calibration result for direction finding application by plotting the MUSIC spatial spectrum on the simulated test data that contains three signals impinging from , and .
Figure 4 shows the comparison of MUSIC spatial spectra before calibration, after calibration, and ideal spectrum for joint gain/phase and mutual coupling calibration. It is clear that before calibration the estimation of DOA is wrong, the spectrum is not acuity, and SNR is low. After calibration, the estimation of DOA is correct, the spectrum is acuity and very similar to the nominal MUSIC spectrum.
To demonstrate the statistical efficiency of the proposed calibration, we perform the following Monte Carlo experiments. The 10-element uniform linear array described above is used with the same parameter setting. The data set for the calibration contains only the calibrating source, while that for the direction finding includes three far-field narrow band emitters. The gains and phases as well as coupling matrix coefficients are unchanged. We run 500 realization for each fixed SNR, and repeat for various SNR values ranging from −20 dB to 20 dB in 5 dB step-size. The SNR is calculated as the ratio of the received signal power and the noise power: SNR , where and are the th received signal power and noise power, respectively. The values of the sensor gains and phases were kept constant throughout these simulations. For each realization, we apply MUSIC-based direction finding method to estimate the DOA.
For each SNR, we obtain the average DOA estimation error out of the 500 realizations and use them to compute the estimated root mean square error (RMSE). Figure 5 depicts the RMSE plot versus SNR for the ideal (known gain/phase and mutual coupling coefficients), as well as before and after calibration. It can be seen that the proposed algorithm has significantly improved the accuracy of the DOA estimation through accurate joint estimation of gain/phase and mutual coupling coefficients. In addition, we have also included the Cramer-Rao Bound (CRB) plot of the DOA estimation under the ideal condition. The CRB plot in Figure 5 depicts the statistical lower bound of the DOA estimation under perfect calibration where the coupling matrix is known. While the RMSE plot under ideal case (known coupling matrix) appears to converge to the CRB for high SNR case, the performance of the proposed calibration does not deviate far from the CRB.
Lastly, we verify the calibration method experimentally using the data collected from a high frequency (HF) radio antenna array system. The measurement setup is described in detail in , therefore it will not be repeated here.
We implement the calibration technique to calibrate jointly the array gain/phase and mutual coupling coefficients. Unlike the calibration described in , we only need one calibrating data set from the short-time Fourier transform (STFT) output that consists of single narrowband signal. It is important to note that although it is impossible to conduct a controlled experiment in which only one calibrating signal is observed, it is possible to extract only the narrowband calibrating signal from the received signal using short-time frequency analysis technique. Even in the congested nature of HF band, we have demonstrated that it is possible to extract the calibrating data set.
After estimating jointly the gain/phase and mutual coupling coefficients, we applied the calibration result for plotting the MUSIC spatial spectrum on the data set that contains three narrowband signals. Figure 6 shows the spatial spectrum plots before and after the proposed calibration. Notice that, unlike the spatial spectra shown in [6, Figure 12], the spatial spectra shown here are plotted with respect to where the angle is measublue with respect to the array normal axis.
It can be seen that although the proposed method does not include the array geometry calibration, all three DOAs are estimated effectively as indicated by the peaks in the plot located at −0.9219, 0.1719, and 0.6266, respectively. Furthermore, the estimation results are in agreement with the results obtained using the method in . This can be understood from the fact that the array geometry error is negligible, as shown in array geometry calibration result (see [6, Figure 11(c)]).
In this paper, a calibration technique based on eigen-structure analysis and the beam pattern characteristics has been used to obtain jointly the estimates of gain, phase, and mutual coupling parameters of the smart antenna system. The effectiveness of this method is illustrated by simulation results, as well as experimental verification. Even though the proposed method is developed under the assumption that the calibrating data contains only calibrating source of known direction with additive noise, it does not render the proposed method impractical since this requirement can be easily overcome through prefiltering approach.
Lemma A.1. For any complex vector and any complex diagonal matrix one has where the components of the vector and the matrix are given by
Lemma A.2. For any complex vector and any complex symmetric circulant matrix , one has where the elements of the vector are given by where when is even and when is odd. The matrix is the sum of the following four matrices: where for even and for odd .
Lemma A.3. For any complex vector and any banded complex symmetric Toeplitz matrix , one has where the vector is given by and is the highest superdiagonal that is different from zero. The matrix is given by the sum of the following two matrices
- R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.
- J. Capon, “High-resolution frequency—wavenumber spectrum analysis,” Proceedings of the IEEE, vol. 57, no. 8, pp. 1408–1418, 1969.
- B. Friedlander and A. J. Weiss, “Direction finding in the presence of mutual coupling,” IEEE Transactions on Antennas and Propagation, vol. 39, pp. 273–284, 1991.
- R. B. Ertel, Z. Hu, and J. H. Reed, “Antenna array hardware amplitude and phase compensation using baseband antenna array outputs,” in Proceedings of the IEEE 49th Vehicular Technology Conference (VTC '99), pp. 1759–1763, September 1999.
- K. R. Dandekar, H. Ling, and G. Xu, “Smart antenna array calibration procedure including amplitude and phase mismatch and mutual coupling effects,” in Proceedings of the IEEE International Conference on Personal Wireless Communications, pp. 293–297, December 2000.
- B. P. Ng, J. P. Lie, M. H. Er, and A. Feng, “A practical simple geometry and gain/phase calibration technique for antenna array processing,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 7, pp. 1963–1972, 2009.
- W. Keizer, “Fast and accurate array calibration using a synthetic array approach,” IEEE Transactions on Antennas and Propagation, vol. 59, pp. 4115–4122, 2011.
- H.-S. Lui, H. T. Hui, and M. S. Leong, “A note on the mutual-coupling problems in transmitting and receiving antenna arrays,” IEEE Antennas and Propagation Magazine, vol. 51, no. 5, pp. 171–176, 2009.
- K. M. Pasala and E. M. Friel, “Mutual coupling effects and their reduction in wideband direction of arrival estimation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 30, no. 4, pp. 1116–1122, 1994.
- H. T. Hui, “A practical approach to compensate for the mutual coupling effect in an adaptive dipole array,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 5, pp. 1262–1269, 2004.
- C. K. E. Lau, R. S. Adve, and T. K. Sarkar, “Minimum norm mutual coupling compensation with applications in direction of arrival estimation,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 8, pp. 2034–2041, 2004.
- B. Friedlander, “A sensitivity analysis of the MUSIC algorithm,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 38, no. 10, pp. 1740–1751, 1990.
- W. Kahn and H. Kurss, “Minimum-scattering antennas,” IEEE Transactions on Antennas and Propagation, vol. 13, pp. 671–675, 1965.
- D. H. Johnson and D. E. Dudgeon, Array Signal Processing: Concepts and Techniques, Simon & Schuster, 1992.
- H. L. Van Trees, Optimum Array Processing, Part IV of Detection, Estimation, and Modulation Theory, Wiley-Interscience, New York, NY, USA, 2002.
- I. A. S. Milton Abramowitz, Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables, Dover Publications, New York, NY, USA, 1965.
- B. H. Wang, H. T. Hui, and M. S. Leong, “Decoupled 2D direction of arrival estimation using compact uniform circular arrays in the presence of elevation-dependent mutual coupling,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 3, pp. 747–755, 2010.
- R. Goossens and H. Rogier, “A hybrid UCA-RARE/Root-MUSIC approach for 2-D direction of arrival estimation in uniform circular arrays in the presence of mutual coupling,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 3, pp. 841–849, 2007.