Next Article in Journal
Prototyping a GNSS-Based Passive Radar for UAVs: An Instrument to Classify the Water Content Feature of Lands
Previous Article in Journal
The Mine Locomotive Wireless Network Strategy Based on Successive Interference Cancellation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Valued Covariance Vector Sparsity-Inducing DOA Estimation for Monostatic MIMO Radar

College of Automation, Harbin Engineering University, No. 145 Nantong Street, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Submission received: 21 July 2015 / Revised: 22 September 2015 / Accepted: 3 November 2015 / Published: 10 November 2015
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, a real-valued covariance vector sparsity-inducing method for direction of arrival (DOA) estimation is proposed in monostatic multiple-input multiple-output (MIMO) radar. Exploiting the special configuration of monostatic MIMO radar, low-dimensional real-valued received data can be obtained by using the reduced-dimensional transformation and unitary transformation technique. Then, based on the Khatri–Rao product, a real-valued sparse representation framework of the covariance vector is formulated to estimate DOA. Compared to the existing sparsity-inducing DOA estimation methods, the proposed method provides better angle estimation performance and lower computational complexity. Simulation results verify the effectiveness and advantage of the proposed method.

1. Introduction

Recently, multiple-input multiple-output (MIMO) radar has drawn considerable attention and become a hot research topic in the field of radar due to its advantages [1,2,3]. Based on the structure of transmit and receive arrays, MIMO radar can be divided into two classes: one is called statistical MIMO radar [2], and the other is called colocated MIMO radar [3] (including bistatic MIMO radar and monostatic MIMO radar). The transmit and receive arrays of the bistatic MIMO radar are separated from each other, so it can obtain the spatial gain from different direction channels. In colocated MIMO radar, the transmit and receive arrays are close to each other, which can improve the spatial resolution by forming a large aperture of the virtual array. In this paper, we focus on the monostatic MIMO radar, which is one kind of colocated MIMO radar.
Direction of arrival (DOA) estimation of multiple targets from the received data is an important aspect of array signal processing [4,5,6] and MIMO radar applications in practice [7,8,9,10,11]. Some subspace-based methods have been proposed for DOA estimation in MIMO radar, such as the multiple signal classification (MUSIC) algorithm [8], the estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithm [9] and tensor analysis-based algorithms [10,11]. In addition, the reduced dimensional technique-based methods [12,13,14] are proposed for DOA estimation in monostatic MIMO radar by using the special structure of the virtual array. However, the above mentioned subspace-based methods have angle estimation performance degradation under lower SNR, closely spaced targets and limited snapshots. Recently, the sparse representation technique for signal recover has attracted much attention, and it has been applied to DOA estimation in array signal processing, such as l 1 -norm singular value decomposition ( l 1 -SVD) [15], l 1 -norm sparse representation of array covariance vectors ( l 1 -SRACV) [16], covariance matrix sparse representation (CMSR) [17] and the real-valued l 1 -SVD method [18]. The simulation results have verified that compared to conventional subspace-based methods, these sparse representation-based methods have a better ability to adapt to many challenging scenarios, such as coherent targets, low SNR and limited snapshots. Furthermore, the l 1 -norm-based sparse representation frameworks are proposed for angle estimation in MIMO radar [19,20]. However, they focus on the sparse representation of the received data in MIMO radar, which involves a high computational burden due to the multiple measurement vector (MMV) problem and complex-valued processing.
In this paper, we propose a real-valued covariance vector sparsity-inducing DOA estimation method for monostatic MIMO radar. The proposed method exploits the sparsity of the real-valued covariance vector for DOA estimation, which is different from the methods in [19,20]. In the proposed method, the reduced dimensional transformation and unitary transformation techniques are used to transform the received data into a low-dimensional real-valued one. Then, based on the Khatri–Rao product, a real-valued sparse representation framework of the covariance vector is formulated for DOA estimation when multiple snapshots are available. Due to the fact that the proposed method only involves real-valued processing and the single measurement vector (SMV) problem, it has lower computational complexity than existing sparsity-inducing DOA estimation methods. Furthermore, the proposed method has superior angle estimation performance when the number of snapshots is reasonable. The simulation results are used to verify the performance of the proposed method.
The rest of the paper is organized as follows. Section 2 introduces the monostatic MIMO radar signal model and the l 1 -SVD algorithm for monostatic MIMO radar. Section 3 presents a real-valued sparse representation framework of the covariance vector for DOA estimation in monostatic MIMO radar. Section 4 presents some related remarks for this paper. Several numerical simulations are carried out in Section 5 to evaluate the performance of the proposed method. Finally, Section 6 concludes this paper.
Notation: ( · ) H , ( · ) T , ( · ) - 1 , ( · ) * and ( · ) + denote conjugate-transpose, transpose, inverse, conjugate and pseudo-inverse, respectively. ⊗ and ⊙ denote the Kronecker product operator and the Hadamard product operator, respectively. I K denotes a K × K -dimensional unit matrix, and E [ · ] denotes the expectation operator. A ( l 2 ) denotes a column vector whose q-th element equals the l 2 norm of the q-th row of A . | | · | | 1 and | | · | | F denote the l 1 norm and the Frobenius norm, respectively.

2. The Monostatic MIMO Radar Signal Model and the l 1 -SVD Algorithm

We consider a narrowband monostatic MIMO radar system equipped with M transmit antennas and N receive antennas. Both the transmit and receive arrays are half-wavelength-spaced uniform linear arrays (ULAs), and all of the elements are omnidirectional. At the transmit array, M transmit antennas emit M orthogonal narrow-band waveforms, which have identical bandwidth and center frequency. Let g ( t ) = [ ϕ 1 ( t ) , , ϕ M ( t ) ] T C M × 1 be the orthogonal narrow-band transmitted waveforms, which satisfy:
0 T q ϕ i ( t ) ϕ j * ( t ) = 1 , i = j 0 , i j
where T q is the pulse duration. It is assumed that there exists P uncorrelated targets located in the far-field region, and the DOA of the p-th target with respect to the transmit array and receive array is denoted as θ p . Then, the N × 1 received data can be written as [7]:
x ¯ ( t , τ ) = p = 1 P β p ( τ ) e j 2 π f p ( τ ) a r ( θ p ) a t T ( θ p ) g ( t ) + n ¯ ( t , τ )
where τ is the slow time index, i.e., the pulse number (denoted as a snapshot in this paper), and β p ( τ ) and f p ( τ ) denote the reflection coefficient and Doppler frequency, respectively. n ¯ ( t , τ ) C N × M is the noise matrix. The transmit steering vector a t ( θ p ) and the receive steering vector a r ( θ p ) are shown as:
a t ( θ p ) = [ e - j M - 1 2 sin θ p , , e - j 1 2 sin θ p , e j 1 2 sin θ p , , e - j M - 1 2 sin θ p ] T , M is even [ e - j M - 1 2 sin θ p , , 1 , , e - j M - 1 2 sin θ p ] T , M is odd
and:
a r ( θ p ) = [ e - j N - 1 2 sin θ p , , e - j 1 2 sin θ p , e j 1 2 sin θ p , , e - j N - 1 2 sin θ p ] T , N is even [ e - j N - 1 2 sin θ p , , 1 , , e - j N - 1 2 sin θ p ] T , N is odd
Exploiting the orthogonality property of the transmitted waveforms, the M × 1 virtual data vectors can be obtained by matched filtering x ¯ ( t , τ ) with each of the orthogonal waveforms ϕ m ( t ) , m = 1 , 2 , . . . . , M , which is given as follows:
x ¯ m ( τ ) = 0 T q x ¯ ( t , τ ) ϕ m * ( t ) d t = p = 1 P β p ( τ ) e j 2 π f p ( τ ) a t m ( θ p ) a r ( θ p ) + n ¯ m ( τ )
where a t m ( θ p ) is the m-th element of the transmit steering vector a t ( θ p ) , and n ¯ m ( τ ) = 0 T q n ¯ ( t , τ ) ϕ m * ( t ) d t is the noise term. Then, stacking the individual vector components into one column vector, an M N × 1 virtual vector can be obtained, which is written as [7,8,9]:
x ( τ ) = [ x ¯ 1 T ( τ ) , x ¯ 2 T ( τ ) , , x ¯ M ( τ ) T ] T = A s ( τ ) + n ( τ )
where A = [ a t ( θ 1 ) a r ( θ 1 ) , , a t ( θ P ) a r ( θ P ) ] C N M × P , s ( τ ) = [ s 1 ( τ ) , , s P ( τ ) ] T C P × 1 is the signal vector with s p ( τ ) = β p ( τ ) e j 2 π f p ( τ ) . n ¯ ( τ ) = [ n 1 T ( τ ) , , n ¯ M T ( τ ) ] T is the noise vector. When the number of snapshots is J, the Equation (4) can be expressed as:
X = A S + N
where X = [ x ( τ 1 ) , x ( τ 2 ) , x ( τ J ) ] C M N × J is the received data matrix. S = [ s ( τ 1 ) , s ( τ 2 ) , s ( τ J ) ] and N = [ n ( τ 1 ) , n ( τ 2 ) , , n ( τ J ) ] are the signal data matrix and noise matrix, respectively. Some statistical assumptions on the target signals and noise are made as follows:
(i)
The target signals are spatially uncorrelated, temporally white and zero-mean, i.e., the Doppler frequencies f p ( τ ) ( p = 1 , 2 , . . . , P ) are different, and the reflection coefficients β p ( τ k ) ( k = 1 , 2 , . . . , J ) obey the Swerling II model.
(ii)
The noise matrix N is zero-mean, complex circular Gaussian and with the variance R n = δ 2 I M N .
(iii)
The noise is statistically independent of all targets.
Based on the above statistical assumptions, the work is to estimate the DOA from the received data X in this paper. Firstly, we introduce the l 1 -SVD algorithm for DOA estimation. Due to the fact that the DOAs of targets are sparse when the discretized sampling grid number of all spatial spaces is much larger than the target number, the l 1 -SVD algorithm [14] can be applied to DOA estimation in monostatic MIMO radar. After using the SVD technique, the received data can be written as:
X S V = A S S V + N S V
where X S V = X V s , S S V = S V s and N S V = N V s . V s C J × P is composed with the singular vectors corresponding to the P largest singular values of X . Then, Equation (6) can be sparsely represented by using the sampling grid Θ = [ θ ^ 1 , θ ^ 2 , , θ ^ L ] , where L P . The transmit and receive complete dictionaries are expressed as A t θ ^ = [ a t ( θ ^ 1 ) , , a t ( θ ^ L ) ] and A r θ ^ = [ a r ( θ ^ 1 ) , , a t ( θ ^ L ) ] , respectively. Thus, the complete dictionary can be constructed as A θ ^ = A t θ ^ A r θ ^ . Then, under the framework of the sparse representation, Equation (6) can be formulated as:
X S V = A θ ^ S S V θ ^ + N S V
where S S V θ ^ C L × P and S S V have the same row support, i.e., the matrix S S V θ ^ is sparse. In order to estimate S S V θ ^ , Equation (7) can be formulated as the l 1 -norm minimization problem, which is shown as follows:
min ( S S V θ ^ ) ( l 2 ) 1 , s . t . X S V - A θ ^ S S V θ ^ 2 η ¯
where η ¯ is the regularization parameter. After obtaining the sparse matrix S S V θ ^ , the DOAs can be estimated by finding non-zero rows of S S V θ ^ . According to Equation (8), the procedure of recovering S S V θ ^ involves the MMV problem and complex-valued processing. Thus, the l 1 -norm-based sparse representation methods lead to high computational complexity. On the other hand, the spatial resolution depends on the aperture of the MIMO radar in Equation (8), which has limited ability to estimate the DOA when the targets are located closely. In the next section, we propose a real-valued sparse representation framework of the covariance vector for DOA estimation to solve the drawbacks mentioned above.

3. Real-Valued Sparse Representation Framework of the Covariance Vector for DOA Estimation

When applying the l 1 -norm-based sparse representation methods to estimate the DOA in MIMO radar, the MMV problem and complex-value processing are involved [19,20]. Thus, the computational complexity can be significantly reduced if the MMV problem can be converted into the SMV problem, and only real-valued processing is required in the l 1 -norm minimization problem. In addition, the aperture of MIMO radar can be enlarged if the covariance matrix can be vectorized. These motivate us to propose a real-valued sparse representation framework of the covariance vector for DOA estimation in monostatic MIMO radar.

3.1. Reduced Dimension and Real-Value Transformation

According to the special configuration of monostatic MIMO radar, there are only Q = M + N - 1 distinct elements. Therefore, the steering vector a t ( θ ) a r ( θ ) can be expressed as:
a t ( θ ) a r ( θ ) = G M N × Q b ( θ )
where G M N × Q and b ( θ ) are the reduced dimensional transformation matrix and the steering vector, respectively, and they can be written as:
G M N × Q = [ J 0 T , J 1 T , , J M - 1 T ] T
b ( θ ) = [ e - j Q - 1 2 sin θ , , e - j 1 2 sin θ , e j 1 2 sin θ , , e - j Q - 1 2 sin θ p ] T , Q is even [ e - j Q - 1 2 sin θ p , , 1 , , e - j Q - 1 2 sin θ p ] T , Q is odd
where J m = [ 0 N × m , I N , 0 N × ( M - m - 1 ) ] , m = 0 , 1 , . . . , M - 1 . Based on Equation (10), we define a matrix as F = G H G , which can be written as:
F = diag [ 1 , 2 , . . . , min ( M , N ) , . . . , min ( M , N ) | M - N | + 1 , . . . . , 2 , 1 ]
In order to hold the characteristic of the noise matrix N , the reduce-dimensional matrix is defined as W = F - ( 1 / 2 ) G H [14], which satisfies W W H = I M + N - 1 . Then, multiplying the received signal X by W , we have:
Y = F ( 1 / 2 ) B S + W N = B ^ S + W N
where B = [ b ( θ 1 ) , b ( θ 2 ) , . . . , b ( θ P ) ] C Q × P and B ^ = F ( 1 / 2 ) B . Then, the received data in Equation (13) correspond to the effective steering vector b ( θ ) , which can be used for DOA estimation. The sparse representation framework for DOA estimation in [19] is based on received data, and the l 1 -norm minimization problem involves the MMV problem and complex-valued processing. In order to avoid this issue, a real-valued sparse representation framework of the covariance vector with the SMV problem is proposed as follows. Based on Equation (13), we note that the steering matrix becomes B ^ , which satisfies:
Π M + N - 1 B ^ * = B ^
where Π K denotes the K × K exchange matrix with ones on its anti-diagonal and zeros elsewhere. It can be concluded from Equation (14) that the linear array is a centro-symmetric array after the reduced-dimensional transformation. Therefore, the unitary transformation can be adopted to convert the complex-valued steering matrix into a real-valued one. Following the convention in [21], the unitary transformation matrix is defined as follows:
U 2 K = 1 2 I K j I K Π K - j Π K
and:
U 2 K + 1 = 1 2 I K 0 j I K 0 T 2 0 T Π K 0 - j Π K
Then, multiplying the received data Y by the unitary transformation matrix U Q , a new received data matrix can be expressed as:
Y ^ = U Q H Y = U Q H B ^ S + U Q H W N
According to Equation (17), after the unitary transformation, the complex-valued steering matrix is converted into real-valued steering matrix B ¯ = U Q H B ^ . Then, the data matrix Y ¯ can be divided into real and imaginary parts, respectively, which are shown as:
Re ( Y ^ ) = B ¯ Re ( S ) + Re ( U Q H W N )
Im ( Y ^ ) = B ¯ Im ( S ) + Im ( U Q H W N )
where Re ( · ) and Im ( · ) denote the real part and the imaginary part, respectively. Using Equations (18) and (19), a real-valued augmentation data matrix can be formulated, which is expressed as:
Y ¯ = [ Re ( Y ^ ) , Im ( Y ^ ) ] = B ¯ S ¯ + N ¯
where S ¯ = [ Re ( S ) , Im ( S ) ] R P × 2 J and N ¯ = [ Re ( U Q H W N ) , Im ( U Q H W N ) ] R Q × 2 J . With the statistical Assumptions (i), (ii) and (iii), the covariance matrix of S ¯ is R S = E [ S ¯ S ¯ H ] = diag ( d ) with diagonal elements d = [ d 1 2 , . . . . , d P 2 ] T being the signal power vector.

3.2. Real-Valued Sparse Representation Framework for DOA Estimation

According to Equation (20), the real-valued covariance matrix of Y ¯ can be expressed as:
R Y ¯ = E [ Y ¯ Y ¯ H ] = B ¯ R S B ¯ H + σ 2 I Q = B ¯ diag ( d ) B ¯ H + σ 2 I Q
Based on the Khatri–Rao product, the vectorization of the covariance matrix provides a new way for array augmentation [22]. Thus, exploiting the vectorization (vec) operator on R Y ¯ , we have:
y = vec ( R Y ¯ ) = ( B ¯ B ¯ ) d + σ 2 vec ( I Q )
According to Equation (22), the data vector y can be considered as a new signal model corresponding to a virtual array output with a single snapshot, and the effective array response vector of the virtual array can be regarded as b ( θ ) b ( θ ) . According to the received data in Equation (5), the effective array response vector of the received data is b ( θ ) . Thus, the proposed signal model in Equation (22) significantly increases the circs of freedom (DOF) of the virtual array and shows an enlarged aperture. Moreover, the sparse representation framework for Equation (22) is only involved with the SMV problem, which can remarkably reduce the computational complexity compared to the MMV problem in Equation (5).
Based on the proposed signal model in Equation (22), a real-valued sparse representation scheme of the covariance vector for DOA estimation is formulated. Let θ ^ 1 , θ ^ 2 , , θ ^ L be the discretized sampling grid of the spatial domain of interest. The number of potential DOAs will be much greater than the number of targets, i.e., L P . Then, the one-dimensional complete dictionary for DOA estimation can be constructed as B ̲ θ ^ = B ¯ θ ^ B ¯ θ ^ , where B ¯ θ ^ = U Q H F ( 1 / 2 ) B θ ^ with B θ ^ = [ b ( θ ^ 1 ) , . . . , b ( θ ^ L ) ] . Under the sparse representation framework, Equation (22) can be represented with the complete dictionary as:
y = B ̲ θ ^ d θ ^ + σ 2 vec ( I Q )
where d θ ^ C L × 1 is a P-sparse vector, which has P non-zero elements. Thus, the DOAs can be estimated by detecting the non-zeros elements of d θ ^ .
In practice, the unknown covariance matrix can be estimated by R ^ Y ¯ = ( 1 / J ) Y ¯ Y ¯ H . Then, we have y ^ = vec ( R ^ Y ¯ ) and y ^ - y = y , where y is the estimation error. Then, combining the sparse representation framework in Equation (23) and the estimation error y , a real-valued l 1 norm minimization problem with SMV is formulated as:
min d θ ^ 1 s . t . y ^ - B ̲ θ ^ d θ ^ - σ 2 vec ( I Q ) 2 2 η
where η is a regularization parameter, which sets the amount of error and plays an important role in the final DOA estimation performance. According to the l 1 norm minimization problem in Equation (24), the selection of parameter η depends on the distribution of y .
Lemma 1: After using the reduced dimensional transformation and unitary transformation, the estimation error y satisfies y A s N ( 0 , 1 J R Y ¯ T R Y ¯ ) , where A s N ( μ , σ 2 ) represents the asymptotic normal distribution with mean μ and variance σ 2 .
Proof: According to the definition of W and U Q , both of them are orthogonal matrices. Thus, [ U Q H B ^ S ] i , j and [ U Q H W N ] i , j have a complex Gaussian distribution with zero mean, because the orthogonal invariance property of the Gaussian random matrix makes the distribution impervious to multiplication by orthogonal matrices, where [ · ] i , j denotes the ( i , j ) -th element of a matrix. Furthermore, due to the fact that the real part and imaginary part of [ U Q H B ^ S ] i , j and [ U Q H W N ] i , j are real a Gaussian distribution with zero mean, [ B ¯ S ¯ ] i , j and [ N ¯ ] i , j satisfy the Gaussian distribution with zero mean. Thus, according to [23], it can be concluded that the estimation error y = y ^ - y satisfies:
y = vec ( R ^ Y ¯ - R Y ¯ ) A s N ( 0 , 1 J R Y ¯ T R Y ¯ )
According to Equation (25), it is noticed that y is not the asymptotic standard normal distribution, and the parameter η cannot be calculated in an easy way. However, it is obvious that a weighted matrix W ¯ - 1 2 = 1 J R Y ¯ - T 2 R Y ¯ - 1 2 can be formulated to make the estimation error y satisfy the asymptotic standard normal distribution; then, we further have:
W ¯ - 1 2 y 2 2 A s χ 2 ( 0 , I Q 2 )
where A s χ 2 ( 0 , I Q 2 ) denotes the asymptotic chi-square distribution with Q 2 circs of freedom. Combining Equation (24) and Equation (26), the weighted real-valued l 1 norm minimization problem with SMV is formulated as:
min d θ ^ 1 s . t . W ^ - 1 2 ( y ^ - σ ^ 2 vec ( I Q ) ) - ( W ^ - 1 2 B ̲ θ ^ ) d θ ^ 2 η 1
where W ^ - 1 2 is the estimation of W ¯ - 1 2 and can be calculated as W ^ - 1 2 = 1 J R ^ Y ¯ - T 2 R ^ Y ¯ - 1 2 . According to Equation (26), the parameter η 1 can be selected with a high probability 1 - ξ , where ξ is a small value. Usually, it is enough to set ξ = 0 . 001 for calculating the value of η 1 , which can be implemented via the function η 1 = c h i 2 i n v ( 1 - ξ , Q 2 ) in MATLAB software. σ ^ r 2 is the estimation of σ ^ r 2 and can be estimated by the average of the Q - P smallest eigenvalue or the minimum eigenvalue of R ^ Y ¯ . Equation (27) can be calculated by SOC (second order cone) programming software packages, such as Sedumi [24] and CVX [25]. Then, a spatial spectrum can be obtained by plotting d θ ^ solved from Equation (27), and the DOA can be estimated by finding the P largest values.

4. Related Remarks

Remark 1: In order to obtain the computational complexity of the sparsity-inducing DOA estimation methods, we firstly analyze the computational complexity of solving the SOC programming, such as in Equations (8) and (27). Referring to the conclusion in [15], for the observing matrix X S V C M N × P , the computational complexity of solving SOC programming in Equation (8) is O ( L 3 P 3 ) flops when L M N , where flops is the abbreviation of floating-point operations. For the the observing vector y ^ C Q 2 × 1 in the proposed method, the computational complexity of solving the SOC programming in Equation (27) is O ( L 3 ) due to L Q 2 . Thus, the computational complexity analysis of the proposed method is given as follows. As shown in the steps of the proposed method, the main computational complexity of the proposed method contains three parts: (1) the reduced dimensional transformation in Equation (13) requires O ( Q M N J ) flops; (2) the unitary transformation in Equation (17) and calculating real-valued covariance matrix in Equation (21) require O ( 2 Q 2 J ) flops; (3) solving the SOC programming in Equation (27) with real-valued processing requires O ( 1 4 L 3 ) flops. Thus, the computational complexity of the proposed method is O ( Q M N J + 2 Q 2 J + 1 4 L 3 ) flops. For the complex-valued l 1 -SVD algorithm [15] with the received data in MIMO radar, the computational complexity focus on the SVD of the received data in Equation (5) and solving the SOC programming in Equation (8), which requires O ( M 3 N 3 + L 3 P 3 ) flops. For the real-valued (RV) l 1 -SVD algorithm, it only requires the additions to transform the complex received data into a real one compared to the l 1 -SVD algorithm. Thus, it requires O ( M N J + 1 4 M 3 N 3 + 1 4 L 3 P 3 ) . Owing to L M , N , J , Q , the proposed method has lower computational complexity than the RV l 1 -SVD and l 1 -SVD algorithms.
Remark 2: Based on the principle of the proposed method mentioned above, it can be concluded that the proposed method has better angle estimation performance and higher angular resolution than conventional l 1 -norm based algorithms when the covariance matrix can be estimated efficiently, i.e, the number of snapshots is reasonable. This is because exploiting the covariance matrix vectorization technique in Equation (22), the virtual aperture corresponding to y has been enlarged remarkably, i.e., the angular resolution can be improved.
Remark 3: The maximal number of identified targets is an important aspect that can be considered for the sparse representation-based algorithms. Then, according to [21], any set of 2 Q - 1 columns of the complete dictionary B ̲ θ ^ is independent, which leads to Spark [ B ̲ θ ^ ] = 2 Q , where Spark [ Ω θ ¯ ] denotes the smallest integer of columns of the complete dictionary B ̲ θ ^ that are linearly dependent. Then, the maximal number of identified targets is Spark [ B ̲ θ ^ ] / 2 - 1 = Q - 1 in the proposed method [26]. Thus, though the virtual aperture is enlarged, the proposed method cannot handle the underdetermined DOA estimation case.
Remark 4: According to Equation (22), if the target signals are coherent, the signal power vector d exists with zero elements, and the vector y is meaningless. Thus, the proposed method is not suitable for the coherent targets. Additionally, as shown in Equations (21) and (22), the proposed method is based on the assumption that R S is a diagonal matrix. The matrix R S , in fact, does not satisfy the diagonal matrix when the number of snapshots is limited. Therefore, the angle estimation performance of the proposed method degrades remarkably with lower snapshots.

5. Simulation Results

In this section, we carry out some simulations to evaluate the proposed method. We also compare the proposed methods with the l 1 -SVD algorithm [15], the real-valued l 1 -SVD algorithm (denoted as RV l 1 -SVD) [18] and the Cram e ´ r–Rao bound (CRB) [27]. In most of the simulations, a monostatic MIMO radar with M = N = 5 colocated antennas for the transmit and receive arrays is used. Both of them are half-wavelength-spaced ULAs. The signal-to-noise ratio (SNR) is defined as 10log ( σ s 2 / σ n 2 ) , where σ s 2 and σ n 2 denote the signal and noise power, respectively. The root-mean-square-error (RMSE) of the DOA estimation obtained from 300 Monte Carlo runs is defined as:
RMSE = 1 300 P i = 1 300 p = 1 P ( θ ^ i , p - θ p ) 2
where θ ^ i , p denotes the estimation of θ p at the i-th trial. For all methods in the simulations, the number of targets is assumed to be known or estimated by the minimum description length (MDL) criterion, and the confidence interval is set to 0.999 for the l 1 -SVD algorithm, the RV l 1 -SVD algorithm and the proposed method (i.e., ξ = 0 . 001 ). The spatial grid is uniform from the range - 90 to 90 with 0 . 1 interval in all simulations.
Figure 1 shows the spectra of these three algorithm, where P = 3 uncorrelated targets are assumed to be with the DOAs as θ 1 = - 10 . 13 , θ 2 = 0 . 13 and θ 3 = 20 . 13 . The SNR and snapshots are fixed with 0 dB and 400, respectively. As can be seen from Figure 1, compared to the RV l 1 -SVD algorithm and the l 1 -SVD algorithm, the proposed method has a lower sidelobe, which means that the proposed method has higher angular resolution.
Figure 1. The spectra of l 1 -SVD, real-valued (RV) l 1 -SVD and the proposed method for uncorrelated targets.
Figure 1. The spectra of l 1 -SVD, real-valued (RV) l 1 -SVD and the proposed method for uncorrelated targets.
Sensors 15 28271 g001
Figure 2 and Figure 3 show the angle estimation performance of these three algorithms versus the SNR and snapshots, respectively, where P = 3 uncorrelated targets are assumed to be with the DOAs as θ 1 = - 10 . 13 , θ 2 = 0 . 13 and θ 3 = 20 . 13 . In Figure 2, we keep the snapshots fixed at 400, and the SNR is varied from −10 dB to 20 dB; in Figure 3, the SNR is fixed at 0 dB, and the number of snapshots is varied from 50 to 600. It can be seen from Figure 2 that the RV l 1 -SVD algorithm provides better angle estimation performance than the l 1 -SVD algorithm at the low SNR region and comparable angle estimation performance at the high SNR region. Compared to the RV l 1 -SVD and l 1 -SVD algorithms, the proposed method outperforms both of them at all SNR regions. The reason is that the virtual aperture is enlarged in the proposed method. From Figure 3, it can be concluded that the proposed method provides better angle estimation performance than both the RV l 1 -SVD and l 1 -SVD algorithms when the number of snapshots satisfies J > 100 . However, the angle estimation performance is deteriorated in the adaptation to much lower snapshots. This is because the covariance vectors become heavily perturbed when the snapshot number is small. Thus, the proposed method provides better angle estimation performance with a sufficiently large sample size.
Figure 2. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 2. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Sensors 15 28271 g002
Figure 3. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the snapshots number varies from 20 to 600.
Figure 3. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the snapshots number varies from 20 to 600.
Sensors 15 28271 g003
Figure 4 shows the target resolution probability of these three algorithms versus angle separation, where the SNR and snapshot number are fixed at 0 dB and 400, respectively. Two uncorrelated targets are considered with the DOAs as θ 1 = 0 . 13 and θ 2 = 0 . 13 + θ , where θ varies from 2 to 20 . In this simulation, two targets can be considered to be resolved if there are at least two peaks in the spatial spectrum, and it satisfies m a x i = 1 , 2 | θ ^ i - θ i | ¯ θ / 2 , where ¯ θ / 2 = | θ 2 - θ 1 | and θ ^ i is the estimation of θ i . As can be seen from Figure 4, compared to the RV l 1 -SVD and l 1 -SVD algorithms, the proposed method shows the highest target resolution probability for closely-spaced targets. The reason is that the virtual aperture is remarkably enlarged in the proposed method, i.e., the angular resolution is improved.
Figure 4. Target resolution probability of l 1 -SVD, RV l 1 -SVD and the proposed method for two uncorrelated targets when the angle separation varies from 2 to 20 .
Figure 4. Target resolution probability of l 1 -SVD, RV l 1 -SVD and the proposed method for two uncorrelated targets when the angle separation varies from 2 to 20 .
Sensors 15 28271 g004
Figure 5 shows the target resolution probability of these three algorithms versus SNR, where the snapshots number is fixed at 400. Two uncorrelated closely-spaced targets are considered with the DOAs as θ 1 = 0 . 52 and θ 2 = 4 . 52 . The definition of target resolution is the same as Figure 4. It can be seen from Figure 5 that all methods exhibit a 100% correct target resolution at the high SNR region. As the SNR decreases, the probability of target resolution starts dropping for each method at a certain point, which is known as the SNR threshold. From Figure 5, it can be seen that the proposed method has a lower SNR threshold compared to the RV l 1 -SVD and l 1 -SVD algorithms, i.e., the proposed method shows the best capability in resolving closely-spaced targets.
Figure 6 shows the RMSE of these three algorithms versus the number of targets, where the SNR = 0 dB and the number of snapshots is fixed at 400. The DOA of the pth target is - 30 + ( p - 1 ) 10 . From Figure 6, it can be seen that the proposed method provides better angle estimation performance than the RV l 1 -SVD and l 1 -SVD algorithms with different numbers of targets. This is because the proposed method can enlarge the aperture of the virtual array. Figure 7 shows the RMSE of the proposed method versus different transmit and receive elements, where P = 3 uncorrelated targets are assumed to be with the DOAs as θ 1 = - 10 . 13 , θ 2 = 0 . 13 and θ 3 = 20 . 13 , and the number of the snapshots is fixed as 400. As shown in Figure 7, the proposed method has improved angle estimation performance when the number of elements in the transmit or receive array increases. This is because of the diversity gain.
Figure 5. Target resolution probability of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 5. Target resolution probability of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Sensors 15 28271 g005
Figure 6. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method versus the number of targets.
Figure 6. RMSE of l 1 -SVD, RV l 1 -SVD and the proposed method versus the number of targets.
Sensors 15 28271 g006
The angle estimation performance of the proposed method has been evaluated with the foregoing simulations. Here, the computational complexity of l 1 -SVD, RV l 1 -SVD and the proposed method is evaluated by using the time in count (TIC) and time out count (TOC) instruction in MATLAB software. The number of targets P is assumed to be known beforehand. The SNRs of all targets are set at 10 dB, and the number of snapshots is 400. The DOAs of targets satisfy θ i + 1 - θ i = 10 ( i = 1 , 2 , . . . , P - 1 ) . For each target number, the computation time of these algorithms is averaged over 300 trials. The computation time comparison of these algorithms is shown in Table 1. As can be seen from Table 1, the computation time of the l 1 -SVD and RV l 1 -SVD algorithms increases markedly when the number of targets increases, but the computation time of the proposed method is similar. On the other hand, the computation time of the proposed method is lower than both the l 1 -SVD and RV l 1 -SVD algorithms.
Figure 7. RMSE of the proposed method versus different elements when the SNR varies from −10 dB to 20 dB.
Figure 7. RMSE of the proposed method versus different elements when the SNR varies from −10 dB to 20 dB.
Sensors 15 28271 g007
Table 1. Computation time comparison of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the targets number varies from 3 to 6.
Table 1. Computation time comparison of l 1 -SVD, RV l 1 -SVD and the proposed method for uncorrelated targets when the targets number varies from 3 to 6.
Target NumberAverage Computation Time (s)
l 1 -SVDRV l 1 -SVDProposed Method
P = 33.73923.28852.4312
P = 47.88784.17072.5367
P = 59.62595.45072.7102
P = 611.11615.94642.7992

6. Conclusions

In this paper, we proposed a real-valued sparse representation framework of covariance vector for DOA estimation in MIMO radar. The proposed method exploits the reduced dimension and unitary transformation technique to turn the received data into a low dimensional real-valued one, then a sparse representation framework of the covariance vector is formulated for DOA estimation. Both the theoretical analysis and simulation results verify that the proposed method provides better angle estimation performance and has lower computational complexity than both the l 1 -SVD and RV l 1 -SVD algorithms.

Acknowledgments

This work is supported by the National Natural Science Foundation (61571148), China postdoctoral special funding (2015T80328), the China Postdoctoral Science Foundation Grant (2014M550182), the Heilongjiang Postdoctoral Special Fund (LBH-TZ0410) and the Innovation of Science and Technology Talents in Harbin (2013RFXXJ016).

Author Contributions

All authors contributed extensively to this paper. X.W. and W.W. provided the idea of this work. X.W. wrote and modified the paper. W.W. provided many valuable suggestions for preparing the paper. X.L. and J.L. performed the experiments and analyzed the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fishler, E.; Haimovich, A.; Blum, R.; Chizhik, D.; Cimini, L.; Valenzuela, R. MIMO Radar: An idea whose time has come. In Proceedings of the IEEE Radar Conference, Philadelphia, PA, USA, 26–29 April 2004; pp. 71–78.
  2. Haimovich, A.M.; Blum, R.; Cimini, L. MIMO radar with widely separated antennas. IEEE Signal Process. Mag. 2008, 25, 116–129. [Google Scholar] [CrossRef]
  3. Li, J.; Stoica, P. MIMO radar with colocated antennas. IEEE Signal Process. Mag. 2007, 24, 106–114. [Google Scholar] [CrossRef]
  4. Wang, W.; Ren, S.; Dong, Y.; Wang, H. An efficient algorithm for direction finding against unknown mutual coupling. Sensors 2014, 14, 20064–20077. [Google Scholar] [CrossRef]
  5. Liang, J.; Liu, D. Joint elevation and azimuth direction finding using L-shaped array. IEEE Trans. Antennas Propag. 2010, 58, 2136–2141. [Google Scholar] [CrossRef]
  6. Liang, Y.; Ying, R.; Lu, Z.; Liu, P. Off-grid direction of arrival estimation based on joint spatial sparsity for distributed sparse linear arrays. Sensors 2014, 14, 21981–22000. [Google Scholar] [PubMed]
  7. Yan, H.; Li, J.; Liao, G. Multitarget identification and localization using bistatic MIMO radar systems. EURASIP J. Adv. Signal Process. 2008. [Google Scholar] [CrossRef]
  8. Gao, X.; Zhang, X.; Feng, G.; Wang, Z.; Xu, D. On the MUSIC-derived aprroaches of angle estimation for bistatic MIMO radar. In Proceedings of the 2009 International Conference Wireless Networks and Information Systems, Shanghai, China, 28–29 December 2009; pp. 343–346.
  9. Chen, D.; Chen, B.; Qin, G. Angle estimation using ESPRIT in MIMO radar. Electron. Lett. 2008, 44, 770–771. [Google Scholar] [CrossRef]
  10. Nion, D.; Sidiropoulos, N.D. Tensor algebra and multidimensional harmonic retrieval in signal processing for MIMO radar. IEEE Trans. Signal Process. 2010, 58, 5693–5705. [Google Scholar] [CrossRef]
  11. Wang, X.; Wang, W.; Li, X.; Wang, J. A tensor-Based subspace approach for bistatic MIMO radar in spatial colored noise. Sensors 2014, 14, 3897–3907. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, X.; Xu, D. Low-complexity ESPRIT-based DOA estimation for colocated MIMO radar using reduced-dimension transformation. Electron. Lett. 2011, 47, 283–284. [Google Scholar] [CrossRef]
  13. Zhang, X.; Huang, Y.; Chen, C.; Li, J.; Xu, D. Reduced-complexity Capon for derection of arrival estimation in a monostatic mutiple-input multiple-output radar. IET Radar Sonar Navig. 2012, 8, 796–801. [Google Scholar] [CrossRef]
  14. Wang, W.; Wang, X.; Song, H.; Ma, Y. Conjugate ESPRIT for DOA estimation in monostatic MIMO radar. Signal Process. 2013, 93, 2070–2075. [Google Scholar] [CrossRef]
  15. Malioutov, D.; Cetin, M.; Willsky, A.S. A sparse signal reconstruction perspective for source localization with sensor arrays. IEEE Trans. Signal Process. 2005, 53, 3010–3022. [Google Scholar] [CrossRef]
  16. Yin, J.; Chen, T. Direction-of-arrival estimation using a sparse representation of array covariance vectors. IEEE Trans. Signal Process. 2011, 59, 4489–4493. [Google Scholar] [CrossRef]
  17. Liu, Z.; Huang, Z.; Zhou, Y. Array signal processing via sparsity-inducing representation of the array covariance matrix. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1710–1724. [Google Scholar] [CrossRef]
  18. Dai, J.; Xu, X.; Zhao, D. Direction-of-arrival estimation via real-valued sparse representation. IEEE Antennas Wirel. Propag. Lett. 2013, 12, 376–379. [Google Scholar] [CrossRef]
  19. Wang, X.; Wang, W.; Liu, L.; Li, X.; Wnag, L. A sparse representation scheme for angle estimation in monostatic MIMO radar. Signal Process. 2014, 104, 258–263. [Google Scholar] [CrossRef]
  20. Li, J.; Zhang, X. Sparse representation-based joint angle and Doppler frequency estimation for MIMO radar. Multidimens. Syst. Signal Process. 2013, 26, 179–192. [Google Scholar] [CrossRef]
  21. Yilmazer, N.; Koh, J.; Sarkar, T.K. Utilization of a unitary transform for efficient computation in the matrix pencil method to find the direction of arrival. IEEE Trans. Antennas Propag. 2006, 54, 175–181. [Google Scholar] [CrossRef]
  22. Ma, W.K.; Hsieh, T.H.; Chi, C.Y. DOA estimation of quasi-stationary signals with less sensors than sources and unknown spatial noise covariance: A khatri-rao subspace approach. IEEE Trans. Signal Process. 2010, 58, 2168–2180. [Google Scholar] [CrossRef]
  23. Ottersten, B.; Stoica, P.; Roy, R. Covariance matching estimation techniques for array signal processing applications. Digit. Signal Process. 1998, 8, 185–210. [Google Scholar] [CrossRef]
  24. Sturm, J.S. Using SeDuMi 1.02, a MATLAB Toolbox for Optimization Over Symmetric Cones. Available online: http://fewcal.kub.nl/strum (accessed on 3 November 2015).
  25. Grant, M.; Boyd, S. CVX: MATLAB Software for Disciplined Convex Programming. Available online: http://cvxr.com/cvx (accessed on 3 November 2015).
  26. Davies, M.E.; Eldar, Y.C. Rank awareness in joint sparse recovery. IEEE Trans. Inf. Theory 2012, 58, 1135–1146. [Google Scholar] [CrossRef]
  27. Hassanien, A.; Vorobyov, S.A. Transmit energy focusing for DOA estimation in MIMO radar with colocated antenna. IEEE Trans. Signal Process. 2011, 59, 2669–2682. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Wang, X.; Wang, W.; Li, X.; Liu, J. Real-Valued Covariance Vector Sparsity-Inducing DOA Estimation for Monostatic MIMO Radar. Sensors 2015, 15, 28271-28286. https://doi.org/10.3390/s151128271

AMA Style

Wang X, Wang W, Li X, Liu J. Real-Valued Covariance Vector Sparsity-Inducing DOA Estimation for Monostatic MIMO Radar. Sensors. 2015; 15(11):28271-28286. https://doi.org/10.3390/s151128271

Chicago/Turabian Style

Wang, Xianpeng, Wei Wang, Xin Li, and Jing Liu. 2015. "Real-Valued Covariance Vector Sparsity-Inducing DOA Estimation for Monostatic MIMO Radar" Sensors 15, no. 11: 28271-28286. https://doi.org/10.3390/s151128271

APA Style

Wang, X., Wang, W., Li, X., & Liu, J. (2015). Real-Valued Covariance Vector Sparsity-Inducing DOA Estimation for Monostatic MIMO Radar. Sensors, 15(11), 28271-28286. https://doi.org/10.3390/s151128271

Article Metrics

Back to TopTop