# An Efficient 2-D DOA Estimation for a Cylindrical Conformal Array with Unknown Mutual Coupling.

1. IntroductionDirection of arrival (DOA) estimation is a major research direction in array signal processing area. It has broad application prospects in many occasions, and lots of algorithms have been developed in the past decades. However, most of them are only applicable to planar arrays. Conformal array has a lot of advantages such as good aerodynamic performance, wide surveillance area, and low radar cross-section. Therefore, it has drawn extensive attention in recent years [1]. Compared with the planar array, the curved structure of a conformal array brings new challenges to the problem of DOA estimation. Moreover, due to the limited arrange space of a conformal array, the effect of mutual coupling may become more obvious and consequently degrade the performance of most existing DOA estimation algorithms [2].

Researchers have been focused on this problem, and a large number of algorithms have been proposed to calibrate the imperfect manifold caused by mutual coupling. In [3], mutual coupling effect was calibrated as well as sensor location error and phase/gain uncertainty by an iterative process. In [4], a unified framework was formulated to express the typical array perturbations and a sparse Bayesian method was proposed for array calibration and DOA estimation. However, the convergence of these iterative methods was not theoretical guaranteed and the associated multidimensional search and iterations were time consuming. By taking advantages of the special structure of uniform linear array (ULA) and uniform rectangular array (URA), auxiliary elements were demonstrated to be effective for the autocalibration of mutual coupling [5-7]. This kind of method was further extended to the situation of non-Gaussian noise [8] and coherent signal condition [9] by using fourth-order cumulants (FOC) and spatial smoothing technique, respectively. These auxiliary algorithms can accurately estimate DOAs without any calibration source or iteration procedure. However, there is an aperture loss for these methods and they can only be applied to uniform arrays.

The rank reduction (RARE) algorithm is another class of methods for DOA estimation with unknown mutual coupling [10-13]. Through matrix transformation, the coupling effect is eliminated and a rank reduction matrix is defined to estimate the DOAs. In [12], an iterative process was introduced to further improve the DOA estimation accuracy. In [13], a recursive RARE method for DOA estimation in the presence of unknown mutual coupling was presented. It recursively built up RARE matrices by the angle estimation to overcome the problem of nonunique coupling coefficient estimation and gave a better DOA estimation performance than the original RARE method. Generally, these methods use all the array elements rather than the middle subarray to avoid the array aperture loss and have a better performance than the auxiliary methods, especially for small arrays. Other methods are also proposed to calibrate the coupling effect. In [14], mutual coupling matrix (MCM) was reconstructed to calibrate the mutual coupling effect, in which pattern matrix for the embedded elements and for the isolated elements needed to be known in prior. By using the special structure of the uniform linear array, a Toeplitz matrix was reconstructed with the receiving data to replace the covariance matrix [15]. Based on that, subspace-based methods can be used without any compensation.

However, most of the above methods can only be applied to the linear or planar array. Mutual coupling calibration for a conformal array is more complicated, and few researches have been involved in this problem so far. In [16], the RARE algorithm was adopted to estimate DOAs for the conformal array. In [17], a tensor 2-D DOA estimation was proposed for cylindrical conformal array (CCA) under unknown mutual coupling. However, only the mutual coupling between elements in one ULA is considered, which simplifies the problem. Moreover, both of the above methods involved spectrum searching and lots of matrix computation and the amount of calculation will dramatically increase in a 2-D case.

In this paper, an efficient decoupling method is proposed to estimate the 2-D DOAs and the mutual coupling coefficients for CCA. The elevation is first estimated based on a subarray ESPRIT method. By using the elevation estimate, azimuth can be further estimated based on a 1-D RARE method. After that, mutual coupling coefficients can subsequently be solved. The 2-D DOA estimation is decoupled and divided into two individual 1-D DOA estimations, which can significantly reduce the computational complexity. Numerical experiments are provided to demonstrate the effectiveness of the proposed method.

2. Signal Model

Consider a uniform CCA with the size of M x N, in which each layer is an N element uniform circular array (UCA) and each column is an M element ULA. The array geometry, coordinate system, and element numbering are shown in Figure 1. The space between neighbouring sensors in UCA is [d.sub.xy] and that of ULA is dz. Assume that there are K signals, [s.sub.i](t), i =1,2, ..., K, arriving at the array from directions ([[theta].sub.1], [[phi].sub.1]), ([[theta].sub.2], [[phi].sub.2]), ..., ([[theta].sub.k], [[phi].sub.K]), respectively. The symbols [[theta].sub.i] and [[phi].sub.i] denote the azimuth and elevation of the ith signal. Assuming that the MCM of CCA is C, then the array output vector x(t) can be written as

x(t) = CAs(t) + n(t). (1)

The MN x 1 output signal vector x(t), the K x1 source signal vector s(t), and the MN x 1 noise vector n(t) in (1) are defined as follows.

[mathematical expression not reproducible]. (2)

The array manifold matrix A is defined as

[mathematical expression not reproducible], (3)

where [k.sub.0] - 2[pi]/[lambda] is the wave number, [r.sub.j] - [[[x.sub.j], [y.sub.j], [z.sub.j]].sup.T] is the coordinate vector of the jth element, [v.sub.i] - [[cos [[theta].sub.i] cos [[phi].sub.i], sin [[theta].sub.i] cos [[phi].sub.i], sin [[phi].sub.i]].sup.T] is the unit vector of direction ([[theta].sub.i], [[phi].sub.i], and [g.sub.j]([theta], [phi]) is the element pattern of the jth element defined in the global coordinate system.

It is found that mutual coupling effect decreases along with the increase of space between elements [18]. Without loss of generality, assume that each sensor only suffers from the mutual coupling effect of its eight neighbouring elements, and the mutual coupling effect model is shown in Figure 2. The coefficients are assumed to satisfy the relationship [absolute value of [c.sub.xz] < [absolute value of cx], [absolute value of [c.sub.xz] < [absolute value of cz]. The whole CCA can be seen as a compound ULA composed of several identical UCAs. Then the MCM C of the whole array can be expressed as follows.

[mathematical expression not reproducible], (4)

where [C.sub.1] is the MCM of each UCA subarray and [C.sub.2] is the MCM between two neighbouring UCAs. Due to the uniformly spaced structure, they are symmetric Toeplitz matrices [18], and can be expressed as

[mathematical expression not reproducible]. (5)

3. DOA Estimation for CCA with Unknown Mutual Coupling

In this section, dimensionality reduction is first carried out to decrease the computational complexity. The 2-D DOA estimation for CCA is divided into two 1-D DOA processes, and then azimuth and elevation of signals are decoupled and estimated, respectively.

3.1. Elevation Estimation. According to the mutual coupling model, the top and bottom layers are set to be auxiliary elements and two subarrays are selected to estimate the elevation of signals. They compose of the same number of layers, which is shown in Figure 1. Hence, the outputs of the two subarrays can be written as

[mathematical expression not reproducible], (6)

in which [mathematical expression not reproducible]. In (6), [P.sub.1] and [P.sub.2] are two screening matrices defined as

[mathematical expression not reproducible], (7)

where 0 is an M x M zero matrix and I is an M x M identity matrix.

In the presence of mutual coupling, the real steering vectors of the two subarrays are [P.sub.1] Ca([theta], [phi]) and [P.sub.2]Ca([theta], [phi]). The ideal steering vector of the whole CCA a([theta], [phi]) can be rewritten as

a([[theta].sub.i], [[phi].sub.i]) - [a.sub.z]([[theta].sub.i], [[phi].sub.i]) [cross product] [a.sub.xy] ([[theta].sub.i], [[phi].sub.i]), i - 1, ..., K (8)

The symbol [cross product] denotes the Kronecker product and [a.sub.z]([[theta].sub.i], [[phi].sub.i]) and [a.sub.xy]([[theta].sub.i], [[phi].sub.i]) denote the steering vectors of the linear and circular subarrays without mutual coupling, respectively.

[mathematical expression not reproducible]. (9)

Taking a and [a.sub.xy] as the shorter terms for a([theta], [phi]) and [a.sub.xy]([theta], [phi]), respectively, then the steering vectors of the two subarrays can be expressed as follows.

[mathematical expression not reproducible] (10)

[mathematical expression not reproducible] (11)

By comparing (10) and (11), it can be found that the steering vectors of the two subarrays have similar structures and it can be simplified as the following equation.

[P.sub.2]Ca = [P.sub.1]Ca[[beta].sub.z]. (12)

Substitute (12) into (6) and the following equations can be obtained.

[mathematical expression not reproducible], (13)

where [PHI] = diag {[[beta].sub.z1], ..., [[beta].sub.zK]}. The rotation matrix [PHI] only contains the information of signal elevation, which is not affected by mutual coupling. Therefore, the elevation can be estimated by using ESPRIT algorithm. Define a new signal vector z(t) = [[x.sub.1][(t).sup.T], [[x.sub.2] [(t).sup.T]].sup.T]. The signal subspaces [U.sub.s] can be obtained via the eigen decomposition of the covariance matrix of z(t), and it can subsequently be divided into the up and down halves [U.sub.S1] and [U.sub.S2]. It is easy to derive the following equations.

[U.sub.S1] = [U.sub.S1][psi] [psi] = [T.sup.-1] [PHI]T. (14)

Since T is a full-rank matrix, [PHI] and [psi] have the same eigenvalues. Therefore, the elevations can be estimated by the eigenvalues [[phi].sub.i], i = 1, ..., K of [psi].

[[phi].sub.i] = [sin.sup.-1] arg ([lambda]/2[pi][d.sub.z] arg ([[phi].sub.i])), i =1, ..., K. (15)

3.2. Azimuth and Mutual Coupling Coefficient Estimation. Based on the above elevation estimation, the azimuth can be estimated by utilizing the output of the whole array. The covariance matrix of x(t) is defined as

R = E{x(t)[x.sup.H](t)} = [CAR.sub.S][A.sup.H] [C.sup.H] + [[sigma].sup.2][I.sub.M], (16)

where [R.sub.S] = E{s(t)[s.sup.H](t)} is the signal covariance matrix. The signal subspace [U.sub.s] and the noise subspace [U.sub.N] can be obtained by taking eigen decomposition of the covariance matrix R. The space spans by the steering vectors coincide with the signal subspace and are orthogonal with the noise subspace.

span{CA} = span{[U.sub.s]} [+ or -] [perpendicular to] span{[U.sub.N]}. (17)

Therefore, we have

[U.sup.H.sub.N]Ca([[theta].sub.i], [[phi].sub.i]) = 0, i = 1,2, ... K. (18)

As the special structure of the MCM C in (4), the steering vector of the whole array can be rewritten as

[mathematical expression not reproducible]. (19)

Define the steering vector of the nth UCA as [a.sub.n] = [[beta].sup.n-1.sub.z] [a.sub.xy], n = 1, ..., N. Since the submatrices [C.sub.1] and [C.sub.2] are Toeplitz matrices, the following equation can be derived [10].

[C.sub.i] [a.sub.n] = [T.sub.n] [c.sub.i], i = 1,2, n = 1, ..., N, (20)

in which [mathematical expression not reproducible] is a matrix of size M x 2 and

[mathematical expression not reproducible], (21)

where p = M/2 when M is even and p = (M + 1/2 when M is the odd. Then (19) can be rewritten as

[mathematical expression not reproducible], (22)

in which c = [[1, [c.sub.x], [c.sub.z], [c.sub.xz]].sup.T] is the vector of mutual coupling coefficients. Substituting (22) into (18), we have

[U.sup.H.sub.N] Tc = 0. (23)

Define the matrix Q([theta], [phi]) = [T.sup.H] [U.sub.N] [U.sup.H.sub.N] T, and the estimation of the DOAs can be obtained by the K roots of the following equation.

[c.sup.H] Q([theta], [phi]) c = 0. (24)

Because matrix Q([theta], [phi]) is a Hermite matrix and c [not equal to] 0, (24) holds if and only if matrix Q([theta], [phi]) is singular. Generally, matrix Q([theta], [phi]) is full rank if MN - K [greater than or equal to] 4, while it will be singular when ([theta], [phi]) coincides with one of the true signal directions ([[theta].sub.1], [[phi].sub.1]), i = 1, ..., K. Therefore, the following azimuth estimators can be derived.

[mathematical expression not reproducible], (25)

where det [Q] is the determinant of matrix Q, [[lambda].sub.min] [Q] is the smallest eigenvalue of Q, and [v.sub.min] {Q} is the corresponding eigenvector. It is noteworthy that each azimuth only needs a 1-D spectrum searching with the fixed elevation to get its estimation, which significantly reduces the calculation amount.

Based on the estimation of DOAs, the mutual coupling coefficient vector c can be obtained as follows.

[mathematical expression not reproducible]. (26)

The above derivation process takes advantage of the special structure of the array MCM, and the method developed here cannot be applied to arbitrary conformal arrays. In spite of that, it is still applicable to the arrays that have similar column structures, which are commonly used in practice.

4. Simulation Results

In this section, several simulations are carried out to illustrate the performance of the proposed method. An M x N uniformly spaced cylindrical array is considered in the following experiments. The observation range of azimuth and elevation of signals are both confined in [0[degrees], 90[degrees]], and the spectrum scanning step is 0.1[degrees]. The space between two neighbouring elements in a UCA is [d.sub.xy] = 0.4[lambda], and the distance between two adjacent UCAs is [d.sub.z] =0.5[lambda]. The two narrowband uncorrelated signals impinging on the array from far-field directions are ([[theta].sub.1], [[phi].sub.1]) - (20[degrees], 30[degrees]) and ([[theta].sub.2], [[phi].sub.2]) - (40[degrees], 50[degrees]). The mutual coupling coefficients are assumed to be [c.sub.1] [[i,0.5 + 0.2i].sup.T] and [c.sub.2] - [[0.4 + 0.1i,0.22 - 0.14i].sup.T]. Each element is assumed to have the following pattern in its local coordinate system.

[mathematical expression not reproducible], (27)

where [theta]' and [phi]' are the azimuth and elevation in the local coordinate system for each element. Element pattern [g.sub.i]([theta], [phi]) in the above analysis is a rotation of g([theta]', [phi]').

For the DOA estimation method, the most timeconsuming process is the spectrum searching, especially in 2-D condition. Under the prescribed condition, the proposed method needs to calculate 901 x 2= 1802 points of the spatial spectrum, while 2-D RARE needs to calculate 901 x 901 = 811,801 points of the spatial spectrum. Although the amount of calculation will increase along with the number of signals, the proposed method still greatly reduces computational complexity.

In the first simulation, the computational complexity is analyzed by comparing the CPU time of the proposed algorithm and the RARE method. Comparisons are carried out for different sizes of arrays. Assume SNR= 10 dB and the number of snapshots is 500. With the Intel (R)Core (TM) i5 CPU 3.2GHz and 4GB RAM, the elapsed time of getting the DOA estimation results for both the proposed method and the RARE method is shown in Table 1. It is shown that the proposed method has a much smaller computing time than RARE and does not increase with the size of the array obviously.

In the second simulation, statistical experiments are carried out to analyze the estimation accuracy of the proposed method. 2-D RARE and MUSIC are also implemented as comparison methods. A uniform CCA with the size of 16 x 10 is used here. 100 Monte Carlo trails are performed for each simulation and the DOA estimation accuracy is measured by the RMSE defined as follows.

[RMSE.sub.[theta]] = [square root of 1/100K [100.summation over (i=1)] [K.summation over (k=1)] [parallel] [[??].sub.k] - [[theta].sub.k] [parallel][sup.2], (28)

where [mathematical expression not reproducible] is the estimation of [[theta].sub.k] - ([[theta].sub.k], [[phi].sub.k]) in the ith Monte Carlo experiment. Figure 3 shows the RMSE of the DOA estimation versus SNR with fixed 500 snapshots. Figure 4 shows the RMSE of the DOA estimation versus the number of snapshots with fixed SNR = 10 dB. These figures clearly indicate that the proposed algorithm can well compensate the mutual coupling effect and achieve satisfactory performance for the DOA estimation. The RMSE decreases with the increase of SNR and the number of snapshots as expected. Since there is a little aperture loss in the process of elevation estimation, the estimation accuracy of the proposed method is slightly lower than that of 2-D RARE. However, it is tolerable for this little performance loss with the distinct improvement in computational complexity.

Based on the estimation of DOAs, the mutual coupling coefficients can be obtained. To demonstrate its estimation accuracy, the RMSE of the mutual coupling coefficients is similarly defined as follows.

[RMSE.sub.[theta]] = [square root of 1/100L [100.summation over (i=1)] [parallel] [[??].sub.i] - c [parallel][sup.2], (29)

where c, is the estimation of c in the ith Monte Carlo experiment and L is the number of the mutual coupling coefficient.

The third simulation considers the same scenario as the second one. By using the above DOA estimation results, mutual coupling coefficients are subsequently estimated. Figure 5 shows the RMSE of the estimated coefficients versus the SNR with fixed 500 snapshots. Figure 6 shows the RMSE of the estimated coefficients versus the number of snapshots with fixed SNR=10dB. Simulation results show that the estimation errors of mutual coupling coefficients for the two methods are very close. It can also be seen that higher estimation accuracy can be achieved with the larger SNR or the number of snapshots.

5. Conclusion

In this paper, an efficient 2-D DOA estimation method is proposed for the cylindrical conformal array in the presence of unknown mutual coupling. Through dimension reduction, the 2-D DOA estimation problem is divided into two 1-D DOA estimations for the cylindrical conformal array. The elevation and azimuth can be decoupled and estimated successively, and the mutual coupling coefficients can be obtained based on the DOA estimates. The proposed method avoids the complicated calculation of the 2-D spatial spectrum and can significantly reduce the amount of computation. Simulation results demonstrate the superior performance of the proposed method in terms of accuracy and computational efficiency.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest. Acknowledgments

This work is supported by the Natural Science Foundation of Anhui Province of China (no. 1608085MF134).

References

[1] X. Lan, Y. Li, and Y. Chen, "A novel 2D DOA estimation via tensor modeling for cylindrical conformal array," in 2017 9th International Conference on Wireless Communications and Signal Processing (WCSP), pp. 1-6, Nanjing, China, October 2017.

[2] Y. Yu, H. Lui, C. Niow, and H. T. Hui, "Improved DOA estimations using the receiving mutual impedances for mutual coupling compensation: an experimental study," IEEE Transactions on Wireless Communications, vol. 10, no. 7, pp. 2228-2233, 2011.

[3] A. Weiss and B. Friedlander, "Direction finding in the presence of mutual coupling," IEEE Transactions on Antennas and Propagation, vol. 39, no. 3, pp. 273-284, 1991.

[4] Z. Liu and Y. Zhou, "A unified framework and sparse Bayesian perspective for direction-of-arrival estimation in the presence of array imperfections," IEEE Transactions on Signal Processing, vol. 61, no. 15, pp. 3786-3798, 2013.

[5] 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.

[6] Z. Ye and C. Liu, "2-D DOA estimation in the presence of mutual coupling," IEEE Transactions on Antennas and Propagation, vol. 56, no. 10, pp. 3150-3158, 2008.

[7] Z. Ye, J. Dai, X. Xu, and X. Wu, "DOA estimation for uniform linear array with mutual coupling," IEEE Transactions on Aerospace and Electronic Systems, vol. 45, no. 1, pp. 280-288, 2009.

[8] Z. Ye, C. Liu, and Y. Zhang, "DOA estimation based on fourth-order cumulants with unknown mutual coupling," Signal Processing, vol. 89, no. 9, pp. 1839-1843, 2009.

[9] J. Dai and Z. Ye, "Spatial smoothing for DOA estimation of coherent signals in the presence of unknown mutual coupling," Signal Processing, vol. 5, no. 4, p. 418, 2011.

[10] C. Qi, Y. Wang, Y. Zhang, and H. Chen, "DOA estimation and self-calibration algorithm for uniform circular array," Electronics Letters, vol. 41, no. 20, pp. 1092-1094, 2005.

[11] B. Liao and S. Chan, "A cumulant-based method for direction finding in uniform linear arrays with mutual coupling," IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 1717-1720, 2014.

[12] B. Liao, Z. Zhang, and S. Chan, "DOA estimation and tracking of ULAs with mutual coupling," IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 1, pp. 891-905,2012.

[13] J. Dai, X. Bao, N. Hu, C. Chang, and W. Xu, "A recursive RARE algorithm for DOA estimation with unknown mutual coupling," IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 1593-1596, 2014.

[14] Q. Huang, H. Zhou, J. Bao, and X. Shi, "Mutual coupling calibration for microstrip antenna arrays via element pattern reconstruction method," IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 51-54, 2014.

[15] W. Mao, G. Li, X. Xie, and Q. Yu, "DOA estimation of coherent signals based on direct data domain under unknown mutual coupling," IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 1525-1528, 2014.

[16] K. Yang, Z. Zhao, W. Yang, and Z. Nie, "Direction of arrival estimation on cylindrical conformal array using RARE," Journal of Systems Engineering and Electronics, vol. 22, no. 5, pp. 767-772, 2011.

[17] X. Lan, L. Wang, Y. Wang, C. Choi, and D. Choi, "Tensor 2-D DOA estimation for a cylindrical conformal antenna array in a massive MIMO system under unknown mutual coupling," IEEE Access, vol. 6, pp. 7864-7871, 2018.

[18] T. Svantesson, "Modeling and estimation of mutual coupling in a uniform linear array of dipoles," in 1999 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings. ICASSP99 (Cat. No. 99CH36258), vol. 5, pp. 2961-2964, Phoenix, AZ, USA, March 1999.

https://doi.org/10.1155/2018/4015980

Chao Liu and Shunian Yin

Department of Communication Engineering, Hefei University of Technology, Hefei, Anhui 230601, China

Correspondence should be addressed to Chao Liu; disneyl@hfut.edu.cn

Received 2 February 2018; Revised 27 June 2018; Accepted 5 July 2018; Published 4 October 2018

Academic Editor: Renato Cicchetti

Caption: Figure 1: Geometry and element numbering of CCA.

Caption: Figure 2: Mutual coupling effect model.

Caption: Figure 3: RMSE of DOA estimation versus SNR with 500 snapshots.

Caption: Figure 4: RMSE of DOA estimation versus the number of snapshots with SNR= 10 dB.

Caption: Figure 5: RMSE of the coupling coefficients versus SNR with 500 snapshots.

Caption: Figure 6: RMSE of the coupling coefficients versus the number of snapshots with 10 dB.

Table 1: The CPU time of the proposed algorithm and the RARE method (s). Array size (M x N) RARE Proposed algorithm 16 x 6 1633.3164 0.7974 16 x 8 1861.5044 1.0987 16 x 10 2267.4448 1.4970

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Liu, Chao; Yin, Shunian |

Publication: | International Journal of Antennas and Propagation |

Date: | Jan 1, 2018 |

Words: | 3992 |

Previous Article: | Codesigned Wideband High-Efficiency Filtering SIW Slot Antenna with High Selectivity and Flat Gain Response. |

Next Article: | Design and Characterization of a Flexible Wideband Antenna Using Polydimethylsiloxane Composite Substrate. |

Topics: |