Enhanced 2 D-DOA Estimation for Large Spacing Three-Parallel Uniform Linear Arrays

An enhanced two-dimensional direction of arrival (2D-DOA) estimation algorithm for large spacing three-parallel uniform linear arrays (ULAs) is proposed in this paper. Firstly, we use the propagator method (PM) to get the highly accurate but ambiguous estimation of directional cosine. Then, we use the relationship between the directional cosine to eliminate the ambiguity. This algorithm not only can make use of the elements of the three-parallel ULAs but also can utilize the connection between directional cosine to improve the estimation accuracy. Besides, it has satisfied estimation performance when the elevation angle is between 70° and 90° and it can automatically pair the estimated azimuth and elevation angles. Furthermore, it has low complexity without using any eigen value decomposition (EVD) or singular value decompostion (SVD) to the covariance matrix. Simulation results demonstrate the effectiveness of our proposed algorithm.


Introduction
Two-dimensional DOA estimation has attracted extensive attention due to its wide range of applications [1][2][3][4], and there are many algorithms been proposed for DOA estimation.Among these algorithms, the subspace-based algorithms, such as MUSIC and ESPRIT, have received a lot of attention for its accurate estimation performance [5,6].However, the complexity of subspace-based algorithms is often too large because of the EVD or SVD.In [7], Marcos et al. put forward the wellknown propagator method (PM) algorithm for 1D-DOA estimation.It uses linear partitioning instead of any EVD or SVD to reduce complexity.Then, Chen et al. [8] extend it to the DOA estimation of noncircular signal.References [9][10][11] extend the PM algorithm to 2D direction estimation.But these algorithms have some drawbacks.Reference [9] requires angle search operations.Reference [10] may fail in practical situation when elevation angle is between 70 °and 90 °.Reference [11] has worse estimation accuracy, because some element information is missing when calculating the propagation matrix (PMA).Based on this, Chen et al. [12] put forward an improved 2D angle estimation algorithm for threeparallel ULAs.It can solve all the problems mentioned in the above literature.But, it still has some shortcomings, that is, the array aperture is reduced because of the adoption of three parallel arrays with half wavelength spacing.And the estimation accuracy is reduced.References [13,14] establish the "array of subarrays" idea and "cyclic ambiguity" idea to improve the estimation accuracy of 2D-DOA.According to this, we expanded the spacing between the ULA and proposed an enhanced 2D-DOA estimation algorithm.What needs to be stressed here is that our method is different from the method in [13,14], although our algorithm and the method in [13,14] both use the large array spacing to obtain high accuracy but ambiguous estimation.However, the principle of resolving ambiguity is different.References [13,14] use coarse estimation without ambiguity or based on eigen space to resolve ambiguity.However, our algorithm uses the triangular relation between the three directional cosine to resolve the ambiguity.By doing this, we do not have to limit ourselves to using subarrays or eigen spaces to solve ambiguity problems.That is the main difference between our method and the method in [13,14].And it is also the main innovation of this paper.The flow of our algorithm is as follows: Firstly, we use the method in [12] to obtain the estimation of three directional cosines.Because of the large spacing between array elements, the estimated directional cosine is high precision and ambiguous.Then, we use the triangle relationship between the three directional cosines to eliminate the ambiguity.Then, we can get the true 2D-DOA of targets.Simulation results show that it cannot only avoid the problem in [8][9][10][11] simultaneously but also has better estimation accuracy than algorithm of [12] because of adopting large aperture.And the complexity of the algorithm is comparable to that algorithm of [12].
Notations.Superscripts ⋅ * , ⋅ T , ⋅ † , ⋅ −1 , and ⋅ H denote complex conjugation, transpose, pseudo-inverse, inverse, and conjugate transpose, respectively.I K is K × K identity matrix.diag ⋅ denotes the diagonalization of the entity inside.⋅ denotes take the absolute value of the element.⋅ denotes round the element to the nearest integers less than or equal to the element.⋅ denotes round the element to the nearest integers greater than or equal to the element.E ⋅ denotes the expectation operation.arg ⋅ and Re ⋅ denote the phase and the real part of a complex number separately.

Problem Formulation
As shown in Figure 1, assume that there are three-parallel ULAs, namely, X,Y, and Z. Array X contains N + 1 sensors.Array Y and array Z have N sensors, respectively.The spacing between adjacent elements is d x The distance between array X and array Y is d y , and the distance between array X and array Z is d z λ is the wavelength.It is assumed that P far-field narrowband uncorrelated signals are incident onto the array.The elevation angle and azimuth angle of ith target are θ i and φ i , respectively.Here, we assume that the range of 2D angle is the same with that of [12].That is to say that the range value of θ i is 0, 1/2π and the range value of φ i is −1/2π, 1/2π Then, the output of the three ULAs can be expressed as follows: and Ω z = diag e −j2πd z sin θ 1 cos φ 1 /λ , e −j2πd z sin θ 2 cos φ 2 /λ , … , e −j2πd z sin θ P cos φ P /λ W x t ,W y t , and W z t are assumed to be Gaussian white noise vectors whose mean value is zero and variance is σ 2 Then a new vector W t = X t T Y t T Z t T T , with L snapshots, and where The difference between our model and the model in [12] is that the spacing between the ULA is larger than half wavelength.The large spacing will lead to ambiguity of the directional cosine, but it would not affect the estimation method for them.So the first several steps of our algorithm is same with the algorithm of [12].We abbreviate the first several procedures as follows:

Enhanced 2D-DOA Estimation
Step 1. Compute the covariance matrix of W t by where R W1 ∈ C 3N+1 ×P and R W2 ∈ C 3N+1 × 3N+1−P Compute the propagator matrix (PMA) Until now, we get the paired estimation of αi , βi , and γi which represent cos θ i ,sin θ i cos φ i , and sin θ i sin φ i , respectively.Because d z is larger than half wavelength, so the βi is high accuracy but ambiguous.And all the ambiguous values are listed as follows: Highly Accurate Estimation of 2D-DOA.Before estimating the 2D-DOA, we need to eliminate the ambiguity first.Note that there exists a triangle relationship between these directional cosines: 2 International Journal of Antennas and Propagation And the true directional cosine of targets satisfies this equation.So, we can use this equation to select the true directional cosine from the ambiguous directional cosine.
But because there are square terms in the constraint, the mirror image of the true value also satisfies this equation.If the mirror image is exactly equal to an ambiguous value of the true value, we may make mistakes and treat the mirror image as the true value.Here, we set the range value of θ and φ as 0, 1/2π and −1/2π, 1/2π , separately.So the true value of sin θ i cos φ i is limited in the range of 0, 1 And the mirror image can be avoided.
And the procedure of eliminating ambiguity can be got as follows: Step 1. Get all the ambiguous value of sin θ i cos φ i by Step 2. Find the true value by Lastly, we can get the 2D-DOA estimation for each source by Here, the analysis of the complexity, estimation performance, and some notices of our algorithm are as follows: (1) Our algorithm has an additional ambiguity elimination operation compared with the algorithm of [12].And the complexity of ambiguity elimination is far less than the complexity of constructing the covariance.Thus, we can say that the algorithm has comparable computational complexity with the algorithm of [12].
(2) Because our array aperture is larger than that in [12], so our proposed algorithm can obtain a better estimation accuracy.
(3) Note that here we expand the spacing between the X-axis and Z-axis.Similarly, we can also expand the spacing between the X-axis and Y-axis or simultaneously expand the spacing in the X-axis, Y-axis, and Z-axis.Then, we can get high accuracy but ambiguous estimation of sin θ i sin φ i and cos θ i , respectively.And we can still use (4) to eliminate the ambiguity.But considering the mirror image problem, the range of angle needs to be adjusted accordingly.
(4) Because the first several steps of our algorithm are the same with those of [12], so our algorithm inherits the merits of algorithm of [12].For example, (1) it can take full advantage of the elements of the threeparallel ULAs to estimate the PMA; (2) it has satisfied estimation accuracy in actual mobile elevation angles; (3) it can automatically pair the estimated 2D-DOA; (4) it has low complexity without using any EVD or SVD.

Computer Simulations
Here, we will conduct several simulation experiments to test the performance of our proposed method.The 2D-DOA of two uncorrelated targets are θ 1 , φ 1 = 35 °, 40 °and θ 2 , φ 2 = 80 °, 75 °, separately.N is set to 6.So the total number of array elements is 19.The spacing between X-axis and Z-axis is set to d z = 3d The other spacing is set to

International Journal of Antennas and Propagation
In the first simulation, Figure 2 shows the scattergram of 300 times 2D-DOA estimate results by our algorithm and the algorithm of [12].The snapshot is set to 50.The SNR = 10 dB.As we can see, both the two algorithms can clearly observe the two targets.But our algorithm has better statistical properties.That is because we make better use of the information between ULAs.
In the second simulation, we compare the estimation performance of the two algorithms under different SNR which vary from 0 dB to 30 dB.The snapshot is set to 300.200 times Monte Carlo simulations are conducted under each SNR.As illustrated in Figure 3, the RMSE curve of azimuth estimated by our algorithm is better than that estimated by the algorithm of [12].The elevation estimation performance of our   4 International Journal of Antennas and Propagation algorithm is comparable with that estimated by the algorithm of [12].The reason is that when estimating azimuth angle, the equivalent array aperture of our algorithm is larger than that of [12].While when estimating elevation angle, the equivalent array aperture of our algorithm is equal with that of [12].
Note that if we want to improve the estimation accuracy of elevation angle, we can expand d y instead of expandingd z In the third simulation, we compare the estimation performance of the two algorithms under different snapshots which vary from 100 to 2100.The SNR is set to 20 dB.1000 times Monte Carlo simulations are conducted under each snapshot.As illustrated in Figure 4, the RMSE curve of azimuth estimated by our algorithm is better than that estimated by the algorithm of [12].The elevation estimation performance of the two algorithms is similar.The reason is the same with the previous one.From the second and third simulations, we can also see that our algorithm has satisfied estimation performance when elevation is between 70 °and 90 °.In the fourth simulation, we compare the complexity of our algorithm and that of [12].We choose the runtime as the evaluation criteria.Two kinds of comparison were conducted.In the first comparison, the number of elements in subarray, that is, N, is set to vary from 10 to 100 and the snapshot is set to 300.The result is shown in Figure 5(a).In the second comparison, the snapshot is set to vary from 100 to 2100.The number of elements in subarray is set to 20.And the result is shown in Figure 5(b).From the two figures, we can see that the runtime of our algorithm is similar with that of [12] which has shown the efficiency of our proposed algorithm.
In the fifth simulation, we compare the azimuth estimation performance of the two algorithms with different separation of angle.The 2D-DOA of two targets are θ 1 , φ 1 = 30 °, 60 °+ Δφ 1 and θ 2 , φ 2 = 50 °, 70 °+ Δφ 2 , respectively.Both Δφ 1 and Δφ 2 vary from 0 °to 18 °.The SNR is equal to 20 dB.The snapshot is set to 500.200 times Monte Carlo simulations are conducted under each separation of angle.Here, the "angle separation" denotes Δφ 1 and Δφ 2 which vary from 0 °to 18 °.It is different from the common definition of "angle separation."In fact, the separation angle of the two targets is fixed, what changes is the azimuth angle of the two targets.The purpose of doing this is to make the azimuth angle traverse all angles between 70 °and 88 °to prove the effectiveness of the algorithm in practical application.The result is shown in Figure 6.We can see that the performance of our proposed algorithm is better than that of the algorithm of [12].
In the sixth simulation, we compare the elevation estimation performance of the two algorithms with different separation of angle.Here, we set d x = d z = d and d y = 3d The 2D-DOA of the two targets are θ 1 , φ 1 = 60 °+ Δθ 1 , 30 °and θ 2 , φ 2 = 70 °+ Δθ 2 , 50 °, respectively.Both Δθ 1 and Δθ 2 vary from 0 °to 18 °The other parameters are the same with the fifth simulation.The result is shown in Figure 7.We can see that the performance of our proposed algorithm is better than that of the algorithm of [12].

Conclusions
An enhanced 2D-DOA estimation algorithm with threeparallel large spacing ULAs is proposed in this paper.It cannot only make use of the elements of the three-parallel ULAs to estimate the PMA but it can also utilize the information between three-parallel ULAs to improve the estimation accuracy.Besides, it has satisfied estimation performance when the elevation angle is between 70 °and 90 °and it can automatically pair the estimated 2D-DOA.Furthermore, it has low complexity without using any EVD or SVD to the covariance matrix.Simulations proved the effectiveness of the proposed algorithm.
Target 1 with algorithm of [12] Target 2 with algorithm of [12] Target 1 with proposed method Target 2 with proposed method

3. 1 .
Highly Accurate but Ambiguous Estimation of Directional Cosine.Assuming that d x = d y = d = λ/2,d z = Kd, where K is a positive integer and K > 1

Step 3 .T A 1 ′
W2and get the extended PMA by P e = I H P P Implement EVD onΨ z = P + x1 P z = A 1 Ω z A 1 −1to get eigenvalues βi and eigenvectors A 1 ′ They correspond to the diagonal elements of Ω z and the estimation of A 1 , separately.And P + x1 is the first N rows of P x Step 4. Let B 1 = P x1 A 1 ′, B 2 = P y A 1 ′ Attain αi by extracting the ith diagonal elements of B + 1 B 2 Similarly, let C 1 = P T x1 P T y1 P T z1 T A 1 ′ and C 2 = P T x2 P T y2 P T z2 Attain γi by extracting the ith diagonal elements of C + 1 C 2 And P y1 and P z1 are the first N − 1 rows of P y and P z P y2 and P z2 are the last N − 1 rows of P y and P z P x2 represents the last N rows of P x

Figure 1 :
Figure 1: Sketch map of array structure.

Figure 3 :
Figure 3: RMSE contrast of two algorithms with different SNR.
RMSE of elevation estimation

Figure 4 :
Figure 4: RMSE contrast of two algorithms with different snapshots.
Runtime versus number of snapshots

Figure 5 :
Figure 5: Runtime comparison of two algorithms.