Next Article in Journal
Patterned Array of Poly(ethylene glycol) Silane Monolayer for Label-Free Detection of Dengue
Previous Article in Journal
A Bevel Gear Quality Inspection System Based on Multi-Camera Vision Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Low-Complexity ESPRIT-Based DOA Estimation Method for Co-Prime Linear Arrays

1
College of Information Science and Engineering, Shandong Agricultural University, Tai’an 271018, China
2
College of Communications Engineering, PLA University of Science and Technology, Nanjing 210007, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(9), 1367; https://doi.org/10.3390/s16091367
Submission received: 7 July 2016 / Revised: 18 August 2016 / Accepted: 22 August 2016 / Published: 25 August 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
The problem of direction-of-arrival (DOA) estimation is investigated for co-prime array, where the co-prime array consists of two uniform sparse linear subarrays with extended inter-element spacing. For each sparse subarray, true DOAs are mapped into several equivalent angles impinging on the traditional uniform linear array with half-wavelength spacing. Then, by applying the estimation of signal parameters via rotational invariance technique (ESPRIT), the equivalent DOAs are estimated, and the candidate DOAs are recovered according to the relationship among equivalent and true DOAs. Finally, the true DOAs are estimated by combining the results of the two subarrays. The proposed method achieves a better complexity–performance tradeoff as compared to other existing methods.

1. Introduction

Direction of arrival (DOA) estimation is a crucial problem in various applications, such as radar, sonar, and wireless communications [1]. Various DOA estimation methods have been studied in uniform linear arrays (ULAs), including multiple signal classification (MUSIC) [2] and estimation of signal parameters via rotational invariance technique (ESPRIT) [3]. In [4], a Khatri-Rao product-based real-valued sparse estimation method is proposed. With respect to the random errors of sensor position, a stochastic framework is established to find the probability density function of the DOA-estimates [5]. However, most of the traditional DOA estimation schemes have focused on ULA structure [6], which, in fact, is not an optimal array geometry.
Recently, sparse array geometry has drawn lots of attention due to its high resolution [7,8,9]. An eigenstructure-based direction-finding algorithm is proposed in [7] for sparse uniform Cartesian arrays. In [8], an ESPRIT-based estimation method is proposed for a sparse array with two different sizes of spatial invariances. Moreover, the literature [9] introduces a novel direction-finding algorithm for a multiscale sensor array, which presents multiple scales of spatial invariance. Except for these sparse structures, the nested arrays [10] and co-prime arrays [11,12,13] have also attracted great attention. The nested array suffers from the mutual coupling problem due to some closely located sensors, while the co-prime array has less of a mutual coupling effect since sensor elements are sparsely located. Consequently, we consider the co-prime array structure in this paper. Various DOA estimation methods have been proposed for co-prime arrays [14,15,16]. A projection-like method is proposed in [14] to estimate DOAs by combining the results of the two subarrays. By applying the MUSIC algorithm, a total spectral search-based method (TSS) is proposed in [15], which, however, suffers from the high complexity caused by the spectral search step. By limiting the searching region to a small sector, a partial spectral search-based method (PSS) is proposed in [16]. It is shown in [16] that the PSS method can achieve almost the same estimation accuracy as TSS, but with a substantially reduced complexity. Further, the PSS method is extended to co-prime planar array structures to estimate two dimensional DOAs in [17]. However, since the complexity is mainly caused by spectral search, the works [16,17] still have a heavy computational burden.
To this end, in this paper, we propose a low-complexity DOA estimation method for co-prime linear arrays. For each subarray of the co-prime array, the true DOAs are mapped into their respective equivalent angles impinging on a traditional ULA with half-wavelength inter-element spacing, which can be the obtained by ESPRIT. Then, according to the relationship among true and equivalent DOAs, the candidate DOAs are recovered immediately. Finally, the true DOAs can be uniquely estimated by seeking the common angles recovered by the two subarrays. Simulation results are provided to verify the effectiveness of the proposed method.

2. System Model

As shown in Figure 1, we consider a co-prime linear array consisting of two uniform linear subarrays, which are located in the same line. The first subarray has M 1 equal-spaced omnidirectional sensors with inter-element spacing M 2 λ λ 2 2 , while the second has M 2 equal-spaced omnidirectional sensors with inter-element spacing M 1 λ λ 2 2 . Here M 1 and M 2 are co-prime integers, and λ is the wavelength. The two subarrays share the first sensor, and consequently the co-prime array has M 1 + M 2 - 1 sensors.
Assume K ( K < min M 1 , M 2 ) uncorrelated narrowband signals impinge on the array. The received signal for the ith (i = 1, 2) subarray at time t 1 t T is
x i ( t ) = k = 1 K a i ( θ k ) s k ( t ) = A i s ( t ) + n i ( t )
where a i ( θ k ) = 1 , e - j M i ˜ π sin ( θ k ) , e - j M i - 1 M i ˜ π sin ( θ k ) T ( i + i ˜ = 3 ) is the steering vector for the kth source θ k , A i = a i ( θ 1 ) , , a i ( θ K ) is the steering matrix, s ( t ) = s 1 ( t ) , , s K ( t ) T is the source vector. The components of n i ( t ) are assumed to be independent and identically distributed additive white Gaussian noise with equal variance in each sensor, and are independent from the sources. · T denotes the transpose operation.
According to the property of sinusoid function, for each source θ k in the ith subarray, there exist some particular angles, denoted as θ i , k e q v , that satisfy π sin θ i , k e q v = M i ˜ π sin θ k + 2 n i π , where n i is an integer. From the definition of a i ( θ k ) , the angle θ i , k e q v can be regarded as the equivalent angle for the true θ k impinging on one traditional M i -element ULA with half-wavelength spacing, which generates the same steering vector as θ k . The relation between true θ k and its equivalent θ i , k e q v is thus given as
sin θ i , k e q v = M i ˜ sin θ k + 2 n i
The signal model in Equation (1) amounts to K equivalent sources from angles θ i , k e q v impinging on a traditional ULA simultaneously. The system model of the ith subarray of the co-prime array and its equivalent uniform linear array is shown in Figure 2.

3. Proposed DOA Estimation Method

The equivalent DOA equals to the true one only when M i ˜ = 1 . In the considered co-prime array, since M i ˜ > 1 , there exist multiple equivalent angles θ i , k e q v for each true θ k . The true DOA θ k may be one of the potential equivalent angles recovered by Equation (2), however, which one is the true DOA cannot be determined by each subarray alone. According to the property of coprimeness, the true DOAs can be uniquely estimated by finding the common angles in the two equivalent DOA sets generated by the two subarrays, respectively.
In the following, we aim to find the equivalent DOAs from Equation (1). To avoid the high complexity caused by spectral search, we propose to use the ESPRIT method.

3.1. Proposed Method

The steering matrix A i of the ith subarray can be rewritten as
A i = A i , 1 the last row = the first row A i , 2
where A i , 1 is the first M i - 1 rows of A i and A i , 2 is constructed from the last M i - 1 rows. The relationship between A i , 1 and A i , 2 is
A i , 2 = A i , 1 J
with J = diag e - j π sin θ i , 1 e q v , e - j π sin θ i , 2 e q v , , e - j π sin θ i , K e q v .
The eigen-decomposition of the array covariance matrix R i = E x i ( t ) x i H ( t ) of the ith subarray yields
R i = U s , i Λ s U s , i H + U n , i Λ n U n , i H
Here E · and · H stand for statistical expectation and Hermitian transpose, respectively. U s , i and U n , i denote the signal- and noise-subspace eigenvectors, respectively. Λ s and Λ n contain the corresponding eigenvalues.
The signal subspace eigenvector U s , i can be similarly partitioned as
U s , i = U s , i , 1 the last row = the first row U s , i , 2
where U s , i , 1 = A i , 1 G , U s , i , 2 = A i , 2 G = A i , 1 JG , and G is a full rank K × K matrix. Then, we have
G - 1 JG = U s , i , 1 H U s , i , 1 - 1 U s , i , 1 H U s , i , 2
The matrix U s , i , 1 H U s , i , 1 - 1 U s , i , 1 H U s , i , 2 has the same eigenvalues as J , while J contains the information of the equivalent DOAs. The equivalent DOAs can then be obtained by finding the eigenvalues of the matrix U s , i , 1 H U s , i , 1 - 1 U s , i , 1 H U s , i , 2 . Specifically, assume the K eigenvalues are denoted as d k ( k = 1 , 2 , , K ) , the equivalent DOAs of the ith subarray can be estimated as
θ ^ i , k e q v = arc sin d k d k π π
Let θ ^ i e q v = θ ^ i , 1 e q v , θ ^ i , 2 e q v , , θ ^ i , K e q v . The candidate true DOAs can be recovered by Equation (2); i.e.,
θ ^ i c a n d = arcsin 1 M i ˜ sin θ ^ i e q v - 2 n i
Here the relationship - M i ˜ sin θ ^ i e q v - 2 n i M i ˜ must hold due to the constraint of the sinusoid function. Hence, the set θ ^ i c a n d contains M i ˜ × K candidate DOAs, among which only K angles are true. According to the coprimeness of M 1 and M 2 , the true DOAs are obtained by finding the K common angles of the two candidate sets θ ^ 1 c a n d and θ ^ 2 c a n d . In practice, due to the effect of noise, the true DOAs are estimated by seeking for the nearest angle pairs.
Remark 1. 
The two candidate sets θ ^ 1 c a n d and θ ^ 2 c a n d both contain the K true DOAs. In particular cases, except for true DOAs, they may have some other common angles (e.g., the candidate DOAs of different sources in different subarrays may have common angles). This may generate more than K common angles and cause ambiguity as a result. Without loss of generality, we consider K = 2 sources, denoted as θ 1 and θ 2 . Let θ c both exist in the candidate set of source 1 in the ith subarray and the candidate set of source 2 in the i ˜ th subarray; i.e., except for θ 1 and θ 2 , they generate another common angle. The relations are denoted as
sin θ c = sin θ 1 + 2 m i ˜ M i ˜ sin θ c = sin θ 2 + 2 m i M i
where m i ˜ and m i are integers. Then we have the relation between θ 1 and θ 2 as
sin θ 1 - sin θ 2 = 2 m i M i - 2 m i ˜ M i ˜
It is observed from Equation (11) that only when m i 0 and m i ˜ 0 , the particular θ 1 , θ 2 pair will generate another common angle. For example, let us consider a co-prime linear array with M 1 = 5 and M 2 = 7 and assume θ 1 = 60.0 ° and θ 2 = 48.7 ° . The two DOAs satisfy the relationship (11). The corresponding ambiguous DOAs generated by each subarray are listed in Table 1. As can be seen, for each source, the estimate of the two subarrays can generate one unique common angle—i.e., the true DOA. However, when the two sources ( θ 1 = 60.0 ° and θ 2 = 48.7 ° ) are considered, expect for true DOAs, the estimate of the two subarrays may generate another two common angles (i.e., - 58.0 ° and 27.8 ° ). This violates the uniqueness and causes ambiguity as a result.
However, notice that the incident DOAs are in the continuous range between - 90 ° and 90 ° , and only some particular angle pairs θ 1 , θ 2 that satisfy Equation (11) may generate another common ambiguous angle, except for the true DOAs. Specifically, for arbitrarily given θ 1 , only some discrete angles for θ 2 follow the relationship Equation (11). Therefore, as compared to the continuous angle range, the number of θ 1 , θ 2 pairs that satisfy Equation (11) is rare. Therefore, the probability of the existence of a different common angle θ c approaches zero and it is consequently a small probability event. Even if this small probability event happens, it can be solved by the entire co-prime linear array. The beamforming-related techniques, such as the classical beamforming approach (CBF) [18] and the Capon’s approach (also known as the minimum variance distortionless response, MVDR) [19,20], can be utilized to eliminate the problem. Specifically, the CBF spatial spectrum of the entire co-prime array is
P C B F = a C H θ R C a C θ
and the MVDR spatial spectrum is
P M V D R = 1 a C H θ R C - 1 a C θ
where a C θ denotes the steering vector with respect to θ of the entire co-prime array, and R C is the covariance matrix of the entire co-prime array. By finding the peaks of P C B F in Equation (12) or P M V D R in Equation (13) the true DOAs can be estimated, and other ambiguous angles are eliminated. Thus, the ambiguity can be resolved successfully.
Remark 2. 
In general, there exist K common angles, which are the true DOAs. However, if the incident angles satisfy Equation (11), more than K common angles can be generated. We need to determine the number K ^ ( K ^ K ) of the common angles. Let N 1 and N 2 , respectively, denote the numbers of the relevant angles (including the true and ambiguous angles) for the first and second subarray, and let d i denote the absolute value of angle difference, where 1 i N 1 N 2 . By sorting d i in ascendant order, we have
d 1 d 2 d K d N 1 N 2
For common angles of the two subarrays, the angle differences are small (nearly 0), while they become larger for distinct angles. Define a decision variable as
D n = d n + 1 - d n d n , n = K , K + 1 , , N 1 N 2 - 1
Notice that when both d n and d n + 1 are the differences of common angles or distinct angles, the decision variable will almost keep small. Meanwhile, if d n is the difference of common angles and d n + 1 is the difference of distinct angles, the decision variable will become much larger. Thus, the number of common angles can be estimated as
K ^ = arg max n D n
When K ^ > K , the true K DOAs can be distinguished according to Equations (12) and (13).

3.2. Complexity Analysis

The complexity of different methods is provided in Table 2. The ESPRIT method with M sensors and K sources requires covariance matrix estimation and the eigen-decomposition of both the covariance matrix R i and the matrix U s , i , 1 H U s , i , 1 - 1 U s , i , 1 H U s , i , 2 , and the corresponding complexities are O M 2 T , O M 3 , and O 3 M K 2 + 2 K 3 , respectively. The resulting complexity of ESPRIT is given as O M 2 T + M 3 + 3 M K 2 + 2 K 3 . The complexity of the proposed method is given similarly. As is shown in Table 2, the complexity of the proposed method is lower than that of ESPRIT.

4. Results

In this section, we compare the proposed method with PSS [16] in co-prime arrays with M 1 = 5 and M 2 = 7 . For fair comparison, the results of the traditional ESPRIT are also provided for two types of half-wavelength spacing ULAs with M 1 + M 2 1 = 11 and M 1 M 2 − min M 1 , M 2 + 1 = 31 sensors, respectively. The former has the same sensor number as the considered co-prime array, but with less aperture length, denoted as ESPRIT with the same sensor number (ESPRIT-SSN). The latter has the same aperture length, but with more sensor elements, denoted as ESPRIT with the same aperture length (ESPRIT-SAL). The Cramer-Rao bound (CRB) for the co-prime array geometry is also given as a benchmark [21]. The searching grid for PSS is set as 0.1 ° .
In the first test, we compare the average root mean square error (RMSE) of different methods. Two uncorrelated sources are assumed to impinge on the array from directions 21 ° and 41 ° . Figure 3 plots the RMSE performance versus signal-to-noise ratio (SNR) via 200 Monte Carlo simulations with T = 200 . As is shown, under the same sensor number condition, the proposed method outperforms ESPRIT-SSN greatly. With respect to the same aperture length, the proposed method can achieve almost the same estimation accuracy as ESPRIT-SAL. The proposed method is inferior to PSS at low SNR, and gradually exceeds PSS with the increase of SNR. Regarding the complexity, the ESPRIT-SSN, ESPRIT-SAL, PSS, and proposed method requires O 2 . 57 × 10 4 , O 2 . 22 × 10 5 , O 3 . 17 × 10 4 , and O 1 . 55 × 10 4 complexities, respectively. The average running time for the four methods are 5.74 s, 15.60 s, 12.98 s, and 2.66 s. The complexity of the proposed method is lower than that of others, implying that the proposed method needs less memory. Especially, to achieve the same estimation accuracy, the proposed method only requires about 10% computational complexity as compared to ESPRIT-SAL. Therefore, the proposed method can achieve a better complexity–performance tradeoff. Further, Figure 4 depicts the RMSE performance versus snapshot number with SNR = 5 dB. The performance of the proposed method is almost the same as ESPRIT-SAL and is superior to ESPRIT-SSN, but is slightly weaker than that of PSS. However, the performance gap becomes smaller with the increase of snapshot number.
In the second test, we investigate the resolution probability performance. The two sources are said to be resolvable [22] if both θ ^ 1 - θ 1 and θ ^ 2 - θ 2 are smaller than θ 1 - θ 2 θ 1 - θ 2 2 2 , where θ k and θ ^ k ( k = 1 , 2 ) are the true and estimated DOAs, respectively. We set θ 1 = 20 ° and θ 2 = θ 1 + Δ θ , where Δ θ is the control variable. We fix S N R = 5 dB and T = 200 and plot the resolution probability against Δ θ in Figure 5. As can be seen, the resolution ability is enhanced with the increase of Δ θ . The proposed method provides almost the same resolution performance as ESPRIT-SAL, but with substantially reduced complexity. The proposed method provides the best resolution performance as compared to ESPRIT-SSN and PSS.
To eliminate the ambiguity for some particular angle pairs that satisfy Eqaution (11), we give the ambiguity check results in Table 3. For θ 1 = 48.7 ° and 60.0 ° as in Table 1, the estimate of two subarrays may generate another two common angles except for true DOAs; i.e., - 58.0 ° and 27.8 ° . We can see that the true DOAs are checked successfully by both CBF and MVDR. Notice that the spatial output powers for MVDR at true DOAs are much higher than that of other common DOAs, thus the MVDR-based check approach has superior resolution ability. However, since an inverse operation of the covariance matrix is required, the MVDR-based approach has a heavy computation burden. Therefore, the CBF-based approach is more feasible when real-time processing is required.

5. Conclusions

In this paper, we have proposed a low-complexity ESPRIT-based DOA estimation method for co-prime arrays. We firstly map the true DOAs into several equivalent angles for each source, which are regarded as the DOAs impinging on traditional uniform linear array. We obtain the equivalent DOAs by ESPRIT and then recover the candidate DOAs according to their relationship. Finally, we uniquely estimate the true DOAs by finding the common angles of the candidate DOAs of the two subarrays. Simulation results show that the proposed method achieves an improved complexity–performance tradeoff in terms of estimation accuracy and resolution ability, as compared to other existing methods.

Acknowledgments

This work was supported in part by Shandong Province High School Science & Technology Fund Planning Project (J12LN02), Research Fund for the Doctoral Program of Higher Education of China (20123702120016), and Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period (2011BAD32B02).

Author Contributions

Fenggang Sun contributed in the conception and design of the study. Peng Lan wrote and revised the manuscript. Bin Gao carried out the simulation and participated in the design of the study. Lizhen Chen helped perform the analysis with constructive discussions and helped to draft the manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DOADirection of arrival
ULAUniform linear array
ESPRITEstimation of signal parameters via rotational invariance technique
MUSICMultiple signal classification
TSSTotal spectral search
PSSPartial spectral search
CBFClassical beamforming
MVDRMinimum variance distortionless response
ESPRIT-SSNESPRIT with the same sensor number
ESPRIT-SALESPRIT with the same aperture length
CRBCrame-Rao bound
RMSERoot mean square error
SNRSignal-to-noise ratio

References

  1. Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  2. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  3. 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]
  4. Chen, T.; Wu, H.; Zhao, Z. The real-valued sparse direction of arrival (DOA) estimation based on the Khatri-Rao product. Sensors 2016, 16, 693. [Google Scholar] [CrossRef] [PubMed]
  5. Inghelbrecht, V.; Verhaevert, J.; van Hecke, T.; Rogier, H. The influence of random element displacement on DOA estimates obtained with (Khatri-Rao-) root-MUSIC. Sensors 2014, 14, 21258–21280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ye, Z.; Xu, X. DOA estimation by exploiting the symmetric configuration of uniform linear array. IEEE Trans. Antennas Propag. 2007, 55, 3716–3720. [Google Scholar] [CrossRef]
  7. Zoltowski, M.; Wong, K. Closed-form eigenstructure-based direction finding using arbitrary but identical subarrays on a sparse uniform Cartesian array grid. IEEE Trans. Signal Process. 2000, 48, 2205–2210. [Google Scholar] [CrossRef]
  8. Wong, K.; Zoltowski, M. Direction-finding with sparse rectangular dual-size spatial invariance array. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 1320–1336. [Google Scholar] [CrossRef]
  9. Miron, S.; Song, Y.; Brie, D.; Wong, K. Multilinear direction finding for sensor-array with multiple scales of invariance. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2057–2070. [Google Scholar] [CrossRef]
  10. Pal, P.; Vaidyanathan, P. Nested arrays: A novel approach to array processing with enhanced degrees of freedom. IEEE Trans. Signal Process. 2010, 58, 4167–4181. [Google Scholar] [CrossRef]
  11. Vaidyanathan, P.; Pal, P. Sparse sensing with co-prime samplers and arrays. IEEE Trans. Signal Process. 2011, 59, 573–586. [Google Scholar] [CrossRef]
  12. Shen, Q.; Zhang, Y.; Amin, M.G. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar]
  13. Pal, P.; Vaidyanathan, P.P. Coprime sampling and the MUSIC algorithm. In Proceedings of the Digital Signal Processing Workshop and IEEE Signal Processing Education Workshop (DSP/SPE), Sedona, AZ, USA, 4–7 January 2011; pp. 289–294.
  14. Weng, Z.; Djuric, P. A search-free DOA estimation algorithm for coprime arrays. Digit. Signal Process. 2014, 24, 27–33. [Google Scholar] [CrossRef]
  15. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA estimation with combined MUSIC for coprime array. In Proceedings of the International Conference on Wireless Communications & Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5.
  16. 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]
  17. Wu, Q.; Sun, F.; Lan, P.; Ding, G.; Zhang, X. Two-dimensional direction-of-arrival estimation for co-prime planar arrays: A partial spectral search approach. IEEE Sens. J. 2016, 16, 5660–5670. [Google Scholar] [CrossRef]
  18. Godara, L.C. Application of antenna arrays to mobile communications. II. Beam-forming and direction-of-arrival considerations. Proc. IEEE 1997, 85, 1195–1245. [Google Scholar] [CrossRef]
  19. Capon, J. High resolution frequency-wavenumber spectrum analysis. Proc. IEEE 1969, 57, 1408–1418. [Google Scholar] [CrossRef]
  20. Benesty, J.; Chen, J.; Huang, Y. A generalized MVDR spectrum. IEEE Signal Process. Lett. 2005, 12, 827–830. [Google Scholar] [CrossRef]
  21. 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]
  22. Stoica, P.; Gershman, A.B. Maximum-likelihood DOA estimation by data-supported grid search. IEEE Signal Process. Lett. 1999, 6, 273–275. [Google Scholar] [CrossRef]
Figure 1. The system model of co-prime linear array.
Figure 1. The system model of co-prime linear array.
Sensors 16 01367 g001
Figure 2. The System model of (a) the ith subarray and (b) its corresponding equivalent array of the co-prime linear array. θ k is the true DOA for the kth source, and θ i , k e q v is the equivalent angle.
Figure 2. The System model of (a) the ith subarray and (b) its corresponding equivalent array of the co-prime linear array. θ k is the true DOA for the kth source, and θ i , k e q v is the equivalent angle.
Sensors 16 01367 g002
Figure 3. Root mean square error (RMSE) vs. signal-to-noise ratio (SNR), T = 200 . CRB: Cramer–Rao bound; ESPRIT: estimation of signal parameters via rotational invariance technique; ESPRIT-SAL: ESPRIT with the same aperture length; ESPRIT-SSN: ESPRIT with the same sensor number.
Figure 3. Root mean square error (RMSE) vs. signal-to-noise ratio (SNR), T = 200 . CRB: Cramer–Rao bound; ESPRIT: estimation of signal parameters via rotational invariance technique; ESPRIT-SAL: ESPRIT with the same aperture length; ESPRIT-SSN: ESPRIT with the same sensor number.
Sensors 16 01367 g003
Figure 4. RMSE vs. Snapshot Number, S N R = 5 dB.
Figure 4. RMSE vs. Snapshot Number, S N R = 5 dB.
Sensors 16 01367 g004
Figure 5. Resolution probability versus Δ θ , with θ 1 = 20 ° , θ 2 = θ 1 + Δ θ .
Figure 5. Resolution probability versus Δ θ , with θ 1 = 20 ° , θ 2 = θ 1 + Δ θ .
Sensors 16 01367 g005
Table 1. Two directions of arrival (DOAs) and their corresponding ambiguous DOAs in each subarray.
Table 1. Two directions of arrival (DOAs) and their corresponding ambiguous DOAs in each subarray.
θ 1 = 60.0 ° Subarray 1 with M 1 = 5 58.0,   −34.2,   −16.1,   0.51,   17.1,   35.5,   60.0
Subarray 2 with M 2 = 7 −47.2,   −19.5,   3.8,   27.8,   60.0
θ 2 = 48.7 ° Subarray 1 with M 1 = 5 −74.4,   −42.6,   −23.1,   −6.1,   10.4,   27.8,   48.7
Subarray 2 with M 2 = 7 −58.0,   −26.7,   −2.8,   20.6,   48.7
Table 2. Computational complexity comparison.
Table 2. Computational complexity comparison.
ESPRIT O ( M 2 T + M 3 + 3 M K 2 + 2 K 3 )
PSS [16] O ( ( M 1 2 + M 2 2 ) T + M 1 3 + M 2 3 + J M 2 M 1 ( M 1 - K ) + J M 1 M 2 ( M 2 - K ) )
Proposed O ( ( M 1 2 + M 2 2 ) T + M 1 3 + M 2 3 + 3 ( M 1 + M 2 ) K 2 + 4 K 3 )
Note: M = M 1 + M 2 - 1 , J is the number of sampling points. PSS: partial spectral search.
Table 3. Ambiguity check results of classical beamforming (CBF) and minimum variance distortionless response (MVDR) with T = 200 and S N R = 5 dB.
Table 3. Ambiguity check results of classical beamforming (CBF) and minimum variance distortionless response (MVDR) with T = 200 and S N R = 5 dB.
θ - 58.0 ° 27.8 ° 48.7 ° 60.0 °
P C B F | P M V D R 0.6122|0.09310.5974|0.09261.3082|0.95611.2529|0.9454
True or falsefalse|falsefalse|falsetrue|truetrue|true

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/s16091367

AMA Style

Sun F, Gao B, Chen L, Lan P. A Low-Complexity ESPRIT-Based DOA Estimation Method for Co-Prime Linear Arrays. Sensors. 2016; 16(9):1367. https://doi.org/10.3390/s16091367

Chicago/Turabian Style

Sun, Fenggang, Bin Gao, Lizhen Chen, and Peng Lan. 2016. "A Low-Complexity ESPRIT-Based DOA Estimation Method for Co-Prime Linear Arrays" Sensors 16, no. 9: 1367. https://doi.org/10.3390/s16091367

APA Style

Sun, F., Gao, B., Chen, L., & Lan, P. (2016). A Low-Complexity ESPRIT-Based DOA Estimation Method for Co-Prime Linear Arrays. Sensors, 16(9), 1367. https://doi.org/10.3390/s16091367

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