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
-norm singular value decomposition (
-SVD) [
15],
-norm sparse representation of array covariance vectors (
-SRACV) [
16], covariance matrix sparse representation (CMSR) [
17] and the real-valued
-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
-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
-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: , , , and denote conjugate-transpose, transpose, inverse, conjugate and pseudo-inverse, respectively. ⊗ and ⊙ denote the Kronecker product operator and the Hadamard product operator, respectively. denotes a -dimensional unit matrix, and denotes the expectation operator. denotes a column vector whose q-th element equals the norm of the q-th row of . and denote the norm and the Frobenius norm, respectively.
2. The Monostatic MIMO Radar Signal Model and the -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
be the orthogonal narrow-band transmitted waveforms, which satisfy:
where
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
. Then, the
received data can be written as [
7]:
where
τ is the slow time index,
i.e., the pulse number (denoted as a snapshot in this paper), and
and
denote the reflection coefficient and Doppler frequency, respectively.
is the noise matrix. The transmit steering vector
and the receive steering vector
are shown as:
and:
Exploiting the orthogonality property of the transmitted waveforms, the
virtual data vectors can be obtained by matched filtering
with each of the orthogonal waveforms
, which is given as follows:
where
is the
m-th element of the transmit steering vector
, and
is the noise term. Then, stacking the individual vector components into one column vector, an
virtual vector can be obtained, which is written as [
7,
8,
9]:
where
,
is the signal vector with
.
is the noise vector. When the number of snapshots is
J, the Equation (
4) can be expressed as:
where
is the received data matrix.
and
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 are different, and the reflection coefficients obey the Swerling II model.
- (ii)
The noise matrix is zero-mean, complex circular Gaussian and with the variance .
- (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
in this paper. Firstly, we introduce the
-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
-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:
where
,
and
.
is composed with the singular vectors corresponding to the
P largest singular values of
. Then, Equation (
6) can be sparsely represented by using the sampling grid
, where
. The transmit and receive complete dictionaries are expressed as
and
, respectively. Thus, the complete dictionary can be constructed as
. Then, under the framework of the sparse representation, Equation (
6) can be formulated as:
where
and
have the same row support,
i.e., the matrix
is sparse. In order to estimate
, Equation (
7) can be formulated as the
-norm minimization problem, which is shown as follows:
where
is the regularization parameter. After obtaining the sparse matrix
, the DOAs can be estimated by finding non-zero rows of
. According to Equation (
8), the procedure of recovering
involves the MMV problem and complex-valued processing. Thus, the
-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.
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
, the computational complexity of solving SOC programming in Equation (
8) is
flops when
, where flops is the abbreviation of floating-point operations. For the the observing vector
in the proposed method, the computational complexity of solving the SOC programming in Equation (
27) is
due to
. 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
flops; (2) the unitary transformation in Equation (
17) and calculating real-valued covariance matrix in Equation (
21) require
flops; (3) solving the SOC programming in Equation (
27) with real-valued processing requires
flops. Thus, the computational complexity of the proposed method is
flops. For the complex-valued
-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
flops. For the real-valued (RV)
-SVD algorithm, it only requires the additions to transform the complex received data into a real one compared to the
-SVD algorithm. Thus, it requires
. Owing to
, the proposed method has lower computational complexity than the RV
-SVD and
-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
-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
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
columns of the complete dictionary
is independent, which leads to
, where
denotes the smallest integer of columns of the complete dictionary
that are linearly dependent. Then, the maximal number of identified targets is
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
exists with zero elements, and the vector
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
is a diagonal matrix. The matrix
, 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
-SVD algorithm [
15], the real-valued
-SVD algorithm (denoted as RV
-SVD) [
18] and the Cram
r–Rao bound (CRB) [
27]. In most of the simulations, a monostatic MIMO radar with
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
, where
and
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:
where
denotes the estimation of
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
-SVD algorithm, the RV
-SVD algorithm and the proposed method (
i.e.,
). The spatial grid is uniform from the range
to
with
interval in all simulations.
Figure 1 shows the spectra of these three algorithm, where
uncorrelated targets are assumed to be with the DOAs as
,
and
. The SNR and snapshots are fixed with 0 dB and 400, respectively. As can be seen from
Figure 1, compared to the RV
-SVD algorithm and the
-SVD algorithm, the proposed method has a lower sidelobe, which means that the proposed method has higher angular resolution.
Figure 1.
The spectra of -SVD, real-valued (RV) -SVD and the proposed method for uncorrelated targets.
Figure 1.
The spectra of -SVD, real-valued (RV) -SVD and the proposed method for uncorrelated targets.
Figure 2 and
Figure 3 show the angle estimation performance of these three algorithms
versus the SNR and snapshots, respectively, where
uncorrelated targets are assumed to be with the DOAs as
,
and
. 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
-SVD algorithm provides better angle estimation performance than the
-SVD algorithm at the low SNR region and comparable angle estimation performance at the high SNR region. Compared to the RV
-SVD and
-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
-SVD and
-SVD algorithms when the number of snapshots satisfies
. 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 -SVD, RV -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 2.
RMSE of -SVD, RV -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 3.
RMSE of -SVD, RV -SVD and the proposed method for uncorrelated targets when the snapshots number varies from 20 to 600.
Figure 3.
RMSE of -SVD, RV -SVD and the proposed method for uncorrelated targets when the snapshots number varies from 20 to 600.
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
and
, where
varies from
to
. 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
, where
and
is the estimation of
. As can be seen from
Figure 4, compared to the RV
-SVD and
-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 -SVD, RV -SVD and the proposed method for two uncorrelated targets when the angle separation varies from to .
Figure 4.
Target resolution probability of -SVD, RV -SVD and the proposed method for two uncorrelated targets when the angle separation varies from to .
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
and
. 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
-SVD and
-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
. From
Figure 6, it can be seen that the proposed method provides better angle estimation performance than the RV
-SVD and
-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
uncorrelated targets are assumed to be with the DOAs as
,
and
, 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 -SVD, RV -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 5.
Target resolution probability of -SVD, RV -SVD and the proposed method for uncorrelated targets when the SNR varies from −10 dB to 20 dB.
Figure 6.
RMSE of -SVD, RV -SVD and the proposed method versus the number of targets.
Figure 6.
RMSE of -SVD, RV -SVD and the proposed method versus the number of targets.
The angle estimation performance of the proposed method has been evaluated with the foregoing simulations. Here, the computational complexity of
-SVD, RV
-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
. 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
-SVD and RV
-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
-SVD and RV
-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.
Table 1.
Computation time comparison of -SVD, RV -SVD and the proposed method for uncorrelated targets when the targets number varies from 3 to 6.
Table 1.
Computation time comparison of -SVD, RV -SVD and the proposed method for uncorrelated targets when the targets number varies from 3 to 6.
Target Number | Average Computation Time (s) |
---|
-SVD | RV -SVD | Proposed Method |
---|
P = 3 | 3.7392 | 3.2885 | 2.4312 |
P = 4 | 7.8878 | 4.1707 | 2.5367 |
P = 5 | 9.6259 | 5.4507 | 2.7102 |
P = 6 | 11.1161 | 5.9464 | 2.7992 |