Research Article  Open Access
Blind DirectionofArrival Estimation with Uniform Circular Array in Presence of Mutual Coupling
Abstract
A blind directionofarrival (DOA) estimation algorithm based on the estimation of signal parameters via rotational invariance techniques (ESPRIT) is proposed for a uniform circular array (UCA) when strong electromagnetic mutual coupling is present. First, an updated UCA model with mutual coupling in a discrete Fourier transform (DFT) beam space is deduced, and the new manifold matrix is equal to the product of a centrosymmetric diagonal matrix and a Vandermondestructure matrix. Then we carry out blind DOA estimation through a modified ESPRIT method, thus avoiding the need for spatial angular searching. In addition, two mutual coupling parameter estimation methods are presented after the DOAs have been estimated. Simulation results show that the new algorithm is reliable and effective especially for closely spaced signals.
1. Introduction
Uniform circular array based DOA estimation methods are always attractive because the UCA has a special symmetric structure that provides almost the same resolution ability along the 360° azimuth angle domain. Except for conventional DOA estimation methods such as beam space searching or the Capon method [1], many other methods can be applied, such as the commonly used multiple signal classification (MUSIC) [2] or the wellknown ESPRIT [3] in phase mode space. More effective methods such as UCARBMUSIC and UCAESPRIT [4] have been introduced, and the mapping error reducing method was developed [5].
However, the electromagnetic mutual coupling effect cannot be ignored in a real array. Generally, this effect will severely degrade the performance of the above methods [6]. A classic DOA estimation method based on an iterative search technique was proposed to estimate the DOA and the mutual coupling matrix (MCM) parameters jointly [7], but it has a high computation cost. The rank reduction (RARE) method introduced by [7] was further developed to obtain the DOA estimate for a UCA [8–10], but an angular search is still needed. Moreover, angular ambiguities are challenges for RAREbased blind methods. For example, two closely spaced signals cannot be differentiated according to the UCARARE [8] spectrum because there is a spurious peak at the middle of the two true angles.
In this paper, we design a new modified ESPRIT method for UCA to estimate the azimuth angle when mutual coupling is present. The proposed method is blind to mutual coupling and can completely avoid angular searching. Reliable DOA estimates can be obtained, especially for closely spaced signals. In addition, we will introduce two methods to estimate the MCM parameters once the DOA values are calculated.
2. UCA Model with Mutual Coupling
Suppose element UCA has a radius (Figure 1). All of the antenna elements are identical, and there are farfield narrow signals impinging from which are the parameters to be estimated. The snapshot can be written as is the total sampling number and is the manifold matrix:where . is the wavelength of the signals, and is the manifold vector for the th signal. is the common directivity factor of the antenna elements. Since all the antenna elements are identical, can be normalized as 1, and we use effective SNR which has already included the directivity factor. is the signal sample vector: is the MCM, which is a symmetric circular Toeplitz matrix with the form is the flooring function and is normalized as 1. For ideal dipole antenna array (Figure 1), induced EMF (Electromotive Force) method can give a close form of selfimpedance and mutual impedance and thus MCM can be calculated according to Gupta and Ksienski’s formulation [6].
Suppose the signals are uncorrelated and is white Gaussian noise:
3. Algorithm
3.1. UCA Model in DFT Beam Space
First we introduce the DFT of the MCM for UCA.
Lemma 1 (see [11–13]). If is a circular matrix with its first column vector as and is a Fourier matrixwhere , then is a diagonal matrix with entries as ’s eigenvalues:in which . For the Fourier matrix If is also symmetric, this means thatThen we can rewrite asand we havewhere . From (8) and (12), we get a linear equation which and should satisfy for even :whereTherefore, if we determine the estimate of , then we can obtain the estimate of the mutual coupling parameters by (13). In addition, there are similar equations where is odd.
According to [4], we setwhere is the number of excited phase modes. We define another Fourier matrix , which is different from (6):Then the snapshot in DFT beam space isif we write aswhere (see (11), (12), and (18))In addition, we can obtain [4]with is order firstkind Bessel function. According to its propertywe getwhere is the updated manifold matrix with a Vandermonde structureFinally, we obtain the snapshot in DFT beam spacein which
3.2. DOA Estimation and MCM Calculation
We carry out the eigendecomposition on the covariance matrix and obtain the signal subspace , which consists of eigenvectors corresponding to maximum eigenvalues. Select the first rows of as and the second rows as . Select the first rows of as and the second rows as . Select the first diagonal matrix of as and the second diagonal matrix as . Use the same notations for other matrices , , , and . Then we havewhere is a nonsingular matrix. Thus we getwith
Since we use sampled data, we should replace with and define the object function as
We use the solution to the above equations from [14, 15]. Considerwithwhere and is the Hadamard product. If is the eigenvector of corresponding to its minimum eigenvalue, then we have
Equation (37) indicates that there is a phase ambiguity for , and thus it will introduce a ambiguity to the DOA values through the eigenvalues of . However, this ambiguity can be cleared by comparing the RARE spectrum [9] values on and . We mark the estimated vector as without ambiguity.
Finally, we can determine the DOA estimates through the eigenvalues of . Following is the detailed procedure of the blind method for DOA estimation of :(1)Calculate the sample covariance matrix and do eigendecomposition. Get the estimated signal subspace and noise subspace .(2)Get the updated signal subspace . Select the first rows as and the second rows as :(3)Estimate the vector and fitting matrix according to (34)–(36). Carry out eigendecomposition on and get the DOA estimates .(4)Compare the blind RARE spectral values [9] on and and clear the ambiguity. Output the final DOA estimates.
We label the above mentioned method as “Blindm1half” because we only select rows of . In addition, we can select the first and the second rows from as and . This method is labeled “Blindm1full.” Furthermore, we can select the first and the third rows from as and or select the first and the third rows from as and . We label the two methods as “Blindm2half” and “Blindm2full,” respectively. We will carry out simulations on these four methods in Section 4.
Now we can estimate the MCM parameters once we get .
Method 1. If , then, according to (32) and (13), we getwhere , , and , which means that consists of the first rows and columns of (see (15)). In (40), is a nonzero parameter, but it will be cleared by the normalization operation.
Method 1 can only estimate mutual coupling parameters.
Method 2. In addition, we can obtain a more accurate estimate of by the method introduced in [7]. Considerwhere . If is the eigenvector of corresponding to its minimum eigenvalue, thenWe give a simpler expression of than that in [7, 9]. Considerwhere is a function that shifts by times circularly. If , then it shifts in the counter direction.
4. Simulations
Consider a 16element halfwave dipole antenna UCA with and . Set . The mutual coupling vector is listed in Table 1. The theoretical parameters are calculated according to Gupta and Ksienski’s formulation [6], and small values are treated as . Real parameters should be estimated by array calibration method.

Two signals are impinging from 35° and 45° with the same signal noise ratio (SNR). We apply the classic MUSIC, blind RARE [9], blind RRARE [10], and UCARARE [8] methods to obtain azimuth estimates. The results are shown in Figure 2.
It shows that the MUSIC method and RAREbased blind methods cannot differentiate these two signals because RAREbased methods will introduce spurious estimates. We should obtain the initial estimates from the MUSIC spectrum to start the iterative method [7], but it is difficult to find two different spectral peaks from the spectrum of “MUSIC without MCM.”
We apply the proposed blind method to the above example. The estimated average DOA absolute bias and root mean square error (RMSE) versus SNR are illustrated in Figures 3–6 (, 100 trials). Method in [7] () and the CramérRao bound (CRB) with a known MCM are also presented [16].
This shows that all the four proposed methods can give satisfactory estimates and that the tendency of the RMSE is the same as the CRB with an increase of SNR. Method Blindm1full and method Blindm2full are more effective than the other two methods because more rows of are involved when we calculate the fitting matrix (see (35)). A comparison of simulations versus sampling number will give similar results. Method in [7] gives a biased estimate.
Table 2 lists the MCM parameter estimates based on the two methods introduced in Section 3.
This shows that Method 2 can give more accurate results.
5. Conclusion
In DFT beam space, we utilize a modified ESPRIT algorithm to obtain a reliable DOA estimate when there is a severe mutual coupling effect. The new blind method is efficient because it avoids searching for the spectral peaks. For closely spaced signals, neither the classic MUSIC nor RAREbased methods provide a good estimate, while the proposed new method can produce an accurate estimate. Moreover, these direction estimates can be used for further MCM parameter estimation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work is supported by the Chongqing Academician Fund under Grant cstc2014yykfys90001, by the Project of Chinese Academy of Engineering under Grant 2014XX05, by the National Natural Science Foundation of China under Grant 61501068, by the Project of Equipment PreResearch of the Ministry of Education of China under Grant 62501040217, and by the Fundamental Research Funds for the Central Universities under Grant 106112013CDJZR165502.
References
 J. Capon, “Highresolution frequencywavenumber spectrum analysis,” Proceedings of the IEEE, vol. 57, no. 8, pp. 1408–1418, 1969. View at: Google Scholar
 R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986. View at: Publisher Site  Google Scholar
 R. Roy and T. Kailath, “ESPRITestimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989. View at: Publisher Site  Google Scholar
 C. P. Mathews and M. D. Zoltowski, “Eigenstructure techniques for 2D angle estimation with uniform circular arrays,” IEEE Transactions on Signal Processing, vol. 42, no. 9, pp. 2395–2407, 1994. View at: Publisher Site  Google Scholar
 F. Belloni and V. Koivunen, “Beamspace transform for UCA: error analysis and bias reduction,” IEEE Transactions on Signal Processing, vol. 54, no. 8, pp. 3078–3089, 2006. View at: Publisher Site  Google Scholar
 I. J. Gupta and A. A. Ksienski, “Effect of mutual coupling on the performance of adaptive arrays,” IEEE Transactions on Antennas and Propagation, vol. 31, no. 5, pp. 785–791, 1983. View at: Google Scholar
 B. Friedlander and A. J. Weiss, “Direction finding in the presence of mutual coupling,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 3, pp. 273–284, 1991. View at: Google Scholar
 M. Pesavento and J. F. Böhme, “Direction of arrival estimation in uniform circular arrays composed of directional elements,” in Proceedings of the IEEE Sensor Array and Multichannel Signal Processing Workshop Proceedings, pp. 503–507, Rosslyn, Va, USA, August 2002. View at: Publisher Site  Google Scholar
 M. Lin and L. X. Yang, “Blind calibration and DOA estimation with uniform circular arrays in the presence of mutual coupling,” IEEE Antennas and Wireless Propagation Letters, vol. 5, no. 1, pp. 315–318, 2006. View at: Publisher Site  Google Scholar
 D. Jisheng, B. Xu, H. Nan, C. Chunqi, and X. Weichao, “A recursive RARE algorithm for DOA estimation with unknown mutual coupling,” IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 1593–1596, 2014. View at: Publisher Site  Google Scholar
 Y. Ma, Y. Yang, Z. He, K. Yang, C. Sun, and Y. Wang, “Theoretical and practical solutions for highorder superdirectivity of circular sensor arrays,” IEEE Transactions on Industrial Electronics, vol. 60, no. 1, pp. 203–209, 2013. View at: Publisher Site  Google Scholar
 R. H. F. Chan and X.Q. Jin, An Introduction to Iterative Toeplitz Solvers, SIAM, Philadelphia, Pa, USA, 2007. View at: Publisher Site  MathSciNet
 S. Liu, L. Yang, and S. Yang, “Robust joint calibration of mutual coupling and channel gain/phase inconsistency for uniform circular array,” IEEE Antennas and Wireless Propagation Letters, 2015. View at: Publisher Site  Google Scholar
 C. M. S. See, “Sensor array calibration in the presence of mutual coupling and unknown sensor gains and phases,” Electronics Letters, vol. 30, no. 5, pp. 373–374, 1994. View at: Publisher Site  Google Scholar
 L. Bin and L. Guisheng, “Method for array gain and phase uncertainties calibration based on ISM and ESPRIT,” Journal of Systems Engineering and Electronics, vol. 20, no. 2, pp. 223–228, 2009. View at: Google Scholar
 P. Stoica, E. G. Larsson, and A. B. Gershman, “The stochastic CRB for array processing: a textbook derivation,” IEEE Signal Processing Letters, vol. 8, no. 5, pp. 148–150, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Song Liu 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.