140 IEEE SIGNAL PROCESSING LETTERS, VOL. 21, NO. 2, FEBRUARY 2014 Improved Unitary Root-MUSIC for DOA Estimation Based on Pseudo-Noise Resampling Cheng Qian, Lei Huang, and H. C. So Abstract A novel pseudo-noise resampling (PR) based unitary root-music algorithm for direction-of-arrival (DOA) estimation is derived in this letter. Our solution is able to eliminate the abnormal DOA estimator called outlier and obtain an approximate outlier-free performance in the unitary root-music algorithm. In particular, we utilize a hypothesis test to detect the outlier. Meanwhile, a PR process is applied to form a DOA estimator bank and a corresponding root estimator bank. We propose a distance detection strategy which exploits the information contained in the estimated root estimator to help determine the final DOA estimates when all the DOA estimators fail to pass the reliability test. Furthermore, the proposed method is realized in terms of real-valued computations, leading to an efficient implementation. Simulations show that the improved MUSIC scheme can significantly improve the DOA resolution at low signal-to-noise ratios and small samples. Index Terms Direction-of-arrival (DOA) estimation, pseudonoise resampling (PR), small sample size, unitary root-music. I. INTRODUCTION MANY applications in radar, sonar and wireless communications require direction-of-arrival (DOA) information of multiple sources impinging upon a sensor array. Numerous DOA estimators have been proposed in the literature. Among them, subspace methods such as ESPRIT [1], [2] and MUSIC [3], offer a good compromise between accuracy and computational complexity. However, this kind of algorithms suffers considerable performance degradation when the number of snapshots and/or signal-to-noise ratio (SNR) is small [4], [5], [6]. The most visual embodiment of this phenomenon is that the observed estimation errors rapidly depart from the Cramér-Rao lower bound (CRLB) below a threshold SNR [5]. Gershman and Böhme [6] have proposed an enhanced weighted MUSIC estimator which can improve the DOA estimation performance via a pseudo-noise resampling (PR) of spatial spectrum. Compared with the conventional MUSIC Manuscript received October 17, 2013; revised November 26, 2013; accepted December 06, 2013. Date of publication December 11, 2013; date of current version December 19, 2013. This work was supported in part by a grant from the NSFC/RGC Joint Research Scheme, sponsored by the Research Grants Council of Hong Kong and the National Natural Science Foundation of China (Projects N_CityU 104/11 and 61110229/61161160564), by the National Natural Science under Grants 61222106 and 61171187, and by the Shenzhen Kongqie talent program under Grant KQC201109020061A. The associate editor coordinating the review of this manuscript and approving it for publication was Prof. Muhammad Zubair Ikram. C. Qian and L. Huang are with the Department of Electronic and Information Engineering, Harbin Institute of Technology Shenzhen Graduate School, Shenzhen, Guangdong 518055, China (e-mail: alextoqc@gmail.com; dr.lei.huang@ieee.org). H. C. So is with the Department of Electronic Engineering, City University of Hong Kong, Hong Kong 852, China (e-mail: hcso@ee.cityu.edu.hk). Color versions of one or more of the figures in this paper are available online at http://ieeexplore.ieee.org. Digital Object Identifier 10.1109/LSP.2013.2294676 algorithm, it significantly increases the DOA resolution in low SNR scenario. But it still needs to search the spectrum, requiring huge computational burden which may hinder its practical application. The idea behind the PR technique which is motivated by the success of modern resampling schemes (e.g., bootstrap [7], [8]), is to utilize synthetically generated psuudo-noise to perturb the original noise in such a way that the outlier will be removed [9]. Vasylyshyn [10] has proposed a variant of the PR-based root-music algorithm which combines the PR technique with the conventional beamformer. In [11], a PR technique based unitary ESPRIT algorithm has been developed to mitigate the effect of outliers. In [12], the author utilizes the PR process and local performance test to improve the accuracy of the beamspace ESPRIT algorithm. However, the ESPRIT-like algorithms suffer performance degradation when the number of overlapping sensors is large. In this paper, we derive an improved unitary root-music algorithm that is based on the PR technique to reduce the threshold SNR. Our approach utilizes a PR process to construct a DOA estimator bank and a corresponding root estimator bank. Then a reliability test [6], [10], [11] is applied to the whole DOA estimator bank and only the DOA estimator that has passed the reliability test can be retained. Furthermore, we have proposed a DOA selection strategy to determine the final DOA estimates by replacing the outliers with appropriate candidates in the retained DOA estimators. To combat the case where there is no DOA estimator passing the reliability test, unlike the conventional median average method [10], [11], we propose a distance detection strategy (DDS) which utilizes the information contained in the root estimator bank to determine the final DOA estimates. Due to the real-valued formulation, our proposal enables an efficient implementation when a moderate number of PR processes is given. Simulations show that the PR based unitary root-music scheme outperforms its counterparts in the low SNR and small sample situations. II. PROBLEM FORMULATION AND UNITARY ROOT-MUSIC A. Signal Model Consider a uniform linear array (ULA) with isotropic sensors. There are uncorrelated narrowband source signals impinging on the array from distinct directions in the far field. The observation vector is Here, is the steering matrix, is the source signal vector, is the number of snapshots and the steering vector due to the th source is expressed as (1) (2) 1070-9908 2013 IEEE
QIAN et al.: IMPROVED UNITARY ROOT-MUSIC FOR DOA ESTIMATION 141 where is the transpose, is the carrier wavelength and is the interelement spacing. It is assumed that the noise is a white Gaussian process with zero mean and covariance,where is the identity matrix. Moreover, the noise is uncorrelated with the signal. The covariance matrix of is where, stands for the mathematical expectation and represents the conjugate transpose. B. Unitary Root-MUSIC Algorithm The unitary root-music algorithm [13] utilizes a real-valued covariance matrix given as where represents the real part, is an exchange matrix with ones on its antidiagonal and zeros elsewhere, represents complex conjugate and isasparseunitarymatrix, defined as [11], [13] (5) Here, is an zero vector. Define the eigenvalue decomposition (EVD) of as where, diag and with being the signal eigenvalues, being its corresponding signal eigenvectors, being the noise eigenvectors and diag being a diagonal matrix. Then the unitary root-music polynomial can be expressed as where with being the root of (7). Through finding the roots which are closest to the unit circle, we determine the DOAs: where represents the angle operator. It has been shown in [13] that the forward-backward (FB) and unitary root-music polynomials are identical. This in turn implies that the unitary root-music has the same performance as that of the FB root-music. But the unitary root-music method is much more computationally efficient since it is realized in terms of real-valued computation. However, the unitary root-music will suffer performance degradation, especially in the low SNR and small sample scenarios. This result is due to the fact that the unitary root-music estimator cannot efficiently handle the outliers. To circumvent this issue, we devise a PR based unitary root-music algorithm for computationally efficient and accurate DOA estimation. for for (3) (4) (6) (7) (8) III. PROPOSED ALGORITHM Let us now formulate the PR based unitary root-music algorithm. Our aim is to remove the outliers and recover the approximate outlier-free performance of the unitary root-music by applying the PR technique. To this end, we resample the data matrix several times using synthetically generated pseudo-noise [6], [11]. For each PR process, we run the unitary root-music method to estimate the DOAs, and then utilize a reliability test to remove the outliers that are rejected by the test and only retain the remaining DOAs. A. Pseudo-Noise Resampling Process The key of the proposed method is to test the following hypothesis for each DOA estimator [6], [10], [11], which enables us to separate the normal and abnormal DOA estimators. : All the DOA estimates in a DOA estimator are localized in. Here, is the pre-estimated angular sectors of source localization. A simple way of determining is to use the conventional beamformer [10], [12], [14]. Let where are the coordinates of the highest peaks of the conventional beamformer output, and are the left and right boundaries of the th subinterval, and both of them can be chosen as angular distances between the maximum of the th peak and the left/right neighbor point with 3 db drop, respectively. When the hypothesis is rejected, the data matrix will be resampled times using synthetically generated pseudo-noise. The resampled data matrix is given as where is the pseudo-noise matrix with mean-zero and covariance matrix. Here, is the variance of the pseudo-noise and its value should be comparable with the variance of the measurement noise.itis shown in [6], [10], [11], [12] that we can estimate as where is a user-defined parameter and is the consistent estimate of given by with being the ordered eigenvalues of the real-valued covariance matrix. Here, is the estimated real-valued sample covariance matrix with being the sample covariance. B. DOA Selection Strategy and DDS For each resampling run, we apply the unitary root-music to obtain a DOA estimator which contains DOA estimates and the corresponding root estimator.letthe th estimator be (9) (10) where are the DOA estimates and is the th DOA root corresponding to. Assuming that after PR runs, we have estimators which are used to form the estimator bank (11)
142 IEEE SIGNAL PROCESSING LETTERS, VOL. 21, NO. 2, FEBRUARY 2014 Performing the reliability test to yields the following two subsets TABLE I PSEUDO-CODE OF PROPOSED ALGORITHM (12) where consists of estimators that are successfully accepted by and contains the remaining estimators that are rejected by. In addition, we define (13) to represent the th estimators accepted by. If, i.e., there is at least one estimator in,areliable process to determine the final DOA estimator is to average the successfully resampled estimators. Since the DOA estimates in each estimator are sorted in ascending order, we have where the th DOA in is calculated as (14) If, which means all the resampled estimators fail to pass,weproposearobustddstodeterminethefinal DOA estimates. It is worth mentioning that the DDS exploits the information contained in the root estimator. However, it should be noted that this information is ignored in [10] and only the median average method [6], [11] is used to determine the final DOA estimates from. SpecificstepsoftheDDSaregivenas follows: 1) Divide into subsets with each subset contains roots corresponding to a common DOA. Then the th subset can be expressed as 2) Calculate the modulus of each in (15) 3) Find the maximum and store its index, i.e., (15) (16) (17) where is the maximum value of and is the corresponding index. 4) Substitute into (8) to obtain the final DOA estimates. When, all the DOA estimators are stored in. The PR unitary ESPRIT algorithm [11] discards the whole estimator bank and utilizes the standard unitary ESPRIT to re-estimate the final DOAs. Therefore, it has lost all the information brought by the PR technique. In [10], the median average is utilized to divide the DOA estimator bank into subsets with each subset containing DOA estimates that are corresponding to a common DOA and sorted in descending order. Then the th final DOA estimate can be determined by choosing the average of the two DOA estimates in the middle of the th subset if is even or choosing the median of the th subset if is odd. In [4], a conventional average method that just takes the average of th subsets as the th DOA estimate is suggested. Clearly, in each subset, there always exists a DOA estimate that is closest to the true DOA. However, both the conventional and the median average methods cannot determine such a relatively optimal DOA estimate. As a result, it may not guarantee that its final DOA estimates are closest to the true DOA values. On the contrary, the DDS can efficiently circumvent this problem. It is observed from (15) (17) that for each subset, i.e.,, the DDS always selects the root that is closest to the unit circle as the th final DOA estimate. Thus, all the DOA estimates determined by theddsarealwaystheoneclosesttothetruedoas.theeffectiveness of the DDS will be further verified in Section IV. Note that the computational complexity of the proposed algorithm is higher than the unitary root-music by a factor,where is the probability of outliers, i.e., the probability that the unitary root-music estimator may be rejected by. The additional computational burden improves the performance in return. Due to the low complexity of unitary root-music, the proposed algorithm can be implemented efficiently when a proper is given. IV. SIMULATIONS We compare the proposed method with the root-music, unitary root-music, PR root-music [10] and PR unitary ESPRIT [11] algorithms in terms of root mean square error (RMSE). In order to illustrate the effectiveness of the DDS, we add another DDS based algorithm for comparison, which is obtained via replacing the median average method in the PR root-music with the proposed DDS, i.e., PR root-music with DDS. Meanwhile, we examine their ability of removing outliers as well, namely, the probability of outliers. The CRLB is plotted as a benchmark. In our simulations, two independent narrowband Gaussian signals are assumed to imping upon a ULA with omnidirectional sensors from directions and. The noise is zero-mean white Gaussian process. The SNR is defined as the ratio of the power of all source signals to that of the additive noise at each sensor. Furthermore, the user defined parameter is set to be 1 in all the simulations. 1000 Monte Carlo simulations are carried out to evaluate the RMSE, which is defined as RMSE.Inall experiments, we assume that the source localization sectors and the number of sources are known or estimated by [15], [16]. According to (9), we use the conventional beamformer to pre-estimate in each independent run. In the first simulation, we study the ability of removing outliers under small sample size as a function of SNR. We set the
QIAN et al.: IMPROVED UNITARY ROOT-MUSIC FOR DOA ESTIMATION 143 Fig. 1. Probability of outliers versus SNR. ( ). Fig. 3. RMSE angle error performance versus SNR. ( ). Fig. 4. RMSE angle error performance versus.( db). Fig. 2. Probability of outliers versus SNR. ( ). number of PR runs as. Figs. 1 and 2 show the outlier probability for and, respectively. It is seen that the proposed algorithm provides a considerable improvement to remove the outliers. When, the proposed algorithm has almost removed all the outliers. However, the other DOA estimation algorithms suffer a high probability of outliers no matter is small or large. Moreover, we find that the DDS based PR root-music has achieved almost the same outlier removing performance compared to the median average based PR root-music algorithm [10]. In the second example, we set and vary the signal power such that the input SNR goes from - to 7 db. The number of snapshots is. We examine the performance with different numbers of PR processes. For the PR root-music and PR unitary ESPRIT algorithms, we set.thermsesare depicted in Fig. 3. It can be seen that in low SNR regime, the proposed algorithm outperforms the conventional and unitary root-music schemes no matter how small or large the is. Fig. 3 also implies that a relatively larger can lead to a better performance for the proposed method and the value of is sufficient to obtain a satisfactory performance. When the SNR is larger than 3 db, the proposed and unitary root-music algorithms converge together. This is due to the fact that, at high SNRs, all the DOA estimates are accepted by the hypothesis test, and the proposed one is reduced to be the unitary root- MUSIC scheme. Let us now study the RMSE versus the sample size. In this example, we fix thesnrto-2db.foralltheprbasedalgorithms, we set the number of PR runs to be for fair comparison. It is observed from Fig. 4 that the proposed approach performs the best. When, the proposed PR unitary root-music algorithm achieves a better performance improvement compared to the unitary root-music scheme. From Figs.3and4,wefind the DDS-based PR root-music algorithm always outperforms the median average based one, especially at low SNRs and small samples. V. CONCLUSION A PR based unitary root-music algorithm has been derived for DOA estimation. In contrast to the existing root-music methods, the proposed scheme combines the unitary root-music and PR techniques to form a DOA estimator bank and a corresponding root estimator bank. Meanwhile, a DOA selection strategy with the DDS is proposed to select the reliable DOA estimates. Here, the DDS which exploits the DOA roots information inherent in the unitary root-music algorithm can help to determine the final DOA estimates when all the DOA estimators are rejected by the reliability test. Simulation results verify the effectiveness of the proposed algorithm. REFERENCES [1] R. Roy, A. Paulraj, and T. Kailath, Estimation of signal parameters via rotational invariance techniques-esprit, in Proc. 19th Asilomar Conf. Circuits, Syst. Comput., Monterey, CA, Nov. 1985. [2] C. Qian, L. Huang, and H. C. So, Computationally efficient ES- PRIT algorithm for direction-of-arrival estimation based on Nyström method, Signal Process., vol. 94, no. 1, pp. 74 80, 2014. [3] R. Schmidt, Multiple emitter location and signal parameter estimation, IEEE Trans. Antennas Propag., vol. 34, no. 3, pp. 276 280, 1986.
144 IEEE SIGNAL PROCESSING LETTERS, VOL. 21, NO. 2, FEBRUARY 2014 [4] A. B. Gershman and J. F. Bohme, A pseudo-noise resampling approach to direction of arrival estimation using estimator banks, in Proc. 9th IEEE SP Workshop Statistical Signal Array Process., Portland, OR, USA, 1998, pp. 244 247. [5] J. Thomas, L. Scharf, and D. Tufts, The probability of a subspace swap in the SVD, IEEE Trans. Signal Process., vol. 43, no. 3, pp. 730 736, 1995. [6] A. B. Gershman and J. F. Böhme, Improved DOA estimation via pseudorandom resampling of spatial spectrum, IEEE Signal Process. Lett., vol. 4, no. 2, pp. 54 57, 1997. [7] A. M. Zoubir and B. Boashash, The bootstrap and its application in signal processing, IEEE Signal Process. Mag., vol. 15, no. 1, pp. 56 76, 1998. [8]A.M.Zoubir, Bootstrap Techniques for Signal Processing. Cambridge, U.K.: Cambridge Univ. Press, 2004. [9] A. B. Gershman and J. F. Böhme, A pseudo-noise approach to direction finding, Signal Process., vol. 71, no. 1, pp. 1 13, 1998. [10] V. Vasylyshyn, Removing the outliers in root-music via pseudonoise resampling and conventional beamformer, Signal Process., vol. 93, no. 12, pp. 3423 3429, 2013. [11] A. B. Gershman and M. Haardt, Improving the performance of unitary ESPRIT via pseudo-noise resampling, IEEE Trans. Signal Process., vol. 47, no. 8, pp. 2305 2308, 1999. [12] V. Vasylyshyn, Improved Beamspace ESPRIT-based DOA estimation via pseudo-noise resampling, in IEEE 9th Eur. Radar Conf., Amsterdam, The Netherlands, 2012, pp. 238 241. [13] M. Pesavento, A. B. Gershman, and M. Haardt, Unitary root-music with a real-valued eigendecomposition: A theoretical and experimental performance study, IEEE Trans. Signal Process., vol. 48, no. 5, pp. 1306 1314, 2000. [14] A. B. Gershman, Pseudo-randomly generated estimator banks: A new tool for improving the threshold performance of direction finding, IEEE Trans. Signal Process., vol. 46, no. 5, pp. 1351 1364. [15] M. Wax and T. Kailath, Detection of signals by information theoretic criteria, IEEE Trans. Acoust., Speech Signal Process., vol. 33, no. 2, pp. 387 392, 1985. [16] L. Huang, S. Wu, and X. Li, Reduced-rank MDL method for source enumeration in high-resolution array processing, IEEE Trans. Signal Process., vol. 55, no. 12, pp. 5658 5667, 2007.