Abstract
We introduce an iterative least squares method (ILS) for estimating the 2D-DOA and frequency based on L-shaped array. The ILS iteratively finds direction matrix and delay matrix, then 2D-DOA and frequency can be obtained by the least squares method. Without spectral peak searching and pairing, this algorithm works well and pairs the parameters automatically. Moreover, our algorithm has better performance than conventional ESPRIT algorithm and propagator method. The useful behavior of the proposed algorithm is verified by simulations.
1. Introduction
Antenna arrays have been used in many fields such as radar, sonar, and mobile communications, and so forth, [1–6]. The direction of arrival and frequency estimation of signals impinging on an array of sensors have received considerable attention in the field of array signal processing. For example, these parameters can be applied to locate the mobiles and allocate pilot tones in space division multiple access (SDMA) systems. Furthermore, a precise estimation of these parameters is helpful to attain a better channel estimate and thus enhances the system performance [7]. Uniform linear arrays for estimation of wave arrival have been studied extensively, and they contain maximum likelihood (ML) [8], multiple signal classification (MUSIC) algorithm [9, 10], estimation of signal parameters via rotational invariance techniques (ESPRIT) [11, 12], propagator method (PM) [13], and so forth.
The ML method is often applicable but might be computationally prohibitive. ESPRIT and MUSIC algorithms are based on signal subspace and have better parameter estimation performance. The main advantages of MUSIC/ESPRIT are the high-resolution estimates of direction of arrivals (DOAs) and frequencies, while the computational effort compared to ML method is significantly reduced. MUSIC requires multiple dimensional spectral peak searching, and it is the search that is still computationally expensive. The primary computational advantage of ESPRIT is that it eliminates the search procedure inherent. ESPRIT method requires eigenvalue decomposition (EVD) to the cross spectral matrix or singular value decomposition (SVD) to the received data. Reference [13] presented propagator method to estimate the angle and frequency with uniform linear array. Propagator method has low complexity, but its parameters’ estimation performance is less than ESPRIT algorithm. Compared with uniform linear array, L-shaped array can identify 2D-DOA and is very close to actual situation [14, 15]. And thus we propose a novel iterative-based angle and frequency estimation algorithm with L-shaped array which can achieve better performance than ESPRIT [11] and propagator method [13]. The proposed algorithm can obtain automatically paired 2D-DOA and frequency estimation. This method is an iterative algorithm, which does not need EVD or SVD, and only requires fewer iterations for convergence. The useful behavior of the proposed algorithm is verified by simulations.
The remainder of this paper is structured as follows. Section 2 develops data model. Section 3 deals with algorithmic issues. Section 4 presents simulation results, and Section 5 summarizes our conclusions.
Denote
We denote by the complex conjugation, by the matrix transpose, and by the matrix conjugate transpose. is to construct a diagonal matrix.
2. Data Model
We consider an L-shaped array with sensors at different locations as shown in Figure 1. A uniform linear array containing elements is located in -axis, and the other uniform linear array containing elements is located in -axis. We suppose that there are narrowband uncorrelated signal sources impinge on the L-shaped array with , where are the elevation angle and the azimuth angle of the kth source, respectively.
The received signal of elements in -axis is shown where is the source matrix, is the received additive white Gaussian noise, and , is the -axis direction matrix, and ,.
The received signal of elements in -axis is denoted as where is the -axis direction matrix and , , , is the received white Gaussian noise.
The received signal of the L-shaped array antennas can be denoted as where the direction matrix and noise matrix are shown as follows: and are Vandermonde matrices. The delayed signal for (3) with can be denoted as where .
We assume that channel state information is constant for symbols. The received signal of array antennas can be denoted as where .
Define the delay matrix as Equation (6) is also denoted as where is to extract the mth row of its matrix argument and constructs a diagonal matrix out of it.
3. Iterative Least Squares Method
It can be seen from (8) that if the direction matrices and are obtained, then the estimates of 2D-DOA and frequency follow immediately based on L-shaped array geometries. Hence our main focus in this section is on estimating and .
3.1. ILS Cost Function
The cost function can be constructed via the least squares criterion and given by where denotes the Frobenius norm. Interestingly, if any two subsets of and are fixed, the remaining parameter subset can be easily obtained by minimizing (9).
Firstly, fix and . Recall that , , and , where . Then (9) can be rewritten as Let the gradient of with respect to , be equal to zero, one obtains for and . The above equation yields Let , and , where denotes the Hadamard product. Then (12) can be rewritten in matrix form Hence, the estimate of can be easily computed as Secondly, let and be fixed. Differentiation of (9) with respect to yields the estimate of Thirdly fix , and . Let the gradient of (9) with respect to be equal to zero, we have which gives the estimate of Based on the above analysis, the ILS algorithm is summarized as follows.
Let the initial values and , the elements of and are equal to 1. For iteration index , do the following three steps.
Step 1. Calculate in (14) using and , for .
Step 2. Calculate in (15) using and .
Step 3. Calculate in (17) using and .
with a threshold ; if , then stop; otherwise, go to Step 1.
The basic idea of the iterative least square algorithm is to minimize the cost function with respect to , , and by using the cyclic minimization (CM) technique [16], which monotonically decreases the cost function. Furthermore, we know . Hence, the proposed algorithm is convergent, only 20 iterations are required to achieve convergence for this iterative algorithm (also shown in Figure 2).
We obtain the estimated matrices , , and , where is a permutation matrix and stand for diagonal scaling matrices satisfying . That is to say, the th column of corresponds to th column of and . So our algorithm can estimate 2D-DOA and frequency estimation without extra pairing.
3.2. Frequency and 2D-DOA Estimation
We use iterative least squares method to attain the direction matrix and the delay matrix , and then angle and frequency are estimated according to the least squares principle and Vandermonde characteristic. Define the delay vector as . is the th column of the delay matrix. We get , where angle(·) is to get phase part of a complex number. The estimated frequency vector (the th column of the estimated delay matrix ) is processed through normalization, which also resolves the scale ambiguity, and then normalized sequence is processed to attain according to above processing. Finally, we use leaset squares principle to estimate. Least squares fitting is , where The least squares solution for is Since , and the delay matrix are with Vandermonde characteristic, the direction matrix has a quasi-Vandermonde structure, so we can use the same method with frequency estimation method to estimate DOAs. The estimated receive array steer vector (the kth column of the estimated matrix ) is processed through normalization, which also resolves the scale ambiguity, then direction estimation matrix is divided into two submatrices, the first rows form submatrices .The last rows form sub-matrix , the kth column of the sub-matrix and is and , respectively. Then normalized sequence is processed to attain where , . Least squares fitting is , where where is the estimated frequency through (19) and is the estimated value of . The LS solution to is Similarly, , where in which is the estimated value of . The LS solution to is The 2D-DOAs are estimated via
3.3. Complexity Analysis and CRB
In contrast to ESPRIT algorithms in [11] and PM algorithm in [13], our algorithm has a heavier computational load, in our method, computing in (14) takes a complexity of , the complexity for computing in (15) is ; the complexity for computing in (17) is . Therefore, the complexity per iteration in the ILS algorithm is ; when the ILS algorithm converge within 20 steps, the total computational complexity of the proposed algorithm is about , and ESPRIT algorithm in [11] needs , and PM algorithm in [13] needs .
According to [17], we derive the CRB for 2D-DOA and frequency estimation with L-shaped array where stands for Hadamard product. , , , , , .
The advantages of the proposed algorithm can be summarized as follows. (1) The 2D-DOA and frequency can be paired automatically. (2) The proposed algorithm has better angle and frequency estimation performance than ESPRIT algorithm and PM algorithm.
4. Simulation Results
We present Monte Carlo simulations that are to assess joint 2D-DOA and frequency estimation performance of the proposed algorithm. The number of Monte Carlo trials is 500. There are three signals impinging on L-shaped array at (, 1.2 MHz), (, 1.6 MHz), (, 1.8 MHz), respectively. We consider an L-shaped array with 15 sensors, which is shown in Figure 1. The spacing between the adjacent elements in each uniform linear array is smaller than the half smallest wavelength of the incoming signals. Define root mean squared error (RMSE) of angle: where is the estimated elevation angle, is the perfect elevation angle, is the estimated azimuth angle, and is the perfect azimuth angle.
Define the root mean squared error: (RMSE) of frequency where is the estimated frequency and is the perfect frequency. Note that is the number of snapshots; is the number of sources; is the number of delay outputs for received signal of array antennas; is the number of antennas.
Simulation 1. We first investigate the convergence performance of the proposed algorithm. Define , where is the estimated matrix of th iteration. Figure 2 presents the algorithm convergence performance of the proposed algorithm with , and SNR = 20 dB. Figure 2 shows that the proposed algorithm needs 20 iterations or so to achieve convergence.
Simulation 2. The performance of our proposed algorithm is investigated. , and is in this simulation. Figure 3 shows 2D-DOA estimation of our proposed algorithm at SNR = 20 dB, and Figure 4 presents elevation angle and frequency scatter of our proposed algorithm at SNR = 20 dB. From Figures 3 and 4 we find that our proposed algorithm works well.
Simulation 3. We compare our proposed algorithm with ESPRIT algorithm, MUSIC algorithm, propagator method and CRB. The simulation parameters are retained as Simulation 1. From Figures 5 and 6 we find that our algorithm has much better angle and frequency estimation performance than ESPRIT algorithm and PM algorithm, and it has a very close angle and frequency estimation performance to MUSIC algorithm, but MUSIC algorithm needs three-dimensional spectral peak searching and the complexity is much larger than our proposed algorithm.
Simulation 4. Our proposed algorithm performance with and different values of is investigated. is set to 13, 15, and 17 in this simulation. It is indicated from Figures 7 and 8 that 2D-DOA estimation performance of our algorithm is improved with the number of antennas increasing. When the number of antennas increases, our algorithm has higher receive diversity.
Simulation 5. The performance of our proposed algorithm with and different values of is investigated. is set to 2, 3, 4 in this simulation. It is indicated from Figures 9 and 10 that 2D-DOA and frequency estimation performance of our algorithm degrade with the increasing of the source .
5. Conclusions
In this paper, we develop a novel method for joint 2D-DOA and frequency estimation based on L-shaped array using iterative least squares technique. Without spectral peak searching and pairing, this algorithm works well. Furthermore, our algorithm has much better 2D-DOA and frequency estimation performance than conventional ESPRIT algorithm and PM algorithm, and it has a very close 2D-DOA and frequency estimation performance to MUSIC algorithm. The useful behavior of the proposed algorithm is verified by simulations.
Acknowledgments
This work is supported by National Nature Science Foundation of China (nos. 61179006, 60801052), Jiangsu Planned Projects for Postdoctoral Research Funds (no. 1201039C), Open project of key laboratory of underwater acoustic communication and marine information technology (Xiamen University) and Nanjing University of Aeronautics and Astronautics Research Funding (nos. NP2011036, NZ2012010, kfjj120115, kfjj20110215). The authors wish to thank the anonymous reviewers for their valuable suggestions on improving this paper.