Abstract

In this paper, an improved propagator method (PM) is proposed by using a two-parallel array consisting of two uniform large-spacing linear arrays. Because of the increase of element spacing, the mutual coupling between two sensors can be reduced. Firstly, two matrices containing elevation angle information are obtained by PM. Then, by performing EVD of the product of the two matrices, the elevation angles of incident signals can be estimated without direction ambiguity. At last, the matrix product is used again to obtain the estimations of azimuth angles. Compared with the existed PM algorithms based on conventional uniform two-parallel linear array, the proposed PM algorithm based on the large-spacing linear arrays has higher estimation precision. Many simulation experiments are presented to verify the effect of proposed scheme in reducing the mutual coupling and improving estimation precision.

1. Introduction

Estimating the directions of arrival (DOA) of spatial signals by sensors array has widespread application in wireless communication [1] and multiple input multiple output (MIMO) radar [2]. Two-dimensional (2D) DOA estimation technology can obtain more angle information of spatial signals, and it has higher practical significance than one-dimensional (1D) DOA estimation technology. L-shaped array and parallel array are frequently used array constructions for 2D DOA estimation. The L-shaped array consists of two orthogonal linear arrays, based on which many effective 2D DOA estimation algorithms [3, 4] were proposed.

Being different from the L-shaped array, the parallel array consists of multiple parallel linear arrays and has more diverse array geometries. Many 2D DOA estimation algorithms [510] were designed based on parallel arrays. In [5], a DOA matrix algorithm was proposed to obtain the estimations of the elevation and azimuth angles by using a two-parallel linear array. In [6], the polynomial rooting technology was used for 2D DOA estimation based on a two-parallel linear array. In [7], a rank-reduction algorithm based on a three-parallel linear array has presented for 2D DOA estimation. However, the three algorithms need to perform eigenvalue decomposition (EVD) or singular value decomposition (SVD) of covariance matrix. The propagator method (PM) [813] gains extensive attention for the lower computational complexity because it does not need to perform EVD of covariance matrix. In [8], authors used PM algorithm to estimate 2D DOA by multiple parallel linear arrays. In [9], a modified 2D PM algorithm based on two-parallel linear array was proposed for improving estimation precision. In [10], another modified 2D PM algorithm based on three-parallel linear array was proposed for reducing computational complexity. In [11], PM algorithm was used with the parallel factor analysis (PARAFAC) model to estimate the 2D DOA, where the PARAFAC model is solved by circulative iteration.

But we should notice that all the algorithms [511] are based on uniform array, the interval of adjacent sensors cannot exceed half-wavelength of the incident signal. In fact, mutual coupling inevitably exists between two sensors, particularly for the closed two sensors. Aiming to extend array aperture and reduce mutual coupling, some nonuniform sparse array constructions such as nested array [1417] and coprime array [1820] were proposed. Nested array consists of multiple subarrays with different intervals. Coprime array consists of two uniform linear arrays, and the intervals of two subarrays are relevant to two mutual prime integers. Compared with the uniform array, nonuniform sparse arrays have extendable array aperture and higher degree of freedom. For the sparse arrays [1419], the interval of adjacent sensors can exceed half a wavelength of the incident signal which can reduce the mutual coupling. Despite all these, in order to avoid direction ambiguity, the minimum element spacing of these sparse arrays still cannot exceed half-wavelength of the incident signal. Hence, it is difficult to eliminate mutual coupling completely. In [20], authors proposed an off-grid DOA estimation algorithm based on an unfolded parallel coprime array. The mutual coupling of array can be reduced to a great extent. In addition, for nonuniform sparse arrays, constructing extended covariance matrix is a universal method to exploit the potential virtual sensors. But the computational complexity of the DOA estimation algorithms is closely related to the order of covariance matrix. This method can improve performance of angle estimation, but also can increase computational complexity.

In this paper, we propose a two-parallel linear array consisting of two uniform large-spacing linear arrays. The element spacing is q units for one linear array and q + 1 units for the other linear array, where q is not smaller than 2. Hence, the mutual coupling in one linear array can be eliminated when q is selected properly. For avoiding direction ambiguity, an improved 2D PM algorithm also is proposed to estimate the elevation and azimuth angles. Compared with the PM [810], the proposed algorithm has two obvious advantages: (1) the proposed algorithm can reduce the mutual coupling considerably due to the use of large-spacing linear array; (2) the proposed algorithm has higher estimation precision than PM [810], even if the mutual coupling between sensors is ignored.

2. Array Received Model

Consider that two uniform large-spacing parallel arrays are located on the plane, as shown in Figure 1 The coordinates of the N + 1 sensors on axis are , , where , is the wavelength of incident signal, and is any integer no less than 2. The coordinates of the N sensors for the other subarray are , , …, , . Figure 2 shows the construction of a 17-element large-spacing two-parallel linear array with q = 2, and Figure 3 shows the construction of a traditional 17-element two-parallel linear array [8, 9].

If we assume that mutual coupling exists between two sensors with the interval no further than , we can see that mutual coupling only exists in three pairs of sensors from Figure 2. Obviously, from Figure 3, we can find that the coupling effect exists between any adjacent sensors for the conventional two-parallel linear array.

Suppose that far-field, uncorrelated, narrowband signals received by the proposed array. We define the elevation angle and azimuth angle of the kth signal as and , respectively. Let received vectors by the two large-spacing linear arrays be , and , respectively. If we ignore the effect of mutual coupling, and can be written as [6, 10]where , , , , , , , , and are white Gaussian noise vectors received by the two linear arrays.

Let , then we have

If we assume that mutual coupling exists in array, can be expressed as [19, 20]where is the mutual coupling matrix and is the mutual coupling coefficient between the ith sensor and the jth sensor. Here, we let the ith sensor be the ith sensor of the first array for and be the (i − N−1)th sensor of the second array for .

The total mutual coupling of an array can be reflected in coupling leakage [17] aswhere with , and is F-norm.

According to (4), it is easy to know that the mutual coupling of proposed array is smaller than the traditional parallel array [811] and larger than unfolded parallel coprime array [20]. But just for the structural feature, the unfolded parallel coprime array cannot be used directly in traditional algorithms such as PM algorithm and estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm [21].

3. Description of Improved PM Algorithm

In this section, we introduce the improved PM algorithm based on the proposed array and the algorithm description is based on (2).

Denote the manifold matrix of the whole array as

There must be a propagator matrix satisfying [810]where consists of the first K rows of matrix and consists of the last 2N + 1 − K rows of matrix .

Denote the covariance matrix , and it can be partitioned as , where and .

Then, can be estimated by , and can be estimated by as [810], where T is the number of snapshots.

Construct a block matrix :where is the first N + 1 rows of and is the last N rows of .

According to (7), we can obtain

Combining (7) with (8), we can know

Denote as the first N rows of matrix , as the last N rows of matrix , as the first N − 1 rows of matrix , and as the last N − 1 rows of matrix , then we havewhere is the first N rows of , is the last N rows of , is the first N − 1 rows of and is the last N − 1 rows of .

Denote the matrix , , and , then we can know

According to formula (11), we can get

According to formula (12), we have

Denote two vectors and , where the (s + 1)th component of or is 1, and all the other elements are zeros. Then, we can construct a matrix and a matrix .

Using formulae (13) and (14), we can obtain two matriceswhere , , and .

Denote a matrixwhere consists of the first q(N + 1) rows of .

Combining (15) with (16), we have

Performing EVD of and can obtain the eigenvalues of the two matrices. Let the eigenvalues of and be and , then we can get the estimation of , as

Using the pairing method [8], we can get the matched 2D DOA estimations.

According to (14), we can know that the row number of and are larger than and . Hence, the construction of and also can see the process of adding the virtual sensors of array. For traditional nonuniform array, spatial smoothing algorithm [14] is the widely used method to obtain extended covariance matrix. And DOA is estimated by dealing with the extended covariance matrix in some existing algorithm. Certainly, this process is different from the proposed algorithm.

4. Complexity Analysis

In this section, the complexity of proposed algorithm is compared with PM [810]. Suppose that the two-parallel array consisting of an (L + 1)-element linear array and an L-element linear array for the proposed algorithm and PM [8, 9]. Consider an (L + 1)-element linear array and two L/2-element linear arrays for PM [10]. Assume that , so we only consider the complexity of main step for four PM algorithms. Both the complexity of PM [9] and proposed PM are . The complexity of PM [8] and PM [10] is and , respectively. Figure 4 lists the comparison of complexity versus snapshots with L = 10. Figure 5 lists the comparison of complexity versus the number of sensors with T = 500.

5. Simulation

In this section, we perform several groups of simulation experiments to confirm the outstanding effectiveness of proposed algorithm in reducing coupling and improving estimation precision. In [9], authors have proved that the accuracy of PM [8] is lower than PM [9], and complexity of PM [8] is higher than PM [9]. Hence, we only compare the proposed algorithm with PM [9, 10]. In the first three experiments, assume that a 21-element two-parallel array is used for the PM [9] and the proposed PM algorithm, and a 21-element three-parallel array are used for PM [10]. Suppose that the elevation and azimuth angles of three signals are and . Denote and as the estimations of and in the jth experiment, respectively. For the sake of fair comparison, only accurate pairing estimation results are given in the experiment. We denote the root mean square error (RMSE) of 2D DOA estimation aswhere J = 1000 is the number of repeated experiments.

Firstly, we compare the estimation precision of three PM algorithms in the absence of mutual coupling. Figure 6 lists the comparison result of RMSE for three PM algorithms and CRLB [9] versus SNR with 500 snapshots. Figure 7 lists the comparison result of RMSE for three PM algorithms and CRLB versus the number of snapshots with 5 dB SNR. According to the results shown in the two figures, we can see that the estimation precision of proposed PM and PM [9] is far higher than the PM [10], and the estimation precision of proposed PM is slightly higher than PM [9] under high SNR. But the proposed PM algorithm shows significant advantage when the SNR is lower, which can be seen from Figure 6.

Secondly, we compare the performance of the proposed method for different q in the absence of mutual coupling. Figure 8 lists the comparison result of RMSE versus SNR with 500 snapshots. Figure 9 lists the comparison result of RMSE for different q versus the number of snapshots with 5 dB SNR. The results in the two figures indicate that the performance of the proposed method can be improved slightly as the growth of q. Combining the results of Figures 8 and 9, we can know that the estimation precision can be improved as the extension of aperture under the same number of sensors. Since the row number of , is larger than , , it is reasonable to find that the proposed algorithm can improve the estimation precision.

Thirdly, we compare the estimation performance of the PM algorithm [9] and proposed algorithm in the appearance of mutual coupling. We only suppose that the mutual coupling effect only exist between two sensors with the interval no further than . Let be the mutual coupling coefficient between two sensors with interval d, and be the self-couple coefficient. Fix SNR at 15 dB and the number of snapshots at 500. Figures 10 and 11 display 100 estimation results of proposed algorithm and PM algorithm [9], respectively. It is clear to see that the estimation error of proposed PM is smaller than the PM [9]. The results also can prove the effectiveness of the proposed array in reducing mutual coupling.

At last, we test the estimation performance of proposed algorithm for more signals with small interval.

Assume that the proposed array consists of a 25-element array and a 26-element array. Suppose that the elevation and azimuth angles of fifteen signals are and . Fix SNR at 20 dB and the number of snapshots at 1000. Figure 12 displays 100 estimation results of the proposed algorithm. The result in this figure can show that the proposed algorithm can distinguish many signals with small interval, but it needs a large number of sensors.

6. Conclusion

In this paper, we have proposed an improved 2D PM algorithm with a two-parallel array consisting of two uniform linear arrays, the spacing of which is larger than half a wavelength of the incident signal. The proposed algorithm has higher estimation precision than many existed 2D PM algorithms. Due to the large-spacing of proposed array, mutual coupling can be reduced. The proposed algorithm has better performance than the traditional algorithm under mutual coupling. But because no decoupling algorithm is used for the proposed algorithm, the error of the DOA estimate is still significant under mutual coupling. In future research, according to the proposed two-parallel array, we will work for presenting the decoupling algorithm like [22, 23].

Data Availability

In this paper, we use simulation to test the performance of proposed algorithm. Hence, all data are generated by simulation. The simulation data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (51877015 and U1831117), the Cooperation Agreement Foundation by the Department of Science and Technology of Guizhou Province of China (LH[2017]7320 and LH[2017]7321), the Innovation Group Major Research Program Funded by Guizhou Provincial Education Department (KY[2016]051), the Foundation of Top-notch Talents by Education Department of Guizhou Province of China (KY[2018]075), and PhD Research Startup Foundation of Tongren University (trxyDH1710).