Abstract

Direction of arrival (DOA) estimation is an essential problem in the radar systems. In this paper, the problem of DOA estimation is addressed in the multiple-input and multiple-output (MIMO) radar system for the fast-moving targets. A virtual aperture is provided by orthogonal waveforms in the MIMO radar to improve the DOA estimation performance. Different from the existing methods, we consider the DOA estimation method with only one snapshot for the fast-moving targets and achieve the super-resolution estimation from the snapshot. Based on a least absolute shrinkage and selection operator (LASSO), a denoise method is formulated to obtain a sparse approximation to the received signals, where the sparsity is measured by a new type of atomic norm for the MIMO radar system. However, the denoise problem cannot be solved efficiently. Then, by deriving the dual norm of the new atomic norm, a semidefinite matrix is constructed from the denoise problem to formulate a semidefinite problem with the dual optimization problem. Finally, the DOA is estimated by peak-searching the spatial spectrum. Simulation results show that the proposed method achieves better performance of the DOA estimation in the MIMO radar system with only one snapshot.

1. Introduction

The orthogonal waveforms are transmitted by multiple-input and multiple-output (MIMO) radar systems [16], so the performance of target detection, estimation, and tracking can be improved by the waveform diversity. Generally, the MIMO radar system can be classified into two types: (1) the colocated MIMO radars [79], where the distance between antennas is comparable with the wavelength, and the waveform diversity is exploited to improve the target estimation performance with the large virtual aperture; (2) the distributed MIMO radar [1013], where the spacing between antennas is large, and the different view-angles from the antennas to targets are used to improve the detection performance with the diversity of radar cross-section (RCS). Since the direction of arrival (DOA) estimation problem is addressed in this paper, the colocated MIMO radar is adopted to improve the DOA estimation performance with larger virtual aperture than that in traditional phased array [14].

Many papers have investigated the DOA estimation problem [1518]. Traditionally, the discrete Fourier transform (DFT) [19] is used to estimate the DOA, where the received signals are sampled by the antennas in the spatial domain, and then the DOA estimation is equal to a corresponding frequency estimation in the transformed domain. Therefore, the frequency (DOA) in the spatial domain is obtained by the DFT methods, but the resolution of DFT method is limited by Rayleigh criterion. The methods that can break through the Rayleigh criterion are called super-resolution methods. Multiple signal classification (MUSIC) method [2022], Root-MUSIC [23], and the estimation of signal parameters via rotational invariant techniques (ESPRIT) method [2426] are three most essential super-resolution methods. The noise subspace and signal subspace are obtained in the MUSIC and ESPRIT methods to estimate the DOA, respectively. A TOD-MUSIC algorithm is proposed in [27] to estimate the DOA in the scenario with low signal-to-noise ratio (SNR) with diversity bistatic MIMO radar.

However, the subspaces are obtained from the estimated covariance matrix of the received signals, so the multiple measurements are needed in MUSIC and ESPRIT methods to achieve a reasonable estimation of the covariance matrix. In the MIMO radar system, the multiple measurements are obtained from the multiple pulses, where the measurements are obtained from the output of the pulse compression. In the scenario with fast-moving targets, the multiple measurements are impossible, so the subspace-based methods cannot be used. For the fast-moving targets, the DOA of targets is changing fast and exceeds the limit of the pulse repetition frequency (PRF) [28], so the statistic information, such as the covariance matrix of received signals, cannot be estimated well from the nonstationary pulse signals, and only the single pulse can be adopted to estimate the DOA. Therefore, the fast-moving targets are defined as the ones with nonstationary DOA between adjacent pulses.

The compressed sensing- (CS-) based methods have been proposed [29, 30] to improve the DOA estimation performance with fewer measurements, where the target sparsity is exploited in the spatial domain [3135]. Therefore, the DOA estimation for the sparse targets is transformed into the sparse reconstruction problem in the CS-based methods [36, 37, 37], and the sparse-based methods can be used in the scenario with fast-moving targets. In [38], a reweighted norm minimization method with fast iterations is proposed for the DOA estimation in MIMO radar. A fast sparse DOA estimation algorithm for both the white and colored Gaussian noises is proposed in [39] in the scenario with multiple measurement vectors. Xie et al. [40] developed a covariance-vector sparsity-aware estimator to estimate the DOA from the MIMO radar.

Usually, the CS-based methods can be classified into the greedy methods and the norm-based methods: (1) in the greedy methods, such as orthogonal matching pursuits (OMP) [41], stagewise OMP (StOMP), and CoSaMP [42], iterations are used to reconstruct the sparse signals; (2) in the norm-based method, the norm minimization problem is transformed into a norm minimization problem, which can be solved efficiently with the convex optimization tools. Additionally, sparse Bayesian learning- (SBL-) based methods with the prior assumption of sparse signals are also propped [43], such as SBL method and OGSBI method [44], which can achieve the excellent performance with relatively high computational complexity. Moreover, different from the traditional methods of the sparse reconstruction using the discretized dictionary matrix, the atomic norm-based method as a new type of the norm-based methods realizes the sparse reconstruction in the continuous domain [4547]. In [48], a semidefinite programming-based method is proposed for the norm optimization over infinite dictionaries. However, the atomic norm has not studied well in the MIMO radar system for the DOA estimation.

In this paper, the DOA estimation problem in the MIMO radar system is addressed, where a larger virtual aperture is provided by MIMO radar than that in phased array. To estimate the DOA of fast-moving targets, we consider the DOA estimation problem in the scenario with only one measurement from the MIMO radar system. Then, by exploiting both the structure of the MIMO radar system and the target sparsity in the spatial domain, a new type of atomic norm is proposed to obtain the trade-off between the target sparsity and the reconstructed signals. A least absolute shrinkage and selection operator (LASSO) with the novel atomic norm is formulated to denoise the received signals, but the LASSO cannot be solved efficiently, so a semidefinite matrix is constructed theoretically from the dual norm of the new atomic norm. Therefore, a dual optimization problem is obtained and transformed into a semidefinite problem, which can be solved efficiently. With the denoised signals, the DOA is estimated by peak-searching the spatial spectrum.

The remainder of this paper is organized as follows. The MIMO radar system is shown in Section 2. A continue domain super-resolution method for fast-moving targets is proposed in Section 3. Section 4 shows the simulation results, and Section 5 concludes the paper.

Notations. , and denote the norm and the norm, respectively. denotes the dual norm. denotes the real part of a complex value. denotes the Hermitian transpose of a complex matrix/vector.

2. MIMO Radar System with Fast-Moving Targets

In this paper, we consider the DOA estimation problem in the MIMO radar with fast-moving targets, and we assume that the number of transmitting antennas is M and that of receiving antenna is N. As shown in Figure 1, the orthogonal waveforms are transmitted in the MIMO radar system, and the waveform in the m-th transmitting antenna is denoted as (). Since the orthogonal waveforms are transmitted, we havewhere T is the pulse duration. The MIMO radar system considered in this paper is a type of colocated MIMO radar. Both the transmitting and receiving antennas are equipped in the same system, so the DOA and direction of departure (DOD) are the same. To simplify the analysis, the transmitting and receiving antennas are both uniform linear arrays (ULA). Therefore, in the n-th receiving antenna, the received signal can be expressed aswhere the number of targets is K, denotes the scattering coefficient of the k-th target, d is the distance between adjacent antennas, λ is the wavelength, is the DOA of the k-th target, and is the additive Gaussian noise.

In the n-th receiving antenna, the matched filters for orthogonal waveforms are adopted to distinguish the waveforms transmitted by different antennas. For the received signals in the MIMO radar system, if we use the signals directly, the SNR can be lower than 10 dB. However, we can use the pulse compression (i.e., matched filters) to increase the SNR of the received signals. For the m-th waveform , the corresponding matched filter is , where is a conjugate function of . Therefore, in the n-th receiving antenna, the signal passing the m-th matched filter and sampled at time T can be given as

We define the steering vectors with the direction θ for the transmitter and receiver, respectively, as follows:where the m-th entry of is and the n-th entry of is . Therefore, the received signal passing the matched filter can be simplified aswhere denotes the additive white Gaussian noise. Then, we can collect all the signals in the n-th receiving antenna after the matched filter as follows:where the noise vector is defined as . Collect the signals from all the receiving antennas, and we can get the vector with all the information aswhere we define , and denotes the Kronecker product. Additionally, we define the matrix with steering vectors asand we define the vector for target scattering coefficients as

Finally, from the system model (7), we try to estimate the DOA of targets from the received signals with only one snapshot. Usually, in the existing methods with super-resolution, the multiple measurements are used, for example, the MUSIC or ESPRIT methods; the multiple measurements are used to estimate the covariance matrix, and the DOAs are estimated from the noise or signal subspaces, which are obtained from the covariance matrix. In this paper, we focus in the scenario with only one measurement, so the proposed DOA estimation method can be used to track the fast-moving targets and can get the target information from just one pulse.

3. Continue Domain Super-Resolution Estimation Method

In this section, we will propose a super-resolution method for the DOA estimation in the continuous domain, named as continue domain super-resolution (CDSR) method. To exploit the target sparsity in the spatial domain, we will formulate a system model from (7) by introducing the sparsity. In the traditional methods to exploit the target sparsity, the compressed sensing- (CS-) based methods have been proposed [49], where the dictionary matrix is formulated by discretizing the spatial domain. Then, the DOA information is obtained from the dictionary matrix. However, the off-grid problem will be introduced during discretizing the spatial domain [50, 51], where the targets cannot be precise at the discretized angles. Different from the traditional methods, we will propose a sparse-based method in the continuous domain. First, the preliminary for the CDSR will be given as the background knowledge.

3.1. Preliminary of Dual Norm

For a norm ( is a vector with the entries being complex), the dual norm of is defined as

Therefore, for the traditional LASSO problem [52], we havewhere is an unknown vector and is a known vector. is a measurement matrix. The LASSO problem tries to estimate the vector from the known vector and measurement matrix . The first term describes the accuracy and the second term is for the reconstruction sparsity. The parameter λ is used to control the tradeoff between the reconstruction accuracy and the sparsity of .

Then, the optimization problem in (11) can be rewritten asas

To solve this optimization problem, the corresponding Lagrangian can be expressed asaswhere is the vector of Lagrangian parameter, and denotes the real part of a complex value. Therefore, the dual optimization problem of (12) can be obtained as

Additionally, we have

Let , and we can obtain

Moreover, from the definition of dual norm, we can obtain [53]where denotes the dual norm of norm. is an indicator function, which is defined as

Therefore, substitute (16) and (17) in (14), and we can simplify the dual optimization problem as follows:which can be rewritten as

The original LASSO problem is transformed into the dual optimization problem using the dual norm, which can be solved more efficiently.

3.2. CDSR Method

We will define a new type of atomic norm, which can exploit the target sparsity in the spatial domain and also use the MIMO advantage. The new type of atomic norm [5457] for MIMO radar is named as MIMO Atomic Norm (AN-MIMO), which is defined aswhere we use the subscript “AN-MIMO, 0” to show the decomposition for the signal with minimum number of terms. The AN-MIMO shows that we can use the minimum number of terms formulated by to describe the signal . However, as show in (21), we must determine the number of terms, i.e., the parameter P, to obtain the decomposition, which is nondeterministic polynomial (NP) hard and cannot be solved easily.

Inspired by the norm used in CS theory, we can relax the original AN-MIMO by the following equation:where we use the footnote “AN-MIMO” to show that the relaxation equation is different from the original one in (21). The AN-MIMO is described by the sum of absolute values and can be solved by the following semidefinite programming- (SDP-) based method. The number of terms P can be obtained by searching the number of peaks in .

With the AN-MIMO, the received signal can be approximated by a signal with the sparse consideration. We can formulate the following optimization (denoising) problem aswhich is a type of LASSO problem, but we use the AN-MIMO to replace the traditional norm. As shown in (22), the AN-MIMO is defined as the sum of the absolute values , which are the coefficients in the atomic decomposition. In the norm minimization problem with the vector , the norm is defined as the absolute sum of the entries in . Therefore, the AN-MIMO norm is similar to the norm, and both use the sum of the absolute values. Then, the method used in the norm minimization problem can be introduced in the AN-MIMO problem. Similar to the LASSO problem in Section 3.1, the denoising problem in (24) can be simplified aswhere denotes the dual norm of AN-MIMO .

Proposition 1. The optimization problem in (23)with the dual norm of AN-MIMO can be solved by the SDP method.

where denotes the approximation of received signal with the sparse constraint. The matrix is a semidefinite matrix and is composed of the block matrix . For the diagonal block matrices (), the sum of diagonal entries is , and the sum of subdiagonal entries is 0. However, for the off-diagonal block matrices (), both the sum of the diagonal entries and that of the subdiagonal entries are 0.

In the optimization problem (25), the parameters , , and are optimized and solved by minimizing . The optimal is denoted as the denoised signal from . The optimization problem in (25) is a type of SDP problem and can be solved by the interior-point method. According to Section 6.4 in [58], the computation complexity is using the interior-point method to solve the SDP problem (25).

After obtaining the denoised signal , the DOA can be estimated by the spatial spectrum, which is calculated by .

3.3. Proof of Proposition 1

With the definition of dual norm in (10), the dual norm of AN-MIMO can be expressed as

The equation (a) is obtained from the definition of AN-MIMO, and the equation (b) is obtained by choosing the DOA to maximize the correlation between and . The equation (c) is obtained by letting

Therefore, based on (26), the constraint in the dual problem (24) can be simplified as

This inequality can be further simplified as a SDP problem. Before getting the SDP problem, we first introduce the Schur complement as follows.

Lemma 2. Given a block matrix , the matrix is a semidefinite matrix if and only if the following conditions are satisfied: the block is a semidefinite matrix and is also a semidefinite matrix.
With the Schur complement, we can formulate a semidefinite matrix ( and ), then we haveHence, from the definition of semidefinite matrix, we can find that for any vector , we have , which can be rewritten asBy choosing , (30) can be simplified aswhere the left side is equal to the square of left side in (28).

Therefore, if we can let , the constraint (28) can be satisfied. Then, the dual optimization problem (24) can be simplified as

This optimization problem is a type of SDP problem, but cannot be solved directly, so we need to further simplify the expressions. Let , and we need to formulate the block matrix to satisfy

We can express the left side as

Therefore, we can formulate the matrix as a following block matrixwhere ( and ). For all blocks, has the following characteristics:where denotes the entry of at the -th row and -th column. When the characteristics in (36) is satisfied, we can find that the result in (34) is 1. Therefore, the condition in (33) is also satisfied.

Finally, the optimization problem in (32) is simplified as

4. Simulation Results

We will show the simulation results in this section, where the simulation parameters are shown in Table 1. We consider the DOA estimation problem in the scenario with the MIMO radar system, where the transmitter and the receiver are colocated. Without additional explanation, the default parameters are the same with Table 1. All the simulation results are obtained from a Work Station with 64G Byte RAM and Intel i7 CPU (3.40 GHz). The MATLAB code for the proposed method is available online: https://drive.google.com/drive/folders/1_iFF1XEQhKOl7-6Y9uKcy8kK6h6tRbtD?usp=sharing.

In this section, three state-of-art methods are compared with the proposed method:(i)Direct Method. In this method, the correlation between the received signal and the steering vector is calculated as to get the spatial spectrum for the DOA estimation.(ii)OMP Method. Orthogonal matching pursuits (OMP) method [59] is a typical method for sparse reconstruction. Since our method is a sparse-based DOA estimation method, this comparison is needed.(iii)SBL Method. Sparse Bayesian learning (SBL) method [44] is a Bayesian-based method for the sparse reconstruction. Usually, this method can achieve better reconstruction performance, but with relatively high computational complexity.

First, we show the estimated spatial spectrum for the DOA estimation in Figure 2, where the number of receiving antennas is . The estimated DOAs using different methods are given in Table 2. The parameter λ is chosen as . As shown in this figure, the DOAs of targets are, respectively, and When we use the proposed method, we can get the estimated DOAs as and To describe the DOA estimation performance, we use the root mean square error (RMSE) to measure the errors, and the RMSE is defined aswhere denotes the number of Monte Carlo simulations and in this paper. is the estimated DOA during the -th Monte Carlo simulation, and is the ground-truth DOA. We can calculate the RMSE (in degree) using the proposed 0.1300. When the direct method is used to estimate the DOA, the RMSE is 0.4508. The RMSEs of OMP and SBL methods are, respectively, 0.4082 and 0.3985. We can find that the proposed method achieves better DOA estimation performance.

Second, in the scenario with less receiving antennas, we give the simulation results for the DOA estimation in Figure 3, where the number of receiving antennas is . As shown in this figure, we can find that the spatial spectrum is worse than the scenario with receiving antennas (Figure 2). The corresponding DOA estimation results are given in Table 3. The RMSE of the proposed method is 0.7009 in degree. The RMSEs of existing methods including the direct method, OMP method, and SBL method are 1.0756, 1.0591, and 0.7623, respectively. Therefore, in the scenario with less receiving antennas, the proposed method also achieves better performance of the DOA estimation.

Then, Figures 4 and 5 are, respectively, the DOA estimation performance with different SNRs of the received signals, where the number of receiving antennas is in Figure 4 and that in Figure 5 is . From these figures, we can find that the proposed method can achieve better performance of DOA estimation in the scenarios with relative high SNR of received signals. Additionally, with , the proposed method has better performance at SNR>25 dB, but the better performance can be achieved at SNR>15 dB with . Therefore, with more antennas, the performance improvement is significant at higher SNR. The corresponding Cramér-Rao lower bound (CRLB) for the DOA estimation can be obtained from [60, 61]. In both Figures 4 and 5, the CRLBs for the DOA estimation are shown. As shown in the figures, the proposed method can approach the CRLB in the scenario with higher SNR ( dB).

Finally, the DOA estimation performance with more transmitting antennas is shown in Figure 6, where 10 transmitting antennas are adopted. As shown in this figure, we can find that the proposed method achieves better DOA estimation performance at SNR>20 dB. Additionally, the performance improvement is more significant in the scenario with higher SNR of the received signals. Compared with Figure 5 using only 5 transmitting antennas, the DOA estimation performance is improved with more antennas. Therefore, the proposed method can achieve better performance of DOA estimation in the scenario with higher SNR, especially using fewer antennas. In Table 4, the computational time is given for the different methods. The computational time of the proposed method is 0.6037 s. The direct and OMP methods have lower computational time than the proposed method. The SBL method has higher computational complexity with the computational time being 1.0702 s. Therefore, the proposed method can achieve better performance for the DOA estimation with the acceptable computational complexity.

5. Conclusion

The DOA estimation problem in the MIMO radar system with fast-moving targets has been addressed, and the new type of atomic norm has been formulated to measure the target sparsity in the spatial domain. The received signals have been denoised by the LASSO-based model. Then, the semidefinite matrix has been constructed to transform the atomic norm minimization problem into the SDP problem, which can be solved efficiently. Simulation results show that the proposed method achieves better DOA estimation performance in the MIMO radar with fast-moving targets. The future work will focus on the radar system optimization to further improve the DOA estimation performance.

Data Availability

The simulation data used to support the findings of this study are included within the article. Additionally, the Matlab codes for the simulations are available online: https://drive.google.com/drive/folders/1_iFF1XEQhKOl7-6Y9uKcy8kK6h6tRbtD?usp=sharing.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (Grant no. 61976113).