Noncircular Sources-Based Sparse Representation Algorithm for Direction of Arrival Estimation in MIMO Radar with Mutual Coupling

In this paper, a reweighted sparse representation algorithm based on noncircular sources is proposed, and the problem of the direction of arrival (DOA) estimation for multiple-input multiple-output (MIMO) radar with mutual coupling is addressed. Making full use of the special structure of banded symmetric Toeplitz mutual coupling matrices (MCM), the proposed algorithm firstly eliminates the effect of mutual coupling by linear transformation. Then, a reduced dimensional transformation is exploited to reduce the computational complexity of the proposed algorithm. Furthermore, by utilizing the noncircular feature of signals, the new extended received data matrix is formulated to enlarge the array aperture. Finally, based on the new received data, a reweighted matrix is constructed, and the proposed method further designs the joint reweighted sparse representation scheme to achieve the DOA estimation by solving the l1-norm constraint minimization problem. The proposed method enlarges the array aperture due to the application of signal noncircularity, and in the presence of mutual coupling, the proposed algorithm provides higher resolution and better angle estimation performance than ESPRIT-like, l1-SVD and l1-SRDML (sparse representation deterministic maximum likelihood) algorithms. Numerical experiment results verify the effectiveness and advantages of the proposed method.


Introduction
With orthogonal transmitted waveforms, multiple-input multiple-output (MIMO) radar has drawn increasing attention in the field of wireless communications.Compared with the conventional phased-array radar, MIMO radar owns a number of advantages, such as higher resolution and better parameter identifiability [1].In addition, MIMO radar can be classified into the following two types: statistical MIMO radar and colocated MIMO radar.Colocated MIMO radar can achieve more degrees of freedom and higher spatial resolution, because its closely-spaced antennas form a virtual array with a large aperture.Colocated MIMO radar includes the monostatic one and the bistatic one.In bistatic MIMO radar, the transmit array and the receive array are separated away, while they are close to each other in monostatic MIMO radar.In this paper, the angle estimation problem is investigated in monostatic MIMO radar.
Parameter estimation is an important aspect in wireless communications and sensor array signal processing [2].In order to estimate the direction of arrival (DOA), a large number of subspace-based algorithms have been proposed [3][4][5].Recently, the emerging sparse representation (SR) and compressed sensing (CS) have attracted more and more attention in various fields [6][7][8], such as sensor array signal processing, imaging and image processing, and so forth.Furthermore, it has been verified that SR-based DOA estimation algorithms have remarkable advantages over other methods [9], for instance they adapt better to challenging circumstances and provide higher angle resolution.A few SR-based algorithms, such as l 1 -SVD (singular value decomposition) [9] and l 1 -SRACV (array covariance vectors) [10], have been proposed via sparse signal recovery.In [11], based on the array covariance vector, the revised real-valued (RV) l 1 -SVD approach has better DOA estimation performance and lower computation burden than l 1 -SVD [9] and RV l 1 -SVD [12].
However, in practical situations, the sensor array manifold is often affected by the errors of unknown mutual coupling.When they are taken into account, the DOA estimation performance of the above-mentioned methods is seriously degraded.With mutual coupling, the subspace-based MUSIC-like (multiple signal classification-like) [13] algorithm is proposed.Besides, a ESPRIT-like (estimation of signal parameters via rotational invariance techniques-like) algorithm is presented in [14], which owns a low calculation burden.In [15], the SR-based revised l 1 -SVD algorithm eliminates the errors of mutual coupling, and it achieves angle estimation by solving the problem of the l 1 -norm constraint minimization.In [16], in the presence of mutual coupling, by introducing the refined deterministic maximum likelihood (DML) procedure, a covariance vector-based sparse representation algorithm is proposed with outstanding angle estimation performance.We call the method in [16] l 1 -SRDML (sparse representation deterministic maximum likelihood).
In practical communication systems and radar applications, the complex noncircular sources are widely used, such as binary phase shift keying (BPSK), multiple amplitude shift keying (MASK) and unbalanced quadrature phase shift keying (UQPSK) modulated signals [17].Noncircular signals can be used to enlarge the array aperture without extra antennas, which contributes to the improvement of angle estimation performance [18].Exploiting the signal noncircularity, some subspace-based methods have been developed for angle estimation in MIMO radar [17,19,20].They have verified that compared with traditional subspace-based algorithms, the methods that are based on noncircular signals can achieve better DOA estimation performance.However, these subspace-based approaches do not take the errors of unknown mutual coupling into account.Moreover, they perform poorly when experiencing challenging circumstances, such as fewer snapshots, which can be solved by sparse representation-based methods.For DOA estimation in MIMO radar with mutual coupling, to the best of our knowledge, up to now, there are no references about SR-based methods that exploit the signal noncircularity to extend the virtual array.Consequently, utilizing noncircular sources and the sparse representation framework, we aim to eliminate the effect of mutual coupling and then achieve a better DOA estimation.
In this paper, by using the noncircular signals, we propose a sparse representation-based DOA estimation algorithm in MIMO radar with mutual coupling.The new method firstly eliminates the effect of mutual coupling by extracting the information of mutual coupling coefficients into a diagonal matrix.Secondly, based on the reduced dimensional transformation and the characteristics of signal noncircularity, the new extended received data matrix is formulated to enlarge the array aperture.Finally, a reweighted matrix is constructed, based on which a joint reweighted l 1 -norm minimization sparse representation framework is designed to obtain the DOAs.In the presence of mutual coupling, the proposed algorithm performs better than the conventional algorithms on account of the error elimination and the application of noncircular signals.
The rest of this paper is organized as follows.In Section 2, the MIMO radar system model with mutual coupling and noncircular signals is described.In Section 3, the implementation process of the proposed method is described in detail, including mutual coupling elimination, noncircular signal-based extended matrix construction and the joint reweighted sparse representation-based DOA estimation scheme.In Section 4, we give some related remarks and discussions regarding the parameter setting and the computational complexity of the proposed method.Then, in Section 5, simulations are implemented to verify the efficiency and the advantages of the proposed method.Finally, our conclusions are summarized in Section 6.

MIMO Radar System Model with Mutual Coupling
Consider a narrowband monostatic MIMO radar system, shown in Figure 1, whose transmit and receive arrays are both half-wavelength d spaced uniform linear arrays (ULAs).In the transmit array, M antennas transmit M orthogonal noncircular waveforms, such as BPSK modulated signals.In the receive array, N antennas are impinged by P uncorrelated targets that can be deemed as point scatterers in the far-field.Different from the mutually-separated arrays in bistatic MIMO radar, the transmitter and the receiver are closely located in monostatic radar, as shown in Figure 2. Since the transmit array and the receive array are collocated at the same site in monostatic MIMO radar systems [17], the echo signal travels the same path as the transmitted signal.Therefore, for the p-th target, the directions of departure of the transmitted signal and arrival of the echo from the target are the same [21], i.e., θ DOD = θ DOA , both of which are denoted as DOA θ p , p = 1, 2, . . .P. When mutual coupling is taken into account in both the transmit array and the receive array, at the receiver, the output of the matched filters at a certain snapshot is given by [16]: where x(t) ∈ C MN×1 is the received data vector, s(t) = [s 1 (t), s 2 (t), . . ., s P (t)] T ∈ C P×1 is the complex-valued noncircular signal vector and n(t) ∈ C MN×1 is the additional Gaussian noise vector with zero-mean and covariance matrix σ 2 I MN , in which σ 2 is the noise power.Meanwhile, the detailed expression of Ã is [14]: where Ã ∈ C MN×P , a t (θ p ) = [1, e jπ sin(θ p ) , e jπ2 sin(θ p ) , . . ., e jπ(M−1) sin(θ p ) ] T is the transmit steering vector and a r (θ p ) = [1, e jπ sin(θ p ) , e jπ2 sin(θ p ) , . . ., e jπ(N−1) sin(θ p ) ] T is the receive steering vector, p = 1, 2, . . .P. C t and C r are mutual coupling matrices of the transmit array and the receive array, respectively.C t and C r can be modeled as banded symmetric Toeplitz matrices, which are expressed as [22][23][24]: where i = t, r, c ij is the non-zero mutual coupling coefficient, j = 0, 1, ..., k.In Figure 1, the coupling effects of the m-th and the n-th elements in the transmit and the receive arrays are described, respectively.β (m,q) is the coupling contribution from the q-th to the m-th elements in transmit array, and β(n,q ) is the contribution from the q -th to the n-th elements in the receive array, 1 ) and β(n,q ) contain the information of c t|q−m| and c r|q −n| , respectively.Mutual coupling coefficients are factors related to the distance between the two antennas [23], and for the k + 1 non-zero coefficients, they satisfy 0 Besides, for the ULA coupling model, just a few mutual coupling coefficients are non-zero ones [22][23][24], and without loss of generality, 2k < min{M, N} is assumed.Then, by collecting J snapshots, the received data in MIMO radar with mutual coupling is represented as follows: where

Noncircular Signals
For complex random variables and vectors, circularity is an important property [18].Making full use of the properties of noncircular signals, the array aperture can be enlarged.The invariance specialityof rotation can be used to judge if a complex random sequence is noncircular [18].More specifically, if stationary complex signal sequence q is circular, it accords with E[ q qH ] = 0 and E[ q qT ] = 0.By comparison, if the sequence is noncircular, then E[ q qH ] = 0 and E[ q qT ] = 0, which can be utilized to increase the number of effective array elements.In addition, for zero-mean stationary complex signal sequence q, it holds the property that: where 0 ≤ ρ ≤ 1 and ψ are noncircular rate and noncircular phase angle, respectively, and ρ depends on the signal modulation type.In this paper, we consider the signals with maximum noncircular rate, namely, ρ = 1 in s(t), such as BPSK and MASK modulated signals that are widely used in MIMO radar systems.It has been pointed out that the noncircular signal vector s(t) in Equation ( 1) can be rewritten as [17,25]: where  [26].This indicates that the relationship between ϕ p and ψ p can be written as ϕ p = ψ p /2.With the application of signal noncircularity, X in Equation ( 4) is expressed as: where S r = [s r (t 1 ), . . ., s r (t J )] is a real-valued matrix.

Mutual Coupling Elimination
To successfully implement the noncircular source-based sparse representation method, we firstly eliminate the effect of mutual coupling.Exploiting the structure of mutual coupling matrices C t and C r , we define two selection matrices as follows: where Γ 1 ∈ C (M−2k)×M , and Γ 2 ∈ C (N−2k)×N .Γ 1 and Γ 2 can be used to choose the central rows of the mutual coupling matrices and then transform the transmit-receive mutual coupling matrix C into a diagonal one.Thus, let z p = e jπ sin(θ p ) , by multiplying the selection matrix Γ 1 on the left side of the transmit steering vector with mutual coupling; the transformation of the i-th element in Γ 1 C t a t (θ p ) ∈ C (M−2k)×1 can be derived as: where ω tp is a scalar that just contains the information of θ p and mutual coupling in the transmit array.

Noncircular Signal-Based Extended Matrix Construction
Based on the new received data matrix X, a reduced dimensional transformation is introduced in the following.For the p-th target, the transmit-receive steering vector ã(θ p ) can be expressed as: where p = 1, 2, . . ., P. Many repeated terms are contained in ã(θ p ), and the non-repeated elements can make up another vector b(θ The relationship between ã(θ p ) and b(θ p ) can be derived as ã(θ p ) = Gb(θ p ), in which G is represented as: where G∈ C M N×( M+ N−1) , and Thus, based on G, a reduced dimensional transformation matrix can be constructed as ) G H , then X in Equation ( 12) turns into: where Y and Nd = R N are the new reduced dimensional data matrix and the new noise matrix, respectively. 1)can be directly calculated as: By the transformation, the dimension of the data matrix is reduced from M N × J to ( M + N − 1) × J.In addition, as RR H = I M+ N−1 , the transformation does not bring in the additional spatial colored noises.In order to enlarge the array aperture by exploiting the noncircular signals, the new received data matrix Y∈ C ( M+ N−1)×J can be extended as Z = [Y T , Y H ] T ∈ C 2( M+ N−1)×J .Since F, D and S r are real-valued matrices, Z is expressed as:

Joint Reweighted Sparse Representation-Based DOA Estimation Scheme
In order to successfully apply the sparse representation theory to estimate the DOAs, we bring in the SVD technique of Z, i.e., Z = UΛV H , where Λ = diag(ω 1 , ω 2 , . . ., ω M+ N−1 ), and are singular values.Let V s be composed of the P vectors in V that correspond to (ω 1 , ω 2 , . . ., ω P ).With V s , the dimension of the observation matrix can be further reduced [9], which is conducive to the signal reconstruction, that is: where where T θ ∈ C L×P and T θ ∈ C L×P hold the same row support with T and T, respectively.For the two models in Equation ( 19), if they are solved independently by conventional sparse representation methods that are based on the l 1 -norm penalty [9], the signal noncircularity cannot be utilized to extend the array aperture.In order to make the best of signal noncircularity, we design a joint , then the sparsities of T θ and T θ are combined.In r(k), T(l 2 ) θ (k) and T(l 2 ) θ (k) are defined as T(l 2 ) θ (k) = T θ (k, :) 2 and T(l 2 ) θ (k) = T θ (k, :) 2 , respectively.More specifically, T(l 2 ) θ (k) is equal to the l 2 -norm of the k-th row in T θ ; similarly, T(l 2 ) θ (k) is equal to the l 2 -norm of the k-th row in T θ .Let r represent the solution vector of the sparse signal reconstruction, r(k) ≥ r(k), k = 1, 2, . . ., L. As a result, based on the l 1 -norm constrained minimization, a joint sparse representation scheme is designed as follows: where η 1 and η 2 are regularization parameters that set the error amount.To make the l 1 -norm penalty better approximate the l 0 -norm penalty, we introduce the reweighted matrix in the following.Let U n be the noise subspace and composed of the ( M + N − 1 − P) vectors in U n that correspond to (ω P+1 , ω P+2 , . . ., ω M+ N−1 ).According to the subspace principle [13], for the true target θ p , [(ω tp ω rp e jϕ p Fb(θ p )) T , (ω tp ω rp e −jϕ p Fb(θ p ) * ) T ] T in Z is orthogonal to the noise subspace U n , that is: where: Since Q(θ p , ϕ p , ω tp , ω rp ) → 0, ω tp , ω rp and e jϕ p are all non-zero scalars and the rank of for the true DOA θ p , the rank reduction of Q(θ p ) will take place, that is det[ Q(θ p )] → 0. Consequently, a reweighted vector is designed as: for i = 1, 2, . . ., L. Thus, the reweighted matrix is: where W r ∈ C L×L is a diagonal matrix.Large weights in W r can punish the entries that are more likely to be zeros in w, while small weights reserve the larger entries.Hence, W r can enhance the solution of the sparse recovery problem in Equation (20).Then, for DOA estimation in the presence of mutual coupling, the noncircular source-based joint reweighted sparse representation framework transforms into: Effective SOC (second order cone) programming software packages, such as SeDuMi [27] and CVX [28], can be used to solve Equation (25).As a result, by plotting r, the DOA estimation is achieved.

Related Remarks
Remark 1.In the problem of the sparse signal reconstruction that is based on the l 1 -norm minimization, the chosen values of the regularization parameters η 1 and η 2 are important for the estimation accuracy.In the proposed algorithm, noise matrix N in Equation ( 4) is complex white Gaussian with zero mean.By mutual coupling elimination and reduced dimensional transformation, the corresponding selection matrix Γ and transformation matrix R transform N into Nd .According to the designed matrix structure of Γ and R, ΓΓ H = I M N and RR H = I M+ N−1 can be obtained.
Based on the invariance speciality of linear transformation in an asymptotic normal distribution [29], if n is asymptotic normal with zero-mean and covariance matrix σ 2 I MN , RΓn complies with a zero-mean asymptotic normal distribution with covariance matrix RΓ(σ 2 I MN )(RΓ) H = σ 2 I M+ N−1 .Therefore, Nd = [ nd1 , nd2 , . . ., ndJ ] in Equation ( 15) is a zero-mean complex Gaussian white noise matrix.
In [9], it has been verified that if the noise nd is independent and identically distributed (i.i.d.) Gaussian, then Nds 2 F has approximately a chi-square distribution with ( M + N − 1)P degrees of freedom upon normalization by the variance of nd , where Nds = Nd V s with V s being the signal subspace of Z [9].As a result, η 1 and η 2 can be chosen as the upper limit value of Nds 2 F and Ñds 2 F upon σ 2 , with a high probability 1 − ε confidence interval, and ε = 0.001 is enough.In addition, σ 2 can be estimated by the average of the squares of 2( M + N − 1)P smallest singular values of Z.By using MATLAB software, the function chi2inv(1 − ε, ( M + N − 1)P)σ 2 can be used to compute the regularization parameters.
Remark 2. The computational complexity of the proposed algorithm mainly focuses on eliminating the effect of mutual coupling, constructing the reweighted matrix and obtaining the sparse solution from Equation (25), ] L} and O{ L3 P} calculation burden, respectively, where L is the number of the discretized sampling grids.Therefore, the total computational complexity of the proposed algorithm is O{ M N MN J + [4( M + N − 1)(2( M + N − 1) − P + 1)] L + L3 P}.Although the constructions of the extended data matrix and the reweighted matrix increase the computational complexity, they enlarge the virtual array aperture and enhance the solution of the sparse reconstruction.As a result, the angle estimation performance of the proposed method is improved.Remark 3. The prior knowledge of the number of emitting sources is important.This is because in the absence of a perfect knowledge of the target number P, either adding spurious sources or missing actual signals may lead to the deviations of the signal and the noise subspaces, which further cause the deviations of extended observation matrix Z v in Equation ( 18) and reweighted matrix W r in Equation (24).Therefore, the selected value of P affects the performance of the proposed algorithm.Fortunately, when P is unknown in practical applications, some effective methods can be used to estimate it, such as the Akaike information criterion (AIC) and the minimum description length (MDL) [30].In the proposed algorithm, the number of emitting sources P is assumed to be known.Remark 4. After performing the mutual coupling elimination and the reduced dimensional transformation, the received data matrix X in Equation ( 4) turns into Y in Equation (15), and the number of effective antenna elements reduces from (M + N − 1) to ( M + N − 1).Then, by using the real-valued signals that satisfy S r = S * r extracted from the noncircular signals, the extended matrix Z with the dimension 2( M + N − 1) × J is constructed in Equation (17).It implies that the extended data matrix Z corresponds to a virtual array, whose virtual antenna number is twice that of the array corresponding to Y.As the exploitation of the noncircularity of emitted signals enlarges the effective array aperture, the more diversity gain can be obtained in the proposed algorithm.In addition, the reweighted matrix enhances the sparse solution.Thus, the proposed algorithm considerably improves DOA estimation performance in the presence of mutual coupling.Remark 5. Mutual coupling coefficients are related to the distance between two antennas.For uniform linear arrays (ULAs) of the considered MIMO radar system in the proposed algorithm, the mutual coupling coefficients between two antenna elements satisfy c t(q,m) = c t(m,q) = c t|q−m| for the transmit array and c r(q ,n) = c r(n,q ) = c r|q −n| for the receive array [22], in which 1 ≤ q, m ≤ M and 1 ≤ q , n ≤ N. Consequently, all of the effects of mutual coupling in ULA can be modeled as a banded symmetric Toeplitz matrix in Equation (3) [23], based on which selection matrices are constructed in Equation ( 8) to eliminate the mutual coupling.Mutual coupling matrix in Equation ( 3) is modeled under the circumstance of symmetric coupling in ULA.Hence, the proposed algorithm is not suitable for asymmetric coupling.

Simulation Results
In this section, using the ESPRIT-like [14] algorithm, the l 1 -SVD [15] algorithm and the l 1 -SRDML [16] algorithm for comparison, some simulation experiments are implemented to demonstrate the efficiency and the advantages of the proposed algorithm.The signal-to-noise ratio (SNR) is defined as SNR = 10log 10 (||CAS|| 2 F /||N|| 2 F ). Let θi p be the estimation of the true DOA θ p for the i-th Monte Carlo trial; the root mean square error (RMSE) of angle estimation is defined as: (1/P) ∑ P p=1 where Q is the total number of the Monte Carlo trials, which is Q = 200 in the simulations.Consider a narrowband monostatic MIMO radar system, whose arrays are both half-wavelength spaced ULAs with the effects of mutual coupling, and the numbers of transmit and receive antennas are M and N, respectively.Let K = k + 1 be the nonzero mutual coupling coefficient number; two different mutual coupling cases of the transmitter and the receiver are considered in the following simulations: The number of uncorrelated targets P is assumed to be known.Besides, the confidence interval is set to 0.99, and the discretized grids are uniform with 0.1 • sampling from −90 • to 90 • , for the proposed method, as well as l 1 -SVD and l 1 -SRDML algorithms.
When performing the experiments, we solve the signal reconstruction problem of the noncircular source-based joint reweighted sparse representation framework in Equation (25).The sparse solution vector r is obtained by using the SOC programming software packages.When there are true DOAs of targets located at certain points of the complete dictionary, the corresponding elements in r are non-zero, and the rest are zero.Therefore, the true DOA locations in dictionary correspond to the largest P values of |r|.In the simulations, the functional form 10log 10 [|r|/max(|r|)] is used to plot the spatial spectrum of angle estimation.Then, by finding P peaks in the spectrum, DOA estimation is achieved.
Figure 3 depicts the spatial spectrum of the proposed method for the number of targets being P = 3 and P = 5, respectively, where M = N = 6, J = 200, SNR is fixed at SNR = 0 dB and mutual coupling Case (1) is considered.For P = 3, the DOAs of the uncorrelated targets are θ 3, it can be observed that the spatial spectrum peaks of the proposed method are sharp, and the sidelobe suppression is low.In addition, for different target numbers P, the spectrum maintains the accuracy of DOA estimation.This indicates that the proposed method is effective to perform DOA estimation accurately.26) with Q trials.As can be seen in Figure 4, l 1 -SRDML owns a lower RMSE than the ESPRIT-like and l 1 -SVD algorithms, because the deterministic maximum likelihood procedure improves the estimation accuracy.In addition, the proposed method provides the best DOA estimation performance in all SNR regions, for the reason that the noncircularity-based joint reweighted sparse representation scheme enlarges the virtual array aperture and enhances the sparse solution.5, it is clear that the performance of l 1 -SVD is inferior to the other algorithms, and the proposed method has the lowest RMSE in all analyzed methods.This means that when the number of nonzero mutual coupling coefficients increases, the proposed method still has superior DOA estimation performance.Figure 6 shows the RMSE of DOA estimation versus snapshots of different methods with mutual coupling Case (1), where M = N = 6, SNR = 0 dB, there are three uncorrelated targets with DOAs being θ 1 = −11.5 • , θ 2 = 0 • , θ 3 = 11.5 • and the number of snapshots varies from J = 50 to J = 550.It can be seen from Figure 6 that the proposed method performs well in the case of fewer snapshots, and with obvious superiority over the ESPRIT-like, l 1 -SVD and l 1 -SRDML algorithms, the proposed method performs the best in all pf the ranges of snapshots on account of the application of the reweighted matrix and the noncircularity of emitted signals., respectively.They can be regarded as successfully detected when all of their absolute DOA estimation errors are within 0.1 • .It can be observed that when SNR is high enough, all methods provide 100% target resolution probability.Owing to the increased number of virtual array elements, more diversity gain is obtained in the proposed algorithm.Consequently, the proposed method provides higher resolution probability than the ESPRIT-like, l 1 -SVD and l 1 -SRDML algorithms in all SNR regions, as shown in the simulation results of Figure 7.

Conclusions
In this paper, we have proposed the noncircular source-based sparse representation algorithm for DOA estimation in MIMO radar with mutual coupling.On the basis of the error matrix structure and the properties of noncircular signals, the proposed method eliminates the effect of mutual coupling and obtains the extended observation matrix.Then, the joint reweighted sparse representation scheme is constructed to achieve the DOA estimation.The computational complexity of the proposed method has been analyzed, and the simulation results have verified that in the presence of mutual coupling, the proposed algorithm provides higher resolution and better angle estimation performance than the ESPRIT-like, l 1 -SVD and l 1 -SRDML algorithms.

Figure 1 .
Figure 1.Configuration of MIMO radar in the presence of mutual coupling.

Figure 3 .
Figure 3. Spatial spectrum of the proposed method with target number P = 3 and P = 5.

Figure 4
Figure4shows the RMSE of DOA estimation versus SNR in different methods with mutual coupling Case(1), where M = N = 6, J = 200 and three targets withθ 1 = −10 • , θ 2 = 0 • , θ 3 = 10 • are considered.Additionally, the values of RMSE are computed by Equation (26) with Q trials.As can be seen in Figure4, l 1 -SRDML owns a lower RMSE than the ESPRIT-like and l 1 -SVD algorithms, because the deterministic maximum likelihood procedure improves the estimation accuracy.In addition, the proposed method provides the best DOA estimation performance in all SNR regions, for the reason that the noncircularity-based joint reweighted sparse representation scheme enlarges the virtual array aperture and enhances the sparse solution.

Figure 5
Figure 5  illustrates the RMSE of DOA estimation versus SNR in different methods with mutual coupling Case(2), where M = N = 7, J = 200, and there are three uncorrelated targets located atθ 1 = −15 • , θ 2 = 0 • , θ 3 = 15 • .From Figure5, it is clear that the performance of l 1 -SVD is inferior to the other algorithms, and the proposed method has the lowest RMSE in all analyzed methods.This means

Figure 5 .
Figure 5. RMSE versus SNR in the ESPRIT-like, l 1 -SVD and l 1 -SRDML methods for nonzero mutual coupling coefficient number K = 3.

Figure 7
Figure7demonstrates the target resolution probability of different methods versus SNR with mutual coupling Case(1), where M = N = 6, J = 100 and the DOAs of three targets are θ 1 = −20 • , θ 2 = 0 • and θ 3 = 20 • , respectively.They can be regarded as successfully detected when all of their absolute DOA estimation errors are within 0.1 • .It can be observed that when SNR is high enough, all methods provide 100% target resolution probability.Owing to the increased number of virtual array

Figure 7 .
Figure 7. Target resolution probability versus SNR in the ESPRIT-like, l 1 -SVD and l 1 -SRDML methods with snapshot number J = 100.