Weighted Block Sparse Recovery Algorithm for High Resolution DOA Estimation with Unknown Mutual Coupling

Based on weighted block sparse recovery, a high resolution direction-of-arrival (DOA) estimation algorithm is proposed for data with unknown mutual coupling. In our proposed method, a new block representation model based on the array covariance vectors is firstly formulated to avoid the influence of unknown mutual coupling by utilizing the inherent structure of the steering vector. Then a weighted l1-norm penalty algorithm is proposed to recover the block sparse matrix, in which the weighted matrix is constructed based on the principle of a novel Capon space spectrum function for increasing the sparsity of solution. Finally, the DOAs can be obtained from the position of the non-zero blocks of the recovered sparse matrix. Due to the use of the whole received data of array and the enhanced sparsity of solution, the proposed method effectively avoids the loss of the array aperture to achieve a better estimation performance in the environment of unknown mutual coupling in terms of both spatial resolution and accuracy. Simulation experiments show the proposed method achieves better performance than other existing algorithms to minimize the effects of unknown mutual coupling.


Introduction
The main research objectives of array signal processing include angle estimation, or known as direction-of-arrival (DOA) estimation, and adaptive beamforming, often called spatial filtering or adaptive array processing [1].In recent decades, there has been great interest in the study of high resolution DOA estimation which has become increasingly important to applications in electronic countermeasures, medical diagnosis, radar, and communication, etc. [2,3].Meanwhile, inspired by the development of multiple-input multiple-output (MIMO) communication technology, MIMO radar has recently been proposed to obtain more degrees of freedom and high resolution for DOA estimation [4].It has become a hot spot in the fields of autopilot, target location, and medicine [5][6][7].One problem we must consider is that with a fixed array size, as the number of antennas is increased (or the distance between array elements is decreased), this inevitably brings about mutual coupling effect in the received data.In this article, we mainly study DOA estimation with a uniform linear array (ULA) to deal with the effect of unknown mutual coupling.
At present, the existing DOA estimation algorithms for the array are categorized into subspace-based methods and sparse signal recovery (SSR) methods.In the 1980s, the multiple signals classification (MUSIC) algorithm was proposed to start the use of the subspace technique for DOA estimation [8].Then, the algorithm of estimating signal parameters via rotational invariant techniques (ESPRIT) algorithm, was reported [9].The above traditional DOA estimation algorithms are generally known as subspace-based algorithms.Because they are mainly based on eigenvalue decomposition (EVD) or singular value decomposition (SVD) of covariance matrix for DOA estimation, the performance is reduced in the case of low signal to noise ratio (SNR) or a limited number of snapshots [10][11][12][13].In order to deal with the problems associated with traditional DOA estimation methods, the SSR algorithms, such as l 1 -SVD algorithm [14], sparse Bayesian learning (SBL) algorithm [15,16], l 1 -sparse representation of array covariance vector (SRACV) algorithm [17], and their derivative algorithms [18] were proposed in the past few years.The research of SSR algorithms has been mainly focused on two core procedures: the sparse representation model and the reconstruction algorithm.Extensive studies have shown that the performance of SSR algorithms is superior compared to the traditional subspace-based DOA estimation algorithms in the environments of low SNR or a small number of snapshots [19][20][21][22].
It is known that as the distance between array elements decreases, there may exist mutual coupling in the actual array antennas because space electromagnetic fields interact with each other [23].In the above-mentioned DOA estimation algorithms, the steering vector depends on the array structure without considering unknown mutual coupling.Hence, the performance of traditional subspace-based algorithms and SSR algorithms will be severely reduced or even invalid under the condition of unknown mutual coupling [24][25][26].To solve the DOA estimation problem with unknown mutual coupling, extensive research efforts have been made for performance improvement on the basis of traditional DOA estimation algorithms [27][28][29][30].The banded complex symmetric Toeplitz structure of the mutual coupling matrix (MCM) was used to remove the influence of unknown mutual coupling with the help of auxiliary arrays to enable the ESPRIT algorithm to be directly used for DOA estimation [30].By exploiting the banded complex symmetric Toeplitz special structure of MCM in [31], the l 1 -SVD based algorithm was proposed under the condition of unknown mutual coupling.However, both of the above methods in [30,31] have lost the array aperture.In order to make full use of the received data, a new method for avoiding the effect of unknown mutual coupling was reported to be used in the MUSIC algorithm for DOA estimation [32].Although there is no loss of array aperture, it also has the same limitation as that of subspace-based algorithms.For the SSR method, dealing with the unknown mutual coupling, an effective over-complete dictionary is proposed to avoid the influence of unknown mutual coupling [33].Since this method takes advantage of the information of the entire received data, the DOA estimation performance is improved.However, this method is not robust since the DOA estimation performance is based on the length of the received data.Meanwhile, the performance of this method is also limited by the use of l 1 -norm, which is only an approximation of l 0 -norm.
In this paper, a high resolution DOA estimation method based on a weighted sparse recovery algorithm is proposed to deal with the effect of unknown mutual coupling.Firstly, a new block representation model is constructed by using an over-complete dictionary obtained with the parameterized steering matrix to avoid the influence of unknown mutual coupling.The block sparse matrix can be recovered by proposing the weighted l 1 -SRACV algorithm, in which the weighted matrix is obtained by the principle of a novel Capon space spectrum function based on the parameterized steering vector for enhancing the sparsity of solution.Finally, the DOAs can be achieved by the position of the non-zero blocks in the recovered block sparse matrix.Our proposed method not only avoids the loss of the array aperture but also makes full use of the noise statistics.Compared with the l 1 -norm constraint model, the proposed method can obtain a more sparse solution by using a weighted l 1 -norm penalty algorithm.Computer simulations prove the efficiency and robustness of the proposed method.
Notation: [•] mathematical expectation of the matrix or vector.A −1 denotes the inverse of a matrix and B l 2 means that the mth element of a column vector equal to the l 2 -norm of the mth row of matrix B.

Date Model
Suppose N narrowband far-field source signals {s n (t)} N n=1 , impinging on a uniform linear array (ULA) with M array elements, where t is the time index.The spacing between two adjacent array elements is generally assumed to be d ≤ λ 2 , where λ indicates the signal wavelength.Because the signals come from different and unknown directions, the DOA is denoted as • to 90 • .If the first array element is regarded as a reference point and there is no unknown mutual coupling, the ideal array output can be expressed as where denotes array manifold matrix, and its nth column vector a(θ In practice, when the adjacent array elements are relatively close, the unknown mutual coupling appears due to the interaction of space electromagnetic field.In order to achieve more accurate estimation of DOAs, the effect of unknown mutual coupling has to be properly dealt with.Because the structure of the array steering vector is changed by the unknown mutual coupling, the actual array steering vector takes the following form: where D is the mutual coupling matrix (MCM) whose element magnitude is related to the distance between array sensors.In this way, the effect of mutual coupling between two closely located array sensors is considered.Based on the above analysis in ULAs, a banded symmetric Toeplitz matrix [10] is adopted to represent the mutual coupling matrix (MCM) in ULA.So D is a banded symmetric Toeplitz matrix, whose specific form is: where represents the value of the mutual coupling coefficient between two sensors.When the distance between two array elements is less than or equal to K inter-sensor spacing, the effect of unknown mutual coupling between them cannot be overlooked.Then the array output expression with unknown mutual coupling is converted to: which shows that the steering vector has been changed greatly due to the existence of unknown mutual coupling.According to [32], a series of transformations are carried out for the steering vectors with unknown mutual coupling.Then the transformed steering vector can be formulated as the following: where, where T(θ) is a constant for each DOA, and the magnitude of its value is related to the mutual coupling.
In this paper, we suppose T(θ) = 0, and where v(θ) is a (2K − 1) × 1 column vector with the Kth element being one, and, where Hence, Equation ( 4) can be further transformed into: where, According to Equation (11), it is known that there are no mutual coupling coefficients in matrix A J , which can be regarded as a new array manifold matrix, similar to A, and where and Q = 2K − 1. Equation (11) only considers the signal reception model at time t.If we consider the total received data of T snapshots, Equation ( 11) can be rewritten as: where Y = [y(t 1 ), y(t 2 ), In this section, the l 1 -SRACV algorithm is given for DOA estimation in ULA.Through the study of the new data model in Equation (11), it is found that the block diagonal matrix Λ can be combined with the unknown signal s(t) to form a new signal vector s (t).Hence, Equation ( 11) can be further rewritten as: where s (t) = Λs(t) is a NQ × 1 new signal vector, the nth block matrix of s (t) is from (Qn − Q + 1)th to (Qn)th rows of s (t) corresponding to the nth element of signal vector s(t).Therefore, the unknown mutual coupling coefficients are added to the part of unknown signals for avoiding the influence of mutual coupling on DOA estimation.By using Equation ( 15), the covariance matrix of array output can be expressed as: where R s = E[s (t)s H (t)] is the covariance matrix of the new signal s (t).In addition, σ 2 is the power of noise and σ 2 I M denotes the noise covariance matrix.At the same time, it can also be shown that R s A J H and s (t) have the same block sparse structure [17].However, in practice, only limited samples are available.Therefore, the covariance matrix is estimated by using the T available snapshots, i.e., In order to exploit the view of SSR based on Equation ( 17), an over-complete dictionary is firstly represents all possible signal directions and G is the total number of all possible DOAs.According to the knowledge of linear algebra, each column vector of covariance matrix R can be represented linearly by the over-complete dictionary.Then, the mth column of covariance matrix R can be represented as: where y m (t) represents the received data of the mth array sensor at time t.
] is a given over-complete dictionary.And {∆( θ g )} G g=1 is the M × Q block matrix, which has the same structure as ∆(θ) in Equation (7).Suppose that the grid points are sufficiently dense; hence, the (Qg − Q + 1)th to (Qg)th rows of ideal b m are non-zero, the position of the non-zero blocks in b m corresponds to the N true DOAs.Meanwhile, e m (m = 1, • • • , M) is a M × 1 column vector, the mth element of vector e m is one, and the rest are zero.However, Equation (18) only considers one column of the covariance matrix.For a covariance matrix, Equation ( 18) can be written as the form of a matrix, and the sparse representation model of the covariance matrix can be obtained as: where From the above analysis, it is obviously shown that the sparsity of vector b l 2 can be used to describe the joint block-sparsity of matrix B. Hence, the DOAs can be achieved by detecting the position of non-zero elements of sparse vector b l 2 .
In sparse representation theory, the constraint recovery model that accurately describes the sparsity of signal is the l 0 -norm constraint.It is known that l 0 -norm optimization is a computationally expensive NP-hard problem, and the l 1 -norm penalty optimization is a convex approximation of l 0 -norm optimization.For reducing the computational complexity, we transform the DOA's estimation problem into a convex optimization.Therefore, the sparse vector b l 2 can be recovered by the following constrained optimization problem: where η represents the regularization parameter, which is used to balance the relationship between data fitting and sparsity of the solution in the l 1 -norm penalty optimization problem.The value of regularization parameter η is related to noise.When the recovered matrix is obtained, the DOAs can be achieved by the position of the non-zero blocks.

A Reweighted l 1 -SRACV Algorithm for DOA Estimation
Since the l 1 -norm is a convex approximation of l 0 -norm, there exists a difference between the recoveries obtained by these two operations.In the case of the l 0 -norm constraint, the solution for any vector only has a value of zero and one without differentiating the contributions from large and small coefficients to the objective function.To achieve better recovery performance in the l 1 -norm constraint, more constraints should be imposed on the large coefficients to objective function for ensuring the sparse solution of the entire cost function.To avoid this problem, a weighted norm constraint is introduced.By reweighting the l 1 -norm constraint model, both large and small coefficients in the reconstructed signal can be equally constrained as much as possible.That is, the smaller weight penalties are used for large coefficients and the larger weight penalties are used for small coefficients.This adaptive adjustment mechanism enables the reweighted l 1 -norm constraint optimization to better approximate the l 0 -norm penalty minimization.In the following section, the basic principle of Capon algorithm is used to construct a weighted matrix for the reweighted l 1 -norm constraint optimization The traditional Capon algorithm is firstly introduced, which is based on the steering vector of the signal and the covariance matrix of the source signal.We first find the minimum value of the following function, and then determine the DOAs with the minimum value for the following function: where a(θ) denotes the steering vector of source signal and R −1 represents the inverse of the covariance matrix of the sample.Therefore, f(θ) will tend to be zero when θ corresponds to the true DOAs.Following the ideal in [25], it is known that the relationship between the new over-complete dictionary matrix A J ( θ g ) and the covariance matrix can be used to construct the weighted matrix.Then the over-complete dictionary can be expressed as follows, Suppose that A J1 ( θ g ) is made up of the N block matrices corresponding to real DOAs, and A J2 ( θ g ) is composed with G − N residual block steering matrices.Now a new function is introduced in the following form: where ŵg denotes the value of the determinant corresponding to the different possible signals.Then the weights can be formulated as follows: where ŵl 2 g denotes the l 2 -norm of ŵg .And we define ].In Equation ( 23), the weights ŵg are more likely to be zero when {θ g } G g=1 corresponds to the true DOAs if the number of samples is sufficiently large.Ŵ1 is assumed to correspond to the true DOAs.Therefore, the weights in Ŵ1 satisfies Ŵ1 → 0 , and the weights in vector Ŵ1 are smaller than that in vector Ŵ2 .Finally, the weighted matrix is defined as: where W is a G × G diagonal matrix.When the weighted matrix is applied to the above l 1 -norm minimization problem, the small weights in W reserve the large coefficients and the large weights punish the small coefficients who tend to be zero.Such a constraint effectively guarantees the fidelity of signal reconstruction.Through the above-mentioned study, the weighted matrix can enforce the sparsity of the solution.The new reweighted l 1 -norm minimization problem is as follows: where the regularization parameter η sets the upper limit of noise power that is to be permitted, which will be introduced in detail in Remark 1. Finally, the above optimization problem can be solved by the second order cone (SOC) programming packages such as CVX.After obtaining the sparse vector b l 2 , DOA estimation of real signals can be obtained by the angle position corresponding to the position of the non-zero values in the sparse vector b l 2 .
Remark 1.In the reweighted l 1 -norm minimization problem, the regularization parameter η plays an important role in recovering the sparse vector b l 2 accurately.The regularization parameter η can also compromise the sparseness of fitting error and resolution.Through the introduction of [14], we find that the selection of parameters is related to the probability distribution of noise covariance.Several studies have shown that the noise approximately follows a χ 2 distribution with M × N degrees of freedom [17].The upper limit value of noise power can be computed as the regularization parameter η with a high probability of 99%.The regularization parameter η can be calculated by chi2inv(p, M × N) with the MATLAB.Furthermore, the noise power, that is related to the singular value of the covariance matrix, can be calculated by averaging the square of the M − N smallest singular values.
Remark 2. The main computational complexity of the proposed method is for the construction of the weighted matrix and the solution of the convex optimization problem in Equation (26).The construction of the weighted matrix W requires O{M 2 T + GQ 4 M(M − N)(M + Q)} flops, and solving the convex optimization problem requires O{(MQG) 3 } flops.Therefore, the proposed method requires more computation than Dai's method [31] and Wang's method [33].Nevertheless, the advantages of the proposed method outweigh the disadvantages.The proposed method can work well in the environment of unknown mutual coupling, and the performance of our proposed method is much better than other methods in terms of accuracy and resolution.

Simulation Results
In this section, the DOA estimation performance of the proposed method is demonstrated by our simulation experiments.The methods in [31,33], referred to as Dai's method and Wang's method, respectively, and the Cramer-Rao bound (CRB) [17] are chosen to compare with our proposed method.In Dai's method, a selection matrix is used to remove the influence of an unknown mutual coupling coefficient and in Wang's method, a novel data model is constructed to avoid the influence of unknown mutual coupling by parameterizing the steering vector.Furthermore, the DOA estimation performance is evaluated by the root mean square error (RMSE) defined as: where θ l,i denotes the estimated value of θ l obtained by the lth Monte Carlo trial, L = 200 is the number of Monte Carlo trials, and N is the total number of source signals.
In the following simulation experiments, we consider a ULA equipped with M = 8 sensors that are separated by a half-wavelength.We assume that there are two uncorrelated narrowband far-field signals whose DOAs are θ 1 = −6 • and θ 2 = 10 • .The mutual coupling coefficients are set as d 1 = 0.6864 − j0.1776, d 2 = 0.20 − j0.0896, and K = 3.The range from −90 • to 90 • in the spatial domain is discretized with a 0.05 • grade resolution for all methods.
Figure 1 represents the spatial spectrum of the proposed method, Wang's method and Dai's method, for two uncorrelated source signals with M = 8, T = 300 and SNR = 10 dB.From Figure 1, it can clearly be seen that our proposed method achieves a spectrum that has a lower side-lobe and sharper peaks than those achieved by the other two methods.Therefore, our method has superiority with regard to the spatial resolution.It can be seen clearly from Figure 1 that our proposed method obtains two sharp peaks located very close to the true DOAs; therefore, the proposed method is more accurate than Wang's method and Dai's method.It can be seen that our proposed method achieves more accurate estimation with a lower RMSE than other two methods.In addition, the RMSE value achieved by our proposed method is closer to CRB.The main reason is that the weighted matrix is used to ensure that the recovery solution is sparser than other methods.Meanwhile, the proposed method takes full advantage of statistical characteristics.However, Wang's method has a loss of the array aperture and the performance of Dai's method is influenced by the data length used for DOA estimation.Therefore, these two methods have a higher RMSE than the proposed method.   .It is seen that as the number of snapshots increases, the RMSE values of all methods are reduced and the performance is improved.However, it is evident that the estimation performance of our proposed method is an improvement compared to those achieved by Dai's method and Wang's method; its RMSE value is closest to CRB.The performance of the proposed method is much higher than that achieved by the other two methods in the entire range of the number of snapshots.We also observe that the performance of Wang's method is slightly better than that of Dai's method.
Figures 4 and 5 show the DOA estimation performance of the proposed method with different numbers of antennas versus different SNRs or number of snapshots, respectively.From both Figures 4 and 5, we observe that when SNR or the number of snapshots are fixed, the DOA estimation performance of the proposed method is improved through increasing the number of the array elements because a greater diversity gain is achieved as the number of antennas is increased.However, it should be noted that the computational complexity is also increased with the increasing number of antennas for better estimation performance.Figure 3 shows the RMSE versus the number of snapshots achieved by different methods for M = 8 and SNR = 10 dB.It is seen that as the number of snapshots increases, the RMSE values of all methods are reduced and the performance is improved.However, it is evident that the estimation performance of our proposed method is an improvement compared to those achieved by Dai's method and Wang's method; its RMSE value is closest to CRB.The performance of the proposed method is much higher than that achieved by the other two methods in the entire range of the number of snapshots.We also observe that the performance of Wang's method is slightly better than that of Dai's method.and 5, we observe that when SNR or the number of snapshots are fixed, the DOA estimation performance of the proposed method is improved through increasing the number of the array elements because a greater diversity gain is achieved as the number of antennas is increased.However, it should be noted that the computational complexity is also increased with the increasing number of antennas for better estimation performance.Figures 6 and 7 reveal the probability of successful detection versus SNR and the number of snapshots for different methods, respectively.In these two simulation experiments, the successful detection rate is defined as the error between the DOAs estimation and the real DOAs that satisfies θn − θ n ≤ 0.5 • , where θn is the estimated value of θ n .It can be seen from Figures 6 and 7, as the number of snapshots or SNR increases, the probability of detection also increases.When the SNR or the number of snapshots is greater than a certain value, the detection probability of the proposed method reaches 100%.When the SNR is low or the number of snapshots is small, the detection probability of the other two methods is relatively lower.Therefore, the detection probability of the proposed method is much higher than that achieved by the other two methods in the entire range of SNRs or the number of snapshots.
detection rate is defined as the error between the DOAs estimation and the real DOAs that satisfies ˆ0.5 n n θ θ − ≤ °, where ˆn θ is the estimated value of n θ .It can be seen from Figures 6 and 7, as the number of snapshots or SNR increases, the probability of detection also increases.When the SNR or the number of snapshots is greater than a certain value, the detection probability of the proposed method reaches 100% .When the SNR is low or the number of snapshots is small, the detection probability of the other two methods is relatively lower.Therefore, the detection probability of the proposed method is much higher than that achieved by the other two methods in the entire range of SNRs or the number of snapshots.

Conclusions
In this article, a weighted block sparse recovery algorithm for high resolution DOA estimation under unknown mutual coupling is proposed.In the proposed method, a block sparse representation model based on the array covariance vectors is formulated by parameterizing the steering vector without the influence of mutual coupling, and then the weighted matrix is constructed by the principle of a novel Capon space spectrum estimation algorithm based on the detection rate is defined as the error between the DOAs estimation and the real DOAs that satisfies ˆ0.5 n n θ θ − ≤ °, where ˆn θ is the estimated value of n θ .It can be seen from Figures 6 and 7, as the number of snapshots or SNR increases, the probability of detection also increases.When the SNR or the number of snapshots is greater than a certain value, the detection probability of the proposed method reaches 100% .When the SNR is low or the number of snapshots is small, the detection probability of the other two methods is relatively lower.Therefore, the detection probability of the proposed method is much higher than that achieved by the other two methods in the entire range of SNRs or the number of snapshots.

Conclusions
In this article, a weighted block sparse recovery algorithm for high resolution DOA estimation under unknown mutual coupling is proposed.In the proposed method, a block sparse representation model based on the array covariance vectors is formulated by parameterizing the steering vector without the influence of mutual coupling, and then the weighted matrix is constructed by the principle of a novel Capon space spectrum estimation algorithm based on the

Conclusions
In this article, a weighted block sparse recovery algorithm for high resolution DOA estimation under unknown mutual coupling is proposed.In the proposed method, a block sparse representation model based on the array covariance vectors is formulated by parameterizing the steering vector without the influence of mutual coupling, and then the weighted matrix is constructed by the principle of a novel Capon space spectrum estimation algorithm based on the parameterized steering vector.Finally, the DOA is estimated from the position of the non-zero blocks of the recovered block sparse matrix.Due to the fact that our proposed method takes advantage of the whole array aperture and enhances the sparsity of solution by the reweighted l 1 -norm minimization scheme, the proposed method not only works well but also has a higher spatial resolution.The simulation results have demonstrated that the performance of our proposed method is better than Wang's method and Dai's method, both in resolution and accuracy.

2 =[b l 2 1 , b l 2 2 ,denotes the gth element of b l 2 ,
B = [b 1 , b 2 , • • • , b M ].Through the above analysis, it can be obtained that all the ideal {b m } M m=1 have the same block sparse structure.In other words, non-zero block vectors of each ideal {b m } M m=1 should appear in the same position of matrix B. Here, for convenience, a new column vector b l • • • , b l 2 G ] which is equal to the l 2 -norm of the (Qg − Q + 1)th to (Qg)th rows of matrix B. Therefore, b l 2 is a sparse vector of N and satisfies b

13 Figure 1 .
Figure 1.The spatial spectrum of the proposed method, Wang's method, and Dai's method.Figure 1.The spatial spectrum of the proposed method, Wang's method, and Dai's method.

Figure 1 .
Figure 1.The spatial spectrum of the proposed method, Wang's method, and Dai's method.Figure 1.The spatial spectrum of the proposed method, Wang's method, and Dai's method.

Figure 2
Figure 2 depicts the RMSE versus SNR for different methods with M = 8 and T = 300.It can be seen that our proposed method achieves more accurate estimation with a lower RMSE than other two

Figure 1 .
The spatial spectrum of the proposed method, Wang's method, and Dai's method.

Figure 2 .
Figure 2. The root mean square error (RMSE) versus signal to noise ratio (SNR) for different methods ( 8 M = and 300 T =).

Figure 3
Figure 3 shows the RMSE versus the number of snapshots achieved by different methods for 8 M = and SNR B 10d =.It is seen that as the number of snapshots increases, the RMSE values of all methods are reduced and the performance is improved.However, it is evident that the estimation performance of our proposed method is an improvement compared to those achieved by Dai's method and Wang's method; its RMSE value is closest to CRB.The performance of the proposed method is much higher than that achieved by the other two methods in the entire range of the number of snapshots.We also observe that the performance of Wang's method is slightly better than that of Dai's method.Figures4 and 5show the DOA estimation performance of the proposed method with different numbers of antennas versus different SNRs or number of snapshots, respectively.From both Figures4 and 5, we observe that when SNR or the number of snapshots are fixed, the DOA estimation performance of the proposed method is improved through increasing the number of the array elements because a greater diversity gain is achieved as the number of antennas is increased.However, it should be noted that the computational complexity is also increased with the increasing number of antennas for better estimation performance.

Figure 2 .
Figure 2. The root mean square error (RMSE) versus signal to noise ratio (SNR) for different methods (M = 8 and T = 300).

Figures 4
Figures4 and 5show the DOA estimation performance of the proposed method with different numbers of antennas versus different SNRs or number of snapshots, respectively.From both Figures4 and 5, we observe that when SNR or the number of snapshots are fixed, the DOA estimation performance of the proposed method is improved through increasing the number of the array elements because a greater diversity gain is achieved as the number of antennas is increased.However, it should be noted that the computational complexity is also increased with the increasing number of antennas for better estimation performance.
and diag{•} denote the transpose, conjugate-transpose, conjugate, and diagonal matrix, respectively.I M represents an M × M dimensional unit matrix and • p represents the p-norm of a matrix or vector.det[•]isthe determinant of a matrix and E[•] represents the T , [•] H , [•] *