International Journal of Antennas and Propagation
Volume 2015 (2015), Article ID 681251, 6 pages
http://dx.doi.org/10.1155/2015/681251

Research Article

Utilizing Principal Singular Vectors for 2D DOA Estimation in Single Snapshot Case with Uniform Rectangular Array

Yuntao Wu,¹ Xiaobing Pei,² and Hing Cheung So³

¹Key Lab for Intelligent Robot in Hubei Province, Wuhan Institute of Technology, Wuhan 430073, China
²School of Software, Huazhong University of Science and Technology, Wuhan 430074, China
³Department of Electronic Engineering, City University of Hong Kong, Kowloon Tong, Hong Kong

Received 18 February 2015; Accepted 4 May 2015

Academic Editor: Atsushi Mase

Copyright © 2015 Yuntao Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The problem of azimuth and elevation directions of arrival (DOAs) estimation using a uniform rectangular array (URA) in single snapshot case is addressed in this paper. Using the principal singular vectors of the observed data matrix, an iterative procedure based on the linear prediction property, and weighted least squares is proposed for finding the DOAs with lower computational complexity. Furthermore, the azimuth and elevation parameters are automatically paired. Computer simulations are included to demonstrate the effectiveness of the proposed algorithm.

1. Introduction

Two-dimensional (2D) DOA estimation of multiple narrowband signals is an active research topic in array signal processing due to its wide applications in radar, sonar, radio astronomy, and mobile communications and so forth. In recent years, many high-resolution subspace-based algorithms such as 2D MUSIC- [1] or 2D ESPRIT-based [2, 3] as well as the modified unitary ESPRIT algorithms [4] have been proposed to jointly estimate the 2D parameters with various array geometries. These algorithms are effective in certain conditions; however, high-dimension data stacking and eigen-decomposition of a constructed covariance of array output are needed, corresponding to high computational load. To reduce the computational complexity of eigendecomposition, fast algorithms have been reported [5, 6]. In general, a direction-finding algorithm in real application must provide a real-time solution that is computationally efficient and uses only a few array snapshots. In some real scenarios, such as fast time-variant channel in mobile communications and sonar signal processing due to physical constraints, only one or a few snapshots are available. Therefore, the problem of parameter estimation in the case of small number of samples has been addressed in the literature [7–9], where the advanced concept of compressed sensing (CS) is applied. In particular, in the worst case such as automotive radar systems [10], only a single snapshot is available for parameter estimation of multiple spatial sources. That is to say, the problem of single snapshot DOA estimation is important in certain application and corresponding DOA estimation algorithms in single-snapshot case have been recently proposed [11–17].

In this paper, an algorithm of joint estimation of the azimuth and elevation angles using a uniform rectangular array (URA) is proposed for the case of single snapshot. The main idea is based on the principal-singular-vector utilization for model analysis (PUMA) [18, 19]. With the use of the PUMA technique, not only the estimation performance in single snapshot case is improved but the computational complexity of the proposed method is also largely reduced. Moreover, the 2D parameters are automatically paired. Simulation results show that the proposed method provides better performance but with largely reduced computational load compared with the existing 2D ESPRIT method.

The rest of the paper is organized as follows. The problem formulation is given in Section 1. In Section 2, the proposed 2D DOA estimator for single snapshot case with a URA array is devised. Simulation results are included in Section 3 to evaluate the performance of the proposed method by comparing with the 2D ESPRIT method as well as the corresponding Cramer-Rao lower bound (CRLB). Finally, conclusions are drawn in Section 4.

2. Problem Formulation

Figure 1 shows the planar URA configuration for receiving the incoming signals. The array is composed of omnidirectional antennas with interelement spacing and we assume that narrowband signal sources imping on this array from the 2D spatial directions , . The sensor located at is the reference sensor.

Figure 1: Uniform rectangular array geometry.

The single-snapshot output of the sensor at can be expressed as

where denotes the wavefront of the th signal source, denotes the wavelength of the impinging wavefront, and is the white Gaussian noise with power received at the th sensor.

Our task is to estimate the parameter pairs , , from the received array output.

3. Proposed Method

In this section, the PUMA algorithm [18] is exploited to estimate the 2D DOA parameters of multiple narrowband sources with URA in the case of single snapshot. Define the spatial frequencies as and ; we first express (1) in matrix form as

where , , and .

The has the form of

It is easy to find that can be factorized as

where

Here, denotes the transpose operator and is a diagonal matrix with vector as its main diagonal. It is seen from (8)-(9) that the elements of and satisfy the linear prediction (LP) relations of and , where stands for the th element of .

On the other hand, can be decomposed using singular value decomposition (SVD) as

where and contain the corresponding left and right singular vectors, respectively, is the diagonal matrix of singular values sorted in descending order, and stands for the conjugate transpose.

From the decomposition in (4)–(10), we have , , and thus the best rank- approximation of according to (10), denoted by , is

where

are the corresponding signal subspace components.

Since and have the same subspace, we have

where is an unknown nonsingular matrix with dimension .

For each , we have the following LP property:

where are the LP coefficients. The spatial frequencies are related to the following polynomial:

whose roots are , . According to (14), the LP error vector can be constructed by

where

Let be the symmetric weighting matrix. The weighted least squares (WLS) estimate of , denoted by , is

Defining a Toeplitz matrix , , where is a zero matrix with dimension and , the optimal is derived as [19]

where is the Kronecker product. As (19) depends on the unknown , we follow [18] to estimate the in an iterative manner and the estimation procedure is as follows:

(1)Set , where is the identity matrix with dimension .(2)Calculate using (18).(3)Compute an updated version of   using (19) with .(4)Repeat Steps (2)-(3) until a stopping criterion is reached.

Substituting in (15) and solving for the roots, denoted by , , we have the spatial frequency estimate :

where denotes the phase angle of . In order to obtain the 2D parameter pairing in an automatic manner, another estimation method for the parameters , , is presented as follows.

From (3)-(4), we have

where is the estimate of which is constructed using and

From (21), the least squares estimate of is

where denotes the pseudoinverse. Noting that the elements of satisfy the same LP property as in , we extract from to construct the equations:

where and , respectively. Following [18], the WLS estimate of is computed as

where the optimum weighting matrix has the form

where

Similar to the iterative estimation of in finding , we begin with in the iterations between (25) and (26) to obtain , . Finally, the spatial frequencies , , are estimated as

Note that and are automatically paired in the whole procedure.

Finally, the estimated azimuth angle and elevation angle are computed as

4. Simulation Results

Computer simulations have been conducted to evaluate the 2D DOA estimation performance of the proposed scheme in the presence of white Gaussian noise by comparing with the 2D ESPRIT method [3] and CRLB. The number of iterations in the proposed algorithm is and the initial parameter estimates are provided by the ESPRIT algorithm [3]. Note that the larger values for have been tried but no significant improvement is observed. The number of sensors in the URA is 180 where and . Assume that two narrowband signals imping on this received array from the directional angles and , respectively, while the additive noise is white Gaussian process. All results provided are averages of 200 independent runs.

Figures 2–5 plot the root mean square error (RMSE) performance of the 2D DOA estimates versus SNR. It is seen that the proposed method is superior to the 2D ESPRIT-based method at all signal-to-noise ratio (SNR) conditions. It is because the latter approach is based on the splitting the measurement space into signal subspace and noise subspace and generally gives a higher threshold SNR value. Nevertheless, all the RMSEs of the proposed method are 2 dB above the corresponding CRLBs at all SNRs. The average computational times of the proposed method and the 2D ESPRIT method are 0.0399 s and 0.2071 s, respectively, indicating that the former is much more computationally efficient.

Figure 2: Root mean square error for versus SNR.

Figure 3: Root mean square error for versus SNR.

Figure 4: Root mean square error for versus SNR.

Figure 5: Root mean square error for versus SNR.

5. Conclusion

A 2D DOA estimation algorithm with a URA in single snapshot case is proposed using the PUMA method. Compared with the 2D ESPRIT method, the proposed method can obtain higher estimation accuracy and its performance is close to the CRLBs at higher SNR. Moreover, the proposed method shows the advantage of lower computational complexity over the 2D ESPRIT method, and the 2D DOA parameters are automatically paired.

Conflict of Interests

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

Acknowledgments

The work described in this paper was jointly supported by a grant from the National Natural Science Foundation of China (Project no. 61172156), the Program for New Century Excellent Talents in University (NCET-13-0940), the Natural Science Foundation of Hubei Province (no. 2014CFB791), and the Research Plan Project of Hubei Provincial Department of Education (no. T201206).

References

Y.-M. Chen, J.-H. Lee, and C.-C. Yeh, “Two-dimensional angle-of-arrival estimation for uniform planar arrays with sensor position errors,” IEE Proceedings Part F: Radar and Signal Processing, vol. 140, no. 1, pp. 37–42, 1993.

R. Roy and T. Kailath, “ESPRIT-estimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989.

T. Kuroda, N. Kikuma, and N. Inagaki, “DOA estimation and pairing method in 2D-ESPRIT using triangular antenna array,” Electronics and Communications in Japan, Part I: Communications, vol. 86, no. 6, pp. 59–68, 2003.

M. D. Zoltowski, M. Haardt, and C. P. Mathews, “Closed-form 2-D angle estimation with rectangular arrays in element space or beamspace via unitary ESPRIT,” IEEE Transactions on Signal Processing, vol. 44, no. 2, pp. 316–328, 1996.

C.-H. Lin, W.-H. Fang, K.-H. Wu, and J.-D. Lin, “Fast algorithm for joint azimuth and elevation angles, and frequency estimation via hierarchical space-time decomposition,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07), vol. 2, pp. II-1061–II-1064, Honolulu, Hawaii, USA, April 2007.

Y. Wu, G. Liao, and H. C. So, “A fast algorithm for 2-D direction-of-arrival estimation,” Signal Processing, vol. 83, no. 8, pp. 1827–1831, 2003.

X. Li, X. C. Ma, S. F. Yan, and C. H. Hou, “Single snapshot DOA estimation by compressive sampling,” Applied Acoustics, vol. 74, no. 7, pp. 926–930, 2013.

M. Zhou, X. F. Zhang, X. F. Qiu, and C. H. Wang, “Two-dimensional direction of arrival (DOA) estimation for rectangular array via compressive sensing trilinear model,” International Journal of Antennas and Propagation. In press

W. Xiu-hong, M. Xing-peng, and Z. Nai-tong, “Single-snapshot DOA estimation based on compressed sensing in PCR systems,” in Proceedings of the IEEE Radar Conference, pp. 528–531, Cincinnati, Ohio, USA, May 2014.

C. J. Huang, C. W. Dai, T. Y. Tsai, W. H. Chung, and T. S. Lee, “A closed-form phasecomparison ML DOA estimation for automotive radar with one single snapshot,” IEICE Electron Express, vol. 10, no. 7, pp. 1–7, 2013.

Z. Hui, G. Lindong, and W. Yuexian, “Two-dimension direction finding using single-snapshot data,” in Proceedings of the International Conference on Communication Technology (ICCT '06), pp. 1–4, November 2006.

Q. S. Ren and A. J. Willis, “Extending MUSIC to single snapshot and on-line direction finding applications,” in Proceedings of the IET International Conference on Radar, pp. 783–787, October 1997.

A. C. Chang and C. C. Shen, “Orthogonal projection DOA estimation with a single snapshot,” IEICE Transactions on Communications, vol. E96-B, no. 5, pp. 1215–1217, 2013.

H. Chen and X. Zhang, “Two-dimensional DOA estimation of coherent sources for acoustic vector-sensor array using a single snapshot,” Wireless Personal Communications, vol. 72, no. 1, pp. 1–13, 2013.

B. M. Radich and K. M. Buckley, “Single-snapshot DOA estimation and source number detection,” IEEE Signal Processing Letters, vol. 4, no. 4, pp. 109–111, 1997.

P. Häcker and B. Yang, “Single snapshot DOA estimation,” Advances in Radio Science, vol. 8, pp. 251–256, 2010.

S. Fortunati, R. Grasso, F. Gini, M. S. Greco, and K. LePage, “Single-snapshot DOA estimation by using Compressed Sensing,” EURASIP Journal on Advances in Signal Processing, vol. 2014, article 120, 2014.

H. C. So, F. K. Chan, W. H. Lau, and C.-F. Chan, “An efficient approach for two-dimensional parameter estimation of a single-tone,” IEEE Transactions on Signal Processing, vol. 58, no. 4, pp. 1999–2009, 2010.

F. K. W. Chan, H. C. So, and W. Sun, “Subspace approach for two-dimensional parameter estimation of multiple damped sinusoids,” Signal Processing, vol. 92, no. 9, pp. 2172–2179, 2012.