Active Vibration Control in Mechanical SystemsView this Special Issue
Research Article | Open Access
2D-DOA Estimation for Cylindrical Array with Mutual Coupling
Most conventional direction-of-arrival (DOA) estimation algorithms are affected by the effect of mutual coupling, which make the performance of DOA estimation degrade. In this paper, a novel DOA estimation algorithm for conformal array in the presence of unknown mutual coupling is proposed. The special mutual coupling matrix (MCM) is applied to eliminate the effect of mutual coupling. With suitable array design, the decoupling between polarization parameter and angle information is accomplished. The two-demission DOA (2D-DOA) estimation is finally achieved based on estimation of signal parameters via rotational invariance techniques (ESPRIT). The proposed algorithm can be extended to conical conformal array as well. Two parameter pairing methods are illustrated for cylindrical and conical conformal array, respectively. The computer simulation verifies the effectiveness of the proposed algorithm.
DOA estimation is a major part of array signal processing which is applied in various areas, such as smart antenna, wireless communication system, radar, sonar, and seismology [1–5]. The conventional algorithms such as the multiple signal classification (MUSIC)  and estimation of signal parameters via rotational invariance techniques (ESPRIT)  possess super-resolution ability. Recently, sparse signal representation and compressive sensing (CS) methods such as -SVD  and JLZA-DOA  are applied to estimate DOA [10–16]. However, the DOA estimation performance will degrade severely due to the effect of mutual coupling.
In order to solve this problem, many algorithms have been proposed to deal with it. By adding instrumental elements, the reconstitution MCM can be used to eliminated the mutual coupling effect . This method circumvented the difficulty of dealing with coherent signals in 2D-DOA estimations based on maximum likelihood (ML) method. However, the algorithm suffers from tremendous computational complexity. Conformal array is the array mounted on the curvature carrier . It has many advantages such as reduction of aerodynamic drag, reduction of the weight of array, space saving, and reduction of radar cross-section (RCS) . Several algorithms have been proposed for parameter estimation of conformal array recently. The joint polarization parameter and DOA estimation algorithm were proposed based on ESPRIT . The four-order cumulant and ESPRIT were used for DOA estimation, which suffer from high computational complexity . The subarray divided technique and MUSIC were adopted for DOA estimation for cylindrical conformal array [22, 23]. The parallel factor (PARAFAC) analysis is used in  for 2D-DOA estimation. The space-time matrix and PARAFAC are used in  for joint frequency and 2D-DOA estimations. The state space and propagator method were used for joint frequency and 2D-DOA estimations .
In this paper, a novel DOA estimation algorithm for conformal array in the presence of unknown mutual coupling is proposed. The selection matrix is applied to obtain the special MCM, which can eliminate the effect of mutual coupling. Then the decoupling between polarization parameter and angle information is accomplished with good array design. Two rotation invariance matrices are constructed for 2D-DOA estimation based on ESPRIT. Then the proposed algorithm is extended to conical conformal array. Two parameter pairing methods are illustrated for cylindrical conformal array, respectively. The proposed algorithm possesses good estimation performance according to the computer simulation results. Compared with ESPRIT algorithm, the computational complexity is not increasing, which is suitable for real-time application in cylindrical array system. In the future work, we will focus on the application of the proposed algorithm [27–29].
The organization of this paper is structured as follows. Section 2 introduces the snapshot data model for conformal array with unknown mutual coupling. Section 3 illustrates the mutual coupling elimination mechanism. Section 4 contains the core contributions of this paper, which give the DOA estimation algorithm in the presence of mutual coupling using ESPRIT. Section 5 presents the simulation results. Section 6 summarizes our conclusions.
2. The Snapshot Data Model
We construct the mathematic model of the conformal array using the model which is proposed in . As shown in Figure 1(a), a narrowband far-field signal impinges on the conformal array. The elevation and azimuth are and , respectively. Thus the steering vector can be expressed as where () represents the distance vector between the mth element and the origin point, is the direction vector of the incident signal, and is wavelength of the incident signal. As shown in Figure 1(b), is the pattern of the mth element and and are the components of , which are projected onto two orthogonal directions of electric field of the incident signal, respectively. and are polarization parameters. and denote the mth element’s response to the incident signal and electric field vector, respectively. is the angle between and . denoting the transpose of matrix .
Assume that uncorrelated narrowband plane waves with unknown DOAs impinge on the conformal array; the array output can be represented as where and denote the pattern matrix and the manifold matrix, respectively. “” stands for the Hadamard product. and are the diagonal matrices whose diagonal entries are and , respectively. and are polarization parameters of the th signal. denotes the signal vector. is additive Gaussian white noise (AGWN) with zero mean and the covariance matrix is denotes conjugate transpose of matrix . is the identical matrix.
The matrix is the MCM of the conformal array. Generally, the effect of mutual coupling among elements is ineluctable in practice. According to the reciprocity principle, the MCM is a symmetrical matrix.
Consider a scenario in which snapshots are accumulated; we have where is the signal matrix and is the noise matrix.
For ULA, the mutual coupling between adjacent elements is almost the same, and the magnitude of the coupling coefficients decays very fast along with the increasing element spacing. In essence, mutual coupling coefficients between two far apart elements can be approximated to zero. For simplicity, only the mutual coupling among elements which belongs to the identical ULA is considered in this paper. The mutual coupling degree of freedom is assumed to be , which means that for the ith element, the coupling comes from the th to th elements. Thus, the MCM is modeled as a banded symmetric Toeplitz matrix, whose first row is satisfying . Then, MCM is represented as
Due to the varying curvature of the carrier, the pattern of each element is defined in the local coordinate, which means that the transformation from the global coordinate to the local coordinate should be done . For the presence of the unknown mutual coupling, we must inhibit its effect. Furthermore, the polarization parameters couple with the signal parameters is another problem that we need to deal with.
3. The Elimination of Effect of Unknown Mutual Coupling
As shown in Figure 2, the neighboring elements spacing mounted on the same intersecting surface of the cylinder is . The two adjacent intersecting surfaces spacing are . The radius of the cylinder is .
Based on the snapshot data model constructed above, we will give an elimination method of mutual coupling using the instrumental elements. In order to illustrate the mechanism how the effect of mutual coupling is eliminated, ~ elements are adopted as an example. The coordinate of each element is defined as , . represents the position vector of the th element. , and represent the unit vector of -axis, -axis, and -axis, respectively. can be expressed as
Because the elements are arranged on the same generatrix, we have and , respectively. Also, the elements possess the identical element's response . The manifold matrix can be expressed as , and the steering vector is where and
In order to eliminate the effect of mutual coupling, the first and last elements of the -elements ULA are selected as instrumental elements . At the first glance, it is an unreasonable idea to abandon some information that the array receives. However, the effect of mutual coupling is eliminated by this sacrifice. Define a selection matrix ; the output of the center array is expressed as  where the MCM of the center array is given by
After adding the instrumental elements, an important relationship between MCM and the steering vector is represented as where is the ideal steering vector of the center array ( elements) and is a scalar only concerned with the mutual coupling coefficients and the direction of the incident signal . is defined as Thus, the covariance matrix of the centre array is given by where
Assume that the power of the th incident signal is . Equation (16) can be written as where
4. DOA Estimation Based On ESPRIT
We can deal with the coupling between the polarization parameters and the signal parameters after the mutual coupling is eliminated. Define two
Choosing , , , and elements constitutes array 1, array 2, array 3, and array 4, respectively. The distance vector between array 1 and 2 is and . The distance vector between arrays 1 and 3 is , which is shown in Figure 3. The elements' patterns of the same generatrix are identical. With special array design, both the eliminations of effect of mutual coupling and the decoupling between polarization parameter and angle information can be realized.
(a) The distance vector
(b) The distance vector
Define a novel array output vector ; the covariance matrix is represented as where and are expressed as Regardless of the noise the signal subspace is obtained by taking the eigenvalue decomposition (EVD) of . Then is divided into two submatrices (the top and the bottom), and we have a full rank matrix satisfying We can estimate by The diagonal matrix constructed by the eigenvalue of must be equal to . The columns of matrix construct the eigenvectors of . Then array 3 and array 4 can construct a group array. Thus the output of the centre array 3 and array 4 is given by respectively, where represents the MCM and is the manifold matrix of array. The group array output is defined as .
The covariance matrix can be expressed as where and () represent the elevation and azimuth of the distance vector in the global coordinate, respectively.
The subspaces of array 3 and array 4 are defined as and , respectively. is similarly defined as , which can be solved by There exists an invertible matrix , satisfying where can be obtained by taking the EVD of matrix and the eigenvalue of it is defined as . There may exist phase ambiguous caused by the positive and negative of and . Because and are real numbers, this problem can be solved by squaring the eigenvalue : where is the angle of complex. Based on (25) and (31), the elevation and azimuth can be obtained by
Then the DOA estimation in the presence of mutual coupling is accomplished. However, when , the parameter pairing between and is a big problem. The eigenvectors of the identical source are strongly correlated, and thus we can construct the sequencing matrix . The parameter pairing can be accomplished by adjusting the order of eigenvectors.
and are eigenvectors which correspond to and , respectively. We have the equation as follows: According to the coordinate of the maximal entry in the columns (or rows) of matrices , we adjust the order of eigenvectors. Then the parameter pairing is completed.
5. Simulation Results
For the sake of simplicity, only the cylindrical conformal array is utilized to simulate the experiments as follows: the array structure is shown in Figure 2. The number of elements is 16; that is, . The lowest order circular patch model [21, 30] is used in the simulation, and the element pattern is given by where and stand for the elevation and azimuth of the th incident signal in the th local coordinate, respectively.
The RMSE against SNR is plotted at first, which is shown in Figure 4(a). Two spatially uncorrelated narrowband incident signals impinge on the cylindrical conformal array with the elevation and azimuth of and , respectively. The measurements are corrupted by spatially and temporally uncorrelated additive Gaussian white noise AGWN. Without loss of generality, , ; , . The number of mutual coupling coefficients is , and , . The number of snapshots is 200. 100 Monte Carlo trials are simulated under the parameters above and the following experiments. It is shown in Figure 4(a) that RMSE of NKMC-ESPRIT is extremely larger than that of ESPRIT and the proposed algorithm. The unknown mutual coupling leads to the degradation of the estimation performance. In this scenario, NKMC-ESPRIT could not obtain the actual DOAs of the incident signals. The RMSE of the proposed algorithm is slightly larger than that of ESPRIT. This is caused by the aperture loss of the proposed algorithm; only 12 elements are applied to estimate DOA, which are less than that of ESPRIT (16 elements).
(a) The RMSE against SNR
(b) The successful probability against SNR
Figure 4(b) shows the successful probability against SNR of the three different algorithms. The successful probability of NKMC-ESPRIT approximates to zero, which implies that NKMC-ESPRIT could not estimate DOA accurately. The result is identical with the large RMSE of NKMC-ESPRIT mentioned above. The successful probability of ESPRIT and the proposed algorithm increase with the SNR increases. When SNR is larger than 12 dB, the successful probability of ESPRIT is almost 100%. The proposed algorithm needs SNR to be larger than 14 dB. Without the effect of mutual coupling, ESPRIT outperforms the proposed algorithm. However, the proposed algorithm possesses good estimation performance with unknown mutual coupling.
The RMSE against snapshot number is shown in Figure 5(a). SNR is 0 dB. Other simulation conditions are identical with the first simulation. It is shown in Figure 5(a) that RMSE of NKMC-ESPRIT is stable with the varying snapshot number. It means that NKMC-ESPRIT could not be used for DOA estimation because of the effect of mutual coupling. RMSE of proposed algorithm is slightly larger than that of ESPRIT. The aperture loss caused by the elimination mechanism of mutual coupling leads to this result. However, with snapshot number increases, RMSE of ESPRIT and the proposed algorithm are not varied too much. This result shows that RMSE is not affected by the varying snapshot number.
(a) The RMSE against snapshot number
(b) The successful probability against snapshot number
The successful probability against snapshot number is shown in Figure 5(b). The successful probability of NKMC-ESPRIT approximates to zero with varying snapshot number. When snapshot number is larger than 500, the successful probability of ESPRIT is almost 100%. The proposed algorithm needs 1000 snapshot number to achieve 95% successful probability. The low SNR mainly leads to this large requirement of snapshot number. This simulation gives different result compared with RMSE simulation, which shows that the successful probability is improved with the snapshot number increases.
In this paper, a novel DOA estimation algorithm for conformal array in the presence of unknown mutual coupling is proposed. The special MCM is applied to eliminate the effect of mutual coupling. With good array design, the decoupling between polarization parameter and angle information is accomplished. The 2D-DOA estimation is finally achieved based on ESPRIT algorithm. The proposed algorithm can also be extended to conical conformal array. The computer simulation verifies the effectiveness of the proposed algorithm.
Conflict of Interests
The paper is approved by all authors for publication. The authors declare that there is no conflict of interests regarding the publication of this paper.
One of the authors would like to declare on behalf of other coauthors that the work described is original research that has not been published previously.
- H. Krim and M. Viberg, “Two decades of array signal processing research: the parametric approach,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996.
- D. H. Johnson and D. E. Dudgeon, Array Signal Processing—Concepts and Techniques, Prentice Hall, Englewood Cliffs, NJ, USA, 1993.
- P. Stocia and R. Moses, Introduction to Spectral Analysis, Prentice Hall, Upper Saddle River, NJ, USA, 2nd edition, 2004.
- S. Yin, S. Ding, A. Haghani, H. Hao, and P. Zhang, “A comparison study of basic datadriven fault diagnosis and process monitoring methods on the benchmark Tennessee Eastman process,” Journal of Process Control, vol. 22, no. 9, pp. 1567–1581, 2012.
- S. Yin, H. Luo, and S. Ding, “Real-time implementation of fault-tolerant control systems with performance optimization,” IEEE Transactions on Industrial Electronics, vol. 64, no. 5, pp. 2402–2411, 2014.
- R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.
- R. Roy and T. Kailath, “ESPRIT-Estimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989.
- D. Malioutov, M. Çetin, and A. S. Willsky, “A sparse signal reconstruction perspective for source localization with sensor arrays,” IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 3010–3022, 2005.
- M. M. Hyder and K. Mahata, “Direction-of-arrival estimation using a mixed norm approximation,” IEEE Transactions on Signal Processing, vol. 58, no. 9, pp. 4646–4655, 2010.
- P. Stoica, P. Babu, and J. Li, “SPICE: a sparse covariance-based estimation method for array processing,” IEEE Transactions on Signal Processing, vol. 59, no. 2, pp. 629–638, 2011.
- N. Nu, Z. Ye, and M. Bao, “DOA estimation for sparse array via sparse signal reconstruction,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 2, pp. 760–773, 2013.
- Z. M. Liu, Z. T. Huang, and Y. Y. Zhou, “Array signal processing via sparsity-inducing representation of the array covariance matrix,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 3, pp. 1710–1724, 2013.
- S. Yin, G. Wang, and H. Karimi, “Data-driven design of robust fault detection system for wind turbines,” Mechatronics, 2013.
- A. C. Gurbuz, V. Cevher, and J. H. McClellan, “Bearing estimation via spatial sparsity using compressive sensing,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 2, pp. 1358–1369, 2012.
- E. T. Northardt, I. Bilik, and Y. I. Abramovich, “Spatial compressive sensing for direction-of-arrival estimation with bias mitigation via expected likelihood,” IEEE Transaction on Signal Processing, vol. 61, no. 5, pp. 1183–1195, 2013.
- M. Carlin, P. Rocca, G. Oliveri, F. Viani, and A. Massa, “Direction of arrival estimation through Bayesian compressive sensing strategies,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 7, pp. 3828–3838, 2013.
- Z. Ye and C. Liu, “On the resiliency of MUSIC direction finding against antenna sensor coupling,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 2, pp. 371–380, 2008.
- L. Josefsson and P. Persson, Conformal Array Antenna Theory and Design, IEEE Press, Piscataway, NJ, USA, 2006.
- W. T. Li, X. W. Shi, Y. Q. Hei, S. F. Liu, and J. Zhu, “A hybrid optimization algorithm and its application for conformal array pattern synthesis,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 10, pp. 3401–3406, 2010.
- L. Zou, J. Lasenby, and Z. He, “Direction and polarization estimation using polarized cylindrical conformal arrays,” IET Signal Processing, vol. 6, no. 5, pp. 395–403, 2012.
- Z.-S. Qi, Y. Guo, and B.-H. Wang, “Blind direction-of-arrival estimation algorithm for conformal array antenna with respect to polarisation diversity,” IET Microwaves, Antennas and Propagation, vol. 5, no. 4, pp. 433–442, 2011.
- P. Yang, F. Yang, and Z. Nie, “DOA estimation using MUSIC algorithm on a cylinderical conformal array,” in Proceedings of the IEEE Antennas and Propagation Society International Symposium, pp. 5299–5302, June 2007.
- P. Yang, F. Yang, and Z. P. Nie, “DOA estimation with sub-array divided technique and interporlated esprit algorithm on a cylindrical conformal array antenna,” Progress in Electromagnetics Research, vol. 103, pp. 201–216, 2010.
- L. Wan, W. Si, L. Liu, Z. Tian, and N. Feng, “High accuracy 2D-DOA estimation for conformal array using PARAFAC,” International Journal of Antennas and Propagation, vol. 2014, Article ID 394707, 14 pages, 2014.
- L. Wan, W. Si, L. Liu, and Z. Tian, “Joint estimation of 2D-DOA and frequency based on space-time matrix and conformal array,” The Scientific World Journal, vol. 2013, Article ID 463828, 10 pages, 2013.
- W. J. Si, L. T. Wan, L. T. Liu, and Z. X. Tian, “Fast estimation of frequency and 2-D DOAs for cylindrical conformal array antenna using state-space and propagator method,” Progress in Electromagnic Research, vol. 137, pp. 51–71, 2013.
- L. Liu, Y. Jiang, L. Wan, and Z. X. Tian, “Beamforming of joint polarization-space matched filtering for conformal array,” The Scientific World Journal, vol. 2013, Article ID 589675, 10 pages, 2013.
- S. Yin, S. Ding, A. Haghani, and H. Hao, “Data-driven monitoring for stochastic systems and its application on batch process,” International Journal of Systems Science, vol. 44, no. 7, pp. 1366–1376, 2013.
- S. Yin, X. Yang, and H. Karimi, “Data-driven adaptive observer for fault diagnosis,” Mathematical Problems in Engineering, vol. 2012, Article ID 832836, 21 pages, 2012.
- B.-H. Wang, Y. Guo, Y.-L. Wang, and Y.-Z. Lin, “Frequency-invariant pattern synthesis of conformal array antenna with low cross-polarisation,” IET Microwaves, Antennas and Propagation, vol. 2, no. 5, pp. 442–450, 2008.
Copyright © 2014 Hao Feng 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.