International Journal of Antennas and Propagation

International Journal of Antennas and Propagation / 2017 / Article
Special Issue

Recent Advances in Array Antenna and Array Signal Processing for Radar

View this Special Issue

Research Article | Open Access

Volume 2017 |Article ID 4890203 | 5 pages | https://doi.org/10.1155/2017/4890203

Estimating Direction of Arrival by Using Two-Dimensional State-Space Balance Method

Academic Editor: Amerigo Capria
Received24 May 2017
Revised20 Sep 2017
Accepted04 Oct 2017
Published02 Nov 2017

Abstract

A study of a two-dimensional state-space balance (2D SSB) method for estimating direction of arrival (DOA) for uniform rectangular array (URA) is presented in this letter. The comprehensive utilization of controllability and observability matrices and automatic pairing technique are considered in this method by using the single snapshot. Therefore, the DOAs of elevation angle and azimuth angle can pair automatically and acquire better estimation performance compared with 2D matrix pencil method or unitary matrix pencil method. In addition, the proposed method can handle correlated signals directly without preprocessing. Simulation is conducted to verify the effectiveness of the proposed method.

1. Introduction

In array signal processing, target direction of arrival (DOA) estimation has been widely applied in radar, sonar, wireless indoor positioning [13] and so on. Since the uniform rectangular array (URA) can provide the information of azimuth and elevation angles, a lot of 2D DOA methods have been developed, such as multiple signal classification (MUSIC) algorithm [46], root-MUSIC [7], estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm [8, 9], and rank reduce (RARE) algorithm [10]. These conventional high resolution DOA parameter estimation techniques are mainly based on the statistical data by using the covariance matrix under the uncorrelated conditions. In the nonstationary environment, these DOA estimation methods performed poorly without preprocessing, but the matrix pencil (MP) method [11] got good performance by processing the received data directly based on the single snapshot. The unitary matrix pencil (UMP) method [12] can reduce the computation compared with the MP method by using a unitary matrix transformation, which can convert the complex operations to the real operations. The accuracy of both methods is nearly the same. The MP or UMP methods can also be regarded as a state-space method via observability or controllability processing. However, the method that using both observability and controllability matrices was not considered to estimate the DOAs in the spatial domain which may be acquire a good performance. On the other hand, the additional search of the angle pairing is essential to the MP/UMP method by maximizing the cost function, which exploits the orthogonal property between signal subspace and noise subspace. However, the computational complexity of the MP/UMP method is high, and the angle pairing is always mismatched in the low SNR condition. Therefore, many automatic pairing techniques were presented to solve this problem, such as simultaneous diagonalization [13, 14], simultaneous Schur decomposition [15], and eigenvalue pairing [16].

This paper presents a two-dimensional state-space balance (2D SSB) method that is applied to the DOA estimation in order to get better DOA estimation than the 2D MP/UMP method. On one hand, the comprehensive utilization of controllability and observability matrices has been considered to improve direction abilities. On the other hand, this method can obtain the elevation angle and azimuth angle without pair matching process via automatic pairing technique.

2. Signal Model

The uniform rectangular array (URA) in Figure 1 consists of identical and omnidirectional sensors with interelement spacing of , where is the carrier wavelength.

Suppose far-field narrowband signals impinge on the URA. denotes the th incident signal, is the th elevation angle, and is the th azimuth angle, respectively. The received data at the feeding point of the omnidirectional antennas without noise can be expressed as [17]where and are direction cosines mapping with azimuth and elevation angles and is the gain of the th incident signal.

3. 2D State-Space Balance Method

Equation (1) can be rewritten in a matrix form akin to that of a state-space impulse response aswhere and are diagonal matrices so that they can commute under state-space system equivalence transformations [18]. denotes the square matrix. is the gain vector. Diagonal matrices and contain the DOA information about and . So we will isolate and solve and by using the 2D SSB method. Firstly, Hankel matrix is formed by using the th row of the received data array as follows: where is the window pencil parameter [17]. We define Hankel matrices and as follows:where and are selection matrices. is identity matrix. Then, Hankel matrices and are used to structure the Hankel block matrices ,  , and as follows: where is also the window pencil parameter. Substituting (2) and (4) into (5), Hankel matrices ,  , and can be rewritten as where is the observability matrix, is the controllability matrix, and denotes the transpose.

Considering the singular value decomposition (SVD) of matrix , (6) can be expressed asThe noise effect should be cut off in order to obtain the high accuracy DOA estimation. Hence, observability matrix and controllability matrix can be acquired aswhere is the number of signals, which can be adopted by Akaike information criterion (AIC) [19] and minimum description length (MDL) [20]. Then, through the comprehensive utilization of controllability matrix and observability matrix ,   and can be obtained aswhere is pseudoinverse. For the implementation simplicity, the eigenvalue pairing [16] technique is applied. Therefore, diagonal matrices and can be obtained through eigendecomposition of matrices and where can diagonalize both and in a similarity transformation. Finally, pairs can automatically be paired from diagonal matrices and . Therefore, we can get the final resultsThen the solving steps can be written as follows

Step 1. Construct the Hankel matrices ,  , and from the received data .

Step 2. Use the singular value decomposition of Hankel matrix to acquire the left singular vector , right singular vector , and singular value .

Step 3. Construct the controllability matrix and observability matrix .

Step 4. Use (11) to obtain matrices and ;

Step 5. Use eigenvalue pairing technique to get the diagonal matrices and .

Step 6. Use (13) to estimate final results .

4. Simulation

Firstly, two far-field narrowband signals impinge on the URA with and , where the source positions are relatively dispersed. The distance between any two elements of the URA is half a wavelength. The window pencil parameters are and . The root mean square errors (RMSEs) of azimuth and elevation angles against with SNR for the URA are illustrated in Figure 2, where 500 Monte Carlo experiments are tested for each data. The RMSEs for the estimates of source DOAs and are defined as where is the number of Monte Carlo experiments. As can be seen from Figure 2, the RMSEs of DOAs and decrease as SNR increases, and the 2D SSB method has a better performance than 2D MP/UMP method [17] at the same SNR, which validates the proposed 2D SSB method.

Secondly, the challenging scenario of closely positioned sources impinging on the same URA is taken into consideration with and . Figure 3 shows the results. The 2D SSB method still has a better performance than the 2D MP/UMP method while the performance of 2D MP method deteriorated sharply even the SNR is very high. In addition, the SNR threshold of 2D SSB is much lower than the 2D UMP method. However, in order to acquire a good DOA estimation, a higher SNR threshold is necessary compared with the scenario in the first simulation. All in all, 2D SSB performs more robustness and effectiveness than 2D MP/UMP method.

Thirdly, same scenario with correlated signals and uncorrelated signals is also taken into consideration in Figures 4 and 5. In order to describe the correlated signals and uncorrelated signals, multiple snapshot case is considered. 2D SSB can easily deal with multiple snapshot case by using the method in [21]. The number of snapshots is set to be 4. As illustrated in Figures 4 and 5, 2D SSB also performs better.

In the last simulation, 9 incident signals impinge on the same URA. The 2D DOA estimation results for 300 Monte Carlo simulations are shown in Figures 6 and 7 with SNR = −5 dB and SNR = 5 dB, and the estimate values are quite close to the real values indicated by the cross mark. Obviously, the results demonstrate that the angles of azimuth and elevation are paired correctly and effectively. In addition, the proposed method can compute the same azimuths when the elevations are different or the same elevations when the azimuths are different.

5. Conclusion

A 2D SSB method is proposed in this paper to estimate the angles of azimuth and elevation. This algorithm indicates high performance by using both observability and controllability matrices. Additionally, this method can provide automatically paired 2D direction cosine estimates without matching pair and deal with the correlated signals directly. Simulations demonstrate the effectiveness of our method.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grants nos. 61671035 and 61501012).

References

  1. A. Gaber and A. Omar, “A study of wireless indoor positioning based on joint TDOA and DOA estimation using 2-D matrix pencil algorithms and IEEE 802.11ac,” IEEE Transactions on Wireless Communications, vol. 14, no. 5, pp. 2440–2454, 2015. View at: Publisher Site | Google Scholar
  2. H. Krim and M. Viberg, “Two decades of array signal processing research,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996. View at: Publisher Site | Google Scholar
  3. M. Y. Cao, S. A. Vorobyov, and A. Hassanien, “Transmit array interpolation for DOA estimation via tensor decomposition in 2-D MIMO radar,” IEEE Transactions on Signal Processing, vol. 65, no. 19, pp. 5225–5239, 2017. View at: Google Scholar | MathSciNet
  4. R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986. View at: Publisher Site | Google Scholar
  5. F.-G. Yan, Z.-K. Chen, M.-J. Sun, Y. Shen, and M. Jin, “Two-dimensional direction-of-arrivals estimation based on one-dimensional search using rank deficiency principle,” International Journal of Antennas and Propagation, vol. 2015, Article ID 127621, pp. 1–8, 2015. View at: Publisher Site | Google Scholar
  6. Y. Dong, C. Dong, W. Liu, H. Chen, and G. Zhao, “2-D DOA Estimation for L-Shaped Array With Array Aperture and Snapshots Extension Techniques,” IEEE Signal Processing Letters, vol. 24, no. 4, pp. 495–499, 2017. View at: Publisher Site | Google Scholar
  7. A. Ghobadzadeh and R. Adve, “Low-complexity 2D root-MUSIC pairing for an L-shaped array,” in Proceedings of the IEEE Radar Conference (RadarConf '17), pp. 0957–0962, Seattle, Wash, USA, May 2017. View at: Publisher Site | Google Scholar
  8. R. Roy and T. Kailath, “ESPRIT-estimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Signal Processing, vol. 37, no. 7, pp. 984–995, 1989. View at: Publisher Site | Google Scholar
  9. Y. Fayad, C. Wang, and Q. Cao, “Improved ESPRIT method used for 2D-DOA estimation in tracking radars,” in Proceedings of the 13th European Radar Conference (EuRAD '16), pp. 113–116, London, UK, October 2016. View at: Google Scholar
  10. X. Lan, Y. Li, and E. Wang, “A RARE Algorithm for 2D DOA Estimation Based on Nested Array in Massive MIMO System,” IEEE Access, vol. 4, pp. 3806–3814, 2016. View at: Publisher Site | Google Scholar
  11. T. K. Sarkar and O. Pereira, “Using the matrix pencil method to estimate the parameters of a sum of complex exponentials,” IEEE Antennas and Propagation Magazine, vol. 37, no. 1, pp. 48–55, 1995. View at: Publisher Site | Google Scholar
  12. N. Yilmazer, J. Koh, and T. K. Sarkar, “Utilization of a unitary transform for efficient computation in the matrix pencil method to find the direction of arrival,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 1, pp. 175–181, 2006. View at: Publisher Site | Google Scholar | MathSciNet
  13. J.-F. Cardoso and A. Souloumiac, “Jacobi angles for simultaneous diagonalization,” SIAM Journal on Matrix Analysis and Applications, vol. 17, no. 1, pp. 161–164, 1996. View at: Publisher Site | Google Scholar | MathSciNet
  14. Y.-Y. Dong, C.-X. Dong, Z.-B. Shen, and G.-Q. Zhao, “Conjugate augmented spatial temporal technique for 2-D DOA estimation with L-shaped array,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1622–1625, 2015. View at: Publisher Site | Google Scholar
  15. M. Haardt and J. A. Nossek, “Simultaneous schur decomposition of several nonsymmetric matrices to achieve,” IEEE Transactions on Signal Processing, vol. 46, no. 1, pp. 161–169, 1998. View at: Publisher Site | Google Scholar
  16. M. L. Burrows, “Two-dimensional ESPRIT with tracking for radar imaging and feature extraction,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 2, pp. 524–532, 2004. View at: Publisher Site | Google Scholar
  17. N. Yilmazer and T. K. Sarkar, “2-D unitary matrix pencil method for efficient direction of arrival estimation,” Digital Signal Processing, vol. 16, no. 6, pp. 767–781, 2006. View at: Publisher Site | Google Scholar
  18. J. H. David, State-Space Approaches to Ultra-Wideband Doppler Processing, Worcester Polytechnic Institute, Worcester, Mass, USA, 2007.
  19. M. Wax and T. Kailath, “Detection of signals by information theoretic criteria,” IEEE Transactions on Acoustics, Speech and Signal Processing, vol. 33, no. 2, pp. 387–392, 1985. View at: Publisher Site | Google Scholar | MathSciNet
  20. M. Wax and I. Ziskind, “Detection of the number of coherent signals by the MDL principle,” IEEE Transactions on Signal Processing, vol. 37, no. 8, pp. 1190–1196, 1989. View at: Publisher Site | Google Scholar
  21. N. Yilmazer, S. Ari, and T. K. Sarkar, “Multiple snapshot direct data domain approach and ESPRIT method for direction of arrival estimation,” Digital Signal Processing, vol. 18, no. 4, pp. 561–567, 2008. View at: Publisher Site | Google Scholar

Copyright © 2017 Jun Wang 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.


More related articles

659 Views | 279 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.