Two-Dimensional DOA Estimation of MIMO Radar Coherent Source Based on Toeplitz Matrix Set Reconstruction

In order to realize the high-precision direction of arrival (DOA) estimation of the coherent source of two-dimensional multipleinput and multiple-output (MIMO) radar, a solution is given by combining Toeplitz matrix set reconstruction. MIMO radar obtains a larger aperture with fewer arrays. Traditional two-dimensional reconstruction Toeplitz-like algorithms use part of the information in the construction of two correlation matrices or covariance matrices to construct the Toeplitz matrix when performing two-dimensional coherent source DOA estimation, which makes the information utilization incomplete and requires additional denoising processing. To solve the above problems, this paper proposes an improved Toeplitz matrix set reconstruction algorithm based on the two-dimensional reconstruction Toeplitz class algorithm. -e complete array element receiving signal vector is used to construct two Toeplitz matrix sets containing complete information, and then their conjugate transposes. Multiply and sum to correct the matrix to obtain a full-rank matrix, so as to achieve the purpose of decoherence and combine the traditional ESPRIT algorithm to perform two one-dimensional reconstruction processing through rotation invariance and then perform angle matching to achieve two-dimensional coherent signal angle estimation, while avoiding additional denoising processing. Finally, the simulation results of the cross array and the L-shaped array verify the effectiveness of the algorithm in this paper and further extend it to the two-dimensional MIMO radar array model and compare it with the traditional ESPRIT-like algorithm and the REC-FBSS-ESPRIT algorithm. In comparison, the algorithm in this paper has better performance under the conditions of successful resolution, DOA estimation accuracy, and low signal-to-noise ratio.


Introduction
e multiple-input and multiple-output (MIMO) radar system is a new radar system proposed in recent years. MIMO technology has brought a new breakthrough to the application performance of the radar system. Compared with the traditional radar, MIMO radar has potential advantages in parameter target estimation, target parameter detection, estimation performance analysis, space-time adaptive processing suppression, radar interference waveform design, and so forth. e direction of arrival (DOA) [1] estimation problem is one of the important research directions in the signal processing of sensor arrays [2,3], and it is widely used in radar [4][5][6]. High-precision DOA estimation based on subspace has become the focus of research [7][8][9]; for example, multiple signal classification (MUSIC) and estimation of signal parameters via rotational invariance techniques (ESPRIT) [10] can provide higher resolution to estimate the direction of arrival of uncorrelated and partially related signals [11][12][13]. But in practice, there are a large number of coherent sources due to multipath propagation and cofrequency interference [14]. When the signals are coherent, the rank of the covariance matrix accepted by the array will be deficient, so that the signal will be diffused to the noise, and the performance of the DOA estimation will be degraded [15,16]. e methods of decoherence include spatial smoothing algorithms, such as forward smoothing, backward smoothing, and forward and backward smoothing, but they are usually only applicable to isometric uniform linear arrays, and the corrected dimension is lower than the original matrix dimension; that is, decoherence is in exchange for lower degrees of freedom [6,14].
Another problem to solve rank deficit is to construct an algorithm based on vector and matrix reconstruction [14].
is algorithm rearranges the elements by covariance processing of the vector of the received signal and then constructs a Toeplitz matrix so that its rank is only the direction of arrival is related without being affected by signal correlation, so as to achieve the purpose of decoherence. Reference [17] proposed the ESPRIT-like algorithm by constructing a special antenna array model to reconstruct a Toeplitz matrix with the covariance matrix of the received signal matrix and using rotation invariance to achieve decoherence, but this method only works for the center array of the covariance matrix. When the center row of the covariance matrix is selected for reconstruction, the noise term is only a scalar quantity; otherwise, the denoising process still needs to be used to estimate the DOA using the ESPRIT algorithm. Because only one line of covariance is used, its information utilization is incomplete and affects DOA estimation performance. Reference [6] extended [17] to the two-dimensional DOA estimation of coherent sources. Reference [18] proposed the REC-FBSS-ESPRIT algorithm for reconstructing a covariant matrix, using the forward and backward spatial smoothing method as a preprocessing method to solve the problem of rank deficiency, and constructed an L-shaped matrix for two-dimensional DOA estimation, but the signal is still affected by noise and requires denoising. Reference [14] proposed an improved algorithm based on Toeplitz matrix reconstruction, by constructing a Toeplitz matrix set on the signal matrix, processing the Toeplitz matrix and the signal matrix, and then multiplying and summing by the conjugate transposed matrix. e matrix is corrected to obtain a full-rank matrix to achieve decoherence, and no additional denoising is required, but it is only suitable for one-dimensional DOA estimation.
Based on the literature [14], this paper proposes an improved Toeplitz matrix set reconstruction algorithm based on the two-dimensional reconstruction of the Toeplitz-like algorithms and uses cross arrays and L-shaped arrays for the two-dimensional DOA estimation [19], further extending to MIMO radar two-dimensional DOA estimation, using the characteristics of MIMO radar, with fewer array antennas to obtain higher DOA measurement accuracy. e algorithm uses two even linear arrays of mutually orthogonal array elements, through two permutation matrices and the relationship between the signal subspace and the direction vector, using rotation invariance, through two one-dimensional processing and then angle matching. us, DOA estimation of two-dimensional coherent signals is realized. Finally, simulation experiments verify the universality of the proposed algorithm for two-dimensional DOA estimation. Compared with the traditional ESPRIT-like algorithm and the REC-FBSS-ESPRIT algorithm in [18], the successful resolution probability and DOA estimation accuracy of the algorithm in this paper have better performance in the case of low signal-to-noise ratio, does not require spectrum peak search, and the amount of calculation is small.

Cross Array.
A cross-shaped array as shown in Figure 1 supposes the intersection of the cross-shaped array is the origin of the coordinates, the array elements are evenly distributed on the x-axis and y-axis, and the number of array elements on the x-axis and y-axis is N � 2M + 1. e total number of array elements is 4M + 1. e array element spacing is d � (λ/2), where λ is the wavelength. Assuming that P signals are incident on the antenna, the two-dimensional arrival angles of the ith signal are θi, ϕi (i � 1, 2, . . ., P), and θi and ϕi are the azimuth angle and the pitch angle, respectively. en, the received signals of x-axis and y-axis are where , a x (θ 2 ), . . . , a x (θ p )] and A y � [a y (θ 1 ), a y (θ 2 ), . . . , a y (θ p )] are the x-axis and y-axis direction matrices, respectively; the dimension is (2M + 1) × P, and a x (θ p ) � [e − j2πM d sin θ p sin ϕ p /λ, . . . , 1, . . . , e − j2πM d sin θ p sin ϕ p /λ ], a y (θ p ) � [e − j2πM d sin θ p sin ϕ p /λ , . . . , 1, . . . , e − j2πM d sin θ p sin ϕ p /λ ]; s � [s 1 , s 2 , . . . , s p ] is P × 1 dimension incident signal vector; P signals can be independent, correlated, or coherent; N x and N y are the additive Gaussian noise of the x-axis and y-axis receiving models, respectively; the dimensions are 2M + 1; and they are not related to the source.
First, construct the Toeplitz matrices and of the received signal. As can be seen from [14], Toeplitz matrices B x (t) and B y (t) can be expressed as follows: In other words, the rank of S is independent of the coherence between signals, and the decorrelation can be achieved. B Nx (t) and B Ny (t) are Toeplitz matrices composed of noise vectors. Figure 2: let the coordinate origin of the L-shaped array be the reference array element. e array elements are evenly distributed on the x-axis and the y-axis. e numbers of array elements on the x-axis and the y-axis are both N � 2M + 1 and the total number of array elements is 4M + 1. e array element spacing is d � (λ/2), where λ is the wavelength. Assuming that P signals are incident on the antenna, the two-dimensional arrival angles of the ith signal are θi, ϕi (i � 1, 2, . . ., P), and θi and ϕi are the azimuth angle and the pitch angle, respectively.

L-Shaped Array. L-shaped array is shown in
Similarly, the received signals of the x-axis and y-axis are In the same way, where Ax and Ay are the direction matrix of the x-axis and y-axis, respectively, the dimension is . , e (j2πM d cos θ p sin ϕ p /λ) ]; and S and N are the same as the cross array.
e Toeplitz matrices B x (t) and B y (t) are constructed as follows: Among them, C 1 and C 2 are constants; since S is a diagonal matrix, the correlation matrix is not affected by the constants. If S � CS, then B x (t) and B y (t) are the same.

MIMO Array.
e single-static MIMO array is shown in Figure 3; set the number of transmitting array elements as 1 × M and the number of receiving array elements as N × 1;  that is, the total number of array elements is M + N. Due to the nature of the MIMO array model, it can be virtualized. e number of arrays is M × N, and the array elements are uniformly distributed on the x-axis and y-axis, respectively, and the total number of array elements is M × N. e element spacing is d � (λ/2), where λ is the wavelength. Assuming that P signals are incident on the antenna, the twodimensional arrival angles of the ith signal are θi, ϕi (i � 1, 2, . . ., P), and θi and ϕi are the azimuth and elevation angles, respectively. en, the output is where θ k , ϕ k are the elevation and azimuth angles corresponding to the kth target; n(l) is the MN × 1 Gaussian white noise vector, the mean is 0, and the covariance is σ2I MN ; and where a ty (θ k , ϕ k ) and a tx (θ k , ϕ k ) are the antenna steering vectors along the y-axis and x-axis of the transmitting array (corresponding to the kth target) and a ry (θ k , ϕ k ) and a rx (θ k , ϕ k ) are the same. So, there are the following relations Correlation calculation of B x (t) and B y (t) with the corresponding received signals X i (t) and Y i (t) of the ith array element: It can be obtained from equation (7) that, for the ESPRITlike algorithm in [17], the method is to use an optional set of Toeplitz matrices R Xi and R Yi to achieve the purpose of decoherence of coherent source signals. If i � 0 is selected as the central array element to construct R Xi and R Yi , then I M+1,i is an identity matrix, and after that the relationship between the signal subspace and the direction vector can be used to directly perform feature decomposition on R Xi and R Yi to find the final DOA information, but when i ≠ 0 is the central array element, I M+1,i is not an identity matrix, so denoising is still required. And ESPRIT-like only constructs one R Xi and R Yi , so the complete information of the array is not fully utilized.
e same is for the MIMO array. R XΣ and R YΣ obviously have a joint diagonalization structure and span the same range of space of A rx and A ry , where R s � E[SS H ] is the signal covariance matrix. For D xk , D yk , C xi , and C yi , see literature [14].
From equations (8) and (9), it can be seen that the reconstructed equivalent covariance matrix is a positive definite matrix with rank P. e rank is not affected by signal correlation, but only related to the number of signals [14]. To further improve accuracy, R XΣ and R YΣ are revised where J is a permutation matrix with an antidiagonal of 1.
Since A rx and A ry are Vandermonde matrices and θi ≠ θ j , A rx and A ry are full ranks, and because M k�0 D k H are diagonal matrices, their feature decomposition can obtain P large eigenvalues and M + 1 − P small eigenvalues, where the eigenvectors corresponding to large eigenvalues constitute the signal subspace Φ s � span v 1 , v 2 , . . . , v p and the eigenvectors corresponding to the small eigenvalues constitute the noise subspace Φ s � span v 1+p , v 2+p , . . . , v M .

Two-Dimensional DOA Estimation of Coherent Sources
, where a 1 and a M+1 are the first and last terms of A rx . Similarly, the signal subspace is where u 1 and u M+1 are the first and last columns of U Sx , respectively. It is easy to get It is a rotation matrix. Since A rx and U Sx are the same signal subspace, U Fx � A Fx T and U Bx � A Bx T, in which T is a full-rank matrix; then, We can derive e − (j2π d sin θ k sin ϕ k /λ) , k � 1, 2, . . . , p, by feature decomposition of Ψ x . Both sin θ k sin ϕ k can be estimated, and cos θ k sin ϕ k can be estimated in the same way. e following two-dimensional DOA angle can be estimated: In summary, the steps of the two-dimensional DOA estimation method based on the Toeplitz matrix set proposed in this paper are shown in Table 1.

Computer Simulation
In order to verify the universality of the proposed algorithm for two-dimensional DOA estimation, the proposed algorithm is compared with the cross-shaped array ESPRIT-like algorithm in [6] and the L-shaped array REC-FBSS-ESPRIT algorithm in [18]. Perform 100 Monte Carlo simulations on each algorithm and define the root mean square error as e DOA of the three sources used in the simulation are (35°, 40°), (75°, 70°), (89°, 80°), and sources 1, 3 are related signals. e numbers of array elements on the x-and y-axis are both 11, and the array element spacing is half a wavelength. Figure 4 shows the DOA estimation performance of each algorithm, where Figures 4(a)-4(e) are the cross-shaped ESPRIT-like algorithm, L-shaped array REC-FBSS-ESPRIT algorithm, MIMO array, and the DOA estimation results of the cross array, L-type array, and MIMO array using the algorithm of this paper when the signal-to-noise ratio is SNR � 10 and the number of snapshots is 1024. From the figure, it can be seen intuitively that the algorithm of this paper is closer to the actual angle and has better two-dimensional DOA estimated performance. Figure 5 shows the performance comparison of azimuth estimation using the cross array, L-shaped array, MMO array, ESPRIT-like algorithm, and REC-FBSS-ESPRIT algorithm using this algorithm and examines the RSME variation of the algorithm under different signal-to-noise ratios (SNR). When the signal-to-noise ratio is 0, the performance of this algorithm is improved by about 40% compared to the ESPRIT-like algorithm and improved by about 25% compared to the REC-FBSS-ESPRIT algorithm. By comparing RSMEs with different signal-to-noise ratios, it is concluded that the proposed algorithm performs better Security and Communication Networks Table 1: Algorithm steps in this paper.
Step 1: construct the received data X and Y from different matrix models and construct the Toeplitz matrix from equations (2) and (4) Step 2: RXi and RYi are obtained from Toeplitz matrix and receiving matrix through formula (7) Step 3: select all RXi and RYi to get RXΣ and RYΣ from formulas (8) and (9) Step 4: R is calculated from the R � (1/2)(RΣ + JR * J) formula Step 5: perform feature decomposition on R to obtain signal subspace US Step 6: get the corresponding UB and UF from US, get the corresponding Ψ, and then perform feature decomposition Step 7: match by equation (13)   than the ESPRIT-like algorithm and REC-FBSS-ESPRIT algorithm at low signal-to-noise ratios. Figure 6 shows the comparison of the resolution probability of the cross array, L-shaped array, MIMO array ESPRIT-like algorithm, and REC-FBSS-ESPRIT algorithm using this algorithm. Here, it is defined that if the estimated θ 1 and θ 2 directions of the two sources satisfy |θ 1 − θ 1 | + |θ 2 + θ 2 | < |θ 1 − θ 2 | , it is said that the two sources are correctly resolved. e resolution probability refers to the ratio of the number of correct resolutions to the total number of experiments. e experimental conditions were 10 experiments per group under the condition that the number of snapshots was 1024 and the signal-to-noise ratio SNR was 20.
e experimental results show that the probability of successful resolution of this algorithm under low azimuth interval is higher than that of ESPRIT-like algorithm and REC-FBSS-ESPRIT algorithm, which shows that the accuracy of DOA estimation of this algorithm is higher than the other two algorithms, and, moreover, the cross array, the L-shaped array, and the MIMO array all have better successful resolution in the case of low azimuth spacing. Figure 7 shows the influence of correlation coefficients on the performance of each algorithm. From Figure 7(a), we can see the successful resolution of the MIMO array using this algorithm when the signal-to-noise ratio SNR is 10 and the number of snapshots is 1024, and with the change of the correlation coefficient, the success resolution is always 100%, and the success rate of the L-shaped array using the algorithm of this paper is only 99% to 100% when the correlation coefficient is 0.9. e success rate of the ESPRIT-like algorithm and the REC-FBSS-ESPRIT algorithm is many times between 99% and 100%. It can be seen from the Figure 7(b) that the two arrays adopting the algorithm of this paper have small changes in RSME when the correlation coefficient changes and slightly rise when the correlation coefficient is close to 1, while the ESPRIT-like algorithm has both changes when the correlation coefficient changes. e square error is also small, but the performance is lower than the algorithm in this paper. e mean square error of the REC-FBSS-ESPRIT algorithm under the change of the correlation coefficient is lower than the ESPRIT-like algorithm but higher than the algorithm in this paper. e experimental results show that the two-dimensional DOA estimation algorithm in this paper is less affected by the correlation coefficient, and the estimation performance is higher than the ESPRIT-like and REC-FBSS-ESPRIT algorithms.    Table 2 shows the performance comparison of the cross array, L-shaped array, ESPRIT-like algorithm, and REC-FBSS-ESPRIT algorithm under different parameters using this algorithm. Experiments show that the algorithm in this paper has a low signal-to-noise ratio and a low number of snapshots; the mean square error is lower than the ESPRITlike algorithm and the REC-FBSS-ESPRIT algorithm. Table 3 shows the comparison of the time required for different algorithms to perform 100, 500, and 1000 experiments on 3 signals. e results show that the simulation time of the algorithm in this paper is similar to that of the ESPRIT-like algorithm, only slightly slower, but compared to the REC-FBSS-ESPRIT algorithm, it is much faster and has a better real-time performance.

Conclusions
is paper proposes a modified Toeplitz matrix set reconstruction algorithm based on the two-dimensional reconstruction of Toeplitz-like algorithms, which is based on the incomplete information of traditional reconstruction Toeplitz algorithm and the need for denoising.
e Toeplitz matrix set is multiplied and summed with the conjugate transposed matrix to achieve solution coherence, avoiding additional denoising. e simulation experiment gives a comparative analysis of the performance of this algorithm, the ESPRIT algorithm, and REC-FBSS-ESPRIT algorithm and utilizes the special performance of MIMO arrays to simulate more arrays with fewer antennas, universality under the array model.

Data Availability
e datasets used or analysed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.