Abstract

A novel efficient method for two-dimensional (2D) direction-of-arrivals (DOAs) estimation is proposed to reduce the computational complexity of conventional 2D multiple signal classification (2D-MUSIC) algorithm with uniform rectangular arrays (URAs). By introducing two electrical DOAs, the formula of 2D-MUSIC is transformed into a new one-dimensional (1D) quadratic optimal problem. This 1D quadratic optimal problem is further proved equivalent to finding the conditions of noise subspace rank deficiency (NSRD), which can be solved by an efficient 1D spectral search, leading to a novel NSRD-MUSIC estimator accordingly. Unlike 2D-MUSIC with exhaustive 2D search, the proposed technique requires only an efficient 1D one. Compared with the estimation of signal parameter via rotation invariance techniques (ESPRIT), NSRD-MUSIC has a significantly improved accuracy. Moreover, the new algorithm requires no pair matching. Numerical simulations are conducted to verify the efficiency of the new estimator.

1. Introduction

Estimation of direction-of-arrival (DOA) of multiple narrow-band signals using sensor arrays is of great interest in many fields such as radar, sonar, and wireless communications [1, 2]. Over several decades, this topic has been extensively addressed and numerical algorithms have been developed. Among those approaches, the multiple signal classification (MUSIC) [3] algorithm firstly exploits the orthogonality between the signal and the noise subspaces to achieve a so-called superresolution DOA estimate. Because those two subspaces are perfectly orthogonal to each other when the signal-to-noise ratio (SNR) is sufficiently high, the MUSIC algorithm can resolve two sources as closely spaced as possible in theory. Another outstanding advantage of the MUSIC algorithm over the other superresolution methods is that the former can be used with arbitrary array geometries [4]. However, the computational complexity of MUSIC is prohibitively expensive for real-time applications, owing to an involved matrix decomposition step for the signal or/and the noise subspace estimation and a tremendous spectral search step for final DOA estimates [5].

In practical applications, it is two-dimensional (2D) DOAs (2D DOAs) (i.e., elevation and azimuth angles) that people are usually most concerned about. Estimators such as MUSIC and estimation of signal parameter via rotation invariance technique (ESPRIT) [6] have been studied in one-dimensional (1D) situations and further extended to 2D scenarios [7]. Due to the increase in the dimensionality of the 2D DOA estimation problem, the complexity of DOA estimation process is severely affected by the array geometry [8], and further the pair matching (i.e., association or alignment) of the estimated elevation and azimuth angles is usually required [9]. Since, for arbitrary array configurations, the elevation and azimuth angles are extremely coupled with each other in the array response manifold, the heavy 2D spectral search cannot be avoided. To reduce the complexity, researchers are obliged to focus on reducing the computational burden of the 2D DOA estimation problem with various specified array structures.

Planar sensor arrays composed of two or more uniform linear arrays (ULAs) with simple geometry configurations have received considerable attention. Exploiting specified geometries such as L-shaped [10, 11], Y-shaped [12], and Z-shaped arrays [13], 2D DOAs of multiple incident signals can be estimated with reduced computational burden by applying most 1D subspace-based estimation methods including MUSIC and ESPRIT. Special array configurations allow utilizing the shift invariance and the partition of array response vector of two overlapping subarrays of each ULA such that the signal or/and noise subspaces can be estimated without eigenvalue decomposition (EVD) or singular valued decomposition (SVD). Nevertheless, spectral search is still required for those methods, and, as a matter of fact, the complexity of spectral search is substantially heavier than that of subspace decomposition, especially for 2D DOAs estimates [14]. Moreover, because the total arrays are usually divided into several subarrays in those methods, the array aperture is in fact reduced and DOA estimation accuracy may be sacrificed [15].

Uniform rectangular array (URA) is another frequently discussed array structure. A geometric formulation of the ESPRIT algorithm exploiting subarrays of URA is presented to reparameterize the weighted subspace fitting (WSF) algorithm in [16], leading naturally to an extension to the 2D DOA estimation problem. With high dimensional signal processing, the 2D subspace-based algorithms in [17] can provide precise estimates at the expense of high computational complexity exponentially increasing with the size of the rectangular arrays. To reduce the computational load, a tree structure algorithm by first performing 1D spatial smoothing and then 1D MUSIC several times successively is proposed in [18] to estimate the azimuth and elevation angles independently. Despite the computational efficiency, almost all of the aforementioned techniques realize 2D DOA estimates at the cost of a reduction in array aperture.

Recently, an outstanding reduced-dimensional MUSIC (RD-MUSIC) algorithm is proposed to avoid the high computational cost within 2D-MUSIC for direction of departure (DOD) and DOA estimation in multiinput multioutput (MIMO) radar [19]. Using DOA information embedded in the velocity sensors, RD-MUSIC starts 1D MUSIC searches for the DOD and DOA in succession without parameter pairing nor 2D search. In fact, the formula of RD-MUSIC is somewhat similar to that of 2D DOA estimates except that the noise projection matrix in the former is of reduced order and is invertible while that of the latter is not. Following this idea, we present in this paper a new efficient technique for 2D DOAs estimation with URAs. With two introduced electrical angles, we show that the formula of 2D-MUSIC can be transformed into a new 1D quadratic optimal problem. Further analysis demonstrates that this 1D quadratic optimal problem is equivalent to finding the conditions of noise subspace rank deficiency (NSRD), which can be solved by an efficient 1D spectral search. Consequently, a novel NSRD-MUSIC estimator at hand is derived. Unlike 2D-MUSIC with exhaustive 2D search, the proposed technique requires only an efficient 1D one. Compared with ESPRIT, NSRD-MUSIC has a significantly improved Root Mean Square Error (RMSE) performance.

The outline of this paper is as follows. The narrow-band signal model and the conventional 2D-MUSIC algorithm are introduced in Section 2. Two transformed electrical angles reformulating the 2D-MUSIC algorithm into a new 1D optimum problem is then considered, and the proposed method is addressed in detail in Section 3. The complexity of our method is analyzed in Section 4 and simulation results are conducted and discussed to validate the effectiveness of new method in Section 5.

2. Signal Model and Standard 2D-MUSIC

Assume that there are narrow-band signals with unknown 2D DOAs simultaneously incident on URA of sensors, where and are the numbers of elements in -direction and -direction, respectively. The corresponding element spacing symbols are and , as depicted in Figure 1. The azimuth angle is defined as the one between the wave direction and the -axis while the elevation angle is defined as the one between the -axis and the projection of wave direction onto the plane. Note that and . It is assumed that there are snapshots available; array output vector at snapshot is given bywhereis the matrix of the signal direction vectors andis the steering vector. In addition, and is the vector of source waveforms; is the vector of white sensor noise; is the center wavelength; is the matrix transpose.

The estimate of the array output covariance matrix (AOCM)is given bywhere is signal covariance matrix; is identity matrix; is the noise variance; is the mathematical expectation; and is the matrix Hermitian transpose.

2.1. Subspace Decomposition and 2D-MUSIC Algorithm

The EVDs of the theoretical AOCM (4) and the practical estimated AOCM (5) can be defined in a standard way:wherewith the subscripts and standing for the signal and noise subspace, respectively.

Based on the orthogonality between and , conventional 2D-MUSIC algorithm [3] suggests searching the following function:over the total 2D angular field-of-view with a fine grid. peaks of indicate source DOAs. One of the advantages of the 2D-MUSIC algorithm over other subspace-based methods is its nondependence on array configurations [4]. However, the complexity of this spectral search step is typically substantially high since, for each point, the product has to be computed.

3. The Proposed Algorithm

The 2D steering vector given in (3) represents the nominal array response. Following the notation in [20], can be rewritten aswhere denotes the Kronecker product and the elements of vectors and are given bywhere denote the th element of the embraced vector. Electrical angles and are related to the physical azimuth and elevation angles by and the corresponding inverse relations for obtaining azimuth and elevation angles are

Let us defineThen the spectral search involved in the 2D-MUSIC algorithm is equivalent to the following optimization problem:Observing that , (13) can be rewritten as whereTherefore, the optimization problem in (14) can be reconstructed as follows:

To solve the above quadratic optimization, it is suggested in [19] to use the following constraint: to eliminate the trivial solution , leading to a so-called RD-MUSIC estimator with the following search function: where is the th element of the embraced matrix. The electrical angle can be estimated by an efficient one-dimensional spectral search over for largest peaks of . However, the constraint in fact only warrants that the first element of equals one. We now propose a new NSRD-MUSIC method which has a similar computational complexity as compared to the RD-MUSIC algorithm.

Because , it is clear that holds if and only if matrix is of rank deficiency; that is,Noting that [21] for or equivalently for is a rank-deficient matrix if and only if is the electrical angle of source DOAs. Therefore, we haveBased on the above analysis, electrical angle that related the source DOAs, that is, , can be estimated by searchingor equivalently by searchingover for largest peaks of . Here, and stand for the minimum eigenvalue and the determinant of the embraced matrix, respectively. angles relating to largest peaks of should be the electrical angles .

With , we can compute steering vectors corresponding to as follows [19]:Combing (10) and (27) leads towhere stands for phase angles for each element of the embrace vector. We now exploit the least squares (LS) principle to estimate . DefineThen we have The LS solution for estimating with is given by Using (29) and (31), is estimated by

Now, we have estimated electrical angles and . Since, for each , we can compute a corresponding by using (27)–(32), no further pair matching is required for the proposed algorithm. Hence, we can directly insert , into (12) to finally obtain the estimates of the azimuth and elevation angles .

4. Complexity Analysis and Summary

Comparison of primary complexity among three algorithms including 2D-MUSIC [3], RD-MUSIC [19], and the proposed NSRD-MUSIC is shown in Table 1, where , , and stand for the search numbers along the directions of , , and , respectively. The common term gives the complexity of EVD on [21]. Because 2D-MUSIC involves a 2D spectral search step, it has to compute spectral points. Since RD-MUSIC and the proposed method involve a 1D search step, the two techniques only need to compute spectral points. Generally, we have . Therefore, NRSD-MUSIC has a much lower complexity than 2D-MUSIC.

Detailed steps for implementing the proposed NSRD-MUSIC algorithm are summarized in Algorithm 1.

Require: : snapshots of received datum.
(1) Initiation: .
(2) for  ; ;   do
(3) ;
(4) end for
(5) Compute .
(6) for each   do
(7) Compute ;
(8) Calculate by (25) or by (26);
(9) end for
(10) Find relating to the largest peaks of ;
(11) Insert into (27), (28) and compute by (32);
(12) Insert into (12) to obtain .
(13) return  : estimated azimuth and elevation angles.

5. Simulation Results

Simulations with 500 independent Monte Carlo methods have been conducted on UCR composed of sensors with half wavelength apart along both - and -axis to assess the performance of the proposed estimator. It is assumed that there are sources at , , , , and , that are implying on the array from far-field. The RMSEs for the estimates of source incident angles and are defined as respectively, where and are the true 2D DOAs while and represent the estimated values of the elevation and azimuth angles of the th trial.

Figures 2 and 3 plot the spectrums of three functions including 2D-MUSIC, RD-MUSIC [19], and the proposed estimator. It is seen from the figures that the proposed method estimates correctly, which generates corresponding peaks at , , and as expected. This indicates that can be estimated by the proposed method with 1D spectral search efficiently. Therefore, the proposed method has a much lower computational complexity as compared to 2D-MUSIC.

Figures 4 and 5 compare the estimation RMSEs of by four algorithms including 2D-MUSIC, 2D-ESPRIT, RD-MUSIC, and the proposed method. In Figure 4, the results are plotted against the SNR, where the number of snapshots is set as . In Figure 5, the results are given as functions of the number of snapshots, where SNR = 10 dB. Note that all the estimation accuracies in both Figures 4 and 5 are respected to the source at , . It can be concluded form Figure 4 that the proposed method provides a significantly improved accuracy performance as compared to 2D-ESPRIT over the total range of SNR = −9 dB to SNR = 15 dB. This is because the array aperture is in fact reduced in 2D-ESPRIT while that is not in the proposed method. It can be also concluded form Figure 4 that our method has a very close accuracy to RD-MUSIC and 2D-MUSIC, especially for . On the other hand, we can see clearly from Figure 4 that the proposed estimator provides an asymptotically similar accuracy performance to RD-MUSIC and 2D-MUSIC, which is much better than 2D-ESPRIT.

To further compare the performance of the new method with performances of 2D-MUSIC, 2D-ESPRIT, and RD-MUSIC, Figures 6 and 7 give the estimation RMSEs of by different algorithms. The results in Figure 4 are plotted against the SNR while those in Figure 5 are plotted as functions of the number of snapshots. In Figure 4, the number of snapshots is set as and in Figure 5 it is fixed as SNR = 10 dB. Similarly, all the results in Figures 6 and 7 are respected to the source at , . It can be concluded again form Figures 6 and 7 that our method significantly outperforms 2D-ESPRIT, which provides an asymptotically similar accuracy performance to 2D-MUSIC and RD-MUSIC.

Finally, we evaluate the performance of the proposed method against the correlation between the incident signals. We consider the case of three signals , , and , where and are correlated with each other, and the correlation coefficient is defined as [22] where and . The 2D DOAs of the three signals are set as , , , , and , . The RMSEs of DOA estimation versus source correlation with respect to the second signal are shown in Figures 8 and 9, where the number of snapshots is , , and .

It can be seen clearly form Figures 8 and 9 that the accuracies of the RMSEs of all the three algorithms increase as increases. This implies that the accuracies of all the three algorithms decrease when correlated source exists. However, the proposed method has a similar performance to the RD-MUSIC algorithm, and it performs closely to the standard 2D-MUSIC over the wide range to . Although the RMSEs of the new method are a little bigger than those of 2D-MUSIC, the complexity of the proposed method is much lower than that of 2D-MUSIC. Therefore, the new estimator makes an efficient trade-off between estimation accuracy and computational complexity.

6. Conclusion

We have proposed a novel efficient NSRD-MUSIC algorithm in this paper for 2D direction-of-arrivals (DOAs) estimation. The new method exploits the specific geometries of URAs to estimate 2D DOAs with efficient 1D spectral search, which costs a significantly reduced computational complexity as compared to 2D-MUSIC. Simulations indicate that NSRD-MUSIC which requires no pair matching has a significantly improved accuracy as compared to ESPRIT.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (61501142), Shandong Provincial Natural Science Foundation (ZR2014FQ003), China Postdoctoral Science Foundation (2015M571414), and the Fundamental Research Funds for the Central Universities (HIT.NSRIF.2016102).