Compressive Sensing for High-Resolution Direction-of-Arrival Estimation via Iterative Optimization on Sensing Matrix

A novel compressive sensing(CS-) based direction-of-arrival (DOA) estimation algorithm is proposed to solve the performance degradation of the CS-basedDOAestimation in the presence of sensingmatrixmismatching. Firstly, a DOA sparse sensingmodel is set up in the presence of sensing matrix mismatching. Secondly, combining the Dantzig selector (DS) algorithm and least-absolute shrinkage and selection operator (LASSO) algorithm, a CS-basedDOA estimation algorithmwhich performs iterative optimization alternatively on target angle information vector and sensing matrix mismatching error vector is proposed. The simulation result indicates that the proposed algorithm possesses higher angle resolution and estimation accuracy compared with conventional CSbased DOA estimation algorithms.


Introduction
The strong scatter centers of targets in area of interest only occupy finite angle resolution cells, and the echo signal of targets is sparse, so compressive sensing (CS) theory is widely studied in direction-of-arrival (DOA) estimation applications [1][2][3][4][5].In [1], a CS-based DOA estimation algorithm of multiple input and multiple output (MIMO) radar is proposed, which makes use of the sparsity of radar echo signals to perform compressive sampling on array receipt signals in time-domain.In [2], an array with element randomly distributed is adopted to perform compressive sampling on space-domain signal, reducing the number of receiving frontend channels of the array.However, both [1,2] treat the overcomplete based matrixes as the redundant dictionaries, obtained from the angle interval of uniform quantization area of interest, which cannot ensure that the corresponding sensing matrix meets the restricted isometry property (RIP) [3].And then, [4] proves the RIP for MIMO radar application.Reference [5] uses random Gauss matrix to perform compressive sampling on space-domain signal and adopts regularized multivectors focal undetermined system solver (RMFOCUSS) algorithm to achieve high-resolution estimation.However, the computation complexity of RMFO-CUSS algorithm increases dramatically with the increase of snapshots.
In addition, the estimation performance degrades seriously in the presence of sensing matrix mismatching in the above algorithms [6,7].The authors in [8] investigate the CS-based DOA estimation in the presence of sensing model mismatching errors, proving that the performance of CSbased DOA estimation algorithm degrades dramatically in that case.References [9][10][11] present a DOA estimation model under sensing model mismatching and then use Bayesian method to realize DOA estimation.Reference [12] proposes a joint least-absolute shrinkage and selection operator (LASSO) algorithm to achieve DOA estimation in the presence of mismatching.
In this paper, a new CS-based DOA estimation algorithm is proposed to decrease the effect of sensing matrix mismatching and achieve high resolution on DOA estimation.Firstly, a DOA sparse sensing model is set up in the presence of sensing matrix mismatching.Secondly, combining the Dantzig selector (DS) algorithm [13] and least-absolute shrinkage and selection operator (LASSO) algorithm [14], a CS-based DOA estimation algorithm which performs iterative optimization 2 International Journal of Antennas and Propagation alternatively on target angle information vector and sensing matrix mismatching error vector is proposed to achieve highresolution DOA estimates.
In practice, targets in area of interest only occupy finite angle resolution cells.So ‖s()‖ 0 =  ≪  and ‖ ⋅ ‖ 0 denotes  0 norm.Thus the output signal of the array x() is  sparse signal, A( θ) is the sparsity-based matrix, and  is the sparsity of target angle information vector.

DOA Estimation Model under Sensing Model Mismatching.
Assume that  targets' angle information vector of the array in area of interest  = [ 1 ,  2 , . . .,   ], and   ∉ θ.That is to say, the th target's angle information mismatches the angle resolution vector defined before, which is called mismatching between sensing matrix and target angle information.According to CS theory, sensing model mismatching will lead to the angle information vector failing to represent target angle precisely, decreasing the estimation accuracy of target angles through conventional CS-based DOA estimation method [8].
Thus the steering vector matrix when sensing matrix mismatches target angle information can be rewritten as where , . . .,   ]  , and ( Therefore, taking no account of the approximation error of measurement noise, (2) can be rewritten as According to CS theory, we can recover the target angle information vector s() by sampling the receipt signal with only finite array elements.So we extract  elements from  elements of the array, and let Ψ be the corresponding line-extraction matrix; thus the output of the array after line extraction could be represented as By observing (7), we can conclude that sampling of spacedomain signals can be regarded as measurement matrix Ψ performing random projection measurements on echo signal x().In addition, sensing matrix  is the product result of matrix Ψ whose elements are randomly distributed and sparsity-based matrix Φ which can be treated as Fourier transform matrix of space-domain signal.Therefore,  meets the RIP condition with great probability, thus ensuring the effectiveness and robustness of using compressive sensing reconstruction algorithm to perform DOA estimation.

Derivation of the Proposed Algorithm. By comparing (2)
and ( 6), we can find that the influences of measurement noise and sensing matrix mismatching error on DOA estimation can be summed up to "additive" disturbance and "productive" disturbance.So far, conventional CS-based DOA estimation algorithms only have constraints on "additive" disturbance but fail to take the influence of "productive" disturbance on the accuracy of target angle information estimation into consideration.Therefore, conventional CS-based DOA estimation algorithms cannot effectively reduce sensing matrix mismatching error when angle resolution vector, previously defined, fails to precisely represent target.That is to say, conventional CS-based DOA algorithms are unable to ensure DOA estimation's effectiveness and robustness when there are sensing matrix mismatching errors.
To solve these problems, a novel CS-based DOA estimation algorithm suitable for the situation when sensing matrix mismatches target angle information is proposed in this paper.The proposed algorithm combines DS algorithm and LASSO algorithm to achieve a high-resolution DOA estimation result by performing iterative optimization alternatively on target angle information vector and sensing matrix mismatching error vector.
First assume the vector of sensing matrix mismatching error  = 0 in lack of prior information.According to CS theory, the optimization problem expressed in (7) can be solved by working out  1 norm optimization under the circumstance of noise, consequently obtaining target angle information vector in space-domain: where constant  is relevant to noise variance.This optimization problem can be perfectly solved by DS algorithm.Take the estimation value of target angle information obtained by solving ( 8) to (7); we can get According to the property of vector, compiling ( 9), it can be achieved that where Ω = diag(ŝ()).
Hence, (11) can be retreated as a CS optimization problem using sensing matrix mismatching error as the sparse signal, and this CS optimization problem can be denoted as The CS optimization problem in ( 12) can be perfectly solved by LASSO algorithm.
Take the vector of sensing matrix mismatching error  worked out from ( 12) to (8); estimation value of target angle information ŝ() can be resolved.Repeat the process mentioned above until the difference of two target angle information vector's norm is less than the certain predefined threshold.That is, Stop iteration at this moment; then [ŝ()] (+1) that we get is the target angle estimation information, where Δ is the certain predefined threshold.

Simulation
In this part, numerical simulations are presented to examine the performance of the proposed method.Consider  = 10 array elements spaced randomly in [0, 30].
In the first examples, the angle resolution cell is set to be  = 1 ∘ , and the snapshots of the echo signal  = 200.The angles of input signals are assumed to be [31.3∘ , 37.6 ∘ , 45.2 ∘ ]. Figure 1 illustrates the spatial spectrum, when the signal-tonoise ratio (SNR) is set to 0 dB.It can be seen that both MUSIC algorithm and CAPON algorithm fail to achieve high-resolution estimation on target angle information.The CS-based DOA algorithm based on RMFOCUSS is endowed with better angle resolution compared to conventional DOA estimation algorithms while its estimation accuracy still suffers from the impacts of mismatching between sensing matrix and target angle information.In contrast, the proposed algorithm remarkably increases the DOA estimation accuracy by performing calibration on sensing matrix.
In the second examples, we consider the root-meansquare error (RMSE) of different DOA estimation algorithms versus different SNR.As shown in Figure 2, we can observe that the proposed algorithm possesses better estimation performance and achieves high resolution on DOA estimation, for the reason that it successfully calibrates sensing matrix mismatching error and impairs the effect from system sensing matrix mismatching in low SNR.
In the third experiment, the simulations on randomly generated DOAs are examined.The proposed method (Δ = 0.1) Method in [11] The proposed method (Δ = 0.01) method in [11], but when the certain predefined threshold Δ = 0.01, the proposed algorithm can achieve higher estimation accuracy compared with other CS-DOA methods which deal with off-grid targets.
In the last experiment, the different angle resolution cells are considered to examine the ability of the proposed method to represent the true signals.The angle resolution cells are selected as  = [1 ∘ , 3 ∘ , 5 ∘ ].The other parameters stay consistent.The RMSE versus different angle resolution cells is depicted in Figure 4, which demonstrate that the performance of the proposed method increases with the decrease of the angle resolution cell.

Conclusion
In this paper, a novel CS-based DOA estimation algorithm is proposed to solve the problem that the CS-based DOA estimation performance deteriorates in the presence of sensing matrix mismatching.The proposed algorithm reduces the estimation error of target angle information through calibrating sensing matrix.The algorithm proposed in this paper is characterized by great value in practical applications, since it improves the performance of CS-based DOA estimation algorithm and achieves high resolution on DOA estimation.

Figure 1 :
Figure 1: Angle resolution of estimated spatial spectrum.
Consider that the directions of the three signals are uniformly generated within direction intervals [20 ∘ , 40 ∘ ], and the other parameters stay consistent.The RMSE of different DOA estimation algorithms versus different SNR is plotted in Figure 3.It is seen from the figure that when the certain predefined threshold Δ = 0.1, the estimation accuracy of the proposed algorithm is less than the International Journal of Antennas and Propagation

Figure 2 :
Figure 2: RMSE of the DOA estimates versus input SNR.

Figure 3 :
Figure 3: RMSE of the DOA estimates versus input SNR.

Figure 4 :
Figure 4: RMSE of the DOA estimates versus input SNR with different angle resolution cells.
]  is steering vector of the th receipt signal,    =  sin   /, and  is the distance between the array elements. is the wavelength of carrier wave, e() = [ 1 (),  2 (), . . .,   ()]  is the array noise vector, and   () is the signal plural envelope.