Abstract

This paper considers the long-time coherent detection problem for maneuvering targets with jerk motion. A novel method based on product-scaled integrated cubic phase function (PSICPF) is proposed. The main strategy of PSICPF is to estimate target’s motion parameters along the slow time for each range frequency cell. In order to eliminate the coupling terms between range frequency and slow time, the scaled nonuniform fast Fourier transform (SNUFFT) is newly defined in the integrated cubic phase function (ICPF). Then, the product operation is employed to coherently synthesize the estimation results, improve the antinoise performance, and suppress the cross terms. Finally, coherent integration is achieved via keystone transform (KT) and fold factor searching. Analysis demonstrates that the SNUFFT has the same computational complexity with nonuniform fast Fourier transform (NUFFT), and thus the PSICPF could be efficiently implemented via complex multiplications, the fast Fourier transform (FFT), and NUFFT. Detailed comparisons with other representative methods in computational cost, motion parameter estimation performance, and detection ability indicate that the PSICPF could achieve a good balance between the computational cost and detection ability. Simulations and real data processing results are presented to verify the effectiveness of the proposed method.

1. Introduction

The illegal flight activities of unmanned aerial vehicle (UAV) have caused great danger to air safety [1]. Therefore, radar maneuvering target detection attracts increasing attention in recent years. In order to detect such low radar cross section (RCS) targets, a long-time coherent integration is always required [28]. Unfortunately, the problems that come with it include not only the linear range migration (LRM) caused by the target’s high speed [9, 10] but also the range curvature and Doppler frequency migration (DFM) caused by the acceleration and jerk [1113]. These unfavorable effects seriously deteriorate detection performance of the conventional integration method, e.g., moving target detection (MTD). Consequently, how to effectively detect maneuvering targets is becoming a hot topic in the field of radar signal processing.

In order to coherently detect the targets, many successful methods have been proposed. In summary, these works could be divided into three categories according to the target’s maneuverability:(a)The first category mainly focuses on the LRM caused by the high speed. Representative methods include the keystone transform (KT) [1416], Radon–Fourier transform (RFT) [1719], scaled inverse Fourier transform (SCIFT) [20], frequency-domain deramp-keystone transform (FDDKT) [21], axis rotation MTD (AR-MTD) [10], and modified location rotation transform (MLRT) [22]. Nevertheless, due to ignoring the acceleration and jerk, the signal energy will be dispersed in the Doppler domain.(b)The second category further considers the uniform radial acceleration in the motion. Several typical algorithms could achieve satisfactory integration performance, such as the KT-Dechirp method [23], KT and Lv’s distribution (KT-LVD) [24], improved axis rotation and fractional Fourier transform (IAR-FRFT) [25], KT and linear canonical transform (KT-LCT) [26], MLRT-LVD [27], minimum entropy and Radon transform [28], Radon-LVD (RLVD) [29, 30], RFT-KT [31], frequency-domain second-order phase difference (FD-SoPD) [32], segmented KT and Doppler Lv’s transform (SKT-DLVT) [33], frequency autocorrelation function and and Lv’s distribution (FAF-LVD) [34], and RFT-FRFT [35]. However, for highly maneuvering targets, jerk motion usually exists and defocuses the Doppler spectrum.(c)The last category considers the jerk and thus has wider applications. The generalized RFT (GRFT) is proposed as an optimal estimator and detector by multidimensional parameter searching [36]. Despite its perfect detection performance, the serious blind speed side lobes (BSSLs) and huge computational burden prevent its applications. The subband dual-frequency conjugate (SDFC) and Radon-chirp rate-quadratic chirp rate (RCR-QCR), i.e., SRQ method proposed in [37], avoid the BSSLs of GRFT and could estimate motion parameters at low SNR. However, its computational complexity is heavier than GRFT and thus is not suitable for real-time processing. The adjacent cross-correlation function (ACCF) is proposed in [38, 39], which corrects the range migration and DFM by the rank reduction operation. Although the ACCF-based methods have low computational cost, their poor antinoise performance becomes a shortcoming that cannot be ignored. The method based on time reversing transform, second-order KT, and LVD (TRT-SKT-LVD) is proposed in [40]. Similar to the ACCF, the TRT operation obtains high efficiency at the cost of sacrificing the detection performance. In [41], two methods, i.e., KT and generalized dechirp process (KT-GDP) and KT and cubic phase function (KT-CPF), are presented to remove DFM and realize coherent integration. However, the incoherent integration in fold factor searching is sensitive to noise, which limits the detection performance.

This paper extends our preliminary works in [32, 34]. We propose a method based on product-scaled integrated cubic phase function (PSICPF) to achieve coherent detection for maneuvering targets with jerk motion. Different from conventional methods, the main strategy of PSICPF is to estimate motion parameters in the range frequency and slow time domain. The scaling property and product operation are introduced to realize decoupling and energy accumulation. Furthermore, the proposed method could be efficiently implemented via complex multiplications, fast Fourier transform (FFT), and nonuniform FFT (NUFFT) [42]. Detailed comparisons with other representative algorithms indicate that the proposed method could achieve a balance between detection performance and computational cost.

The remainder of this paper is organized as follows: in Section 2, the signal model for a maneuvering target with jerk motion is established. The proposed method is illustrated in detail in Section 3. In Section 4, the properties are analyzed. Experiments with the synthetic and raw data are performed in Section 5 to validate the effectiveness of the proposed method. Finally, Section 6 gives the conclusions.

2. Signal Model and Problem Formulation

Suppose the radar transmits a linear frequency-modulated (LFM) signal, which could be expressed aswhereis the rectangular window function and is the fast time. and indicate, respectively, the pulse width and carrier frequency. is the chirp rate, with the bandwidth of radar signal.

In most applications, a three-order polynomial could be used to approximate the instantaneous slant range between the maneuvering target and radar; i.e.,where , , , and denote, respectively, the initial slant range, radial velocity, radial acceleration, and jerk of the target. denotes the slow time variable, represents the pulse repetition time, and is the number of integrated pulses.

Ignoring the influence of noise, the received signal after pulse compression can be written as [34]where and denote the signal amplitude and the light speed, respectively.

Inserting (3) into (4) yieldswhere is the wavelength.

As can be seen in equation (5), the first term indicates the signal envelope, which changes nonlinearly with the slow time. When the offset during the observation time exceeds the range resolution , range migration would occur. The exponential term indicates the Doppler frequency of the target, which is defined as [34]

Due to the acceleration and jerk, the Doppler frequency changes nonlinearly with the slow time, which forms a quadratic frequency-modulated (QFM) signal. When the offset exceeds a Doppler resolution, the DFM would occur and seriously defocus the target energy in the Doppler domain.

Performing the Fourier transform (FT) on equation (5) along the -axis, we obtain the compressed signal in the slow time-range frequency domain; i.e.,where is the amplitude of the spectrum and denotes the range frequency variable.

It is obvious in equation (7) that the coupling between and is the essential cause of range migration. In order to detect the maneuvering target in low SNR background, the coupling terms and DFM should be eliminated correctly by an effective detection method.

3. Principle of the Proposed Method

In conventional methods, the range migration is firstly corrected in the slow-range time domain, and then the DFM is estimated and eliminated. However, in our proposed method, the PSICPF is derived in range frequency-slow time domain to estimate the acceleration and jerk. Then, the KT and fold factor searching are employed to achieve coherent integration.

3.1. Parameter Estimation via PSICPF

In equation (7), the signal along the slow time could be expressed aswhere is the scaling factor dependent on the range frequency and , , , and denote the constant phase, centroid frequency, chirp rate, and quadratic chirp rate, respectively.

The cubic phase function (CPF) [43] is defined as a bilinear autocorrelation operation over the slow time; i.e.,where is the lag variable.

In [44], the NUFFT is performed with respect to to accumulate the signal energy; i.e.,where is the frequency variable corresponding to and is the Dirac delta function.

In equation (10), the signal energy peaks along the inclined line . However, because of the coupling between and , the interception and slope of this line indicate the wrong values. Therefore, we consider the scaling factor and propose the scaled NUFFT (SNUFFT), which could be written aswhere is the scaled frequency variable corresponding to

Different from NUFFT, the peak position is corrected to the right place after SNUFFT. In this case, the parameter estimation result will keep consistent for all range frequency cells, which benefits to the coherent integration of signal energy. At the same time, the SNUFFT could also be efficiently implemented via NUFFT without increasing the computational burden.

In equation (9), the slow time and lag variable are coupled with each other, which is the essential cause of the inclined line. Different from the LVD [30], the first exponential term in equation (11) makes difficulties in the decoupling process. In the integrated CPF (ICPF) [45], the complex modulus and inverse FT (IFT) are performed to eliminate the disturbance of ; i.e.,where is the time variable with respect to .

Finally, the scaled Fourier transform (SFT) [46] and FT are employed to coherently integrate the signal energy in the chirp rate-quadratic chirp rate (CR-QCR) domain:

As can be seen from equation (13), the signal is accumulated into a peak, of which the position is . For different range frequency cells, the peak locates at the same position. In order to improve the antinoise performance of parameter estimation and suppress the side lobes, the product operation could be employed over the range frequency. In this case, the PSICPF is defined, which coherently integrates the signal energy aswhere is the amplitude after product.

By equation (14), we could simultaneously estimate the chirp rate and quadratic chirp rate from the peak position; i.e.,

Finally, the acceleration and jerk of the target are derived as

3.2. Coherent Integration via KT

Referring to equation (7), we may utilize the estimated parameters to construct the following phase compensation function to compensate the range curvature and DFM caused by acceleration and jerk:

Multiplying equation (17) by equation (7), we have

The expression in equation (18) shows that the effects of acceleration and jerk are removed and only the linear coupling term caused by the target velocity exists. However, due to the high speed of target and the low-pulse repetition frequency (PRF) of radar, Doppler ambiguity would often occur. In this case, the velocity of the target can be stated aswhere is the ambiguity integer named the fold factor, is the blind speed with the PRF, denotes the ambiguous velocity which satisfies , and is the remainder function.

Substituting equation (19) into equation (18), we obtainwhere the fact is used.

To achieve coherent integration, many successful methods have been proposed [1416]. In this paper, the KT is employed for convenience, which scales the slow time for each range frequency:

Substituting the scaling formula into (20), we get

By equation (22), the coupling caused by the ambiguous velocity is removed. However, the residual range migration caused by integral multiples of blind speed, i.e., the coupling between and , is still present. Thus, another phase compensation function is established to search the unknown fold factor and then to achieve coherent integration; i.e.,where is the searching fold factor.

Multiplying (23) with (22) and performing IFT along the range frequency, we may obtain

When , the range migration will be eliminated and the FT could be performed with respect to the slow time to achieve coherent integration; i.e.,

The searching fold factor could be estimated by comparing the amplitude of :

Finally, the target velocity is estimated as

In (25), the received signal of a maneuvering target is integrated into a single peak with the peak value . Thereafter, the constant false alarm rate (CFAR) technique could be used for target detection; i.e.,where is the detection threshold. The adaptive threshold is obtained by the reference unit after coherent integration via the integration method, and then the test statistic in the detecting unit is compared with the threshold to confirm the presence or absence of a target. If the test statistic is smaller than the threshold, there will be no moving target or a target is missed, and then it moves on to the next detecting unit. Meanwhile, if the test statistic is larger than the threshold, the target detection is declared.

4. Analysis of the Proposed Method

4.1. Cross Terms Analysis

Similar to (7), the received signal of targets in range frequency domain could be written as

The cubic phase function of (29) iswhereare the self-terms, andare the cross terms.

According to Euler’s formula, i.e.,we obtain

The cross terms for multiple targets at a certain range frequency cell could be obtained as

It can be seen from (34) and (35) that the cosine function and the nonlinear coupling terms would disturb the integration of the cross terms after SNUFFT, complex modulus, IFT, and SFT. Only when , , and , the cosine function could be eliminated. In this case, the cross terms become the self-terms. Thus, the cross terms of PSICPF cannot be accumulated.

4.2. Computational Complexity Analysis

In this part, the computational burden of the proposed method is analyzed in detail. The GRFT, KT-GDP, KT-CPF, and TRT-SKT-LVD, which could also complete coherent detection for maneuvering target, are selected for comparisons. Denote the number of range cells, echo pulses, searching velocity, searching acceleration, searching jerk, and searching fold factor as , , , , , and , respectively.

For GRFT, a four-dimensional parameter searching is needed. The computational complexity can be easily obtained as [36].

For KT-GDP, the KT and fold factor searching are firstly performed to eliminate the LRM, and the computational costs are, respectively, [16] and . The GDP is then employed to estimate the acceleration and jerk, of which the computational cost is at the order of . Therefore, the total computational complexity is about [39].

Compared with KT-GDP, the acceleration and jerk are estimated via CPF, of which the computational cost is about [44]. In this case, the computational complexity of KT-CPF is about .

For TRT-SKT-LVD, to complete the chirp-z-based SKT and LVD, the computational complexities are and , respectively. Therefore, the total computational complexity is about [40].

For the proposed method, the main procedures in estimating acceleration and jerk for each range frequency bin include the bilinear autocorrelation operation , the SNUFFT , the IFT along the scaled frequency , the SFT with respect to the slow time , and the FT to integrate the signal energy . The computational complexity of KT and fold factor are, respectively, and . Thus, the total computational complexity is in the order of

The computational complexities of the above methods are listed in Table 1. In practical applications, , , and are usually much larger than , while is much smaller than . Under the assumption of , , and , Figure 1 intuitively presents the computational complexities of the above methods. It is obvious that the GRFT and the KT-GDP take too much time due to the multidimensional parameter searching. The proposed method has a moderate computation complexity compared with KT-CPF and TRT-SKT-LVD. However, these two methods have poor detection performance under low SNR environment, which will be shown in Section 5.

5. Simulations and Real Data Processing Results

In this section, several numerical experiments are given to demonstrate the effectiveness of the proposed method. The simulation parameters are given in Table 2.

5.1. Coherent Integration for a Maneuvering Target

First of all, a single maneuvering target is used to evaluate the coherent integration ability of the proposed method. The motion parameters of the target are , , , and . Figure 2(a) gives the target trajectory. Due to the high speed of the target, serious range migration could be observed. Figure 2(b) shows the parameter estimation result when . Thanks to the proposed SNUFFT, correct values of acceleration and jerk could be obtained. Then, the PSICPF is achieved by multiplying all the range frequency cells, as presented in Figure 2(c). It is clear that the product operations greatly strengthen the main lobe and suppress the side lobes. From the peak position, we could read that and . After compensating for the range migration and DFM caused by acceleration and jerk, the KT and fold factor searching are performed. From Figure 2(d), we could estimate the fold factor, i.e., . The LRM elimination result is given in Figure 2(e). Finally, coherent integration is realized and a well-focused peak could be seen in the range-Doppler domain, as shown in Figure 2(f). Considering the peak position, we may estimate the target velocity and initial range, i.e., and . This experiment preliminarily demonstrates the integration capability of the proposed method.

5.2. Coherent Integration for Multiple Targets

The integration ability of the proposed method for multiple targets (i.e., Tr1 and Tr2) is evaluated in this part, where the motion parameters are listed in Table 3. In this simulation, Gaussian noise with zero mean is added to the signal after pulse compression. The input SNR is defined aswhere is the variance of Gaussian noise. Results of integration and parameter estimation are shown in Figure 3.

Figure 3(a) gives the result of range compression. Since the input SNR is 0 dB, the target trajectory could be hardly recognized. Figure 3(b) shows the parameter estimation result of PSICPF. As can be seen, two sharp peaks are obtained in the CR-QCR domain. By product operation, the self-terms are strengthened while the cross terms are greatly suppressed. From the peak position, we get and for Tr1 and and for Tr2. After phase compensation with the estimated parameters, the fold factor estimation results are presented in Figures 3(c) and 3(d). For Tr1, we have , and for Tr2, we have . Finally, the coherent integration results are, respectively, given in Figures 3(e) and 3(f). The velocities are estimated as for Tr1 and .

In Figure 3, the proposed method could achieve coherent integration for multiple targets. In order to intuitively evaluate the impact the cross terms, we could define the normalized cross terms aswhere denotes the self-terms.

The simulation results of cross terms for a certain range frequency cell are shown in Figure 4, where the Gaussian noise is ignored. It is evident that the cross terms show significant oscillation characteristics. The maximum ratio of is about 6%. After the product operation of different range frequency cells, the ratio will further decrease. Therefore, we may get the conclusion that the cross terms have little impact on parameter estimation.

5.3. Performance of Target Detection and Parameter Estimation

The detection performance of the proposed method is also evaluated by using the CFAR detector. Complex white Gaussian noise is added to the signal after pulse compression to yield SNRs varying from with −20 dB to 20 dB. 500 independent Monte Carlo trials are performed for each SNR value. The false alarm ratio is set as . The representative GRFT, KT-GDP, KT-CPF, and TRT-SKT-LVD are selected for comparisons. The detection performance is given in Figure 5.

First of all, it is important to state that the GRFT has the optimal detection performance since it is designed based on maximum likelihood. Obviously, with the GRFT as a benchmark, the KT-GDP suffers from 6 dB performance loss due to the incoherent procedure in fold factor searching. The KT-CPF has inferior detection probability at −5 dB than the proposed method since the least squares (LS) method fails at low SNR. The TRT-SKT-LVD is only effective when the SNR is high because the TRT operation loses much signal energy, even though it is beneficial to reducing computational complexity. Therefore, when considering the computational cost, the proposed method is more suitable for practical applications than GRFT.

As the detection performance is determined by the target’s acceleration and jerk estimation accuracy, the parameter estimation performance is also evaluated by the experiment. The root-mean-square error (RMSE) is utilized as the criterion. The RMSEs of initial range, velocity, acceleration, and jerk are given in Figures 6(a)6(d), respectively. Corresponding to the detection performance, the GRFT and KT-GDP have better antinoise performance than the proposed method, and the KT-CPF may get wrong estimation results at low SNR. The TRT-SKT-LVD could not estimate motion parameters and thus is not compared with other methods.

5.4. Real Data Processing Results

In this subsection, we adopt the measured data of a DJI Phantom 3 commercial UAV to demonstrate the proposed method. The data were collected in March 2017 by the National University of Defense Technology, Hunan, China. Figures 7(a) and 7(b) present the experimental scene and frequency-modulated continuous wave (FMCW) radar system, respectively. Radar parameters are given in Table 4. It is worth noting that, in order to obtain the LRM and Doppler ambiguity, we deliberately chose a large bandwidth and a low PRF before experiment.

Figure 7(c) shows the target trajectory after pulse compression. During the coherent time, the UAV moves across 24 range cells and causes serious range migration. Figure 7(d) gives the parameter estimation result at . From this range frequency, we may obtain the target acceleration and jerk. However, to improve the antinoise performance, the PSICPF is proposed by multiplying all the range frequency cell. Figure 7(e) presents the result of PSICPF, where the noise and clutter are greatly suppressed. The motion parameters of the target are and . Then, the KT and fold factor searching are employed to remove the LRM and achieve coherent integration, of which the results are shown in Figures 7(f) and 7(g), respectively. It is easy to see that a well-focused peak is obtained in the range-Doppler domain. The initial range and velocity of the UAV are and . At the same time, Figure 7(h) gives the integration result of MTD. Due to ignoring the LRM and DFM effects, the signal energy is dispersed over multiple ranges and Doppler cells, which creates difficulties for target detection. This experiment verifies the practicality of the proposed method.

6. Conclusions and Future Work

In this paper, we have presented a coherent integration method for maneuvering target with jerk motion based on PSICPF. The main strategy of the proposed method is to coherently estimate motion parameters along the slow time for each range frequency. In order to eliminate the coupling between range frequency and slow time, the scaling property is introduced into NUFFT and then the SNUFFT is proposed. When integrating the signal energy, the product operation is applied to improve the antinoise performance and suppress the side lobes. Finally, the KT combining with fold factor searching is performed to achieve coherent integration. The contributions of the proposed method include the following:(1)The SNUFFT is newly defined to realize decoupling between range frequency and slow time so that the parameter estimation results of different range frequency cells could remain coincident(2)The product operation is employed to coherently synthesize the estimation results, improve the antinoise performance, and suppress the cross terms(3)The proposed method could be implemented via complex multiplications, FFT, and NUFFT, which guarantee the high efficiency(4)The PSICPF achieves a good balance between detection performance and computational complexity

As the bilinear autocorrelation operation in CPF brings about estimation performance loss, linear methods will be considered in the future to further improve the detection ability. The core problem may focus on the efficient implementation.

Data Availability

The radar data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This article was supported by the National Natural Science Foundation of China (no. 61501513).