Direction-of-arrival estimation based on joint sparsity.

We present a DOA estimation algorithm, called Joint-Sparse DOA to address the problem of direction-of-arrival (DOA) estimation using sensor arrays. Firstly, DOA estimation is cast as the joint-sparse recovery problem. Then, norm is approximated by an arctan function to represent joint sparsity and DOA estimation can be obtained by minimizing the approximate norm. Finally, the minimization problem is solved by a quasi-Newton method to estimate DOA. Simulation results show that our algorithm has some advantages over most existing methods: it needs a small number of snapshots to estimate DOA, while the number of sources need not be known a priori. Besides, it improves the resolution, and it can also handle the coherent sources well.

Introduction

Direction-of-Arrival (DOA) estimation using sensor arrays has been an active research area, playing a fundamental role in many applications involving electromagnetic, acoustic, and communication systems [1]. Many classical algorithms are available, and the popular methods include beamforming [2], MUSIC [3], ESPRIT [4] and the maximum likelihood method [5], etc. The beamforming method has low angle resolution and suffers from the Rayleigh resolution limit. MUSIC, ESPRIT and the maximum likelihood method all rely on the statistical properties of the data, and thus, require a sufficiently large number of samples for accurate estimation. Besides, MUSIC and ESPRIT cannot handle strongly coherent sources, while the maximum likelihood method has high computation costs.

The problem of sparse recovery has evolved rapidly recently [6,7] and it has been applied in DOA estimation with array processing. Gorodnitsky et al. [8] used a weighted least-squares algorithm named FOCUSS for DOA estimation, but this algorithm can only be used for single snapshots. Cotter [9] combined multiple measurement vectors (MMV) and matching pursuit (MP) to solve the joint-sparse recovery problem in DOA estimation, but it has low angle resolution. JLZA-DOA is proposed in [10]; it minimizes a mixed L2,0 norm to deal with the joint-sparse recovery problem, and a fixed point method is used for DOA estimation. This algorithm doesn’t satisfy numerical stability, as matrix inversion is inevitable in every iteration. Stoica et al. [11] presented a novel SParse Iterative Covariance-based Estimation approach, abbreviated SPICE. However, this algorithm needs more snapshots to estimate DOA. Wide-band covariance matrix sparse representation (W-CMSR) is proposed in [12] for DOA estimation of wideband signals. So far, the most successful joint-sparse recovery algorithm for DOA estimation is L1-SVD [13,14]. It combines the SVD step of the subspace algorithms with a sparse recovery method based on l2,1 –norm minimization. However, the number of sources needs be known a priori.

In this paper, we present Joint-Sparse DOA estimation, abbreviated as JSDOA, for sensor array DOA estimation. First, DOA estimation is cast as a joint-sparse recovery problem. Then, L2,0 norm is approximated by the arctan function to represent spatial sparsity and DOA estimation can be obtained by minimizing the approximate L2,0 norm. Finally, the minimization problem is solved by a quasi-Newton method to estimate DOA. The proposed algorithm has some advantages over most existing methods: it needs a small number of snapshots to estimate DOA, an the number of sources need not be known a priori. Besides, it improves the probability of resolution, and it can also handle coherent sources well.

The outline of the paper is as follow. In Section 2, the DOA estimation problem is formulated. The new algorithm, called JSDOA, is proposed in Section 3. In section 4, the validity of the proposed algorithm is proved by a number of simulations. Finally, conclusions are presented in Section 5.

Problem Formulation

DOA Estimation Problem

Consider a linear array consisting of M identical sensors and receiving signals from K narrowband signals s1(t), s2(t), ⋯, s(t), which arrive at the array from directions θ̄1, θ̄2, ⋯ θ̄ with respect to the line of array. The received signal y(t) at the m sensor can be written as: where denote steering vectors, n(t) (m = 1, 2 ⋯, M) stand for the additive noise.

Let y(t) = [y1(t), y2(t), ⋯, y(t)], n(t) = [n1(t), ⋯, n(t)], Equation (1) can be written as: where the manifold matrix A (θ̄) consists of the steering vectors :

Therefore, DOA estimation is to find K and θ̄ from T snapshots .

Joint-Sparse Recovery for DOA Estimation Problem

Because sources are sparse in space, DOA estimation with sensor arrays can be expressed as a joint-sparse recovery problem. Let Ω denote the set of possible locations, denotes a grid that covers Ω. We assume that the grid is fine enough such that the true location parameters of the existing sources lie on the grid. Let:

Then the received signal y(t) at the m sensor can be written as: where the n element x(t) of x(t) is nonzero only if θ = θ̄, k ∈ [1,2, ⋯, K] and in that case x(t) = s(t):

Let Φ = [a(θ1), a(θ2), ⋯, a(θ)], (3) can be expressed as: when the number of snapshots is denoted as T, Equation (4) can be written as: where Y = [y(t1), y(t2), ⋯, y(t)], X = [x(t1), x(t2), ⋯, x(t)]. As we know, X is row-sparse and only K rows have nonzero elements, so DOA estimation can be obtained by a joint-sparse recovery problem, which is also called a multiple measurement vectors (MMV) problem. Using l norm to express joint sparsity, the MMV problem can be converted to l norm minimization: where p and q are non-negative, and ‖X‖ is defined as:

Concretely, Eldar and Mishali [15] use p = 2, q = 1. Chen and Huo [16] study for p = 1, q = 1. The above algorithms do not perform very well for ‖X‖ 2,1 and ‖X‖ 1,1 can’t sufficiently reflect joint sparsity. Considering ‖X‖ 2,0 can reflect joint sparity sufficiently, we minimize ‖X‖ 2,0 norm to solve the MMV problem. However, ‖X‖ 2,0 norm minimization can hardly be solved directly. Therefore, in this paper we approximate ‖X‖ 2,0 norm by an arctan function and estimate DOA by solving an approximate ‖X‖ 2,0 norm minimization problem.

Joint-Sparse DOA Estimation Algorithm

In this section, the DOA estimation problem is converted into an approximate ‖X‖ 2,0 norm minimization problem. Firstly, the L2,0 norm is approximated by an arctan function to construct an approximate ‖X‖ 2,0 norm minimization problem. Then this problem is solved by a quasi-Newton method to estimate DOA.

Basic Idea of the Proposed Method

Let ξ = ‖X(i,:)‖2, i = 1, 2, ⋯, M, we have: where ξ = (ξ1, ⋯, ξ). Considering the following arctan function: where δ is a positive parameter and s is a variable parameter. Then f(s) has the following property:

Let

From Equation (10), we have:

From Equations (8) and (10), then:

So DOA can be obtained by solving the approximate ‖X‖2,0 norm minimization: where δ is a small positive constant. Because noise is unknown, we synchronously hope to minimize F(X) and . Then Equation (12) is converted into a multiple objective optimization:

Using the linear weighting method, Equation (13) can be written as: where ‖•‖ denotes the Frobenious norm and the parameter λ will be discussed in Section 3.2. When the parameter δ is very small, the objective function F(X) is highly unsmoothed and contains a lot of local minimization, so its global minimization is not easy. On the other hand, if the parameter δ is larger, the objective function F(X) is smoother and contains less local minimization. In order to obtain global minimization, we select a decreasing sequence for δ, denoted as: where δ is a relatively large value, and δ is a small value. For δ = δ, the solution of Equation (14) is denoted as x, and x is used as the initial value for δ = δ, thus we hope the proposed algorithm can escape from getting trapped into a local minimum and reach the global minimization for a small value δ = δ.

For some fixed value δ = δ, minimization problem Equation (14) is solved by the quasi-Newton method in this paper. One of the most successful quasi-Newton methods is the BFGS algorithm, which is second order convergent and has good numerical stability. Therefore, we use the BFGS algorithm to solve Equation (14) for some fixed value δ = δ.

The conjugate gradient for a matrix variable is defined as: where X and X denote the real part and the imaginary part respectively, X*denotes the conjugate of X. Then the conjugate gradient of L(X) (for its derivation see Appendix I) can be expressed as: where

Solving by the BFGS algorithm, the main iterative steps are as follows: where s = X – X,

Search step length t, satisfy:

Iteration:

BFGS adjustment:

Algorithm Description

Based on the above idea, JSDOA can be described as shown in Table 1.

Table 1.

The main steps of JSDOA.

Algorithm 1: joint-sparse DOA estimation
Input: A, y
Initialization:
(1) Set X0 = AT(AAT)−1Y
(2) Select a decreasing sequence δ = [δ1δ2 ⋯ δJ], set ɛ and parameter λ
Iteration:
(1) for j = 1,2, … J
(2) solving (14) by BFGS algorithm
(2–1) x = xj − 1, H = I (I is unit matrix);
(2–2) while norm (∂Lδj,λ (X)/∂X) > δjɛ;
(2–3) Search step length t˜:t˜=mintL_δ_j,λ(X−tH∂L_δ_j,λ(X)/∂X)
(2–4) Let X = X – t̃H∂Lδj,λ (X)/∂X
(2–5) Update the matrix by (16)
(2–6) end
(3) Let Xj = X
(4) end
Output: X̂ = XJ, spatial spectrum P(θ̄i) = 10 log10 (‖X̂(i,:)‖2)

Remark:

(1) Select parameter λ: in this paper, we select the parameter λ by the α– method [17]. Set f1(X) = F(X), , and minimizing f1(X) and f2(X), we have;

It is easy to know that X((i = 1,2) is 0 and A(AA)−1Y respectively. Then we have:

Set λ = λ2/λ1 and introduce auxiliary parameter α, we can have:

Setting the coefficient matrix f(i,j) = f, from the above equations, the solution is: where e = [1,1], and λ = λ2/λ1.

(2) Select parameter δ: In this paper, we set δ = γδ–1, j = 2, ⋯, J, γ ∈ (0.5,1). Let , we hope parameter δ1 satisfies:

In order to save computation cost, we set . When δ → 0, F(x) → ‖x‖0. However, if δ is too small, F(x) will be sensitive to noise, so δ shouldn’t be too small.

(3) For coherent sources: Many of popular methods, such as MUSIC and ESPRIT require the assumption that sources are uncorrelated, because the source covariance matrix remains nonsingular so long as none of these sources are coherent. However, the algorithm proposed in this paper doesn’t need to compute the source covariance matrix, and joint sparsity is used to estimate DOA, so the proposed algorithm can handle highly coherent sources as well.

(4) Limitation: The algorithm proposed in this paper suffers from two problems. First, there exist no clear-cut guidelines for the selection of δ, especially when knowledge of noise is unknown. Second, the proposed algorithm can’t perform well when the SNR is relatively low.

Simulation Results

In this section, we present several numerical simulation results to illustrate the performance of the proposed algorithm. First, we compare the spatial spectrum of JSDOA to those of beamforming [2], MUSIC [3], and L1-SVD [13] under various snapshot scenartios. Next, we compare the spatial spectrum with more than two sources. Then, the probability of resolution is compared. Finally, we compare the spatial spectrum for coherent sources. In the following numerical simulations, we consider a uniform linear array consisting of M = 8 identical sensors and receiving signals are narrowband signals. The sensors are uniformly placed with a spacing of half a wavelength. The interval for the DOA is Ω = [−90,90], we use a uniform grid to cover Ω with a step of 1°, which means N = 181.

Spatial Spectrum Comparison under Various Snapshots

In this simulation, we compare spatial spectrum of different algorithms under various snapshots. We consider two uncorrelated sources located at 13° and 20° with SNR = 10. For JSDOA, we set , δ = γδ–1, γ = 0.5, J = 7. In Figures 1(a,b), we set T = 10 and T = 5, respectively. It can be seen from Figure 1 that JSDOA and L1-SVD can resolve the two sources as T = 10. However, only JSDOA can resolve the two sources as T = 5.

Figure 1.

Spatial spectrum comparison under various snapshots for two uncorrelated sources: (a) T = 10; and (b) T = 5.

Spatial Spectrum Comparison with More Than Two Sources

In this simulation, we compare the spatial spectrum of different algorithms with more than two uncorrelated sources. Set SNR = 10, T = 20, , δ = γδ–1, γ = 0.5, J = 7. In Figure 2(a), there are three sources located at 10, 20 and 30. In Figure 2(b), there are five sources located at 0, 10, 20, 30 and 40. It is shown from Figure 1 that JSDOA and L1-SVD can resolve the three sources, but only JSDOA can resolve the five sources.

Figure 2.

Spatial spectrum comparison for more than two uncorrelated sources: (a) Three sources; and (b) Five sources.

Probability of Resolution Comparison under Various Conditions

In this simulation, we compare probability of resolution for two uncorrelated sources. JSDOA, MUSIC and L1-SVD are used for comparison and we set T = 5, δ = γδ–1, γ = 0.5, J = 7. For each Δθ or SNR, 200 independent simulation are carried out. In Figure 3(a), one source is fixed at 10°, the other source is located at 10° + Δθ. vary from 0 to 30. In Figure 3(b), two sources are located at 10° and 20°. SNR vary from −5 to 20. It is shown from Figure 3 that JSDOA has the higher probability of resolution than MUSIC and L1-SVD.

Figure 3.

Probability of resolution comparison under various conditions. (a) Against Δθ; and (b) Against input SNR.

Spatial Spectrum Comparison for Coherent Sources

In this simulation, we compare spatial spectrum for completely coherent sources. Set SNR = 10, T = 20, , δ = γδ–1, γ = 0.5, J = 7. In Figure 4(a) two coherent sources, located at 10° and 20°, are considered.

Figure 4.

Spatial spectrum comparison for coherent sources. (a) Two coherent sources; and (b) Three coherent sources

In Figure 4(b) three coherent sources, located at 10° and 20°, are considered. It is shown from Figure 4 that JSDOA and L1-SVD can be used to estimate DOA for coherent sources and that JSDOA has higher resolution than L1-SVD.

Conclusions

In this paper, a joint-sparse recovery algorithm, called JSDOA, is proposed to estimate DOA with sensor arrays. We pose the DOA estimation problem as a joint-sparse recovery problem, which can be solved by minimizing the approximate L2,0 norm. In particular, the L2,0 norm is approximated by an arctan function to represent spatial sparsity and the approximate L2,0 norm minimization problem is solved by a quasi-Newton method. Finally, the proposed algorithm is examined by simulations. Several advantages over existing DOA estimation algorithms were identified. It can perform well with a limited number of snapshots, while the number of sources need not be known a priori. Besides, it improves the resolution, and it can also handle the coherent sources well.