Abstract

This paper proposes a method to address the problem of the joint direction of departure (DOD) and direction of arrival (DOA) estimation with augmented coprime array (CPA) bistatic multiple input multiple output (MIMO) radar using interpolating sensors. At first, we deduce the regular pattern of hole positions in the virtual array and interpolate a small number of sensors to augmented CPA to form a partially contiguous virtual array. Then, we use the diversity smoothing algorithm to reconstruct the Toeplitz matrix to obtain a spatial smoothing matrix. Finally, we combine the RD-MUSIC algorithm with the spatial smoothing matrix to estimate the spatial spectrum and achieve automatic matching of DODs and DOAs for the targets. Simulation results clearly illustrate the superiority of the method.

1. Introduction

Direction of arrival (DOA) estimation and Kalman filter algorithm for target tracking [1] are significant areas of research in array signal processing. It can estimate the angular position of different signals in a certain airspace and plays a key role in radar, sonar, and other target detection fields. Since the last century, DOA estimation has gone through three stages, which are beamforming [2], subspace fitting [36], and compressed sensing [7], and made a great contribution to the development of the DOA estimation algorithm. Bistatic multiple input multiple output (MIMO) radar has the advantages of waveform diversity [8] and spatial separation; it estimates target angle including DOA and direction of departure (DOD), so the joint estimation of DOD and DOA [9, 10] has also become the focus of research. Although the uniform linear array (ULA) can also solve the related problems of DOA estimation, it still has certain limitations. For example, the number of physical array sensors usually cannot be less than the number of signal sources; otherwise, the estimation accuracy will be badly affected, and thus the degree of freedom (DOF) is limited by physical array elements and other factors [11]. In order to solve these problems, a sparse array is introduced by scholars. The difference between the sparse array [12, 13] and the traditional uniform array mainly includes that the sparse array is formed by setting different interelement spacings in the array to form a sparse array structure, and part of the sparse array has a larger element spacing, so it forms longer virtual array to increase DOF. Coprime array (CPA) [14, 15] has been widely studied by scholars because of its larger element spacing. Thus, the CPA has a low mutual coupling effect and high accuracy of target angle estimation. Moreover, the CPA can construct a longer virtual array through its own differential coarray [16] and use its equivalent vector to estimate DOA. However, although the virtual array formed by CPA has a long virtual array aperture [17], there are array holes that make the virtual array discontinuous, which limits the expansion of DOF. Scholars have conducted a series of studies to settle the question.

A method was advanced to fill the virtual array hole by moving the sparse array in [18], and it used the differential coarray of the original array and the differential coarray of the moved array to form a composite array without holes to estimate DOA effectively. But this method has a small estimated error that originated from array movement distance and signal source movement distance, and its theoretical feasibility is insufficient. Interpolated virtual array element was considered in [19] and solved the problem of virtual array discontinuity caused by virtual array holes, but the method will be affected when there are continuous virtual holes, and its DOF will reduce. Some scholars utilize the virtual array element interpolation theory [19] to the joint estimation of DOD and DOA of bistatic radar targets; this method reduces redundancy of the virtual array and improved DOF and estimated target resolution, but this method required more calculations. Combined with the traditional ESPRIT algorithm [20], Li et al. proposed a joint estimation of target DOD and DOA with bistatic CPA MIMO radar based on virtual aperture expansion in [21]. This method achieved better estimation performance than the traditional method, but it did not use discrete virtual array elements. Recently, a method which interpolated array elements in CPA was proposed to solve the virtual array hole problem in [22], but augmented CPA which the method used cannot nicely reflected the advantage of the interpolated array elements to expand the virtual array aperture. We propose a joint DOD and DOA estimation of bistatic MIMO radar for coprime array based on array elements interpolation. The method uses augmented CPAs as the transmit array and receive array and interpolates a small number of sensors to the particular holes in the virtual array to expand the aperture. Then, the method uses the selection matrix to reconstruct the Toeplitz matrix based on diversity smoothing to estimate the DODs and DOAs of the sources.

The remaining sections are as follows: in Section 2, we reduce the math model of the bistatic CPA MIMO radar. We deduce the law of hole position and propose a diversity smoothing algorithm for reconstructing the Toeplitz matrix to estimate the DODs and DOAs of the targets in Section 3. We perform the same simulation experiments with the proposed method and other methods, which clearly illustrates the superiority of the method in Section 4. Finally, we present conclusions for this paper in Section 5.

Notations: we use italicized boldface characters to represent vectors and matrices in this paper. Superscripts and represent transpose and conjugate-transpose, respectively, diag [·] denotes diagonal matrix, and and denote the Kronecker product and the Hadamard product, respectively.

2. Mathematical Model

The conventional bistatic CPA MIMO radar model is presented in Figure 1. Transmit array and receive array in the model are augmented CPA which consists of two uniform linear arrays (ULA). One ULA has 2M1 − 1 sensors in transmit array, and red circles represent sensors of subarray 1; another ULA of transmit array has N1 sensors, and the black circle represents sensors of subarray 2. Similarly, receive array contains two ULAs, which, respectively, have 2M2 − 1 sensors and N2 sensors, where M1 and N1 are two coprime integers and M2 and N2 are two coprime integers. The number of sensors of arrays are, respectively, M = 2M1 + N1 − 1 and N = 2M2 + N2 − 1. The unit interelement spacing of the array is d, which is the half wavelength (λ/2). The sensor positions are given by the following equation:

Suppose there are K uncorrelated signals in the space, the DOD and DOA of signals are given by and . represents the position of the sensors in the transmit array, and denotes the position of the sensors in the receive array. Set the reflection coefficient of signals as , , and L represents the number of snapshots. The received signal after matched filtering is given by the following equation:where n(t) represents a matrix composed of Gaussian white noise and it follows the Gaussian distribution . and are also given by the following equation:where and , respectively, denote manifold matrices, and they are given by the following equation:

Therefore, the covariance matrix of the received signal R is given by the following equation:where represents an MN × MN dimensional identity matrix and denotes the covariance matrix of the received target.where represents the power of the kth target signal.

3. Joint Diversity Smoothing DOD and DOA Estimation Algorithm Based on Interpolated Sensors

3.1. Expansion of Virtual Array by Interpolating Sensors to CPA Holes

In this section, we will interpolate sensors to CPA holes to expand a longer virtual array aperture and then we can obtain a new equivalent vector to execute transmit-receive diversity smoothing. Assume transmit array and receive array are identical augmented CPA, and taking transmit array as an example, sensors distribution is illustrated in Figure 2.

The virtual array is formed by the sum-difference array of CPA, and the position of virtual sensors is given by the following equation:

According to (8), the virtual sensors are only generated by the difference coarray of transmit array and the difference coarray of receive array when DOD and DOA are estimated, respectively. Therefore, corresponds to parameter and corresponds to parameter .

The position of the virtual sensors that are formed by the difference coarray of transmit array is given by the following equation:where the value range of is [−(3M1N1 − M1 − N1), 3M1N1 − M1 − N1].

Similarly, the position of the virtual sensors that are formed by the difference coarray of receive array is given by the following equation:where the value range of is [−(3M2N2 − M2 − N2), 3M2N2 − M2 − N2].

For transmit array, M1 and M2 can determine the distribution of the virtual sensors and the length of the virtual array aperture, and the distribution of virtual sensors that are formed by augmented CPA is shown in Figure 3.

In Figure 3, black filled circles represent virtual sensors and dotted circles represent virtual holes. Although augmented CPA can form more virtual sensors than traditional CPA, some virtual holes are located next to the center point. Therefore, other continuous virtual sensors cannot form a virtual array with a larger aperture.

With the purpose of settling the question, we present a method that utilizes a small number of sensors to interpolate a part of virtual holes. The originally discontinuous virtual array will become continuous to achieve more DOF.

We can know a special kind of holes next to the center point in Figure 3 and call them central virtual holes. Central virtual holes are continuous virtual holes which locate next to the center point. When M1 = 2 and M2 = 3, there are two central virtual holes in the virtual array; when M1 = 3 and M2 = 4, there are four central virtual holes in the virtual array; and when M1 = 3 and M2 = 5, there are four central virtual holes in the virtual array.

Assume we use sensors to interpolate the virtual central holes where next to the center point completely. Then, a new distribution of virtual sensors emerges in Figure 4.

From Figure 4, black filled circles represent virtual sensors, dotted circles represent virtual holes, red filled circles represent sensors, and red circles represent new formed virtual sensors. We found that we only need to interpolate the virtual holes located next to the center point to make the virtual array continuous. For example, when M1 = 3 and M2 = 4, we use two sensors to interpolate the virtual central holes which are located in {1d,2d}. Then, the array can form additional new virtual sensors which are located in {−1d,−2d,±5d}. Most of the original virtual holes are interpolated. Comparing Figure 3 with Figure 4, it can be found that the interpolated virtual array is a continuous virtual array that all the central virtual holes are replaced with sensors or new virtual sensors, and the DOF is 47. Therefore, we deduce the number of central virtual holes and DOF for different augmented CPAs from continuous experiments and find that the number of virtual holes to be interpolated is only related to M1. The derived regular pattern and formulas are shown in Table 1.

In general, we interpolate (M1 − 1) sensors into the virtual holes which located next to the center point, the original virtual array becomes continuous, virtual array apertures become larger to (4M1N1 − 1), and DOF also becomes larger to (4M1N1 − 1).

3.2. Reconstructing Toeplitz Matrix Algorithm Based on Diversity Smoothing

In recent years, scholars have proposed some spatial smoothing algorithms for bistatic radars in [2325]. This paper chooses an algorithm that reconstructs the Toeplitz matrix based on diversity smoothing.

After the sensors shave been interpolated, there are continuous virtual sensors in the virtual array of transmit and continuous virtual sensors in the virtual array of receive, and denotes the sensors’ position in the virtual transmit array and denotes the sensors’ position in the virtual receive array. Therefore, the steering matrices of transmit and receive virtual array are given by the following equation:where and , respectively, denote manifold matrices, and they are given by the following equation:

A new covariance matrix of the received signal can be formed by processing continuous virtual sensors. However, the ideal is hard to get in practice, so is usually estimated by using L available snapshots and it is given by the following equation:where is the new received signal of the virtual array.

The new covariance matrix is now vectorized to obtain a new equivalent vector as follows:where and is a vector whose elements are all zeros except for the row.

In fact, we regard vector as the received signal of a single snapshot in bistatic MIMO radar. According to the principles of spatial smoothing algorithms, the number of smoothing times should be no less than the number of overlapping subarrays to form a full rank covariance matrix. Since the number of virtual sensors is always an odd number, the best smoothing results are achieved when the number of smoothing times is equal to the number of overlapping subarrays. Therefore, we divide the virtual array of transmit into Ms overlapping subarrays to carry out Ms smoothing and divide the virtual array of receive into Ns overlapping subarrays to carry out Ns smoothing. Thus, and .

Next, suppose there are two selection matrices and corresponding to the virtual array of transmit and the virtual array of receive, which are given by the following equations:where and . denote the zero matrix, denote identity matrix, and denote the zero matrix. is similar.

Utilize selection matrices and to form the transmission smoothing and receiving smoothing on the equivalent vector . Then, we can change the values of m and n to achieve the smoothing of overlapping subarrays in the virtual array, m and n represent smoothing times. So subvector is given by the following equation:where , , , , and is the corresponding noise vector. The reconstructed Toeplitz matrix is given by the following equation:where represents the new steering matrix.

Therefore, the covariance matrix of each sub-array is given by the following equation:

The values of m and n are changed to achieve the effect of spatial smoothing in the virtual array and the spatial smoothing algorithm is used to solve the single snapshot problem caused by the vectorized covariance matrix. The spatial smoothing matrix is given by the following equation:

Through comparison, we find that and Toeplitz matrix have the same form, so the spatial smoothing matrix is also given by the following equation:

When performing spatial spectrum estimation, we choose the RD-MUSIC algorithm instead of the ESPRIT algorithm. Although we use the ESPRIT algorithm for spatial spectrum estimation to effectively reduce the complexity of the algorithm, the estimated DODs and DOAs of multiple targets need to be matched manually, whereas the RD-MUSIC algorithm can automatically match the DODs and DOAs of targets.

The RD-MUSIC algorithm is combined with to estimate the spatial spectrum, and the spatial spectrum function is given by the following equation:where denotes signal subspace of the covariance matrix. is also given by the following equation [25]:

Therefore, the spatial spectrum function is also given by the following equation:where . The DOA of the kth signal is given by the following equation:where .

Similarly, The DOD of the kth signal is given by the following equation:

The steps of the proposed Algorithm 1 are as follows.

Input: receive signal: x(t) = As(t)+n(t), t = 1,2, …, L;
Output: DODs and DOAs: ;
Step:
(1) estimate covariance matrix based on continuous virtual array as in (12);
(2) build the selection matrix of the transmit array as in (14) and the selection matrix of receive array as in (15);
(3) rebuild a new Toeplitz matrix as in (17);
(4) form a spatial smoothing matrix as in (18);
(5) utilize the RD-MUSIC algorithm, to estimate DODs and DOAs ;

4. Simulation Results

We assume the distance of array sensors M1 = M2 = 2, N1 = N2 = 3, the number of transmit and receive array sensors M = N = 6, and the original positions of transmit and receive sensors are [−9d, −6d, −3d, 0d, 2d, 4d]. The number of central virtual holes to be interpolated is Hf = 1, and the positions of transmit and receive sensors after interpolating the holes are [−9d, −6d, −3d, 0d, 1d, 2d, 4d]. There are continuous virtual sensors in transmit array and receive array, overlapping subarrays of transmit and receive array, and smooth times are Ms = Ns = 12. The signals in the simulation are incoherent.

Compare the performance in target estimation of four methods which include diversity smoothing algorithm based on interpolating sensors (DSIAS) in this article, interpolation virtual array element algorithm (IVAE) in [19], conventional CPA diversity smoothing algorithm (CDS) in [21], and conventional ULA RD-MUSIC algorithm in [26]. It includes three specific simulation experiments, namely, target number detection and angular resolution. In the simulation of root mean square error (RMSE), we not only compared the above-given methods, we also added a unitary dual-resolution ESPRIT (U-ESPRIT) method [27] to the comparative simulation.

4.1. Number of Detectable Targets

In this part, the signal-to-noise ratio (SNR) is set as 10 dB and the number of snapshots is set as 200. There are 15 signal targets distributed over the range [−70°, 70°], where located at [−70°, −60°, −50°, −40°, −30°, −20°, −10°, 0°, 10°, 20°, 30°, 40°, 50°, 60°,70°] and [−70°, −60°, −50°, −40°, −30°, −20°, −10°, 0°, 10°, 20°, 30°, 40°, 50°, 60°,70°]. Figure 5 shows the results of DOD and DOA joint estimation of targets. The four pictures in Figure 5 represent the algorithms in this article, [18, 20], and [25], respectively. As shown in Figure 5, it shows the contour map of the spatial spectrum peak, red circles represent the true direction of 15 targets and the spectral peak contour represents the estimated direction of 15 targets. The spectral peak contour completely overlaps in the red circle in Figure 5(a), so the DSIAS algorithm can accurately estimate 15 targets. Other methods’ spectral peak contours do not overlap completely in the red circle. The DSIAS algorithm is better than other methods for a number of detectable targets.

4.2. Angular Resolution

In this part, with the purpose of comparing the performance of the four algorithms in angular resolution, assume that the SNR is 10 dB and the number of snapshots is 200. Figure 6 shows situations of the angular resolution comparison of different algorithms. The position of two adjacent targets are and , and the red line denotes the real target's direction. From the simulation results, we can also see that the MIMO radar which uses the DSIAS algorithm can estimate the targets well under the condition of two targets close to each other, and the other three algorithms cannot accurately distinguish two similar targets.

4.3. Root Mean Square Error (RMSE)

We compare the RMSEs of different algorithms. RMSE is a common standard that reflects the accuracy of angle estimation, and the average RMSE is defined by the following equation:where Q denotes Monte Carlo simulation times, K denotes the number of targets, and denotes the joint estimated DOD and DOA of the kth target for the ith Monte Carlo simulation, i = 1,2, …, Q. Estimation number of targets will affect the estimation accuracy. In order to reflect the comprehensiveness of the simulation results, we carried out simulation experiments for two targets and four targets, respectively, and the results of the simulation experiments are shown in Figures 7 and 8.

Figure 7 shows the relationship of RMSE with SNR and the number of snapshots under the condition of detecting two targets, and the targets locate at , . Figure 7(a) depicts the variation of the RMSE curve with SNR, where the number of snapshots is set as 200. The DSIAS algorithm has a higher estimation accuracy than other methods at low SNR but has similar performance with the other two methods [19, 21] at high SNR. Figure 7(b) depicts the variation of the RMSE curve with a number of snapshots, where the SNR is set as 10. The DSIAS algorithm has high estimation accuracy at a different number of snapshots.

As shown in Figure 7, the accuracy of all four algorithms is high when estimating two target angles, and the estimation accuracy of the DSIAS algorithm proposed in this paper is only slightly higher than the others. Figure 8 presents the relationship of RMSE with SNR and the number of snapshots under the condition of detecting four targets, and the targets locate at , , , and . Figure 8(a) depicts the variation of the RMSE curve with SNR, where the number of snapshots is set as 200. In contrast to the previous experiments, the DSIAS algorithm clearly performs better than other methods in estimation accuracy. Figure 8(b) depicts the variation of the RMSE curve with a number of snapshots, where the SNR is 10. It can be found that the DSIAS algorithm estimates more targets with greater accuracy by comparing two experiments which detects a different number of targets.

5. Conclusions

We propose a joint estimation of the DOD and DOA method that interpolate a small number of sensors to a specific location in an augmented CPA virtual array. The method can expand the aperture of the virtual array and maintains the maximum DOF of augmented CPA despite the virtual holes which cannot be exploited in the virtual array. Meanwhile, we reconstructed the Toeplitz matrix based on diversity smoothing to obtain a spatial smoothing matrix. Finally, we combined the RD-MUSIC algorithm with the spatial smoothing matrix to estimate the spatial spectrum and accurately estimate the DOD and DOA of the targets. Simulation results illustrate the proposed method has better performance than other methods.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request and with permission from the funders.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China under Grant No. 62071476.