ISRN Signal ProcessingVolume 2011 (2011), Article ID 751670, 8 pagesdoi:10.5402/2011/751670

Research Article

Fast DOA Estimation Algorithm Based on a Combination of an Orthogonal Projection and Noise Pseudo-Eigenvector Approach

Ching Jer Hung

Department of Electrical Engineering, Kun-Shan University, No. 949, Dawan Rd., Yongkang City, Tainan County 710, Taiwan

Received 3 November 2010; Accepted 5 December 2010

Academic Editors: C.-W. Kok, E. D. Übeyli, and Y. Xiang

Copyright © 2011 Ching Jer Hung. Ching Jer Hung. 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.

This paper presents a new fast direction of arrival (DOA) estimation technique, using both the projection spectrum and the eigenspectrum. First, the rough DOA range is selected using the projection spectrum; then, a linear matrix equation is used to acquire a noise pseudo-eigenvector. Finally, the fine DOA estimation is obtained from an eigenspectrum approach based on the noise pseudo-eigenvector. Without the need to form the covariance matrix from a block of the array data and without a prior knowledge of the number of incoming signals, reduced complexity is achieved, in contrast to conventional subspace-based algorithms. Simulation results show that the proposed algorithm has a good resolution performance and deals well with both uncorrelated and correlated signals. Since the new approach can reduce computational complexity while maintaining better or similar resolution capability, it may provide wider application prospects in real-time DOA estimation when contrasted to other comparable methods.

1. Introduction

The direction of arrival (DOA) has been applied in many fields, including speech, radio, telecommunication, and medical signal processing. In the past few decades, a large number of accurate high-resolution DOA estimation techniques have been proposed. One popular approach is the subspace-based method, used in the well-known MUltiple Signal Classification (MUSIC) algorithm [1] and Estimation of Signal Parameters by Rotational Invariance Techniques (ESPRIT) [2]. However, the conventional subspace-based method has some inherent problems that require the assumption of the uncorrelated incident signals, prior knowledge of the number of incoming signals, a block of data snapshots to form a covariance matrix, and an eigendecomposition computation. We know that the estimations for the covariance matrix and the number of the incoming signals, the computation for eigendecomposition, and the procedure for decorrelating the high correlated signals are computationally intensive and time-consuming. Therefore, the conventional subspace-based will remain severely limited in fast DOA estimation under a nonstationary environment if the above issues are not resolved.

A direct data domain least squares (D³LS) approach [3–7] has recently been proposed and applied effectively in adaptive array processing. Based on the D³LS approach, methods [8, 9] developed a fast DOA estimation algorithm process, called the pseudo-covariance matrix technique, which estimated fast varying signals in two steps. The rough incidence angle ranges are first obtained using the bearing response; then, exact incidence angles are determined by combing the bearing response and the directional spectrum. The bearing response and the directional spectrum are obtained using a pseudo-covariance matrix, based on a single snapshot [8]. The advantage of the approach given in [8] is that the DOA can be estimated on a snapshot-by-snapshot basis, without the need of prior knowledge of the number of incoming signals and without the need of the covariance matrix from a block of the array data. Kim’s algorithm performs well for both uncorrelated and correlated signals coexisting in nonstationary environments. However, it requires the solution of a nonlinear-generalized-eigenvalue equation of a pseudo-covariance matrix, resulting in a high computational burden, and, using this method, the number of resolved signals is at most half of the number of the sensor array. In order to enhance the resolution of Kim’s method, an improved pseudo-covariance matrix scheme was proposed. It was Wen’s method [9]which is called forward-backward pseudo-covariance matrix. Wen’s method processes more signals and obtains a better resolution than the forward-only matrix. However, it uses singular-value decomposition (SVD), requiring a greater computational burden than does Kim’s method.

In this paper, an effective method is proposed for developing a fast DOA estimation algorithm when uncorrelated and correlated signals are present in highly nonstationary environments. The proposed method consists of two approaches. First, we use a suboptimal projection spectrum approach to select the rough DOA ranges of the incoming signals. This is done forming small spectrum peaks, generated by removing the signal from a given direction, instead of solving a more exact projection operator, as in Kim’s and Wen's methods. The whole process reduces algorithm computational cost. Secondly, the solution of a linear-matrix equation, formed by a forward-backward data matrix, is applied as a noise pseudo-eigenvector, which performs a function similar to a noise eigenvector. Then, a fine DOA estimation can be obtained from an eigenspectrum. This approach allows us to reduce computational cost and obtain a good algorithm resolution. In contrast to Kim’s and Wen's algorithms, the proposed algorithm has better or similar capability of resolution but offers lower computational complexity. Moreover, the proposed algorithm deals well with both uncorrelated and correlated signals and processes only a signal data snapshot without the use of the covariance matrix. Therefore, it may provide wider application prospects in real-time DOA estimation.

This paper is organized as follows. In Section 2, we briefly introduce the problem statement and signal model. In Section 3, we present the proposed DOA estimation algorithm. In Section 4, some simulation results and a comparison of the proposed algorithm to MUSIC, ESPRIT, Kim’s, and Wen’s methods are presented. Finally, conclusions are presented in Section 5.

2. Problem Statement and Signal Model

In the section, we prepare mathematical models and representations of signals. In what follows, the operators , , , and denote transpose, conjugate transpose, conjugate, and vector norm, respectively. represents a diagonal matrix with diagonal entries .

Our objective is to estimate the DOA from incident signals coming from a received single snapshot. Let a single data snapshot vector received from a uniform linear array with -omnidirectional sensors, representing a linear combination of narrow-band signals added with zero mean white Gaussian noise process, be given by where is a signal vector, is an noise process vector, and represents a matrix of steering vectors. The th column is the steering vector coming from , . Parameter is the distance between the sensors and is the wavelength of the carrier.

3. Fast DOA Estimation of Proposed Method

In this section, the proposed DOA estimation is derived, determining the rough and fine DOA estimation algorithms, presented in Sections 3.1 and 3.2, respectively.

3.1. Rough DOA Estimation Algorithm

In Kim’s and Wen’s methods, the rough incidence angle ranges are estimated by using the bearing response, requiring the use of nonlinear-generalized-eigenvalue or SVD algorithm in the different look directions. However, such processing may bring high computation complexity against real-time requirement. In the subsection, we present a simple novel DOA estimation algorithm, called the orthogonal projection technique, to estimate the rough DOA ranges from incoming signals using a single snapshot for achieving fast tracking varying signals under highly nonstationary environments. Let us define an orthogonal projection matrix as where denotes each search angle, is search angle step, and is the normalized version of , determined as . Then, use (2) to form the following projected power spectrum: When for , the projected signal is given by

Then, (4) is composed of all incident signals with noise added. However, when for ,

Then, (5) is composed of incident signals with noise added, where the steering vector is removed. The norm of (4) will have much greater value than the norm of (5). Thus, it is easy to show that we can obtain a peak in the projected power spectrum when the steering vector coincides with the source direction, that is, , . By using the projection matrix for each steering vector , the rough DOA incident signals ranges are estimated by searching local peaks from , .

3.2. Fine DOA Estimation Algorithm

In the conventional subspace-based methods, the covariance matrix is decomposed into its constituent eigenvectors and eigenvalues. Recalling that the eigenvectors of the noise eigenvalues are orthogonal to the signal subspace spanned by the eigenvectors of the signal eigenvalues, then MUSIC algorithm uses the orthogonality of these subspaces efficiently to get a high-resolution DOA estimation. However, such processing is time-consuming against real-time applications and may lose performance for correlated signals. This subsection presents a procedure of finding a noise pseudo-eigenvector based on the solution of a linear-matrix equation, formed by a forward-backward data matrix, for real-time requirement and for the purpose of working well in correlated signals.

By using the orthogonal projection method, we state that if has a minimum value at angle , the steering vector should be guaranteed to be too far from the signal directions. In order to acquire a noise pseudo-eigenvector, the steering vector may be used. First, we form two data matrices, , and , by dividing the uniform linear array into high overlapping subarrays of size and . Thus, the sensors form the first subarray of and the sensors form the second subarray of and so on; the sensors form the first subarray of and the sensors form the second subarray of and so on, that is, we rearrange the received single snapshot to form the following two data matrices: where is equal to , is equal to , and is the closest integer to . Finally, we combine the forward data matrix (6) and the backward data matrix (7) to form an forward-backward data matrix asWe know that the resolution capability of the DOA estimation is relative to the number of sensors or to the degree of freedom (DOF) of the solved data matrix [7, 9]. For both the forward-only and backward-only methods, the maximum number of the DOF we can consider is given by . In the proposed algorithm, we combine the forward data matrix and backward data matrix to double the received data and thereby increase the DOF significantly over that of either the forward-only or backward-only method alone. The maximum number of the DOF in the forward-backward method can be increased significantly to without increasing the number of sensors; thus, this is the reason that we set and . Then, the proposed algorithm with about DOF can provide a better resolution performance than Kim’s forward-only method with about DOF, and a resolution performance similar to Wen’s forward-backward method with about DOF. Let denote the th column vector of the data matrix (8). Then, we can write where and can be expressed as where , . Suppose that the first two signals are highly correlated or coherent among the incident signals, that is, , where denotes a complex coefficient describing the phase and gain relationship between the two coherent signals. Under this assumption, we can rewrite the first column vector of the data matrix (8) as then Based on (12), excludes the desired signal , that is, it only includes other desired signals and noise vector . If we can find a weight vector to satisfy , then (12) can be rewritten as Therefore, By the same process, we can then get

Combining (12-1) and (12-2) into matrix notation gives where then Note that (16) skips the scalar because it does not affect the vector direction of weight vector . Weight vector has pattern nulls at the incident angles of the rest of incoming signals, . Therefore, if the steering vector is used in (16), then the corresponding weight vector will be orthogonal to the signal steering vectors of the incoming signals, that is, , where . In this paper, weight vector is a noise pseudo-eigenvector because it performs a function similar to that of a noise eigenvector. Then, fine peaks corresponding to the DOA of the incoming signals can be found in eigenspectrum using Moreover, since an eigenspectrum provides higher resolution DOA estimation than a projection spectrum, steering vector with minimum eigenspectrum via (17) is more likely to be constrained too far from the signal steering vectors than steering vector . The corresponding weight vector is obtained by using in (16), which is more orthogonal to the signal steering vectors than the weight vector . Then, we can obtain a more accurate DOA estimation of the incoming signals than (17) by performing an eigenspectrum as follows: It is known that a linear matrix equation can be solved very efficiently by applying a conjugate gradient (CG) algorithm. The conjugate gradient method can provide a faster rate of convergence than other gradient methods and can be implemented to operate in real-time utilizing a digital signal processing chip [7, 10]. Therefore, the use of the CG algorithm makes our proposed algorithm suitable for real-time implementation. The procedure of our proposed fast DOA estimation is as follows.(1)Obtain the rough DOA ranges of the incoming signals and a steering vector by using the projection method via (3).(2)Use the steering vector in (16) to obtain a noise pseudo-eigenvector .(3)Obtain the fine DOA estimation of the incoming signals by searching for local peak values in the eigenspectrum via (17) within the rough DOA ranges.(4)Use the steering vector in (16) to obtain a noise pseudo-eigenvector .(5)Fine-tune the DOA estimation of the incoming signals in step (3) by searching for local peak values of the eigenspectrum via (18).

The computational cost of solving a nonlinear generalized eigenvalue algorithm is about flops for a matrix. The computational cost of solving an SVD algorithm is about flops for a matrix. The computational cost of solving the CG algorithm is about flops per iteration for a matrix. The size of resolved data matrix is approximately in Kim’s method, in Wen’s method and in our proposed algorithm. The main computational complexity difference between the other methods and the proposed algorithm is the weight vector calculation, which requires about from solving a nonlinear generalized eigenvalue algorithm in Kim’s method, from solving a SVD algorithm in Wen’s method, and in our proposed method, which is based on the CG algorithm. A flop is defined as a complex floating-point multiplication operation, where denotes “the order of ”, is the number of the sensors, denotes the total scanning number in the possible DOA range (i.e., the different looking directions), and is the number of iterations for executing the CG algorithm. We can see that the computational complexity of the DOA estimation for fast varying signals for our proposed method is significantly smaller than those that of Kim’s, and Wen’s methods because in general.

4. Computer Simulation

To compare the performance of the described DOA finding method, MUSIC, ESPRIT, Wen’s, and Kim’s method, in terms of the probability of resolution, the methods were investigated in the narrowband farfield DOA estimations with computer experiments. The probability of resolution is defined as follows. Two signals with DOA and are said to be resolved if the respective estimates and are such that both and are less than [11]. Four different experiments were done with variations in SNR value, the separation between sources, and the number of sensors.

In each experiment, it was assumed that the powers of all the incoming sources were equal. The noise at each array sensor was assumed to be additive white Gaussian noise process with zero mean. We also assumed a uniform linear array with half-wavelength element spacing. We examined the conventional MUSIC and ESPRIT algorithm when the number of snapshots was 100, whereas the proposed algorithm, Kim’s, and Wen’s methods only used a single snapshot. A total of 300 independent test runs were done to get each simulated point, and the DOA scanning was performed over [0°,90°] with a step size of . In these simulations, the solution of a linear matrix equation was solved using the CG method.

In the first experiment, three signals with the same SNR of 20 dB were impinging on the array of 23 sensors. Figures 1(b) and 1(c) show that we can estimate the three uncorrelated signals from [2.3°,9.7°,17.8°] using the proposed algorithm and MUSIC, respectively. Figure 2(b) shows that we can still estimate the two coherent signals from [2.3°,9.7°] and one uncorrelated signal from [17.8°] using the proposed algorithm but Figure 2(c) shows that MUSIC was unable to estimate the DOA of the two coherent signals. Also, Figures 1(a) and 2(a) show that the angles [2.3°,9.7°,17.8°] of the three incident signals were within the estimated possible ranges.

(a)

(b)

(c)

Figure 1: DOA estimation using the proposed algorithm and MUSIC for three uncorrelated signals from [2.3°,9.7°,17.8°] under the SNR of 20 dB and the number of sensors .

(a)

(b)

(c)

Figure 2: DOA estimation using the proposed algorithm and MUSIC for two coherent signals from [2.3°,9.7°] and one uncorrelated signal from [17.8°] under the SNR of 20 dB and the number of sensors .

In the second experiment, we considered three uncorrelated signals from [2.3°,9.7°,17.8°], impinging on the array of 23 sensors. Figure 3 shows the comparison result of the probability of resolution versus various SNR values. As seen, ESPRIT has best resolution, the proposed algorithm and Wen’s method have almost identical resolution, and Kim’s method has worst resolution. When an SNR = −5 dB, the resolution probability of the proposed algorithm, Kim’s and Wen’s methods was less than 50% while it is close to 60% for ESPRIT.

Figure 3: Comparison of results of the probability of resolution from ESPRIT, the proposed algorithm, Kim’s, and Wen’s methods, with for three uncorrelated signals from [2.3°,9.7°,17.8°] having different SNR, where the number of sensors .

In the third experiment, three uncorrelated signals from [2.3°,2.3° + δθ,2.3° + 2δθ], impinging on the array of 23 sensors and SNR of 20 dB, were considered where the δθ denote a small angle spacing. Figure 4 shows the comparison result of the probability of resolution versus the angle spacing of the incoming signals. As seen in that figure, ESPRIT achieves a higher rate of successful source separation in comparison with the proposed algorithm, Kim’s, and Wen’s methods. For an angle spacing , the resolution probability of the proposed algorithm, Kim’s and Wen’s methods was close to 0%. When an angle spacing , the Wen’s method and the proposed algorithm had better performance than Kim’s method.

Figure 4: Comparison of results from ESPRIT, the proposed algorithm, Kim’s, and Wen’s methods, for the probability of resolution versus the angle spacing of the incoming signals for three uncorrelated signals from [2.3°,9.7°,17.8°] having the SNR of 20 dB, where the number of sensors .

In the fourth experiment, we considered three uncorrelated signals from [2.3°,9.7°,17.8°], impinging on the array with SNR of 20 dB. Figure 5 shows the comparison result of the probability of resolution versus the number of sensors for ESPRIT, the proposed algorithm, Kim’s, and Wen’s methods. The simulations showed that the performance could be improved by using more number of sensors for all algorithms, but that more processing time is required.

Figure 5: Comparison of results of results of the probability of resolution from ESPRIT, the proposed algorithm, Kim’s, and Wen’s methods, with three uncorrelated signals from [2.3°,9.7°,17.8°] having the different number of sensors, where the SNR of 20 dB.

As expected, from Figure 1 to Figure 5, the simulations showed that the proposed method outperformed Kim’s methods and performed almost the same as Wen’s method with regards to resolution capabilities because the effective array size in both the proposed algorithm and Wen’s method is 15, but in Kim’s method is 11. If the number of signal sources is known exactly and has enough input data snapshots, of course, the MUSIC or ESPRIT provides us with better performance than the proposed algorithm in cases with uncorrelated signals.

5. Conclusion

A new fast direction of arrival estimation technique is presented based on both a simple projection spectrum and the eigenspectrum. The advantages of this method over MUSIC, ESPRIT, Kim’s, and Wen’s methods are as follows: (1) it has a lower computation loads; (2) it allows arbitrary signal statistics, for example, nonstationary and coherent; (3) it is capable of tracking high varying signals. Computer simulations have shown that the proposed algorithm handles the coherent or uncorrelated signals well, whereas MUSIC does not, as well as better resolution as compared with Kim’s method, or similar resolution to Wen's method. Since this new approach improves performance by reducing computational complexity while maintaining sufficient resolution in highly nonstationary environments, it will have a wider range of prospective applications in real-time DOA estimation than other comparable methods.

References

R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.

A. Paulraj, A. Roy, and T. Kailath, “Estimation of signal parameters via rotational invariance techniques—esprit,” in Proceedings of the 19th Asilomar Conference on Circuits, Systems and Computers, pp. 83–89, November 1985.

W. Choi, T. K. Sarkar, H. Wang, and E. L. Mokole, “Adaptive processing using real weights based on a direct data domain least squares approach,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 1, pp. 182–190, 2006.

K. Kim, T. K. Sarkar, H. Wang, and M. Salazar-Palma, “Direction of arrival estimation based on temporal and spatial processing using a direct data domain (D) approach,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 2, pp. 533–541, 2004.

K. Kim, T. K. Sarkar, and M. S. Palma, “Adaptive processing using a single snapshot for a nonuniformly spaced array in the presence of mutual coupling and near-field scatterers,” IEEE Transactions on Antennas and Propagation, vol. 50, no. 5, pp. 582–590, 2002.

T. K. Sarkar, J. Koh, R. Adve, R. A. Schneible, M. C. Wicks, S. Choi, and M. S. Palma, “A pragmatic approach to adaptive antennas,” IEEE Antennas and Propagation Magazine, vol. 42, no. 2, pp. 39–55, 2000.

T. K. Sarkar, et al., Smart Antennas, IEEE Press, New York, NY, USA; Wiley-Interscience, Hoboken, NJ, USA, 2003.

J. T. Kim, S. H. Moon, D. S. Han, and M. J. Cho, “Fast DOA estimation algorithm using pseudocovariance matrix,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 4, pp. 1346–1351, 2005.

Z. Wen, L.-P. Li, and P. Wei, “Fast direction finding using modified pseudocovariance matrix,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 12, pp. 3914–3918, 2006.

G. H. Golub and C. F. Van Loan, Matrix Computations, Johns Hopkins University Press, Baltimore, Md, USA, 1983.

R. Grover, D. A. Pados, and M. J. Medley, “Subspace direction finding with an auxiliary-vector basis,” IEEE Transactions on Signal Processing, vol. 55, no. 2, pp. 758–763, 2007.