Abstract

In bistatic multiple-input multiple-output (MIMO) radar, range migration and invalidly synthesized virtual array resulting from the serious mismatch of matched filter make it difficult to estimate direction of departure (DOD) and direction of arrival (DOA) of high speed target using the traditional superresolution algorithms. In this study, a method for joint DOD and DOA estimation of high speed target using bistatic MIMO radar is proposed. After multiplying the received signals with the conjugate of the delayed versions of the transmitted signals, Fourier transform (FT) of the multiplied signals over both fast time and slow time is employed. Then, the target components of radar return corresponding to the different transmitted waveforms can be perfectly separated at the receivers by extracting the target frequency-domain data along slow-time frequency dimension when the delay between the transmitted signals and their subsequent returns is timed. By splicing the separated target components distributed along several range cells, the virtual array can be formed, and then DOD and DOA of high speed target can be estimated using the superresolution algorithm with the range migration and the mismatch of matched filter properly removed. Simulation results have proved the validity of the proposed algorithm.

1. Introduction

Bistatic radar with the widely separated transmitter and receiver has the advantages of covert receivers, giving increased immunity to antiradiation missiles and electronic countermeasures, and possible antistealth capabilities [1]. However, the need for synchronization between the transmitter and receiver leads to an increase in complexity. Unlike bistatic phased-array radar, bistatic multiple-input multiple-output (MIMO) radar is characterized by using multiple transmit antennas to simultaneously transmit orthogonal waveforms. Therefore, bistatic MIMO radar will simultaneously illuminate a very broad angular sector instead of a focused angular sector. The virtual array can be formed by matching the received signals with the transmitted signals, and then the direction of arrivals (DOAs) and the direction of departures (DODs) of targets can be estimated by processing the output of virtual array using the existing angle estimation algorithms, that is, the algorithms proposed in [24]. Due to the broad illuminated angular sector and the target location determined by its DOA and DOD in this case, the beam scanning and range information of the target are redundant and then the time and space synchronization are relaxed in bistatic MIMO radar [2, 3].

With the military technology evolution and the increasing exploitation of space resources, high speed maneuvering targets, that is, aircraft, missile, and space debris, pose severe challenges to modern radar in recent years. Joint DOD and DOA estimation methods have been widely investigated in bistatic MIMO radar [28]. But all of them, which are designed for low speed target, will be seriously deteriorated when dealing with the target moving at a high speed. The successful virtual array formation and effective target energy accumulation during the long observation time are vital for angle estimation in bistatic MIMO radar. However, the variety of the phase caused by the large-scale Doppler frequency of high speed target within the repetition interval of radar signal results in the severe distortion of target echo. Then significant mismatch loss would occur in matched filters. That is, the virtual array cannot be effectively formed for high speed target [9]. Meanwhile, during the long observation time, the target with high speed easily goes through several range cells so that the target energy is distributed along multiple range cells [1014]. As a consequence, for high speed target, the angle estimation performances of the existing methods, that is, the proposed methods in [68], are limited with both the mismatch loss in matched filter and the target residence time in a single range cell. Persy and Dipietro [15] have introduced the keystone transform to compensate for the linear range migration during the long integration time for synthetic aperture radar (SAR) ground moving targets imaging without knowing the target motion information. For high speed target, it is highly possible that both Doppler frequency ambiguity and target radial acceleration exist. However, in this situation, the detection performance of the keystone transform would be seriously reduced. At the same time, the data interpolation operator is employed to implement the keystone transform, which results in the large computational load. Dorp [16] introduced the keystone transform into MIMO radar for correcting the envelope migration of high speed target. It is worth pointing out that time division, rather than code division, strategy for transmitting signals is applied in MIMO radar to synthesize the virtual elements one by one at different time. The basic idea of the transmitter time division is to avoid the effect of the large-scale Doppler frequency on the virtual array formation. However, for high speed target, the switching time between the transmitters is normally longer than the inverse of the relative Doppler frequency of the target. Thus, the phase variety of the return signal caused by the target’s Doppler frequency during the virtual array forming time alters the array manifold and then degrades the accuracy of the angle estimation. In [17], a long-time integration method based on Hough transform is proposed to detect the weak target without acceleration. This method is capable of overcoming the problem of the range migration, but it performs poorly under the condition that the virtual array is unsuccessfully formed because of the large-scale Doppler frequency of target, and its computational burden is also very huge. A novel long-time coherent detection method, Radon-Fourier transform (RFT), is proposed in [1820], which realizes the echo spatial-temporal decoupling via joint searching along range and velocity dimensions. RFT may obtain the significant coherent integration gain without the Doppler ambiguity restriction, but, for the high speed maneuvering target detection using bistatic MIMO radar, a huge computational complexity is needed for five-dimensional joint searching of range, velocity, acceleration, DOD, and DOA. In [21], the space-time RFT for wideband digital array radar is proposed for high speed target detection by jointly realizing digital beamforming (DBF), range compression, and long-time coherent integration. However, prior to the spatial and temporal information processing for bistatic MIMO radar, the virtual array should be effectively formed when dealing with the target moving at a high speed.

In bistatic MIMO radar, suffering from the range migration and the invalidly synthesized virtual array when dealing with high speed target, the performances of the existing angle estimation methods will rapidly degrade. This paper makes an effort to provide a method for overcoming the effect of the large-scale Doppler frequency on the virtual array formation and removing range migration during the long-time integration, thus facilitating the angle estimation for high speed target in bistatic MIMO radar. In order to overcome the effect of the large-scale Doppler frequency on the virtual array formation, the two-dimensional (2D) Fourier transform (FT) is employed to transform the multiplied signals, which are obtained by multiplying the received signals with the conjugate of the delayed versions of the transmitted signals, into the fast-time frequency and slow-time frequency domains, and then the target components of the radar return corresponding to the different transmitted waveforms perfectly separated by extracting the target frequency-domain data. By splicing the separated target components that are distributed along several range cells, the virtual array is synthesized with the range migration removed. Thus, the DOD and DOA of high speed target can be obtained from the output of virtual array using the superresolution algorithm. Simulation results demonstrate the effectiveness of the proposed method.

The remainder of this paper is organized as follows. In Section 2, the signal model for high speed target in bistatic MIMO radar is established. In Section 3, the proposed angle estimation algorithm for high speed target is described. Moreover, the simulation results of the proposed algorithm are presented and the performances are investigated in Section 4. Finally, Section 5 concludes the paper.

2. Signal Model

Consider a bistatic MIMO radar system [2] with closely spaced transmit elements and closely spaced receive elements, as shown in Figure 1. Both the transmit array and the receive array are uniform linear arrays and the spacing between adjacent elements of the transmit array and the receive array is denoted by and , respectively. We assume that the range between the target and the transmit array or receive array is much larger than the aperture of the transmit array or receive array. different continuous coded periodic signals with identical bandwidth are employed, which are temporally orthogonal. The vector including orthogonal transmitted baseband coded signals can be expressed as [2] where denotes the vector/matrix transpose, , is the slow time, () is the fast time, is the period of the transmitted signal, and , , denotes the transmitted baseband coded signal of the th transmitter.

We assume that there are high speed targets with different Doppler frequencies located at the same initial range cell. The directions of targets with respect to the transmit array normal (i.e., DODs) are denoted by , respectively, and the directions of targets with respect to the receive array normal (i.e., DOAs) are denoted by , respectively. Thus, denotes the location of the th target. Because the range walk of target within the observation time is much smaller than the target range, DODs and DOAs of targets are assumed to be constant during the observation time. The radial velocities of the th target () with respect to the transmit array and the receive array are and , respectively, and then are the sum of two radial velocities for the th target. The radial accelerations of the th target with respect to the transmit array and the receive array are and , respectively, and then are the sum of two radial accelerations for the th target. The received baseband signal through reflections of targets can be written as [2, 3, 13, 14] where is the received signal vector; denotes the complex amplitude of the reflected signal of the th target involving the reflection coefficients and path losses of the target; and are the steering vectors of the receive array and transmit array, respectively, where denotes the wavelength; denotes the Doppler frequency of the th target, where is the radar carrier frequency; and the noise vector is assumed to be independent and identically distributed, and zero-mean complex white Gaussian distribution with , where is noise variance and is the identity matrix with the size of . The contribution of the target’s acceleration to the envelope migration during the observation time can be neglected and the effect of the target’s velocity on the envelope migration during the fast time can also be ignored. The reason is that the range change caused by the above two factors is much less than the radar range resolution. Under this condition, (2) can be simplified as By (3), in general cases, the envelope migration of high speed target echo caused by the high velocity of the target would exceed one range cell during the integration time, thus affecting the accumulation of target energy. There are three exponential terms in (3). The first and the second are the Doppler terms with respect to the fast time and the slow time induced by the target’s radial velocity. It is highly possible that the variety of the phase caused by the large-scale Doppler frequency of high speed target within the fast time is more than half of a turn of the unit circle; that is, . So the sharp distortion of target echo will occur, which results in the severe mismatch loss in matched filters, and then virtual array cannot be effectively formed. The third is the frequency modulation term because of the target’s radial acceleration ignoring the effect of the target’s radial acceleration on the phase variety during the fast time. Due to the range migration and the unsuccessfully formed virtual array, bistatic MIMO radar is incapable of effectively estimating DOD and DOA of high speed target.

3. Joint DOD and DOA Estimation for High Speed Target

Assume that the number of the range cells that the th high speed target has gone through at the slow time is denoted by , , which is an integer and can be written as where is the range resolution of bistatic radar, is the bandwidth of the transmitted signal, and means the ceil operator of integer. Ignoring the truncated effect of the integer operator on the envelope migration of echo, (3) can be rewritten as

Assume that during the observation time the total number of the across range cells of target is element of , where is an integer corresponding to the maximum possible number of the across range cells and minus denotes that the target moves towards the radar. Multiplying the received signal of the th receive element with the conjugate of the delayed version of the th transmitted signal , , one can construct a signal corresponding to the th searching across range cell as follows: where denotes the set for the variable , that is, when the equality should hold, , is the initial phase of the th target reflected signal when this target is located at the th range cell, and , where denotes complex conjugate. The first term in the right hand side of (6) holds under the condition with , which becomes a sinusoid signal over the fast time. The second term in the right hand side of (6) is given by Following that, a Fourier transform (FT) to the variable is employed to transform the signal into the fast-time frequency domain, which can be written as

Because of the high target speed and low repetition frequency of radar signal, it is highly possible that undersampling will occur [13, 14]. In this situation, the target’s true Doppler frequency can be expressed as [13, 14] where is the ambiguous Doppler frequency and is the fold factor. Substituting (9) into (8), we can obtain It is worth pointing out that in the fifth exponential term is integer times of because of and the integer . This term in (10) becomes . So (10) can be rewritten as

Assume that the reflected signal of the th target located at the th range cell is distributed along the repetition period index . Following that, a Fourier transform is taken to the signal with variable , and with the principle of stationary phase [22] we have where , , and , are the Fourier transform of and , respectively. Thus, Fourier transforms over both fast time and slow time have transformed the multiplied signals into the fast-time frequency and slow-time frequency domains. However, in the slow-time frequency domain, the target’s acceleration broadens the spectrum and the Doppler frequency ambiguity occurs.

Due to the target echoes existing in the initial range cell, that is, , the frequency-domain data   (, ) corresponding to the initial range cell can be exploited to estimate the Doppler frequencies and the ambiguous Doppler frequencies of targets. When and , , will generate the peaks. In order to improve the signal-to-noise ratio (SNR) for detection, the frequency-domain data with all transmitters and receivers of the bistatic MIMO radar is incoherently accumulated. Therefore, and can be estimated as Unfortunately, because the estimated is integer times of the repetition frequency of radar signal, the fold factor in (9) cannot be solved by the estimated and . Then the th target’s velocity can be obtained only by the estimated unambiguous Doppler frequency corresponding to the fast-time frequency domain, which can be expressed as

Using the very short echo time, the Fourier transform over fast-time domain has poor frequency resolution so that the Doppler frequency shifting value caused by target’s acceleration during the observation time is much less than half of the fast-time frequency resolution. Then the target’s energy can be approximately regarded as being concentrated in the same fast-time frequency cell during the observation time. However, due to the high frequency resolution of Fourier transform over slow-time domain, the target’s energy may be distributed along several slow-time frequency cells.

The frequency-domain data of the th target can be extracted from along the slow-time frequency dimension when , which can be written as where .

In fact, can be considered as the separated target component from the th transmitter to the th receiver via the th target located at the th range cell, ,  . The separated target components for the th target with all transmitters and receivers of bistatic MIMO radar can be expressed as where , , is the Kronecker product, and the vector represents a noise vector.

Assume that the th target moves away from radar. Analyzing the results of (13), we find that the indexes of the range cells that the th high speed target has gone through during the long observation time are , respectively. Then the separated target components distributed along all the across range cells for the th target are spliced as follows: where is the output of the formed virtual array with virtual elements at the slow-time index . Using (17), we can form the covariance matrix of for the th target from the spliced data of snapshots, which is given by where denotes the Hermitian transpose. It is worth pointing out that the estimation precision of covariance matrix after splicing the separated target components that are distributed along all the across range cells can be remarkably improved due to the average of the covariance matrixes corresponding to all the across range cells. Thus, the angle estimation precision can also be increased using the more accurate covariance matrix.

Assume that the th target moves towards radar. Analyzing the results of (13), we find that the indexes of the range cells that the th high speed target has gone through during the long observation time are , respectively. Then the separated target components for the th target distributed along all the across range cells are spliced as follows:

After synthesizing the virtual array with the range migration removed, most of the existing superresolution angle estimation algorithms can be applied to obtain DOD and DOA of the high speed target. In this paper, only the ESPRIT algorithm is used for angle estimation because of its simplicity and efficiency. Using (19), we can form the covariance matrix of for the th target from the spliced data of snapshots, which can be written as

Performing the eigendecomposition on the sample covariance matrix [5], we then have where is a scalar representing the largest eigenvalue because of only one existing target in , is a diagonal matrix constructed by the remaining eigenvalues, is the signal subspace composed of the eigenvector corresponding to the largest eigenvalue, and is the noise subspace containing the remaining eigenvectors of . Because the signal subspace spans the same space with the steering matrix, we have , where is a nonzero scalar.

Define , which is row equivalent to [7]. Suppose that is a signal subspace vector formed from by the same row interchange operations as is formed from . Let and be the first and the last rows of , and let and be the first and the last rows of .

Then let us define the average estimators as follows [7]: where and are the th elements in and , respectively, and and are the th elements in and , respectively. Therefore, the DOA and DOD for the th target, , can be written as The DOAs and DODs for the other targets can also be estimated by the same method. In summary, the DOAs and DODs of high speed targets for bistatic MIMO radar can be estimated via the following procedure.

Step 1. The signal corresponding to the th searching across range cell is constructed by multiplying the received signal of the th receive element with the conjugate of the delayed version of the th transmitted signal , .

Step 2. Following that, a Fourier transform (FT) to the variable is employed to transform the signal into the fast-time frequency domain, thus obtaining the signal .

Step 3. A Fourier transform is taken to the signal with variable , and we have the signal .
Fourier transforms over both fast time and slow time shown in Steps 2 and 3 have transformed the multiplied signals into the fast-time frequency and slow-time frequency domains.

Step 4. The frequency-domain data of the th target can be extracted from along the slow-time frequency dimension, which can be considered as the separated target component from the th transmitter to the th receiver via the th target located at the th range cell. Then the separated target components for the th target with all transmitters and receivers of bistatic MIMO radar are expressed as .

Step 5. The separated target components for the th target distributed along all the across range cells are spliced. Thus, the output of the formed virtual array for the th high speed target is expressed as .

Step 6. The DOD and DOA of th high speed target can be obtained from the output of virtual array using the superresolution algorithm, that is, ESPRIT algorithm. The DOAs and DODs for the other targets can also be estimated by repeating Steps 46.

4. Simulation Result

Some numerical examples are presented to illustrate the performance of the proposed method. A bistatic MIMO radar with transmitters and receivers is adopted. Radar carrier frequency , and the spacing between adjacent elements used for both transmit array and receive array is . The Gold sequences are selected as the transmitted baseband signals of bistatic MIMO radar, which are nearly orthogonal binary codes with good autocorrelation and intercorrelation. The duration of a code is assumed to be , the period length of Gold sequence is 1023, and then the repetition period of the transmitted signals is . Suppose that the number of the repetition periods for the transmitted signals during the observation time is 128.

4.1. Results of 2D Fourier Transform of the Multiplied Signals

Assume that three high speed targets located at the same initial range cell and their DODs and DOAs are , , and , respectively. The sums of two radial velocities with respect to the transmit array and the receive array for three targets are given by 7500 m/s, 9000 m/s, and 6500 m/s, respectively. The sums of two radial accelerations with respect to the transmit array and the receive array for three targets are 400 m/s2, 500 m/s2, and 450 m/s2, respectively. The SNRs of three targets are all −25 dB. The received signals are multiplied with the conjugate of the delayed versions of the transmitted signals, and then the Fourier transform (FT) of the multiplied signals over both fast time and slow time is employed.

Figures 2(a)2(f) show the results of 2D Fourier transform of the multiplied signals corresponding to the different searching across range cells, where the peak value indicates the existing target. By comparing the results of Figures 2(a)2(f), we can show that the target’s energy is concentrated in the same fast-time frequency cell when the target is located at the different range cells. The outputs are below noise level and unobservable as in Figures 2(a) and 2(f) with the searching range cell indexes . However, there are three or two peaks that are much higher than the noise background and may easily be detected as in Figures 2(b), 2(c), 2(d), and 2(e) with the searching range cell indexes . Thus, the energy of one target is distributed along the range cells with the indexes , while the other two go through the range cells with the indexes . The estimated target’s velocities are 7625.8 m/s, 9092.3 m/s, and 6452.6 m/s, which are obtained by the estimated unambiguous Doppler frequency corresponding to the fast-time frequency domain. The relative errors of velocity estimation are large because of poor frequency resolution in fast-time frequency dimension. However, the Doppler frequency ambiguity occurs in the slow-time frequency dimension.

4.2. DOD and DOA Estimation for High Speed Target

In this simulation, we evaluate the DOD and DOA estimation performance of the proposed method and the standard ESPRIT method [5] that is employed in bistatic MIMO radar for angle estimation. The simulation parameters for the three targets are the same as those in Section 4.1. Figure 3 shows the obtained result by using the standard ESPRIT method for all three high speed targets with 150 Monte Carlo tests. We can observe that the high speed targets are unsuccessfully localized by the standard ESPRIT method due to the invalidly formed virtual array and the range migration. Figure 4 shows the obtained result by using the proposed method for all three high speed targets over 150 Monte Carlo tests. It is clear that the DOAs and DODs of the three high speed targets are well estimated and automatically paired. The crosses denote the true locations of the three high speed targets in Figures 3 and 4.

4.3. Angle Estimation Performance for High Speed Target versus SNR

In this simulation, the SNRs of three targets identically range from −30 dB to 0 dB, while the other parameters for the three targets are the same as those in Section 4.1. The root mean square error (RMSE) of two-dimensional angle estimation is defined as where is the number of the independent trials and and are the estimates of true DOD and DOA in the th Monte Carlo trial. Figure 5 shows the RMSEs of DOD and DOA estimation versus SNR for target 1, which is located at () with the velocity and . The number of Monte Carlo tests is 200. It is shown that the standard ESPRIT method fails to estimate the angles of the high speed target, while the proposed method has a marked performance improvement over the standard ESPRIT method.

Because only the DOD and DOA estimation problem for high speed target is considered in this paper, the Cramer-Rao bounds (CRBs) for target’s motion parameters of different order [23] are not given. Moreover, the CRBs for target’s DOD and DOA parameters in bistatic MIMO radar are derived under the condition that the virtual array can be successfully formed by matched filters [2], which cannot be used as a proper benchmark for evaluating the performance of the proposed method. The proposed method is provided for overcoming the effect of the large-scale Doppler frequency on the virtual array formation and removing range migration during the long-time integration, thus facilitating the angle estimation for high speed target in bistatic MIMO radar. Therefore, it is a good way to make the angle estimation performance of the standard ESPRIT method when dealing with the low speed target without acceleration serve as a benchmark for performance comparison as shown in Figure 5. It is clearly shown that the standard ESPRIT method has a good performance in the estimation of DOD and DOA of the low speed target due to the effective virtual array formed by conventional matched filters and absence of range walk, while the proposed method exploited to estimate angles of the high speed target performs only a little worse than the standard ESPRIT method used for dealing with the low speed target. Moreover, for the high speed target, the angle estimation performances of the proposed method in the case of target acceleration and the proposed method using the target energy only in the initial range cell are also shown in Figure 5. These results illustrate that whether there exists target acceleration, the proposed method has a robust performance, and the proposed method using the target energy in all the across range cells attains better performance compared with that using the target energy only in the initial range cell. Thus, splicing the separated target components that are distributed along several range cells can improve the angle estimation accuracy.

5. Conclusions

A method for high speed target angle estimation using bistatic MIMO radar is proposed. For the target with high velocity, range migration always occurs and the serious mismatch induced by the target’s large-scale Doppler frequency usually exists in matched filter, which creates difficulties for target parameter estimation. After 2D Fourier transform of the multiplied signals which are obtained by multiplying the received signals with the conjugate of the delayed versions of the transmitted signals, the proposed method perfectly separates the target components of the radar return corresponding to the different transmitted waveforms and then splices the separated target components that are distributed along several range cells to successfully synthesize the virtual array. Thus, the DOD and DOA of high speed target can be estimated by using the superresolution algorithm. The effectiveness of the proposed method has been demonstrated by the numerical experimental results.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant no. 61302188, Grant no. 61372066, Grant no. 61071164, Grant no. 61071163, Grant no. 61271327, and Grant no. 41075115, in part by the Natural Science Foundation of Jiangsu Province of China under Grant BK20131005, and in part by the project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions. The authors wish to thank the anonymous reviewers for their valuable comments and suggestions which greatly improved the paper.