Abstract
Direction of arrival (DOA) estimation problem for multiple-input multiple-output (MIMO) radar with unknown mutual coupling is studied, and an algorithm for the DOA estimation based on root multiple signal classification (MUSIC) is proposed. Firstly, according to the Toeplitz structure of the mutual coupling matrix, output data of some specified sensors are selected to eliminate the influence of the mutual coupling. Then the reduced-dimension transformation is applied to make the computation burden lower as well as obtain a Vandermonde structure of the direction matrix. Finally, Root-MUSIC can be adopted for the angle estimation. The angle estimation performance of the proposed algorithm is better than that of estimation of signal parameters via rotational invariance techniques (ESPRIT)-like algorithm and MUSIC-like algorithm. Furthermore, the proposed algorithm has lower complexity than them. The simulation results verify the effectiveness of the algorithm, and the theoretical estimation error of the algorithm is also derived.
1. Introduction
Multiple-input multiple-output (MIMO) radars have attracted a lot of attention recently for their potential advantages over conventional phased-array radars [1–4]. MIMO radar systems can overcome fading effect, enhance spatial resolution, strengthen parameter identifiability, and improve target detection performance for the additional degrees of freedom [5–7]. Angle estimation is a key issue in MIMO radar and has been studied by a lot of researchers. Estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithms [8–10] can obtain the closed-form estimations via rotational invariance in the subspace. Capon algorithm [11] and multiple signal classification (MUSIC) algorithm [12] both estimate the angles via the peak searches. Polynomial root finding technique can transform the peak searches into root finding problem when the arrays are uniform linear arrays (ULA) [13]. Based on the uniqueness of trilinear decomposition, parallel factor analysis (PARAFAC) algorithms [14–16] can also be adopted for the angle estimation. For monostatic MIMO radar, low-complexity ESPRIT [17] and Capon [18] can achieve better angle estimation performance with lighter computation burden based on the reduced-dimension transformation. However, these methods strongly depend on the array manifold, which is often perturbed by the mutual coupling in practical situations. The mutual coupling will make the above methods have performance degradation or even fail to work, and various methods for mutual coupling compensation have been studied. The Receiving-Mutual-Impedance Methods (RMIM) [19, 20] can obtain the decoupled voltages via the known terminated impedance and the receiving mode antennas. The calibration methods [21–23], which will be studied in this paper, can utilize the special characteristics of the mutual coupling coefficients of uniform linear arrays (ULA). Reference [24] proposed a MUSIC-like algorithm for the angle estimation in the presence of mutual coupling. Through reconstructing the MUSIC function based on the Toeplitz structure of the mutual coupling matrix, the angles can be estimated via the peak searches. To avoid the high-computational peak search, an ESPRIT-like algorithm was proposed in [25] by removing the auxiliary sensors to eliminate the influence of mutual coupling. The rotational invariance can be extracted in the new data for the angle estimation.
Although MUSIC-like algorithm can achieve high resolution angle estimation with unknown mutual coupling, the high computational peak searches will make it ineffective. As a faster method, ESPRIT-like can obtain the close-form solution of the angle estimation, but it only exploits the relationship between the subarrays, which will lead to a performance loss.
Our goal is to develop an effective method for angle estimation in MIMO radar with unknown mutual coupling. So, in this paper, we propose a Root-MUSIC based method, which is fast and can deal with the mutual coupling problem, which often occurrs in practical situation. The nonauxiliary sensors are also used to eliminate the influence of mutual coupling. Then a transformation is utilized to reduce the dimension of the data, meanwhile a direction matrix with Vandermonde structure can be obtained. Finally, Root-MUSIC, which is a root-finding form of the MUSIC peak searching, can be adopted for the angle estimation. The angle estimation performance of the proposed algorithm is better than that of ESPRIT-like algorithm and MUSIC-like algorithm. Furthermore, it may also have lower complexity.
Notation. , , , and denote transpose, conjugate-transpose, inverse, and pseudo-inverse operations, respectively. stands for diagonal matrix whose diagonal element is a vector . is a identity matrix. is the Kronecker product. is to get real part of the complex, and means taking phase operation.
2. Data Model
As shown in Figure 1, the monostatic MIMO radar system is equipped with both uniformlinear arrays (ULA) for the transmit and receive arrays, and the arrays are both located in the -axis with half-wavelength spacing between adjacent antennas, respectively. Assume that there are targets in the - plane, and the output of the matched filters at the receiver can be expressed as [5, 23, 24] where is the direction of arrival (DOA) of the th target with respect to the transmit array normal or the receive array normal, , with and being the Doppler frequency and the reflect coefficient, respectively, is an Gaussian white noise vector of zero mean and covariance matrix , and and are the receive and transmit steering vectors for the th target under the influence of mutual coupling, respectively. They satisfy where and represent the matrices containing the mutual coupling coefficients. With the assumption of ULA, they are Toeplitz matrices [24]where and stand for the mutual coupling coefficients. If the first sensors are the reference points, then . and in (2a) and (2b) are expressed as Vandermonde structureswhere .
By collecting samples, and representing the received data as , the received data can be expressed as where is the direction matrix, , and stands for the noise matrix.
3. Root-MUSIC Based Angle Estimation for MIMO Radar with Unknown Mutual Coupling
3.1. Mutual Coupling Elimination
According to [25], define the selecting matrices asThen,where , , and and are the first and elements of and , respectively. Define and , thenThe new steering vector of the output can be now characterized as According to (9), the steering vector can be expressed as a kronecker product of two Vandermonde vectors without considering the scalar coefficients and . Thus, the mutual coupling can be considered to be eliminated through this procedure.
3.2. Reduced-Dimension Transformation
The length of the new steering vector is , which is too long and will add high computation burden. Furthermore, it will make the later root finding hard to implement for the high-order polynomials. So the reduced-dimension transformation is necessary. According to [17], steering vectors with Vandermonde structures will satisfy where . Then define , Using the reduced-dimension transformation for the receive signal after the mutual coupling elimination, and combining (9) and (10), then where , and . The covariance matrix of is where and . can be regarded as the new direction matrix.
3.3. Root-MUSIC Based Angle Estimation
The covariance matrix in (14) can be decomposed as where denotes a diagonal matrix formed by largest eigenvalues and denotes a diagonal matrix formed by the rest smaller eigen-values (). and represent the signal subspace and noise subspace, respectively, of which stands for the eigenvectors corresponding to the largest eigen-values, and consists of the rest eigenvectors.
According to the MUSIC peak search function [12], the space spanned by the direction matrix is orthogonal to the noise subspace, and the angle can be estimated by searching the following function: where . To avoid the peak searching, let ; then can be expressed as . Equation (16) can be transformed into the root finding problem where . According to [26], the polynomial in (17) can be expressed as where is a column vector which contains the coefficients of the polynomial and the th coefficient is given by the sum of the th diagonal elements of .
The roots which are most close to the unit circle will be the solutions, and they are denoted as . The angle can be estimated via With respect to (14), the covariance matrix can be estimated via
The major steps of the proposed algorithm for DOA estimation in MIMO-radar are shown as follows.(1)Construct the new output via (13) to eliminate the mutual coupling and reduce the data dimension.(2)Estimate the covariance matrix of the data through .(3)Perform eigen-value decomposition of to get the noise subspace .(4)Construct the MUSIC function and use root finding technique to obtain the angle estimation via (18)-(19).
3.4. Complexity Analysis
The proposed algorithm has lower complexity than ESPRIT-like algorithm and MUSIC-like algorithm for it requires no peak searches and has lower dimension. Figure 2 shows the run time of the three algorithms in computer versus the number of antennas, and here we choose for simplify. From Figure 2, we find that our algorithm has lower complexity, and the change trend versus the number of antennas of the proposed algorithm is smaller than those of the other two algorithms.
4. Theoretical Analysis and Cramer-Rao Bound (CRB)
4.1. Theoretical Analysis of the Proposed Algorithm
The decomposition of in (20) can be expressed as where and denote the estimated signal and noise subspaces, respectively.
According to [27], Root-MUSIC and MUSIC peak search method have the same angle estimation error as they both seek for the extreme points of , where . The first-order derivative of on will be equal to zero; that is, where . Taking the first-order Taylor series expansion, then where some small items are neglected. The mean square error of the angle estimation can be shown as According to [28], where with . Thus the theoretical mean error of the angle estimation of the proposed algorithm can be derived as The theoretical result and the simulation result of the proposed algorithm will be compared in Section 5, where it can be shown that they are almost the same, which demonstrate the correctness of the derivation above.
4.2. Advantages of the Proposed Algorithm
The advantages of the proposed algorithm can be summarized as follows.(1)It works well with unknown mutual coupling. According to (9), (13), and (14), it can be shown the mutual coupling can be eliminated regardless of the specific distribution of the mutual coupling (the number of the antennas that are influenced by the mutual coupling should be known, i.e., and [21, 22]).(2)It requires no peak search and has lower complexity, which has been shown in Figure 2.(3)It has better angle estimation performance than MUSIC-like algorithm and ESPRIT-like algorithm. Part of the reason of this advantage can be refered to “Remark 1”, and another reason may be the reduced-dimension transformation. According to (13)-(14), the signal to noise ratios (SNRs) of the data before the reduced-dimension transformation and the data after the reduced-dimension transformation are where . It can be shown that , so the reduced-dimension transformation can obtain SNR gain. Advantage 3 can also be verified in Section 5, where the three algorithms are compared.
Remark 1. Although MUSIC-like algorithm [24] can achieve high resolution angle estimation with unknown mutual coupling, the peak searches with high complexity will make it ineffective. ESPRIT-like algorithm [25] can obtain the close-form solution of the angle estimation with lower complexity, but it only exploits the relationship between the subarrays, which will lead to the performance loss.
4.3. Cramer-Rao Bound (CRB)
According to [29], the CRB for DOA estimation in monostatic MIMO radar with unknown mutual coupling can be derived as where with ; with .
5. Simulation Results
Define root mean square error (RMSE) as , where is the estimate of DOA of the th Monte Carlo trial. We assume that there are targets with the angles being . The reflect coefficients are and the Doppler frequencies are . The mutual coupling coefficients are and for the transmit and receive arrays, respectively.
Figure 3 depicts angle estimation result of the proposed algorithm with , , , and dB. It is shown that the DOAs can be clearly observed.
The angle estimation performance comparison between the proposed algorithm, ESPRIT-like algorithm, and MUSIC-like algorithm is shown in Figure 4, from which we can find that the angle estimation performance of the proposed algorithm is better than that of ESPRIT-like and MUSIC-like. It can also be indicated that the theoretical estimation error and the simulation error of the proposed algorithm are almost the same, which demonstrate the correctness of the derivation in (26).
Figures 5 and 6 present the estimation performance of the proposed algorithm with and variable values of or . It can be indicated that the angle estimation performance of the proposed algorithm is improved with the increase of receive or transmit antenna numbers because of the diversity gain.
Figure 7 presents the angle estimation performance of our scheme with , , and different values of . It has been shown in (26) that the angle estimation error will be reduced when increases.
To further demonstrate the advantages of the proposed algorithm, another set of parameters is chosen to test the algorithms. Assume that the angles are . The reflect coefficients are and the Doppler frequencies are . The mutual coupling coefficients are and for the transmit and receive arrays, respectively. The RMSE comparison of the algorithms is shown in Figure 8, from which we can still find that the angle estimation performance of the proposed algorithm is better than the other two algorithms. It can also be indicated that the theoretical estimation error and the simulation error of the proposed algorithm are still the same.
6. Conclusion
In this paper, a DOA estimation algorithm based on Root-MUSIC has been proposed for monostatic MIMO radar with unknown mutual coupling. The proposed algorithm has low complexity and works fast, requires neither peak searches nor iterations, and can work well with practical unknown mutual coupling. So it can support practical online processing. Furthermore, the proposed algorithm has better angle estimation performance then the ESPRIT-like algorithm and MUSIC-like algorithm, which have higher complexity. Many thanks are due to the reviewer’s comments; our future work is to test the algorithms with real world data, and try to convert the theories into practical applications.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the China NSF Grants (61371169, 61071164), the Jiangsu Planned Projects for Postdoctoral Research Funds (1201039C), the China Postdoctoral Science Foundation (2012M521099, 2013M541661), the Open project of key laboratory of underwater acoustic communication and marine information technology (Xiamen University), the Hubei Key Laboratory of Intelligent Wire1ess Communications (IWC2012002), the Open project of Key Laboratory of modern acoustic of Ministry of Education (Nanjing University), the Aeronautical Science Foundation of China (20120152001), Funding of Jiangsu Innovation Program for Graduate Education (CXZZ13_0165), Funding for Outstanding Doctoral Dissertation in NUAA (BCXJ13-09), the Qing Lan Project, priority academic program development of Jiangsu high education institutions, and the Fundamental Research Funds for the Central Universities (NS2013024, kfjj130114).