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

Research Article

A FPC-ROOT Algorithm for 2D-DOA Estimation in Sparse Array

Wenhao Zeng,¹ Hongtao Li,¹ Xiaohua Zhu,¹ and Chaoyu Wang²

¹Nanjing University of Science and Technology, Nanjing, Jiangsu 210094, China
²China Shipbuilding 724 Research Institutions, Nanjing, Jiangsu 210003, China

Received 24 December 2015; Revised 11 March 2016; Accepted 22 March 2016

Academic Editor: Jeich Mar

Copyright © 2016 Wenhao Zeng 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.

To improve the performance of two-dimensional direction-of-arrival (2D DOA) estimation in sparse array, this paper presents a Fixed Point Continuation Polynomial Roots (FPC-ROOT) algorithm. Firstly, a signal model for DOA estimation is established based on matrix completion and it can be proved that the proposed model meets Null Space Property (NSP). Secondly, left and right singular vectors of received signals matrix are achieved using the matrix completion algorithm. Finally, 2D DOA estimation can be acquired through solving the polynomial roots. The proposed algorithm can achieve high accuracy of 2D DOA estimation in sparse array, without solving autocorrelation matrix of received signals and scanning of two-dimensional spectral peak. Besides, it decreases the number of antennas and lowers computational complexity and meanwhile avoids the angle ambiguity problem. Computer simulations demonstrate that the proposed FPC-ROOT algorithm can obtain the 2D DOA estimation precisely in sparse array.

1. Introduction

Two-dimensional direction-of-arrival (2D DOA) [1–4] estimation plays an important role in array signal processing, and it usually utilizes L-shaped arrays, planar arrays, vector sensors, and so forth, to estimate the two-dimensional target angle. Compared with 2D DOA estimation of other arrays [5, 6], 2D DOA estimation of planar array is asymptotically unbiased where the target position can be determined precisely without the problem of angle ambiguity. Moreover, the phased-array radar in the planar array can utilize waveform diversity which makes it suitable for multiple targets and complex environments with enhanced performance of target recognition and target tracking. However, the planar array requires a large amount of array elements which results in tremendous demands of hardware equipment and increases the complexity and cost of design. Comparatively, the sparse array needs only a quite smaller number of elements, so it is a feasible remedy for the abovementioned problems, whereas, in this case, the problem of substantially rising average side lobes occurs, which accounts for its uncommon applications so far. Paper [5] lowers computational complexity of DOA estimation by introducing a preprocessing transformation matrix but has poor performance in sparse array. Paper [7] achieves high-resolution DOA estimation in sparse array; however, it only works in uniform linear array (ULA).

Matrix completion [8–10] is an extension of compressive sensing (CS) [11, 12] and has been widely applied to image processing, remote sensing, and many other engineering fields [13, 14]. On the basis of low-rank matrix, matrix completion can recover the full matrix from an incomplete set of matrix entries with high probability by solving nuclear-norm optimization. Paper [15] proposed a fixed point and Bregman iterative methods for matrix rank minimization, which is a very fast, robust, and powerful algorithm. Paper [16] introduces a particularly simple yet highly efficient alternating projection algorithm and this algorithm is able to recover matrix from the minimum number of measurements necessary. Paper [17] provides a necessary and sufficient condition that quantifies when this heuristic successfully finds the minimum rank solution of a linear constraint set. The sparse received signal can be recovered to the full received signal by matrix completion, so that matrix completion has been applied to DOA estimation [14, 18].

This study proposes a Fixed Point Continuation Polynomial Roots (FPC-ROOT) algorithm. The proposed algorithm establishes a signal model of matrix completion for 2D DOA estimation based on low-rank property of target in two-dimensional space domain, which turns out to meet the Null Space Property (NSP) and ensure the feasibility of 2D DOA estimation via matrix completion. By contrast to conventional algorithms, the proposed algorithm obtains left and right singular vectors of received signal matrix directly from the output of matrix completion algorithm instead of eigendecomposing autocorrelation matrix of received signal, for the benefit of lower dimensions. Besides, the proposed algorithm can avoid scanning of the 2D spectral peak by solving the polynomial roots, which reduces the computational complexity. In addition, it can estimate target angle accurately with less number of array units by adopting matrix completion.

2. Signal Model in Sparse Array

2.1. The Signal Model of 2D DOA Estimation

Suppose that the uniform rectangular array (URA) [3] is depicted in Figure 1. and are the number of elements in -direction and -direction, respectively. The corresponding element spacing is and .

Figure 1: The uniform rectangular array.

Assume that there are multiple targets in the space. The number of snapshots is and the number of targets is , where the waveform of the th () target is and . The target angles are , where and express the elevation and azimuth angles of the th target, respectively. Suppose the uncorrelated narrow-band signals are received in far field [19]; thus, the received signal of elements in -direction can be written as

where is a diagonal matrix of dimensions, the first elements of are nonzero which correspond to targets, and is the noise vector of the array. is a dimensional steering vector, which can be given by

Assume the wavelength of received signal is ; then, is

Similarly, the received signal of elements in -direction can be represented as

where is the noise vector of the array. is dimensional steering vector, which can be denoted as

where is defined as

In conclusion, the received signal model of full array is demonstrated as

When the power of noise matrix is much smaller than the signal power, we can get that ; in other words, matrix is low-rank [10]. And matrix in (7) is consistent with strong incoherence property [15], so it can be recovered precisely with high probability via matrix completion.

2.2. Matrix Completion

With regard to the sparsity of signals, CS can sample signals at far lower sampling frequency than Nyquist sampling frequency and meanwhile reconstruct original signals precisely. In CS, the target to be recovered is a vector; however, on some practical occasions, it normally refers to a matrix and is sensitive to data missing, data corruption, and so on.

Matrix completion is an extension of CS. CS exploits the sparsity of signals; nevertheless, matrix completion utilizes the sparsity of matrix singular values, namely, the low-rank property, and reconstructs full matrix by solving nuclear-norm optimization.

Suppose is an original low-rank matrix, where , . Matrix completion refers to recovering the whole elements of matrix from its partial elements. Let the available entries of matrix be , where is a random sampled subset of the whole elements; then, the matrix completion can be denoted as [10]

where matrix is the estimated value of matrix ; represents the nuclear norm of the matrix .

2.3. The Sparse Array Model

We sample elements from URA at random. Let the number of sampled elements be with the corresponding element positions unchanged and the unsampled units removed; thus, a new sparse array is set up. To recover matrix precisely via matrix completion, the required number of samples is [10], while , are the numbers of rows and columns of the matrix, and is the rank of the matrix. Suppose the received signal matrix in sparse array is ; then, is related to by

where is the set of element positions in sparse array.

3. The Proposed Algorithm

3.1. The Signal Model for 2D-DOA Estimation Based on Matrix Completion

According to (7), the received signal matrix in URA is low-rank. Regard as the recovered matrix by matrix completion, so it can be deduced that

If a signal model meets NSP, the rank minimization of this model is equivalent to its nuclear-norm minimization [13]. So we introduce the projection as follows:

From (11), we can see that (10) can be expressed as

where is the estimated matrix. The null space of projection is

According to NSP, we are unable to recover a matrix if it belongs to the null space of projection . From (4), we can get that

In (7), the arbitrary elements of and are both nonzero and the first elements in are nonzero, so the diagonal elements of dimensional diagonal matrix are nonzero. Arbitrary elements of are nonzero based on matrix multiplication property. Therefore, it can be concluded that whichever sampling operator is selected, there is always ; in other words, , which meets the Null Space Property.

3.2. The FPC-ROOT Algorithm

Suppose and matrices and are the left and right singular value matrix of . FPC-ROOT algorithm can obtain and directly by matrix completion based on the low-rank property of and realizes 2D DOA estimation by computing the polynomial roots.

Firstly, to solve problem (12), we achieve the left and right singular value matrices and of estimated value in full array by FPC algorithm. And then, we solve the autocorrelation matrix of received signal in full array, so the eigendecomposition of the autocorrelation matrix is

where , is the conjugate transpose matrix of , is the signal subspace, is the noise subspace, and is

And in (15) can be expressed as

It is obvious that the spanned subspaces of and are identical. We solve the conjugate transpose of autocorrelation matrix and then it can be concluded that

where , and similarly is

can be expressed as

So the space spanned by is the same as the space spanned by .

Based on singular value decomposition, it is easy to know that left and right singular value matrices and are related to the eigendecomposition of autocorrelation matrix by

So the left and right singular vectors of received signal matrix can be obtained directly from the output of matrix completion algorithm instead of eigendecomposing autocorrelation matrix of received signal and the computational complexity of the proposed algorithm obtained a corresponding reduction.

Consider a polynomial as

where is the th eigenvector of matrix and . It can be deduced that the corresponding ROOT-MUSIC polynomial is

By solving the polynomial roots, the proposed algorithm can avoid scanning of the 2D spectral peak. Taking use of the orthogonality of signal subspace, the above equation can be converted to

The dimensionality of is smaller than ; hence, the proposed algorithm effectively reduced the computational complexity by reducing the dimension of the subspace matrix. Then, the vector corresponding to the roots of polynomial is

Thus,

Similarly, it can be obtained that

Hence, the estimated angles of targets are

In summary, the proposed FPC-ROOT algorithm in this paper can be programmed as follows.

Fixed Point Continuation Polynomial Roots (FPC-ROOT) Algorithm

 Initialize. Given , the projection is corresponding to element positions in the sparse array:

(1)Solve (12) to get left and right singular value matrices and .(2)Construct the polynomial roots (24).(3)Solve (24) to get

(4)Determine the targets angles

 End.

3.3. Dimensionality Reduction

In normal 2D-DOA estimation, DOA estimation is achieved through autocorrelation matrix of a vector which is transformed from the signal matrix. Suppose the signal matrix is , so the dimension of the transformed vector is ; then, the dimension of the autocorrelation matrix is . In the proposed method, we obtain left and right singular vector directly from the signal matrix, avoiding solving the autocorrelation matrix, greatly reducing the dimension of the target matrix.

4. Simulations

In this section, several simulations for 2D DOA estimation are conducted to demonstrate the feasibility and effectiveness of the proposed algorithm in sparse array. In these experiments, we sample 1200 elements from full array at random to formulate a sparse array. The full array is an URA of elements, where the total number of elements is and the corresponding element spacing in -direction and -direction is .

In the first examples, 2D DOA estimation of sparse array is shown. Let three targets be in the space domain, whose elevation and azimuth angles of 2D-DOA are , , and , respectively, and input signal-to-noise ratio (SNR) is 20 dB. Let the number of snapshots be 50. Suppose the received signals are narrow-band in far field and the signal sources are uncorrelated; then, 100 experiment results of 2D-DOA estimation by FPC-ROOT algorithm are depicted in Figure 2. The proposed algorithm samples 1200 elements to construct a sparse array; in other words, 70% units are removed from full array. From Figure 2, it can be seen that the targets’ angles can be achieved precisely by the proposed algorithm in sparse array.

Figure 2: 2D DOA estimation via FPC-ROOT algorithm.

In the second experiment, recovery errors by matrix completion in sparse array with different SNR are examined. We recover a full matrix from a sparse matrix using FPC algorithm. Suppose the full matrix is and the recovered matrix is ; then, the measurement criterion of recovery errors by matrix completion is . Figure 3 illustrates the variation of recovery errors by matrix completion with different SNR when the numbers of sparse array elements are 900, 1200, and 1500, respectively. It can be deduced easily from this experiment that recovery errors by FPC algorithm are inversely proportional to SNR. The larger the SNR is, the more similar the matrix is to low-rank matrix. Simultaneously, the more elements the sparse array has, the closer the recovered signals are to the full array and the smaller the corresponding recovery errors are.

Figure 3: Recovery errors by matrix completion with different elements and SNR.

In the third experiment, root mean square error (RMSE) of 2D DOA estimation based on matrix completion is discussed. When 500 Monte Carlo simulations are performed, Figure 4 exhibits RMSE with different SNR via various algorithms. When the SNR is high, the recovery of matrix completion has a high precision, so the RMSE of DOA estimation is decreased. And it is evident that the proposed algorithm has higher accuracy for both elevation and azimuth angles estimation than 2D-MUSIC algorithms.

Figure 4: RMSE of 2D DOA estimation based on matrix completion.

In the last experiment, RMSE by FPC-ROOT algorithm with different SNR and elements is demonstrated. Let the number of snapshots be 50 and let 500 Monte Carlo experiments be implemented. Figure 5 shows RMSE by FPC-ROOT algorithm when SNR and the number of elements vary. As the recovery errors by FPC algorithm are inversely proportional to SNR, the RMSE of DOA estimation is decreased at high SNR. From this simulation, it can be concluded that the estimation accuracy by FPC-ROOT algorithm is proportional to the number of elements.

Figure 5: RMSE versus SNR with different elements by FPC-ROOT algorithm.

5. Conclusion

In this paper, a FPC-ROOT algorithm is proposed based on matrix completion, which can achieve high accuracy of 2D DOA estimation with reduced antenna units. The proposed algorithm obtains left and right singular vectors of received signal by the output of matrix completion algorithm directly instead of eigendecomposing the autocorrelation matrix of received signal, for the benefit of lower dimensions. Besides, by computing polynomial roots, the proposed algorithm can avoid the scanning of two-dimensional spectral peak, which cuts down the computational complexity.

Competing Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 61401204), Postdoctoral Science Foundation of Jiangsu Province (Grant no. 1501104C), Technology Research and Development Program of Jiangsu Province (Grant no. BY2015004-03), and the Fundamental Research Funds for the Central Universities (Grant no. 30916011319).

References

J. Li and X. Zhang, “Closed-form blind 2D-DOD and 2D-DOA estimation for MIMO radar with arbitrary arrays,” Wireless Personal Communications, vol. 69, no. 1, pp. 175–186, 2013.

Z. Zheng, G. Li, and Y. Teng, “Simplified estimation of 2D DOA for coherently distributed sources,” Wireless Personal Communications, vol. 62, no. 4, pp. 907–922, 2012.

P. Heidenreich, A. M. Zoubir, and M. Rubsamen, “Joint 2-D DOA estimation and phase calibration for uniform rectangular arrays,” IEEE Transactions on Signal Processing, vol. 60, no. 9, pp. 4683–4693, 2012.

M. Agatonović, Z. Stanković, N. Dončov, L. Sit, B. Milovanović, and T. Zwick, “Application of artificial neural networks for efficient high-resolution 2D DOA estimation,” Radioengineering, vol. 21, no. 4, pp. 1178–1186, 2012.

W. Zhang, W. Liu, J. Wang, and S. Wu, “Computationally efficient 2-D DOA estimation for uniform rectangular arrays,” Multidimensional Systems and Signal Processing, vol. 25, no. 4, pp. 847–857, 2014.

L. Yuan, R. Jiang, and Y. Chen, “Gain and phase autocalibration of large uniform rectangular arrays for underwater 3-D sonar imaging systems,” IEEE Journal of Oceanic Engineering, vol. 39, no. 3, pp. 458–471, 2014.

N. Hu, Z. Ye, X. Xu, and M. Bao, “DOA estimation for sparse array via sparse signal reconstruction,” IEEE Transactions on Aerospace & Electronic Systems, vol. 49, no. 2, pp. 760–773, 2013.

O. Klopp, “Noisy low-rank matrix completion with general sampling distribution,” Bernoulli, vol. 20, no. 1, pp. 282–303, 2014.

H. Ji, C. Liu, Z. Shen, and Y. Xu, “Robust video denoising using low rank matrix completion,” in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR '10), pp. 1791–1798, San Francisco, Calif, USA, June 2010.

E. J. Candès and B. Recht, “Exact matrix completion via convex optimization,” Foundations of Computational Mathematics, vol. 9, no. 6, pp. 717–772, 2009.

D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006.

J. Cheng, Q. Ye, H. Jiang, D. Wang, and C. Wang, “STCDG: an efficient data gathering algorithm based on matrix completion for wireless sensor networks,” IEEE Transactions on Wireless Communications, vol. 12, no. 2, pp. 850–861, 2013.

W. Chen, Q. Tian, J. Liu, and Q. Wang, “Nonlocal low-rank matrix completion for image interpolation using edge detection and neural network,” Signal, Image and Video Processing, vol. 8, no. 4, pp. 657–663, 2014.

D. S. Kalogerias and A. P. Petropulu, “Matrix completion in colocated mimo radar: recoverability, bounds & theoretical guarantees,” IEEE Transactions on Signal Processing, vol. 62, no. 2, pp. 309–321, 2014.

S. Ma, D. Goldfarb, and L. Chen, “Fixed point and Bregman iterative methods for matrix rank minimization,” Mathematical Programming, vol. 128, no. 1-2, pp. 321–353, 2011.

J. Tanner and K. Wei, “Normalized iterative hard thresholding for matrix completion,” SIAM Journal on Scientific Computing, vol. 35, no. 5, pp. S104–S125, 2013.

B. Recht, W. Xu, and B. Hassibi, “Necessary and sufficient conditions for success of the nuclear norm heuristic for rank minimization,” in Proceedings of the 47th IEEE Conference on Decision and Control (CDC '08), pp. 3065–3070, Cancún, Mexico, December 2008.

P. Pal and P. P. Vaidyanathan, “A grid-less approach to underdetermined direction of arrival estimation via low rank matrix denoising,” IEEE Signal Processing Letters, vol. 21, no. 6, pp. 737–741, 2014.

H. Krim and M. Viberg, “Two decades of array signal processing research: the parametric approach,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996.