Novel Diagonal Reloading Based Direction of Arrival Estimation in Unknown Non-Uniform Noise

Nested array can expand the degrees of freedom (DOF) from difference coarray perspective, but suffering from the performance degradation of direction of arrival (DOA) estimation in unknown non-uniform noise. In this paper, a novel diagonal reloading (DR) based DOA estimation algorithm is proposed using a recently developed nested MIMO array. The elements in the main diagonal of the sample covariance matrix are eliminated; next the smallest MN-K eigenvalues of the revised matrix are obtained and averaged to estimate the sum value of the signal power. Further the estimated sum value is filled into the main diagonal of the revised matrix for estimating the signal covariance matrix. In this case, the negative effect of noise is eliminated without losing the useful information of the signal matrix. Besides, the degrees of freedom are expanded obviously, resulting in the performance improvement. Several simulations are conducted to demonstrate the effectiveness of the proposed algorithm.


Introduction
Direction of arrival (DOA) estimation, which has been studied for decades, has been widely used in target localization, wireless communication, and so on.Solving the non-uniform noise problem is inevitable and important for improving the DOA estimation performance, no matter for the DOA estimation algorithm or array structure design.
For the DOA estimation algorithm, many classical methods like the Multiple Signal Classification (MUSIC), Estimation Signal Parameters via Rotational Invariance Techniques (ESPRIT), Maximum Likelihood (ML), and so on have been proposed [1].Also, with the introduction of sparse recovery theory, the DOA estimation has come into a new era and algorithms such as Khatri-Rao MUSIC (KR-MUSIC) and Spatial Smoothing MUSIC (SS-MUSIC) are developed to improve the DOA estimation performance [2].Most of these algorithms are based on the assumption that the received noise is uniform white Gaussian noise.However, in practical application, this assumption can hardly be guaranteed; thus the estimation performance will degrade.
For the array structure design, the structure of sensors has also been well studied to expand the degrees of freedom (DOF).Based on the subspace algorithm, traditional uniform linear arrays with  elements can resolve −1 targets [3].The minimum redundancy (MR) array aims to achieve optimal structure by reducing the redundancy, but it needs complex optimization and has no closed-form expression for sensors positions [4,5].The Multiple Input Multiple Output (MIMO) structure can effectively expand the DOFs by generating a virtual array [6].But it only uses the sum co-array and there is still space to improve the DOFs.The nested array [7] and the coprime array [8] are thus proposed, which brings obvious increase of the DOFs.Taking the two level nested arrays as an example, it can obtain ( 2 ) DOF using only () physical sensors.The aforementioned MR array [9] and coprime array [10] are combined with the MIMO array to generate the difference co-array of the sum coarray (DCSC), which further improves the DOFs greatly.Recently, a novel structure of this type, which is called nested MIMO array, is proposed [11].It not only expands the DOFs obviously, but also has closed-form expressions.However, the non-uniformity problem of the received noise occurs with the repaid development of these new structures.For closely spaced arrays whose inter-element spacing is less than or equal to half the wavelength, this problem is not serious and we can assume the noise obeys white Gaussian distribution.However, as the sensors of sparse arrays are usually widely separated, they receive noise with different power.In some special circumstance such as low altitude, the non-uniformity may be severe and will degrade the performance of classical algorithms.Therefore, it is important to solve the nonuniform noise problem for sparse array.
ML methods are classical algorithms to deal with the nonuniform noise problem in DOA estimation.The ML methods can avoid the eigen-decomposition, and the corresponding noise matrix can be estimated and eliminated to improve the accuracy of DOA estimation.A deterministic non-uniform ML estimator is developed in [12], and the Cramér-Rao bounds (CRB) of both the deterministic and stochastic signal model are given.However, the stochastic ML DOA estimator is not derived in the literature.On the basis of this paper, Chen et al. propose a stochastic ML estimator in non-uniform noise [13].All of the ML methods need iterative optimization, which brings high computation complexity.Though Liao B. et al. improve the method to reduce the computational complexity [14,15], stepwise concentration is still essential and time consuming.Besides, there is a chance that the obtained result is a local optimum but not a global one.
Subspace based algorithms are the dominant methods in DOA estimation.Many methods are proposed to solve the non-uniform noise problem by eliminating the main diagonal elements or averaging these main diagonal elements [15].Based on the matrix completion theory, the estimation of the noise-free covariance matrix is solved by the rank minimization [16], and nuclear norm minimization is employed to transform the non-convex problem to a convex one.This method avoids the iterative optimization; thus it is computationally efficient.In [17], the covariance matrix is divided into some sub-matrices and the noise covariance is estimated based on the relationship between the divided submatrices.However the number of detectable targets is limited.
In recent years, the DOA estimation based on sparse recovery method has aroused a lot of attention [18,19].To solve the non-uniform noise problem in sparse signal recovery, He Z. et al. propose a novel covariance sparsity aware (CSA) DOA estimation method in [20], which directly removes the unknown noise variances by a linear transformation.This method improves the accuracy, but it also causes the loss of effective aperture.
In summary, the current methods have not completely settled the non-uniform noise problem.Based on the above analysis, a novel diagonal reloading (DR) DOA estimation algorithm for nested array is proposed in unknown nonuniform noise.The elements in the main diagonal of the received data covariance matrix are eliminated, and the smallest MN-K eigenvalues of the revised matrix are obtained and averaged to estimate the sum of the signal power.Then the estimated sum is filled into the main diagonal of the revised matrix for estimating the signal covariance matrix.Different from the method of completely removing the received data which contains noise, the proposed method only eliminates the negative effect of noise while keeping the useful information of the signal matrix.
The remainder of this paper is organized as follows.The nested MIMO radar structure and the corresponding signal model are introduced in Section 2. The proposed method is derived in Section 3. Several simulations are conducted in Section 4 to validate the performance of the proposed method and conclusions are drawn in Section 5.

Signal Model for Nested MIMO Radar
Suppose the transmitter and receiver of the MIMO radar are all parallel linear arrays with the sensor numbers  and , respectively.The code length and wavelength of each signal are  ( < ) and , respectively.S = [s 1 , s 2 , . . ., s  ]  is the transmitting signal matrix, among which s  = [s  (1), s  (2), . . ., s  ()]  is the th transmitting signal.[⋅]  denotes the transpose.The number of far field targets is  and their DOAs are   ( = 1, 2, . . ., ).The received signal vector of a single snapshot is where ] are the transmitting and receiving steering vector respectively.  (  ) and   (  ) are functions with respect to   for given coordinates of the sensors.U = [u 1 u 2 . . .u  ] ∈ C × denotes the noise matrix, and The noises are assumed to be nonuniform and uncorrelated with each other, namely, Conducting matched filtering with the  orthogonal transmitting waveforms in each receiver, we obtain that where vec(⋅) is the operation to rearrange matrix into one column and (⋅)  denotes Hermitian transpose.As the transmitting signals are orthogonal with each other, z() can be simplified as wherein A tr = [a tr ( 1 ), a tr ( 2 ), . . ., a tr (  )] is the steering matrix of the sum co-array, and ⊗ denotes the Kronecker product.() = [ 1 ,  2 , . . .,   ]  is the reflection coefficient of the  targets.k = vec(US  /) is the output of noise after matched filtering.Define k  as the element in the th row and th column of  ×  matrix US  /, i.e., As both u  and s  are unique vectors, namely, . ., ,  = 1, 2, . . ., ) is unique, too.Therefore, the noise remains non-uniform after matched filtering.
When the radar transmits  pulses, we have where z(), ( = 1, 2, . . ., ) is the receiving signal corresponding to the th transmitting pulse, that is, As a result, the total received signal can be expressed as where B = [(1), (2), . . ., ()] is the  ×  matrix of the reflection coefficient and V = [k(1), k(2), . . ., k()] is the  ×  matrix of the noise.Since the noise is uncorrelated with each other, the covariance of the array is where [⋅] denotes the mathematical expectation.
where ⊙ denotes the Khatri-Rao product operation, and b = [ 2 1 ,  , where e  ( = 1, . . ., ) is  × 1 column vector of all zeroes except a 1 in the th entry.The obtained difference co-array is Let Ã() = A * tr () ⊙ A tr ().Compared with A(), Ã() expands the effective DOFs by constructing the difference coarray, whose virtual sensors' positions are expressed as where 1 ≤  1 ,  2 ≤ , 1 ≤  1 ,  2 ≤ ,  , and  , are the sensors' positions of transmitting array and receiving array, respectively.Now the question turns to maximize the number of unique difference results in (11).Here we adopt a new structure of nested MIMO array to settle the problem [11].
First, the geometry of the receiver is designed as a two level nested array, whose locations are where {  }  =1 is the location of the nested array.This is a hole free array after difference.The number of DOFs in the difference co-array of the transmitting array is   , which is summarized as The receiving array is constructed as another nested array with a large unit inter-element spacing , whose locations are denoted by {  }  =1 is the location sequence of the nested array.In all, the sum co-array of the MIMO radar can be expressed as The location of the difference co-array is where ,   = 1, 2, . . ., , ,   = 1, 2, . . ., .If  =   , a hole free ULA is obtained with expanded DOFs.If the number of sensors satisfies then the best values of and  are Correspondingly, we obtain the optimal geometry of the transmitting and receiving nested array.For example, if  is even,  =  = /2, the final obtained DOF is We can find that, by this kind of structure, the nested array can provide ( 4 ) DOFs with only  sensors.Figure 1 shows the structure of the nested MIMO radar array with a fourelement transmitting (Tx) array and a four-element receiving (Rx) array.The numbers under the sensors indicate their positions.It should be mentioned that only the nonnegative part of the DCSC is given in this figure.Compared with the former structure, the number of DOFs is greatly expanded using the same number of sensors, and the DOA estimation accuracy will be highly promoted.The application of nested MIMO array effectively improves the virtual aperture; however, it also brings another problem.Due to the large spacing between the receiving array, the received noise of different sensors no more satisfy the uniform assumption.Therefore Q in (8) will turn to a non-uniform diagonal matrix and the DOA estimation performance will degrade significantly.In this paper, we will analyze the negative effect of non-uniform noise and propose a novel diagonal recovery algorithm to solve this problem.

Diagonal Reloading Based DOA Estimation Algorithm
3.1.Negative Effect Analysis of Non-Uniform Noise.In ideal situation, there is no noise.So the matrix Q in (8) shall be equal to 0. As a result, R = A tr R  A  tr .Conducting eigendecomposition with A tr R  A  tr , we have where U 1 consists of eigenvectors corresponding to the  largest eigenvalues and U 2 is consisted of the rest of  −  eigenvectors corresponding to the other  −  smaller eigenvalues.However, noise is inevitable and will degrade the DOA performance.It can be explained from the structure of the covariance.Let a ℎ ≜ [ ℎ ( 1 ),  ℎ ( 2 ), . . .,  ℎ (  )] represent the ℎth row of A tr . ℎ (  ) =   ℎ () , wherein  ℎ (  ) is function with respect to   of the virtual sensors.Therefore | ℎ (  )| 2 = 1, and the elements in the main diagonal of A tr R  A  tr can be expressed as It can be found that [A tr R  A  tr ] ℎ,ℎ equals the sum value of all targets' power.As Q is a diagonal matrix, it will be added to the main diagonal of A tr R  A  tr and only affect the element It can be found that the elements in the main diagonal remain uniform, so it can be seen as the change of signal power and the structure of the main diagonal has not been changed.The matrix U 1 corresponding to the eigenvectors associated with the  largest eigenvalues of R remains the same as before.It is the same with the noise subspace U 2 .However, when the noise is non-uniform, the elements in the main diagonal no longer equal to each other.As Q =  2 I  no longer holds, A tr R  A  tr + Q no longer equals to U(Σ +  2 I  )U  .Correspondingly, the signal subspace U 1 and the noise subspace U 2 will be changed and the performance of DOA estimation will degrade significantly.

Diagonal Reloading (DR) Algorithm.
In [15], the diagonal elements of the covariance matrix R are completely eliminated (as shown in ( 23)).As the covariance of noise is a diagonal matrix, it will be removed from R. However, the main diagonal elements of the signal covariance will also be eliminated: where (A tr R  A  tr ) denotes the operation to construct a diagonal matrix whose main diagonal elements are the same as A tr R  A  tr .The main diagonal contains the power of target signal.If they are replaced with zeroes as (23), it will certainly lose valuable information.An improved method which only eliminates noise and keeps the elements of signal is thus proposed.
Set  = ∑  =1  2  ; there is Therefore, Ȓ can be expressed as We can find that the smallest  −  eigenvalues corresponding to U 2 happen to be −, which means that  can be estimated by the average of  −  smallest eigenvalues.
Based on this principle, we first conduct eigen-decomposition of Ȓ as (25) and get an estimation of , then the elements in main diagonal of Ȓ can be replaced with γI  , namely, where γ is the estimation of .R is an estimate of A tr R  A  tr .If the estimation of  is accurate, namely, γ = , then R = A tr R  A  tr , and the covariance of the signal matrix can be accurately reconstructed.
In conclusion, the proposed method removes the negative influence of noise while keeping the main diagonal elements of signal matrix.

Sparse Signal Recovery.
It can be seen that the effective DOFs have obvious improvement after vectorization, leading to a more precise estimation.Now we can divide the whole space to grids with a series of angles θ1 , θ2 , . . ., θ  .The grid angles are sufficiently dense and each of them represents a possible target direction.Since the actual targets only appear in a few angles and the coming signals are sparse in the whole space, the nonzero elements among the whole grid points correspond to the real DOAs.Set Ã( θ) = [a( θ1 ), a( θ2 ), . . ., a( θ  )] ∈ C ×  as the over-complete basis, where  is given by (19).The covariance can be expressed in a sparse way as where b  = [ 1  2 . . .   ]  is the sparse vector with rank .
where  is the penalty factor and  is the sparse set we tend to solve.In summary, the steps of the proposed methods are as follows: (1) Construct covariance matrix of the received signal and remove the elements in the main diagonal as (23).
(2) Conduct the eigen-decomposition of Ȓ and use the smallest  −  estimated eigenvalues to construct matrix γI  .

Simulation
The first simulation is to validate the ability of MIMO and nested and MIMO nested array in distinguishing underestimated targets.The transmitter array and receiver array both have 4 sensors.We suppose that the targets are of the same power  2 .The noise is modeled as white Gaussian process of power  2  ( = 1, 2, . . ., ).To express the non-uniformity of noise, the worst noise power ratio (WNPR) is defined as WNPR =  2 max / 2 min .As (3) shows, after matched filtering, the array has  ×  = 16 distinct virtual receivers; the corresponding noise power of each virtual receiver is [30, 9,3,6,5,9,7,5,1,8,5,6,8,9,7,2].So in this simulation, the WNPR = 30.The average SNR is defined as Here we set the SNR as 0dB.  2. The MIMO array and the nested array cannot distinguish all the ten targets, but for nested MIMO array, all of the targets can be estimated accurately.
Next, several simulations are conducted under different conditions to examine the effectiveness of the proposed DR method.All of the results are obtained using 500 independent trials.
Simulation 2 compares the DOA estimation performance of different methods versus SNR.The number of targets and their directions are as simulation 1.The WNPR equals 30 and Figure 3 shows all the RMSEs of DOA estimation decrease with the increase of SNR.It can be found that the KR-MUSIC performs the worst, while the proposed method has obvious improvement in accuracy, especially in low SNR region.It is because that the power of noise is equivalent or even stronger than that of signal in low SNR condition.The nonuniformity of noise will obviously affect the results of eigendecomposition.Simulation 3 compares the DOA estimation performance of different methods versus the number of snapshots.The number of targets and their directions are same as simulation 1. SNR = −5 dB and the WNPR equals 30. Figure 4 shows the RMSEs versus the number of snapshots.With the increase of the number of snapshots, the accuracy of all these methods gets improved.Besides the proposed method outperforms the KR-MUSIC and CSA method.Simulation 4 examines the DOA estimation performance of different methods under different levels of non-uniform noise.We assume that if the maximum DOA estimation error does not exceed 1 ∘ , the estimation is successful, and otherwise it fails.It can be seen from Figures 5 and 6 that, with the increase of WNPR, the RMSEs of these methods increase and the success probabilities decrease.In addition, compared with KR-MUSIC and the covariance sparsity aware (CSA) method in [20], the proposed method has better performance against the increase of WNPR.

Conclusion
In this paper, we proposed a novel DR based DOA estimation algorithm to solve the non-uniform noise problem for sparse array.First the elements in the main diagonal of received data covariance matrix were eliminated; then the smallest MN-K eigenvalues of the revised matrix were obtained and averaged to estimate the sum of the signal power.Next the estimated sum was filled to the main diagonal of the revised matrix.In this way, the covariance of signal matrix was recovered, where the negative effect of noise was eliminated without losing the useful information in the main diagonal of signal covariance matrix.Therefore the performance of proposed method was improved.Simulation results revealed that the proposed method not only made fully use of the DOF expansion of nested MIMO radar, but also suppressed the its non-uniform noise effect.Compared with the KR-MSUIC and the CSA method, the proposed method had better accuracy and higher success ratio.Furthermore, it could also be applied to arrays with closely spaced sensors.

Figure 1 :
Figure 1: Structure of the nested MIMO radar array with a four-element transmitting (Tx) array and a four-element receiving (Rx) array ( = 11).

Figure 6 :
Figure 6: Success ratio of different DOA estimation methods versus WNPR.