Next Article in Journal
A MEMS-Enabled Deployable Trace Chemical Sensor Based on Fast Gas-Chromatography and Quartz Enhanced Photoacousic Spectoscopy
Next Article in Special Issue
Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization
Previous Article in Journal
An RSSI-Based Low-Power Vehicle-Approach Detection Technique to Alert a Pedestrian
Previous Article in Special Issue
Joint Angle and Frequency Estimation Using One-Bit Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SPICE-ML Algorithm for Direction-of-Arrival Estimation

1
College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
2
The 14th Research Institute Electronics Technology Group Corporation, Nanjing 210000, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 119; https://doi.org/10.3390/s20010119
Submission received: 29 November 2019 / Revised: 17 December 2019 / Accepted: 20 December 2019 / Published: 24 December 2019

Abstract

:
Sparse iterative covariance-based estimation, an iterative direction-of-arrival approach based on covariance fitting criterion, can simultaneously estimate the angle and power of incident signal. However, the signal power estimated by sparse iterative covariance-based estimation approach is inaccurate, and the estimation performance is limited to direction grid. To solve the problem above, an algorithm combing the sparse iterative covariance-based estimation approach and maximum likelihood estimation is proposed. The signal power estimated by sparse iterative covariance-based estimation approach is corrected by a new iterative process based on the asymptotically minimum variance criterion. In addition, a refinement procedure is derived by minimizing a maximum likelihood function to overcome the estimation accuracy limitation imposed by direction grid. Simulation results verify the effectiveness of the proposed algorithm. Compared with sparse iterative covariance-based estimation approach, the proposed algorithm can achieve more accurate signal power and improved estimation performance.

Graphical Abstract

1. Introduction

Direction-of-arrival (DOA) refers to the process of retrieving the direction information of several electromagnetic waves from the outputs of a number of receiving antennas that form a sensor array. DOA estimation is a major topic in array signal processing and plays a very important role in many applications, e.g., wireless communications, medical imaging, and radar systems [1,2,3],. Moreover, wireless sensor networks (WSNs) are an emerging paradigm in wireless communications [4]. In WSNs one of the most significant challenges is localization in which DOA is a major parameter to estimate [5,6]. The research on DOA estimation is mainly divided into two aspects: (1) array geometry, such as the recently proposed nested array [7] and coprime array [8,9]. (2) DOA estimation algorithm with higher estimation accuracy. In general, the algorithm applied to DOA estimation may be classified into three parts: nonparametric, parametric, and semiparametric [10].
A main example of the nonparametric methods is delay-and-sum (DAS) beamformer [10] in which the received signal from each sensor are weighted and delayed so as to focus on different points in space. It is a data-independent estimation technique which is traditionally adopted due to its low computational burden and high-signal-to-noise ratio (SNR) properties. However, data-independent approaches suffer from leakage problem. The local leakage will reduce the resolution, which makes the DAS beamformer unable to distinguish the incident signals with close frequency components. Global leakage will lead to false alarms. Adaptive Capon beamformer can improve the DAS method, but it is limited to independent signal [11]. Parametric methods, especially subspace methods such as multiple signal classification (MUSIC) [12], estimation of signal parameters via rotational invariance techniques (ESPRIT) [13], where parameters are estimated by studying the subspace of the data covariance. They can provide higher resolution than nonparametric methods, but adimit two shortcomings: (1) they require the knowledge of the number of source, and (2) work well only when the sources are independent. Recently, these two kinds of approaches have been combined, leading to the so-called semiparametric approache. The recently developed iterative adaptive algorithm (IAA) [14] largely eliminates the leakage problem of the beamformer and is robust to coherent sources. However, the IAA algorithm still has resolution limitations, especially at low SNR situations. Stoica et al. have recently proposed a user parameter-free sparse iterative covariance-based estimation (SPICE) approach in [15,16], based on minimizing a covariance fitting criterion. SPICE is able to provide excellent resolution and low sidelobe levels while maintaining robustness to coherent sources. Moreover, the power of the signal can be estimated simultaneously. However, as with most power-based sparse methods, the estimation accuracy of the SPICE algorithm is limited to the direction grid. In addition, the signal power estimated by SPICE is inaccurate, especially in the case of coherent signal. In [17], a method named likelihood-based estimation of sparse parameters (LIKES) is proposed, where the same covariance fitting criterion is adopted as SPICE algorithm. The difference is that the weights in LIKES are adaptive rather than constant in SPICE. Compared with the SPICE algorithm, LIKES can obtain more accurate DOA estimation, but at the expense of computational cost. The DOA estimation accuracy will be limited to direction grid because the spatial domain is also discretized in LIKES algorithm. In [18], an iterative algorithm based on the combination of Capon and SPICE (C-SPICE) is proposed, and a mobile average initialization technology is introduced to realize DOA estimation by using the spatial spectrum information estimated in the previous snapshot in the next snapshot.
DAS has low computational burden and accurate signal power estimation, but low resolution. The resolution of IAA is improved, but still not high enough to distinguish two signals which are close to each other. SPICE has high resolution, but the signal power estimated by SPICE is not accurate, and the DOA estimation accuracy is limited to the direction grid. In order to obtain an algorithm with high resolution, accurate signal power estimation and high estimation accuracy, we propose a SPICE and maximum likelihood (ML) estimation approach (SPICE-ML) in this paper. An iterative correction process is derived through the asymptotically minimum variance (AMV) criterion [19,20] so that the power of the signal estimated are more accurate. Moreover, a maximum likelihood cost function is used to refine the DOA estimation, thereby combating the limit of estimation accuracy caused by the direction grid. Numerical simulations are designed to show the effectiveness and superiority in estimation accuracy of the proposed algorithm.
This paper is organized as follows. The signal model for DOA estimation is introduced in Section 2 for far-field, narrowband sources. The concept of sparse representation is also discussed in Section 2. The SPICE algorithm is introduced in Section 3, including the covariance fitting criterion and updating formulas. Section 4 is the main highlight of this paper in which the SPICE-ML is proposed. Compared with SPICE algorithm, SPICE-ML can receive more accurate signal power estimation and DOA estimation accuracy. Numerical simulations are presented in Section 5, and the conclusion and some future research directions are drawn in Section 6.
Notations: Vectors, matrices are denoted by boldface lowercase letters and boldface uppercase letters, respectively. · T , · * , · 1 and · H represent the transpose, the conjugate, the inverse and the complex conjugate transpose of vectors or matrices. ⊗ is the Kronecker product and tr · represents the trace of a matrix. E · , det · and vec · denote the operator of expectation, determinant and vectorization. · denotes the Euclidean norm for vectors and the Frobenius norm for matrices.

2. Signal Model

Assume N far-field, uncorrelated narrowband signals from directions Ω = θ 1 , θ N impinging on an uniform line array (ULA), the number of ULA is M and the inter-element spacing equals to a half of the signal wavelength. Let θ k k = 1 K denote a direction grid that covers Ω , where K N , θ k 90 90 . Assume the direction gird is fine enough so that the true location parameters of the sources lie on (or, practically, close to) the grid. The signal model is shown in Figure 1.
Then the output of the array can be modeled as:
x t = k = 1 K a θ k s k t + n t = A θ S t + n t , t = 1 , , T
where A θ = a θ 1 , a θ 2 , , a θ K denotes the array manifold matrix and a θ k is the steering vector corresponding to θ k . T denotes the number of snapshots. S t = s 1 t , s 2 t , , s K t T contains the K unknown complex-valued signals, and n t is the additional noise term. We assume that E n t n H t = σ I M , where σ is the power of noise, and I M is a M × M identical matrix. Let us further assume that E S t S H t = P k I M , P K = Diag p 1 , p 2 , p K , where p k represents the unknown signal power at θ k . Moreover the signal S t and noise n t are assumed to be statistically independent. Therefore, we have the covariance matrix of x t
R = A P K A H + σ I M .
In practice, this covariance matrix is usually estimated by the sample covariance matrix
R ^ = 1 T t = 1 T x t x H t ,
Note that only few sources exist in practice, therefore in the signal matrix
S = s 1 1 s 1 2 s 1 T s K 1 s K 2 s K T ,
only a small number of rows are different from zeros. Thus, the DOA estimation problem is then changed to decide from x t which rows of the signal matrix S are non-zero.

3. The SPICE Algorithm

The SPICE algorithm is derived from a robust covariance fitting criterion. It has a sound statistical foundation, it does not require any hyper-parameters, and yet it has global convergence properties. Moreover, the signal and noise power are estimated in a natural manner. Although the signal are assumed to be independent, the SPICE algorithm is still robust to coherent signal [15]. In this paper, we only consider the case in which T > M and the variances of noise are identical (see [15] for the case of T < M and noise with different variances).
Equation (2) can be rewrite as
R = A a u g PA a u g H ,
where
A a u g = a θ 1 a θ 2 a θ K I M = a 1 a 2 a K a K + 1 a K + M ,
P = p 1 0 0 0 p K 0 0 σ 1 0 0 σ M = p 1 0 0 0 0 p K 0 0 p K + 1 0 0 p K + M .
The SPICE algorithm considers the following weighted covariance fitting criterion:
f = R 1 / 2 R ^ R R ^ 1 / 2 2 ,
a simple calculation shows that Equation (8) can be expressed as:
f = tr R 1 R ^ R R ^ 1 R ^ R = tr R 1 R ^ I I R 1 R ^ = tr R 1 R ^ + tr R ^ 1 R 2 M ,
where
tr R ^ 1 R = k = 1 K + M p k a k * R ^ 1 a k .
The minimization of f obtained from Equations (9) and (10) is equivalent to the minimization of the function
g = tr R ^ 1 / 2 R 1 R ^ 1 / 2 + k = 1 K + M a k * R ^ 1 a k p k ,
The minimization problem with respect to p k in Equation (11) is a semidefinite program (SDP) and therefore is a convex problem [21]. However, the calculation of Equation (11) as a SDP is very computationally intensive. It can be seen from Equation (10) that a consistent estimation tr R ^ 1 R = M can be obtained when T tends to be infinite. Hence, the probem of minimizing g can be reformulated as the following constrained minimization problem:
min p k 0 tr R ^ 1 / 2 R 1 R ^ 1 / 2 s . t . k = 1 K + M ω k p k = 1 ,
where
ω k = a k * R ^ 1 a k / M .
The minimization of the objective in Equation (12) can be solved by means of a cyclic algorithm [15], leading to the updated formulas of the SPICE algorithm
p k i + 1 = p k i a k * R 1 i R ^ 1 / 2 ω k 1 / 2 ρ i k = 1 , , K ,
ρ i = k = 1 K ω k 1 / 2 p k i a k * R 1 i R ^ 1 / 2 + γ 1 / 2 σ i R 1 i R ^ 1 / 2 ,
σ i + 1 = σ i R 1 i R ^ 1 / 2 γ 1 / 2 ρ i , γ = k = K + 1 K + M ω k .
where index i denotes the number of iteration. The algorithm can be initialized by means of the DAS method. The cyclic operation of the cyclic algorithm makes the objective function in Formula (12) monotonically decrease, and the minimization process is a convex problem. Therefore, the result of SPICE algorithm has global convergence. In [22], it is proved that the limit point of the iterative process of SPICE is the global solution of the minimization problem (12) under the weak condition where p k 0 > 0 and R ( i ) 0 in each iteration.

4. The SPICE-ML Algorithm

SPICE algorithm has a bias in estimating the signal power, especially in the case of coherent sources (as shown in the simulation part). To solve this problem, an iteration calibration process of signal power is obtained by using AMV criterion. After applying the vectorization operator to the matrix R , we can receive the vector r p
r p = vec R = A ˜ p ,
where A ˜ = A ˜ 1 , a ¯ K + 1 , A ˜ 1 = a ¯ 1 , , a ¯ K , a ¯ k = a k * a k , k = 1 , K . a ¯ K + 1 = vec I M , p = p 1 , , p K . Note that the Gaussian circular asymptotic covariance matrix r ^ = def vec R ^ is give by [20]
C r = R * R .
Suppose that p can be identifiable from r p . According to [19,20], it can be proved that the covariance matrix Cov p Alg is bounded below by the following real symmetric positive definite matrix
Cov p Alg A ˜ d H C r 1 A ˜ d 1 ,
where A ˜ d = d r p / d p , d · / d · is the operation of differential. In addition, we can obtain this lower bound by minimizing the following AMV criterion
p ^ = arg min p f p f p = def r ^ r p H C r 1 r ^ r p .
From Equations (18) and (20), the estimation of p and σ are given by the following iterative s
p ^ k i + 1 = a k H R 1 i R ^ R 1 i a k a k H R 1 i a k 2 + p k i 1 a k H R 1 i a k , k = 1 , K ,
σ ^ i + 1 = tr R 2 i R ^ + σ ^ i tr R 2 i tr R 1 i / tr R 2 i ,
where R i = A P i A H + σ ^ i I M , P i = Diag p ^ 1 i , , p ^ K i . The initialization of p ^ k 0 k = 1 K , σ ^ 0 is provided by the SPICE algorithm. In practice the noise power σ k , k = 1 , M estimated by the SPICE algorithm may be different. [23] shows that the degradation of accuracy comapared with that achieved by imposing σ k 2 σ 2 , k , is not significant. Through the iterative Formula (21), the signal power are continuously corrected.
The estimation accuracy of the SPICE algorithm is limited to the fineness of the direction grid. A coarse grid would lead to the degradation of estimation accuracy and a problem of computational complexity arises if a high-density grid is employed. Therefore, how to choose the direction grid becomes a difficult problem. In order to overcome the limit of estimation accuracy caused by the direction grid, the DOA estimation Ω ^ = θ ^ 1 θ ^ N is refined by iteratively minimizing a stochastic ML cost function.
The stochastic negative log-likelihood function of x t can be expressed as [24]
L θ = ln det R + tr R 1 R ^ .
The covariance matrix of noise is defined as
Q k = R p k a k a k H , k = 1 , , K .
Applying matrix inverse lemma to Formula (24)
R 1 = Q k 1 p k β k b k b k H , k = 1 , , K .
where b k = Q k 1 a k , β k = 1 + p k a k H Q k 1 a k 1 .
Besides,
tr R 1 R ^ = tr Q 1 R ^ p k β k b k H R ^ b k .
Then use the algebraic identity det I + AB = det I + BA and obtain
ln det R = ln det Q k + p k a k a k H = ln 1 + p k a k H Q k 1 a k det Q k = ln det Q k ln β k .
From Equations (23), (26) and (27), we can derive
L θ = ln det Q k + tr Q k 1 R ^ ln β k + p k β k b k H R ^ b k = L θ k l θ k ,
l θ k = ln 1 1 + p k α 1 , k + p k α 2 , k N 1 + p k + α 1 , k ,
where α 1 , k = a k H Q k 1 a k , α 2 , k N = a k H Q k 1 R ^ Q k 1 a k 1 , Q k = R p k a k a k H , k = 1 , , K .
The ML cost function can be decomposed into two parts: L θ k , the marginal likelihood function with parameter θ k excluded, and l θ k concerning θ k . Therefore, the ML cost function (23) with respect to θ k is equivalent to the ML function (29). Although the calculation of multi-dimensional optimization problem is very complex, the minimization problem in Equation (29) can be effectively solved by using the Nelder–Mead algorithm [25]. Moreover, the Nelder–Mead algorithm has already been built into the “fminsearch” function in MATLAB. Suppose that we have obtained the corrective parameters p ^ k i k = 1 K , σ ^ i and θ ^ k 0 k = 1 N , then the refinement result can be achieved by minimizing Formula (29).
The SPICE-ML algorithm is summarized in Table 1.

5. Simulation Results

We evaluated the performance of the proposed SPICE-ML algorithm and compared it with DAS, IAA, and SPICE in this section. We employed a uniform linear array with 12 sensors, and the inter-element spacing was half-wavelength. The direction grid θ k k = 1 K uniformly covered the entire DOA space Ω = 90 90 with a step size of 0.5°. The SNR was defined as
SNR = 10 log 10 p a v g σ dB
where p a v g denotes the average power of the signal. The iteration termination condition is set to p ^ i + 1 p ^ i 2 / p ^ i 2 < 10 3 .
Three signals with power 10 dB, 8 dB, and 5 dB from θ 1 = 45.3 , θ 2 = 38.7 , θ 3 = 30.8 impinged on the array. The additional noise was Gaussian white noise with equal power, and the SNR is set to 20 dB. The number of snapshots was given by T = 50 . Both independent and coherent sources are shown in Figure 2 and Figure 3. For the case of coherent sources, the signals at θ 1 and θ 3 shared the same phases but were independent to the signal at θ 2 .
It can be seen from the simulation results of Figure 2a and Figure 3a, the DAS method failed to seperate the two close sources at θ 1 and θ 2 owing to the smearing effects and limited resolution. The IAA algorithm could significantly reduce the smearing effect so that the sidelobes in Figure 2b and Figure 3b were lower, but the resolution was still not high. The SPICE algorithm was capable of resolving the three sources, but the signal power estimated was not accurate, especially in the case of coherent sources, as shown in Figure 2c and Figure 3c. The SPICE-ML algorithm proposed in this paper further corrected the signal power estimated by the SPICE algorithm, leading to a more accurate signal power estimation. Figure 2d and Figure 3d illustrate the effectiveness and high resolution of the proposed SPICE-ML algorithm.
Next, we evaluate the root mean square error (RMSE) of DOA estimation through Monte-Carlo simulations. The defination of RMSE is:
RMSE = 1 L K l = 1 L k = 1 K ( θ ^ k , l θ k ) 2 ,
where θ ^ k , l is the estimated DOA of the k signal in the l-th Monte-Carlo trial, and L is the total number of Monte-Carlo trials. Assume that three signals are randomly located at 60 60 . Considering both independent and coherent sources, respectively. The number of Monte-Carlo trials is 500. The RMSE curves are shown in Figure 4 and Figure 5. Moreover, we display the values of RMSE at each SNR and snapshot in Table 2, Table 3, Table 4 and Table 5.
From Figure 4a,b, it can be seen that in the case of independent signal the RMSE of all algorithms decreased with the increase of SNR and snapshots, and the RMSE of SPICE-ML proposed in this paper was the lowest, which illustrates that the proposed SPICE-ML approach outperformed the other methods. From the Table 2 and Table 3, specifically, compared with SPICE, SPICE-ML could efficiently reduce the RMSE vs. SNR from 0.147°–0.330° to 0.022°–0.238°, the RMSE vs. snapshots from 0.065°–0.150° to 0.051°–0.113° in the case of independent signal. From Figure 5a,b, we can see the estimation accuracy of the SPICE algorithm declined dramatically in the case of coherent signal compared with that in the case of independent signal. However, source coherence did not degrade the DAS and IAA algorithm. Notably, SPICE-ML still offered the best estimation performance in the case of coherent signal, because the angle estimated by SPICE was further refined by minimizing an ML function. From the Table 4 and Table 5, we can see that compared with SPICE, the RMSE vs. SNR and snapshots were reduced from 0.195°–2.777° to 0.023°–0.244°, from 0.276°–1.950° to 0.053°–0.114°, respectively. We can also observe that there existed plateau effects for DAS, IAA, and SPICE. This is because all these methods estimate DOA by means of selecting one element from a fixed set of direction grid, but there always exists an estimation bias no matter how fine the direction grid is. Theoretically, bias can be reduced by selecting a dense gird. However, a dense grid will lead to large computation cost and is not applicable in practice. On the contrary, the SPICE-ML algorithm would not suffer from the plateau effect because of a refinement DOA estimation procedure based on minimizing a ML cost function.

6. Conclusions

In this paper, we combined the SPICE DOA estimation algorithm and maximum likelihood estimation and proposed a SPICE-ML algorithm. Compared with the SPICE algorithm, the SPICE-ML method derive an iterative correction procedure for signal power estimation based on the AMV criterion, and combat the limitation of the direction grid by minimazing a Maximum Likelihood cost function. The simulation results verify the effectiveness of the proposed algorithm. In addition, the superiority of the SPICE-ML algorithm in DOA estimation accuracy is illustrated compared with the DAS, IAA, and SPICE algorithms. Specifically, comparing with SPICE, SPICE-ML can efficiently reduce the RMSE vs. SNR from 0.147°–0.330° to 0.022°–0.238° in the case of independent signal, from 0.195°–2.777° to 0.023°–0.244° in coherent signal situation, the RMSE vs. snapshots from 0.065°–0.150° to 0.051°–0.113° in the case of independent signal, from 0.276°–1.950° to 0.053°–0.114° in coherent signal situation. In the future, it will be of interest to study the performance of SPICE-ML algorithm in the case of single snapshot. This paper only considers one-dimensional DOA estimation, so another future direction of interest would be the extension of SPICE-ML to two-dimensional DOA estimation.

Author Contributions

L.L. and Y.Z. proposed the original idea. Y.Z. and X.Y. designed the simulations and provided all of the figures and data. Y.Z. wrote the paper under the guidance of L.L., and X.Y. checked the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Natural Science Foundation of China under Grant 61571146, and in part by the Fundamental Research Funds for the Central Universities under Grant HEUCFP201908.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gorcin, A.; Arslan, H. A two-antenna single RF front-end DOA estimation system for wireless communications signals. IEEE Trans. Antennas Propag. 2014, 62, 5321–5333. [Google Scholar] [CrossRef]
  2. Synnevag, J.F.; Austeng, A.; Holm, S. Adaptive beamforming applied to medical ultrasound imaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2007, 54, 1606–1613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yao, B.; Wang, W.; Yin, Q. DOD and DOA estimation in bistatic non-uniform multiple-input multiple-output radar systems. IEEE Commun. Lett. 2012, 16, 1796–1799. [Google Scholar] [CrossRef]
  4. Strumberger, I.; Minovic, M.; Tuba, M.; Bacanin, N. Performance of elephant herding optimization and tree growth algorithm adapted for node localization in wireless sensor networks. Sensors 2019, 19, 2515. [Google Scholar] [CrossRef] [Green Version]
  5. Paul, A.K.; Sato, T. Localization in wireless sensor networks: A survey on algorithms, measurement techniques, applications and challenges. J. Sens. Actua. Netw. 2017, 6, 24. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, S.C.; Li, F.G. A novel positioning system based on coverage area pruning in wireless sensor networks. Sensors 2018, 18, 4469. [Google Scholar] [CrossRef] [Green Version]
  7. Pal, P.; Vaidyanathan, P.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] [Green Version]
  8. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  9. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-Grid direction-of-arrival estimation using coprime array interpolation. IEEE Signal Process Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
  10. Stoica, P. SPECTRAL ANALYSIS OF SIGNALS: The missing data case. Synth. Lect. Signal Process. 2016, 1, 90–108. [Google Scholar]
  11. Liu, H.; Liao, G.; Zhang, J. A robust adaptive Capon beamforming. Signal Process. 2006, 86, 2820–2826. [Google Scholar]
  12. Schmidt, R.; Schmidt, R.O. Multiple emitter location and signal parameters estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  13. 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]
  14. Yardibi, T.; Li, J.; Stoica, P. Source localization and sensing: A nonparametric iterative adaptive approach based onweighted least squares. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 425–443. [Google Scholar] [CrossRef]
  15. Stoica, P.; Babu, P.; Li, J. SPICE: A sparse covariance-based estimation method for array processing. IEEE Trans. Signal Process. 2011, 59, 629–638. [Google Scholar] [CrossRef]
  16. Stoica, P.; Babu, P.; Jian, L. New method of sparse parameter estimation in separable models and its use for spectral analysis of irregularly sampled data. IEEE Trans. Signal Process. 2010, 59, 35–47. [Google Scholar] [CrossRef]
  17. Stoica, P.; Babu, P. SPICE and LIKES: Two hyperparameter-free methods for sparse-parameter estimation. Signal Process. 2012, 92, 1580–1590. [Google Scholar] [CrossRef]
  18. Cai, S.; Shi, X.; Zhu, H. Direction-of-arrival estimation and tracking based on a aequential implementation of C-SPICE with an Off-Grid model. Sensors 2017, 17, 2718. [Google Scholar] [CrossRef] [Green Version]
  19. Abeida, H.; Delmas, J.P. Efficiency of subspace-based DOA estimators. Signal Process. 2007, 87, 2075–2084. [Google Scholar] [CrossRef] [Green Version]
  20. Delmas, J.P. Asymptotically minimum variance second-order estimation for noncircular signals with application to DOA estimation. IEEE Trans. Signal Process. 2004, 52, 1235–1241. [Google Scholar] [CrossRef]
  21. Yurii, N. Lectures on convex optimization. IEEE Trans. Autom. Control 2006, 51, 1859. [Google Scholar]
  22. Yu, Y. Monotonic convergence of a general algorithm for computing optimal designs. Ann. Stat. 2010, 38, 1593–1606. [Google Scholar] [CrossRef]
  23. Stoica, P.; Zachariah, D.; Li, J. Weighted SPICE: A unifying approach for hyperparameter-free sparse estimation. Digital Signal Process. 2014, 33, 1–12. [Google Scholar] [CrossRef] [Green Version]
  24. 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]
  25. Nelder, J.A.; Mead, R. A simplex method for function minimization. Commput J. 1965, 7, 308–313. [Google Scholar] [CrossRef]
Figure 1. Signal model.
Figure 1. Signal model.
Sensors 20 00119 g001
Figure 2. Root mean square error (RMSE) with independent signal.
Figure 2. Root mean square error (RMSE) with independent signal.
Sensors 20 00119 g002
Figure 3. RMSE with coherent signal.
Figure 3. RMSE with coherent signal.
Sensors 20 00119 g003
Figure 4. RMSE with independent signal.
Figure 4. RMSE with independent signal.
Sensors 20 00119 g004
Figure 5. RMSE with coherent signal.
Figure 5. RMSE with coherent signal.
Sensors 20 00119 g005
Table 1. The SPICE-ML algorithm.
Table 1. The SPICE-ML algorithm.
step1. Initialization: p ^ k 0 k = 1 K , σ ^ 0 and θ ^ k 0 k = 1 N are obtained from SPICE algorithm.
repeat
step2. Compute R i , Q k i .
step3. Update p k and σ according to Formula (21) and (22).
step4. Minimizing the ML cost function with respect to θ k and obtain the refined estimated DOA θ ^ k k = 1 N .
end for
until (convergence)
Table 2. RMSE vs. signal to noise ratio (SNR) with independent signal.
Table 2. RMSE vs. signal to noise ratio (SNR) with independent signal.
SNR02468101214161820
DAS0.3390.3150.2850.2740.2710.2620.2580.2600.2530.2580.249
IAA0.3650.2970.2500.2260.1980.1800.1700.1610.1560.1530.149
SPICE0.3300.2630.2230.1920.1770.1660.1580.1520.1490.1480.147
SPICE-ML0.2380.1870.1440.1120.0880.0700.0570.0440.0350.0280.022
Table 3. RMSE vs. snapshots with independent signal.
Table 3. RMSE vs. snapshots with independent signal.
snapshots2030405060708090100
DAS0.3420.2570.2410.2410.2090.1960.1840.1800.178
IAA0.1480.1280.1250.1120.1100.1050.1000.1020.102
SPICE0.1500.1050.0950.0810.0780.0730.0680.0630.065
SPICE-ML0.1130.0890.0820.0690.0640.0610.0550.0510.051
Table 4. RMSE vs. SNR with coherent signal.
Table 4. RMSE vs. SNR with coherent signal.
SNR02468101214161820
DAS0.6250.6170.5930.5890.5890.5900.5880.5820.5790.5810.579
IAA0.3760.3220.2730.2330.2100.1900.1790.1660.1590.1550.150
SPICE2.7770.8110.7770.4070.3600.3120.2750.2520.2260.1990.195
SPICE-ML0.2440.1870.1470.1100.0910.0690.0550.0450.0340.0280.023
Table 5. RMSE vs. snapshots with coherent signal.
Table 5. RMSE vs. snapshots with coherent signal.
snapshots2030405060708090100
DAS0.6370.6260.6170.6130.6170.6060.6030.6090.609
IAA0.1630.1430.1350.1300.1240.1240.1170.1120.117
SPICE1.9500.3640.3210.3110.2970.2850.2850.2920.276
SPICE-ML0.1140.0940.0780.0740.0650.0620.0560.0520.053

Share and Cite

MDPI and ACS Style

Zheng, Y.; Liu, L.; Yang, X. SPICE-ML Algorithm for Direction-of-Arrival Estimation. Sensors 2020, 20, 119. https://doi.org/10.3390/s20010119

AMA Style

Zheng Y, Liu L, Yang X. SPICE-ML Algorithm for Direction-of-Arrival Estimation. Sensors. 2020; 20(1):119. https://doi.org/10.3390/s20010119

Chicago/Turabian Style

Zheng, Yu, Lutao Liu, and Xudong Yang. 2020. "SPICE-ML Algorithm for Direction-of-Arrival Estimation" Sensors 20, no. 1: 119. https://doi.org/10.3390/s20010119

APA Style

Zheng, Y., Liu, L., & Yang, X. (2020). SPICE-ML Algorithm for Direction-of-Arrival Estimation. Sensors, 20(1), 119. https://doi.org/10.3390/s20010119

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