Next Article in Journal
3D Exploration and Navigation with Optimal-RRT Planners for Ground Robots in Indoor Incidents
Previous Article in Journal
Prediction of 30-Day Readmission for COPD Patients Using Accelerometer-Based Activity Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A High-Resolution and Low-Complexity DOA Estimation Method with Unfolded Coprime Linear Arrays

1
School of Information Engineering, Guangdong University of Technology, Guangzhou 510006, China
2
Institute d’Electronique et Télécommunications de Rennes (IETR), Université de Nantes, 44306 Nantes, France
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 218; https://doi.org/10.3390/s20010218
Submission received: 8 November 2019 / Revised: 4 December 2019 / Accepted: 28 December 2019 / Published: 30 December 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
The direction-of-arrivals (DOA) estimation with an unfolded coprime linear array (UCLA) has been investigated because of its large aperture and full degrees of freedom (DOFs). The existing method suffers from low resolution and high computational complexity due to the loss of the uniform property and the step of exhaustive peak searching. In this paper, an improved DOA estimation method for a UCLA is proposed. To exploit the uniform property of the subarrays, the diagonal elements of the two self-covariance matrices are averaged to enhance the accuracy of the estimated covariance matrices and therefore the estimation performance. Besides, instead of the exhaustive peak searching, the polynomial roots finding method is used to reduce the complexity. Compared with the existing method, the proposed method can achieve higher resolution and better estimation performance with lower computational complexity.

1. Introduction

Direction-of-arrival (DOA) estimation is one of the most active research topics in the field of array signal processing, and it has been widely used in radar, sonar, radio astronomy and other fields [1,2,3,4]. Many DOA estimation methods, such as Multiple Signal Classification (MUSIC) [5] and Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) [6], have been well developed for uniform linear arrays (ULAs), in which the inter-element spacing is restricted to the half-wavelength of incoming signals, leading to a possible mutual coupling effect and inferior estimation performance. To solve these problems, coprime linear arrays (CLAs) have been proposed in [7]. Composed of two ULAs with larger inter-element spacing, a CLA can achieve a larger array aperture with less mutual coupling effect, and consequently better effective estimation performance can be obtained.
The research of the DOA estimation with a CLA can be classified into two categories, which are difference-coarray-based methods and subarray-based methods. The difference-coarray-based methods aim to increase the degrees of freedom (DOFs). However, this class of methods requires a great number of snapshots, which makes it computationally complex [8,9,10,11]. In the subarray-based methods, a CLA is treated as two sparse uniform subarrays. Separately dealing with them, the uniform property of the two subarrays can be directly exploited, making low-complexity DOA estimation possible. Besides, the true DOAs can be determined based on the coinciding results of the two subarrays, and the ambiguities caused by the large inter-element spacing can be avoided thanks to the coprime property [12,13,14]. As they are more suitable for practical applications, many subarrays-based methods have been proposed in recent years. In [12], the MUSIC algorithm is performed separately in both subarrays, and the true DOAs are determined by finding the coinciding peaks of the two spectrums. In [13], the above-mentioned method is improved; by taking advantage of the properties of complex exponential functions, the method limits the peak-searching region and reduces the computational complexity. Another method is proposed in [14]. Benefitting from the uniform property of the two subarrays, ESPRIT is employed. Without the step of peak-searching, the complexity is significantly reduced. Besides, the potential matching error problem in [12,13] is fixed with beamforming-based methods.
Unfortunately, dealing with the two subarrays separately, the subarray-based methods have the following problems: (i) the DOF is determined by the subarray with fewer sensor elements, and therefore the number of detectable signals is limited; (ii) the mutual information of the two subarrays is lost, resulting in an inferior estimation performance; and (iii) an additional step is needed to pair the results obtained from the two subarrays, which increases the computational complexity.
In order to solve these problems, a method based on a new geometry of unfolded coprime linear arrays (UCLAs) is proposed in [15]. By rotating a subarray of a CLA 180°, a non-uniform linear array with a larger aperture can be obtained. Instead of treating the two subarrays separately, MUSIC is performed on the outputs of the whole array. Being superior to a CLA and subarray-based methods, the UCLA can achieve full DOFs as well as better estimation performance due to the exploitation of both self and mutual information. Meanwhile, thanks to the coprime property, only the peaks associated with the true DOAs exist in the MUSIC spectrum. However, this method still has some problems: (i) the uniform property of the two subarrays is wasted; and (ii) the step of peak-searching is involved, increasing the computational complexity.
In this paper, an improved DOA estimation method is proposed. Taking advantage of the uniform property of the two subarrays, we average the diagonal elements of the estimated self-covariance matrices of the observation of the two subarrays to make the covariance matrix of the outputs of the whole array partially Toeplitz. Besides, instead of traditional MUSIC, root-MUSIC is used to avoid the step of peak-searching. Compared with the method in [15], the proposed method can achieve higher resolution and better estimation performance with lower computational complexity. Simulation results are provided to show the performance of the proposed method.

2. System Model

Consider a UCLA composed of two uniform linear subarrays, in which the numbers of sensors are M 1 and M 2 , and the inter-element spacings are d 1 = M 2 λ / 2 and d 2 = M 1 λ / 2 , respectively, with M 1 and M 2 being two coprime integers and λ the wavelength of incoming signals. One sensor element is shared by the two subarrays and is set as the reference point. The total number of sensors M = M 1 + M 2 1 . Figure 1 shows the case of M 1 = 5 and M 2 = 7 .
Assume that K ( K is supposed to be known or can be correctly estimated by the Akaike Information Criterion (AIC) or Minimum Description Length (MDL) method [16], and K < M ) uncorrelated, far-field and narrowband signals impinge on the UCLA from directions { θ 1 ,   θ 2 , , θ K } , with θ k ( π / 2 , π / 2 ) and k [ 1 ,   K ] . The signals received at the two subarrays can be respectively written as
x 1 ( t ) = A 1 s ( t ) + n 1 ( t )
x 2 ( t ) = A 2 s ( t ) + n 2 ( t )
where A i = [ a i ( θ 1 ) a i ( θ 2 ) a i ( θ K ) ] , ( i = 1 ,   2 ) is the directional matrix with
a 1 ( θ k ) = [ 1 e j M 2 π sin ( θ k ) e j ( M 1 1 ) M 2 π sin ( θ k ) ] T
a 2 ( θ k ) = [ e j ( M 2 1 ) M 1 π sin ( θ k ) e j ( M 2 2 ) M 1 π sin ( θ k ) 1 ] T
which are the directional vectors of the two subarrays, respectively. s ( t ) K × 1 denotes the source vector, and n i ( t ) M i × 1 the white Gaussian noise vector with zero-mean and covariance matrix σ 2 I M i , which is independent from the source signals. Therefore, the total outputs of the UCLA can be written as
x ( t ) = [ x 1 ( t ) x 2 ( t ) ] = [ A 1 A 2 ] s ( t ) + [ n 1 ( t ) n 2 ( t ) ]

3. Proposed Method

3.1. Self-Covariance Matrix Reconstruction

The covariance matrix of the total output of the UCLA R can be estimated with L snapshots as
R ^ = 1 L t = 1 L x ( t ) x H ( t ) = [ R ^ 11 R ^ 12 R ^ 21 R ^ 22 ]
with R ^ u v = 1 L t = 1 L x u ( t ) x v H ( t ) and u , v { 1 , 2 } .
The self-covariance matrix of the observation of either subarray can be written as
R u u = [ r u ( 0 ) r u ( 1 ) r u ( 1 M u ) r u ( 1 ) r u ( 0 ) r u ( 2 M u ) r u ( M u 1 ) r u ( M u 2 ) r u ( 0 ) ]
where r u ( m ) = k = 1 K P s k e j m M v sin ( θ k ) , u v , and r u ( m ) = r u * ( m ) . P s k denotes the power of the signal coming from direction θ k .
It can be observed that, because both the subarrays are uniform, the self-covariance matrix of the signal received at either subarray is Toeplitz. All the diagonal elements are equal, and this can be uniquely determined by its first column of the self-covariance matrix. To exploit this property, we can average the diagonal elements of the estimated self-covariance matrix R ^ u u by
r ^ u ( m ) = 1 M u m n = 1 M u m R ^ u u ( m + n , n )
According to Equation (7), depending on the two sets of averaged elements { r ^ 1 ( 0 ) , r ^ 1 ( 1 ) , , r ^ 1 ( M 1 1 ) } and { r ^ 2 ( 0 ) , r ^ 2 ( 1 ) , , r ^ 2 ( M 2 1 ) } , two Toeplitz self-covariance matrices R ^ 11 , T and R ^ 22 , T can be constructed. Replacing R ^ 11 and R ^ 22 by R ^ 11 , T and R ^ 22 , T in R ^ , respectively, we can make the covariance matrix of the total outputs of the UCLA partially Toeplitz, which can improve the estimate performance. In contrast with the traditional covariance matrix averaging technique, which is performed on the whole covariance matrix of the received signals, the proposed partial Toeplitz averaging method is particularly well adapted to unfolded coprime linear arrays. The reconstructed partial Toeplitz covariance matrix R ^ T is given as
R ^ T = [ R ^ 11 , T R ^ 12 R ^ 21 R ^ 22 , T ]

3.2. DOA Estimation

After the reconstruction of the covariance matrix, the eigenvalue decomposition of the obtained partial Toeplitz matrix R ^ T can be expressed as
R ^ T = U ^ s Λ ^ s U ^ s H + U ^ n Λ ^ n U ^ n H
where U ^ s contains the eigenvectors spanning the signal subspace and Λ ^ s a diagonal matrix composed of the K largest eigenvalues of R ^ T , the eigenvalue matrix corresponding to U ^ s ; U ^ n contains the eigenvectors spanning the noise subspace and Λ ^ n the eigenvalue matrix corresponding to U ^ n . According to the orthogonality between the signal and noise subspaces, the spectrum function can be written as
P ( θ ) = 1 [ a 1 H ( θ ) a 2 H ( θ ) ] U ^ n U ^ n H [ a 1 T ( θ ) a 2 T ( θ ) ] T
and the DOAs can be found by searching the peaks of the spectrum P ( θ ) .
To reduce the computational complexity, the polynomial root finding method can be used instead of an exhaustive search. Define
p ( z ) = [ p 1 T ( z ) p 2 T ( z ) ] T
with
p 1 ( z ) = [ 1 z M 2 z ( M 1 1 ) M 2 ] T
p 2 ( z ) = [ z ( M 2 1 ) M 1 z ( M 2 1 ) M 1 1 ] T
which are related to the directional vectors of the two subarrays by
a 1 ( θ ) = p 1 ( z = e j π sin ( θ ) )
a 2 ( θ ) = p 2 ( z = e j π sin ( θ ) )
Then, the exhaustive peak search in Equation (11) can be transformed to the root finding of the following polynomial:
p T ( z 1 ) U ^ n U ^ n H p ( z ) = 0
It can be seen that, since the sensor elements are sparsely and non-uniformly located, p ( z ) contains only several discrete powers of z and Equation (17) is not a full polynomial. In order to solve this problem, we define two transformation matrices as
G 1 = [ 0 1 H 1 ] M 1 × ( 2 M 1 M 2 M 1 M 2 + 1 )
G 2 = [ H 2 0 2 ] M 2 × ( 2 M 1 M 2 M 1 M 2 + 1 )
where H 1 is a selection matrix with a dimension of M 1 × [ ( M 1 1 ) M 2 + 1 ] , of which the [ ( i 1 ) M 2 + 1 ] th element of the i th row is one and the other elements are zeros; H 2 is a selection matrix with the dimension of M 2 × [ ( M 2 1 ) M 1 + 1 ] , of which the [ ( i 1 ) M 1 + 1 ] th element of the i th row is one and the other elements are zeros. 0 1 and 0 2 are two zero matrices with the dimension of M 1 × ( M 2 1 ) M 1 and M 2 × ( M 1 1 ) M 2 , respectively.
Then, we get
p ( z ) = [ G 1 G 2 ] p u ( z )
where
p u ( z ) = [ z ( M 2 1 ) M 1 z ( M 2 1 ) M 1 + 1 1 z ( M 1 1 ) M 2 1 z ( M 1 1 ) M 2 ] T
which contains all continuous power of z . Therefore, Equation (17) can be transformed as
p u T ( z 1 ) U n n U n n H p u ( z ) = 0
where
U n n = [ G 1 T G 2 T ] U ^ n
Equation (22) is a full polynomial to which the root finding technique can be directly applied. It is known that if θ i corresponds to the direction of an actual source, z i = e j π sin ( θ i ) would be a root of Equation (22), and | z i | = | e j π sin ( θ i ) | = 1 . However, because of the presence of noise, the roots may not be precisely located on the unit circle. Besides, note that if z i is a root of Equation (22), so is 1 / z i * . Therefore, half of the roots will be inside the unit circle and half will be outside. The DOAs can be decided by the K roots inside and closest to the unit circle as
θ ^ k = arcsin ( arg ( z ^ k ) π )

4. Simulation and Analysis

4.1. Estimation Performance

In the simulations, the UCLA shown in Figure 1 with M 1 = 5 and M 2 = 7 is considered, and the root mean square error ( RMSE ) is used for the performance assessment, which is defined as
RMSE = ( 1 Q K ) q = 1 Q k = 1 K ( θ ^ k , q θ k ) 2
with K the number of source signals, Q the number of Monte Carlo trials, and θ ^ k , q the estimate of the true DOA θ k of the q th Monte Carlo trial. Q = 500 is used in this paper. The Cramér–Rao lower bound (CRB) for the unconditional model is also given as a benchmark [17].
Figure 2 depicts the RMSE performance of Zheng’s method in [15] and the proposed method versus the signal-to-noise ratio ( SNR ) with K = 2 and L = 200 , in both distantly separated angles situation (denoted as general angles in Figure 2), where signals come from { 20 ° , 50 ° } , and a close angle situation, where signals come from { 24 ° , 25 ° } . It can be seen that, in the distantly separated angles situation, the estimation performance is comparable to Zheng’s method; in the close angles situation, due to the reconstruction of the covariance matrix and the exploitation of the uniform property of subarrays, the proposed method has higher resolution and better estimation performance than Zheng’s method.
Figure 3 depicts the RMSE performance of the two methods versus the snapshots number with SNR = 0   dB , in both the distantly separated angles situation (denoted as general angles in Figure 3) and close angles situation. As shown in the figure, in the distantly separated situation, the two methods can achieve similar estimation performance; in the close angles situation, the performance of Zheng’s method decreases greatly for small numbers of snapshots, as MUSIC depends on the accuracy of the estimated covariance matrix. On the contrary, benefiting from the reconstructed partially Toeplitz covariance matrix, the proposed method remains robust and reliable even in the case of a small snapshots number.
To investigate the resolution of the proposed method, two signals are assumed to come from the two close directions θ 1 and θ 2 = θ 1 + θ , respectively, where θ 1 is fixed at 20 ° , and θ is a small and controllable variable. The two signals are said to be successfully resolved if the two following equations are satisfied:
| θ ^ 1 θ 1 | < θ 2
| θ ^ 2 θ 2 | < θ 2
where θ ^ 1 and θ ^ 2 are the estimations of θ 1 and θ 2 , respectively [14,18]. Figure 4 shows the comparison of the resolution probability, which is calculated from the percentage of the success trials among 200 Monte Carlo trials, of Zheng’s method in [15] and the proposed method, with SNR = 0   dB and L = 200 . It can be seen that, benefiting from the partial Toeplitz averaging, the accuracy of the estimated covariance matrix is enhanced, and the proposed method exhibits much better resolution performance.

4.2. Computational Complexity

Based on root-MUSIC, the number of complex multiplications of the proposed method is O ( ( M 1 + M 2 ) 2 L + ( M 1 + M 2 ) 3 ) , which is obviously lower than Zheng’s method, which is O ( ( M 1 + M 2 ) 2 L + ( M 1 + M 2 ) 3 + ( M 1 + M 2 ) ( M 1 + M 2 K ) T ) , where L and T denote the number of snapshots and the times of spectral searching respectively. The computational complexity comparison versus the number of sensors (   M 1 + M 2 1 ) is given in Figure 5, with K = 2 and L = 200 . It can be seen that, without the exhaustive searching, the proposed method has much lower computational complexity with better resolution and estimation performance.

5. Conclusions

In this paper, an improved DOA estimation method with a UCLA is proposed. Exploiting the uniform property of the subarrays, we average the diagonal elements of the two estimated self-covariance matrices to enhance the accuracy of the estimated covariance matrices and the estimation performance. Besides, the polynomial root finding method is utilized instead of exhaustive searching to reduce the computational complexity. The simulation results show that the proposed method can achieve higher resolution and better estimation performance with lower computational complexity.

Author Contributions

The work presented here was carried out in collaboration between all authors. W.H., X.Y. and Y.W. proposed the general idea; W.H. performed the simulations and wrote the paper; X.Y. and Y.W. supervised the research, reviewed the manuscript and provided many valuable and constructive suggestions. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the China Scholarship Council (File No.201708440322) and in part by the 2017-2019 Sino-French CAI Yuanpei Program.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DOADirection-of-arrival
UCLAUnfolded coprime linear array
DOFDegree of freedom
MUSICMultiple Signal Classification
ESPRITEstimation of Signal Parameters via Rotational Invariance Techniques
ULAUniform linear array
CLACoprime linear array
AICAkaike Information Criterion
MDLMinimum Description Length
RMSERoot mean square error
CRBCramér–Rao lower bound
SNRSignal-to-noise ratio

References

  1. Krim, H.; Viberg, M. Two decades of array signal processing research. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  2. Bencheikh, M.L.; Wang, Y.; He, H. Polynomial root finding technique for joint DOA DOD estimation in bistatic MIMO radar. Signal Process. 2010, 90, 2723–2730. [Google Scholar] [CrossRef]
  3. Zhang, D.; Zhang, Y.; Zheng, G.; Feng, C.; Tang, J. ESPRIT-like two-dimensional DOA estimation for monostatic MIMO radar with electromagnetic vector received sensors under the condition of gain and phase uncertainties and mutual coupling. Sensors 2017, 17, 2457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ma, Y.; Miao, C.; Zhao, Y.; Wu, W. An MIMO radar system based on the sparse-array and its frequency migration calibration method. Sensors 2019, 19, 3580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  6. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef] [Green Version]
  7. Vaidyanathan, P.P.; Pal, P. Sparse sensing with co-prime samplers and arrays. IEEE Trans. Signal Process. 2011, 59, 573–586. [Google Scholar] [CrossRef]
  8. Kou, J.; Li, M.; Jiang, C. A robust DOA estimator based on compressive sensing for coprime array in the presence of miscalibrated sensors. Sensors 2019, 19, 3538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Meng, Z.; Zhou, W. Direction-of-arrival estimation in coprime array using the ESPRIT-based method. Sensors 2019, 19, 707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Qin, S.; Zhang, Y.D.; Amin, M.G. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar] [CrossRef]
  11. Liu, C.L.; Vaidyanathan, P.P. Remarks on the spatial smoothing step in coarray MUSIC. IEEE Signal Process. Lett. 2015, 22, 1438–1442. [Google Scholar] [CrossRef]
  12. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA estimation with combined MUSIC for coprime array. In Proceedings of the 2013 International Conference on Wireless Communications and Signal Processing, Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  13. Sun, F.; Lan, P.; Gao, B. Partial spectral search-based DOA estimation method for co-prime linear arrays. Electron. Lett. 2015, 51, 2053–2055. [Google Scholar] [CrossRef]
  14. Sun, F.; Gao, B.; Chen, L.; Lan, P. A low-complexity ESPRIT-based DOA estimation method for co-prime linear arrays. Sensors 2016, 16, 1367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Zheng, W.; Zhang, X.; Gong, P.; Zhai, H. DOA estimation for coprime linear arrays: An ambiguity-free method involving full DOFs. IEEE Commun. Lett. 2018, 22, 562–565. [Google Scholar] [CrossRef]
  16. Wax, M.; Kailath, T. Detection of signals by information theoretic criteria. IEEE Trans. Acoust. Speech Signal Process. 1985, 33, 387–392. [Google Scholar] [CrossRef] [Green Version]
  17. Stoica, P.; Nehorai, A. Performance study of conditional and unconditional direction-of-arrival estimation. IEEE Trans. Acoust. Speech Signal Process. 1990, 38, 1783–1795. [Google Scholar] [CrossRef]
  18. Zhou, C.; Zhou, J. Direction-of-arrival estimation with coarray ESPRIT for coprime array. Sensors 2017, 17, 1779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. System model.
Figure 1. System model.
Sensors 20 00218 g001
Figure 2. Root mean square error ( RMSE ) performance versus signal-to-noise ratio ( SNR ).
Figure 2. Root mean square error ( RMSE ) performance versus signal-to-noise ratio ( SNR ).
Sensors 20 00218 g002
Figure 3. RMSE performance versus snapshots number.
Figure 3. RMSE performance versus snapshots number.
Sensors 20 00218 g003
Figure 4. Resolution probability versus θ .
Figure 4. Resolution probability versus θ .
Sensors 20 00218 g004
Figure 5. Complexity comparison.
Figure 5. Complexity comparison.
Sensors 20 00218 g005

Share and Cite

MDPI and ACS Style

He, W.; Yang, X.; Wang, Y. A High-Resolution and Low-Complexity DOA Estimation Method with Unfolded Coprime Linear Arrays. Sensors 2020, 20, 218. https://doi.org/10.3390/s20010218

AMA Style

He W, Yang X, Wang Y. A High-Resolution and Low-Complexity DOA Estimation Method with Unfolded Coprime Linear Arrays. Sensors. 2020; 20(1):218. https://doi.org/10.3390/s20010218

Chicago/Turabian Style

He, Wei, Xiao Yang, and Yide Wang. 2020. "A High-Resolution and Low-Complexity DOA Estimation Method with Unfolded Coprime Linear Arrays" Sensors 20, no. 1: 218. https://doi.org/10.3390/s20010218

APA Style

He, W., Yang, X., & Wang, Y. (2020). A High-Resolution and Low-Complexity DOA Estimation Method with Unfolded Coprime Linear Arrays. Sensors, 20(1), 218. https://doi.org/10.3390/s20010218

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop