Abstract

A novel algorithm is proposed for two-dimensional direction of arrival (2D-DOA) estimation with uniform rectangular array using reduced-dimension propagator method (RD-PM). The proposed algorithm requires no eigenvalue decomposition of the covariance matrix of the receive data and simplifies two-dimensional global searching in two-dimensional PM (2D-PM) to one-dimensional local searching. The complexity of the proposed algorithm is much lower than that of 2D-PM. The angle estimation performance of the proposed algorithm is better than that of estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithm and conventional PM algorithms, also very close to 2D-PM. The angle estimation error and Cramér-Rao bound (CRB) are derived in this paper. Furthermore, the proposed algorithm can achieve automatically paired 2D-DOA estimation. The simulation results verify the effectiveness of the algorithm.

1. Introduction

Direction-of-arrival (DOA) estimation is a fundamental problem in array signal processing and has been widely used in many fields [15]. Till now, many algorithms have been proposed for DOA estimation with uniform linear array. Among them, multiple signal classification (MUSIC) algorithm [6] and estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm [7] are widely used superresolution methods. The two-dimensional DOA (2D-DOA) estimation problem, which plays an important role in array signal processing, has received more attention. This problem is usually considered with rectangular array, two parallel uniform linear arrays, L-shape array, and so forth. Also, many algorithms have been considered for 2D-DOA estimation, which include 2D-MUSIC algorithm [8], 2D Unitary ESPRIT algorithm [9], modified 2D-ESPRIT algorithm [10], matrix pencil methods [11, 12], maximum likelihood method [13], parallel factor (PARAFAC) algorithm [14], and high order cumulant method [15].

Propagator method, which is known as a low complexity method without eigenvalue decomposition (EVD) of the covariance matrix of the received data, has been proposed for DOA estimation through peak searching [1720]. Due to its low complexity, PM algorithms are widely used for 2D-DOA estimation. In [21], a PM-based DOA estimation algorithm is proposed for two parallel uniform linear arrays via the rotational invariance property of propagator matrix, which requires extra pairing match. Reference [22] presents an efficient 2D-DOA estimation algorithm with two parallel uniform linear arrays using PM. In [23], an improved PM algorithm is proposed for 2D-DOA estimation, which has better angle estimation performance than the algorithms in [21, 22]. The above-mentioned PM algorithms can be extended to the rectangular array for 2D-DOA estimation; however, they only employ the rotational invariance property of propagator matrices, and as a result of that, the accuracy of angle estimation performance is not sufficient in some cases, especially with low signal to noise ratio (SNR). Two-dimensional PM (2D-PM) algorithm through peak searching can be extended for 2D-DOA estimation; however, the high computational complexity caused by 2D peak searching makes it inefficient.

In this paper, we derive a reduced-dimension PM (RD-PM) algorithm, which reduces the high complexity for 2D-DOA estimation with uniform rectangular array compared with 2D-PM algorithm. The proposed algorithm applies the rotational invariance property of propagator matrix to get the initial angle estimation and then employs one-dimensional local searching to get more accurate angle and finally obtains the other angle via least square (LS) method and estimate azimuth and elevation angles. The proposed algorithm has the following advantages: the proposed algorithm has lower computational complexity than 2D-PM algorithm since it requires no EVD of the covariance matrix of the receive data and it only requires one-dimensional local searching; the angle estimation performance of the proposed algorithm is better than that of ESPRIT algorithm and the PM algorithms in [21, 23], also very close to that of 2D-PM algorithm; it can obtain automatically paired parameters estimation.

The remainder of this paper is structured as follows: Section 2 shows the data model for uniform rectangular array, while Section 3 proposes the RD-PM algorithm for 2D-DOA estimation. The angle estimation error is derived as well as CRB in Section 4. In Section 5, the simulation results verify the feasibility and effectiveness of the proposed algorithm, and the conclusions are showed in Section 6.

Notion. , , , and denote transpose, conjugate-transpose, inverse, and pseudoinverse operations, respectively; stands for diagonal matrix whose diagonal element is a vector ; , , and denote a vector of ones, a identity matrix, and a vector of zeros, respectively. , , and are the Kronecker product, the Khatri-Rao product, and the Hadamard product, respectively; is to get the phase angle.

2. Data Model

As illustrated in Figure 1, consider a uniform rectangular array having sensors; and are the numbers of sensors in -axis and -axis, respectively. The distance between the two adjacent elements is in both -axis and -axis. Assume that there are uncorrelated sources, and are the elevation and azimuth angles of the th source.

The received signal of the first subarray can be expressed aswhere ,  , , and is the wavelength. is the additive white Gaussian noise of the first subarray and is the source vector.

Similarly, the received signal of the th subarray in the rectangular array is also denoted as where and is the additive white Gaussian noise of the th subarray. So the integrated received signal of the rectangular array is shown as follows:where stands for the received signal of the th subarray. Equation (3) is also denoted aswhere , , , and .

For the signal model in (4), we collect snapshots to get the covariance matrix as

3. Two-Dimensional DOA Estimation Algorithm

3.1. -PM Algorithm

In this section, we give a brief introduction of 2D-PM algorithm.

Partition the matrix as where and . Since is a nonsingular matrix, is a linear transformation of as [17]where is the propagator matrix. We define the matrix as According to (6) and (7), we have Then, we construct the function of the 2D-PM spatial spectrum as By searching and , the largest peaks of the costing function in (10) correspond to the angles information. However, the high computational cost caused by the 2D peak searching makes this algorithm inefficient.

For estimating the propagator matrix , we partition the covariance matrix aswhere . The propagator matrix meets the minimization problem , and the estimation of is

3.2. Reduced-Dimension Propagator Method for 2D-DOA Estimation
3.2.1. Initial Estimation

Define Then, and can be rewritten asConstruct matrix as Divide into blocks as , and every block contains elements. Construct and ; then we havewhere . Use EVD of and obtain the th eigenvalue as ; then the initial estimation of is

3.2.2. 2D-DOA Estimation

Define Equation (18) can also be denoted aswhere . Equation (19) is the problem of quadratic optimization. We consider the constraint , where . The optimization problem can be reconstructed via linearly constrained minimum variance as [24] The costing function is where is a constant. We have According to (22), we obtain , where is also a constant. As , . Then, can be expressed as Combining (20) and (23), the estimation of is Equation (24) also can be present as We find the largest peaks of the element of through searching , where is a small value. The largest peaks correspond to the estimations of ; then we obtain vectors of via (23). Since , we havewhere . Normalize and use LS method to estimate ; the LS fitting is where , , , and is the estimation of . The solution to is Finally, we get the estimation of as Till now, we have achieved the proposed algorithm for 2D-DOA estimation. The major steps and complexity of the algorithm are shown as follows:(1)construct the covariance matrix via (5);, ;(2)partition the to get ; then construct and obtain the initial estimation of from the EVD of ;, ;(3)find the largest peaks of the element of through searching and obtain the estimation of ;, ;(4)compute via (23), and then can be estimated through the LS fitting; ;(5)estimate via (29).

Remark 1. In this paper, we assume that the target number is a prior known parameter. If not, we can estimate by applying the method in [25, 26].

Remark 2. The proposed algorithm gets the estimation of through searching and then obtains the corresponding and estimates corresponding . Since and are correspondence, the proposed algorithm can obtain automatically paired parameters estimation.

Remark 3. The proposed algorithm requires being different. If the sources have the same or , then the direction matrix is a matrix of rank loss according to the property of Khatri-Rao product [27]. That will lead many algorithms invalid, such as ESPRIT, PMs in [21, 23], and our algorithm.

3.3. Complexity Analysis

The proposed algorithm has much lower complexity than 2D-PM algorithm. The major complexity of the proposed algorithm is , where is the number of steps within the local searching range; 2D-PM algorithm requires , where is the number of steps within the global searching range and , while ESPRIT algorithm costs .

The complexity comparison with different parameters is shown in Figures 2 and 3. From Figures 2 and 3, we can find that the proposed algorithm has much lower complexity than 2D-PM algorithm since the proposed algorithm only requires a one-dimensional local searching.

4. Performance Analysis

This section aims at analyzing the theoretic estimation error of the proposed algorithm and deriving the Cramér-Rao bound (CRB).

4.1. Error Analysis

Assume that the receive data in (1) can be written as where is the noisy data matrix and is the additive perturbation matrix with mean zero and variance statistically independent. Since , we havewhere and . According to (11), we also have The estimation of propagator matrix can be denoted as With a first-order expansion in and , is rewritten as where . We also obtainwhere with .

Since , using first-order expansion we obtain The proposed algorithm has the same cost function as 2D-PM algorithm, which is a reduced-dimension form of 2D-PM algorithm. For simplification, we derive the theoretical error of the proposed algorithm through the cost function of 2D-PM algorithm.

The function from which we search for the minima in order to determine the estimates is Define , , , and . Use to denote the true angles corresponding to the th target and denotes the estimated ones. Then, After the Taylor expansion around (use instead of for simplification),where is a vector of some value on the line segment joining and .

Define Under large sample, is asymptotically zero-mean Gaussian distributed [1719], so the associated asymptotic covariance matrix is given byAccording to (36) and (41) and taking into account that , we can determinewhere , , and .

That is, , , and .

According to (29), the mean square error (MSE) of azimuth and elevation can be derived asThe comparison of theoretical estimation error and simulation results will be shown in Section 5.

4.2. CRB

According to [28], we can derive the CRB for uniform rectangular array where , and is the th column of ; , , , and is the covariance of the noise. is the Hadamard product.

5. Simulation and Analysis

5.1. Simulation Results

We present 1000 Monte Carlo simulations to assess the angle estimation performance of the RD-PM algorithm. Define root mean square error (RMSE):where are the estimates of in the th Monte Carlo trial, respectively. We assume that there are noncoherent sources with , , and ; the spacing between two adjacent elements is . Note that and are the numbers of sensors in -axis and -axis, respectively. is the number of snapshots.

Simulation 1. Figure 4 depicts angle estimation performance of the proposed algorithm with , , , and  dB. The results show that elevation and azimuth angles can be clearly observed.

Simulation 2. Figures 5 and 6 present angle estimation performance comparison among the proposed algorithm, 2D-PM algorithm, the PM algorithm in [21], the PM algorithm in [23], ESPRIT algorithm, RD-MUSIC algorithm, CRB, and the theoretic estimation error of RD-PM. We can see that the estimation performance of the proposed algorithm is better than ESPRIT algorithm and the PM algorithms in [21, 23] and is very close to 2D-PM algorithm and RD-MUSIC algorithm. Compared with 2D-PM algorithm, the proposed algorithm has the same costing function and simplifies two-dimensional global searching to one-dimensional local searching. As a result of that, our algorithm has very close angle estimation performance to 2D-PM algorithm and greatly reduces complexity. Furthermore, Figures 5 and 6 also show that the RMSE of the proposed algorithm is almost the same as the theoretic estimation error.

Simulation 3. Figure 7 presents angle estimation performance of the proposed algorithm with , , , and different values of . We can find that the angle estimation performance is improved when increases. When increases, we get more samples to estimate more accurate propagator matrix , and the angle estimation performance is enhanced.

Simulation 4. Figure 8 shows angle estimation performance of the proposed algorithm with , , , and different values of . It is shown that the angle estimation performance of the proposed algorithm is improved when decreases. More sources enhance the interference, so the angle estimation performance of the proposed algorithm is affected.

Simulation 5. Figures 9 and 10 investigate angle estimation performance of the proposed algorithm with , , and different values of or ,respectively. It is indicated that the performance of the proposed algorithm is improved when or increases. Multiple sensors enhance the aperture of the array as well as diversity gain, so the angle estimation performance is improved.

Simulation 6. Figure 11 depicts angle estimation performance of the proposed algorithm when the two sources are closely spaced with , , , and SNR 20 dB. It is shown in Figure 11 that our algorithm works well for the closely spaced sources.

6. Conclusions

In this paper, we have presented a novel algorithm for 2D-DOA estimation for uniform rectangular array using RD-PM algorithm. The proposed algorithm, which only requires one-dimensional local searching and avoids the EVD of the covariance matrix of the receive data, has a much lower complexity than 2D-PM algorithm. Simulation results show that the angle estimation performance of the proposed algorithm is better than that of ESPRIT algorithm and conventional PM algorithms, also very close to that of 2D-PM algorithm. Furthermore, the proposed algorithm can achieve paired 2D-DOA estimation automatically.

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 China NSF Grants (61371169, 61301108,61471191,61471192, and 61271327), Jiangsu Planned Projects for Postdoctoral Research Funds (1201039C), China Postdoctoral Science Foundation (2012M521099, 2013M541661), Open Project of Key Laboratory of Modern Acoustic of Ministry of Education (Nanjing University), the Aeronautical Science Foundation of China (20120152001), Qing Lan Project, priority academic program development of Jiangsu high education institutions, and the Fundamental Research Funds for the Central Universities (NZ2012010, NS2013024, kfjj130114, kfjj130115).