A recognition algorithm of VGPO jamming based on discrete chirp-Fourier transform

This paper addresses the recognition problem of velocity gate pull-off (VGPO) jamming from the target echo signal for the velocity-based tracking system. The discrete chirp-Fourier transform (DCFT) is studied in this paper to jointly estimate the chirp rates and frequencies of the target and jamming signals. Firstly, the scaling characteristic of the DCFT algorithm is explored. Then, we focus on the quantitative effect of the VGPO jamming signal by analyzing the jointly estimated results in each pulse. The quantitative effect indicates that, as long as the estimated frequency is unchanged, the relationship between the estimated chirp rate and the pulse numbers is similar to the relationship between the frequency offset of VGPO jamming and the time. Finally, by utilizing the joint estimated results in each pulse repetition interval and calculating the mean square to variance ratio (MSVR) of the normalized estimated chirp rate, the VGPO jamming can be recognized. Simulation results show that the jamming signal and the target echo become distinguishable with the proposed feature. Comparing to the existing works, the proposed method can correctly recognize the jamming signal with lower jamming-to-noise ratio (JNR) 5 dB with less data needed, which means it can work effectively in the early stage of interference implementation and shows great potential in practical applications.


Introduction
Electronic countermeasure (ECM) has always been an important research direction in radar [1,2]. As a typical category, the deceptive jamming is often used to deceive a hostile radar system by generating a large number of active false targets in terms of direction, position, and velocity [3]. With the widespread application of the digital radio frequency memory (DRFM) technique, the interference effect of ECM is further reinforced. In a DRFM system, the intercepted radar signal is firstly down-shifted in frequency. And then it is sampled with an analog-to-digital converter and stored in a high-speed digital memory where the samples can be manipulated in amplitude, frequency, and phase to produce (2020) 2020:37 Page 3 of 18 Currently, many methods are available on joint estimation, such as maximumlikelihood [20], adaptive chirplet transform [21], and phase function [22]. However, the theoretical value of the rate of frequency of the jamming signal is far less than the minimum interval in the chirp rate domain in a general way. Consequently, most of these methods will lose efficacy. For the discrete chirp-Fourier transform (DCFT) algorithm [23], it is generalized from the discrete Fourier transform (DFT). Besides matching the multiple frequencies similar to the DFT algorithm, it uses the multiple chirp components to match the chirp rates. Most important of all, it can increase the resolution in the chirp rate domain by introducing a scaling factor, which is termed as the scaling property and will be analyzed in this paper. As a result, it is still applicable in this situation. Nowadays, it has applied in many cases [24,25].
This paper is organized as follows. In Section 2, the received signal models of both real target echo and VGPO jamming are first formulated; then, we introduce the modified DCFT and analyze its scaling property, derive the estimated chirp rate and frequency in each pulse repetition interval (PRI), and further define a new feature parameter. Section 3 shows the simulation results. Section 4 is the concluding remarks.

Methods
In this section, we elaborate on the theory and the method of jamming signal recognition. First, we present the signal models for VGPO jamming and target echo. Second, we introduce the modified DCFT and study its scaling property. Then, the quantitative effect is analyzed when performing the joint estimation with the DCFT algorithm, and it is also discussed in the non-linear scenario. Finally, we define a feature to distinguish the VGPO jamming and the target echo.

Signal model
The whole process that a jammer implements VGPO interference can be divided into three steps: pull-off halting, pulling off, and closing [3]. After the radar signal is intercepted, the jammer retransmits the signal with a Doppler frequency shift inflicted on. To deceive the radar effectively, the power of the jamming signal is always greater than that of the real target echo. Therefore, the tracking loop of the radar would track the false Doppler frequency and be slowly pulled away from the velocity gate of the real target.
Assume the signal transmitted by radar is s(t). Then, the signal reflected by a hypothetical target in the far-field can be modeled as [16] where A r and t r denote the reflection coefficient and the time delay of the hypothetical target, respectively, and f d is its Doppler frequency. Comparing with the real target echo, the VGPO jamming has a pseudo-Doppler frequency shift based on f d . For convenience, we consider the scenario that only one interference exits, and the jamming signal can be expressed as where A J denotes the equivalent reflection coefficient of the interfering source, t J = t r + t d is the jamming propagation delay, and t d is the time delay of jammer relay which is often far less than t r , f dj = f d + f d is the modulated pseudo-Doppler frequency. Here, we define jamming-to-noise ratio (JNR) A 2 J /σ 2 , and signal-to-noise ratio (SNR) A 2 r /σ 2 , where σ 2 is the variance of Gaussian noise n(t), and jamming-to-signal ratio (JSR) is defined as JNR-SNR.
Similar to the RGPO system described in [5,26], a VGPO system should linearly increase the frequency offset to create a series of constant velocity-rate false targets. This scenario is also the one we often use when analyzing VGPO [19,27,28] and is equivalent to the constant acceleration profile in the range domain. Then, the f dj at the mth PRI can be expressed as [28]: where V f denotes the slope of f d with respect to t which is also called acceleration and T r denotes the pulse repetition period.
To enable radar to accurately capture and track the velocity loop of the jamming signal, V f would not be greater than the maximum trackable acceleration of a target to radar [1]. On the other side, V f would not be particularly small for implementing interference effectively and protecting the real target. In a word, V f has a reasonable value range.
Substitute (3) into (2), and we can find that, despite the signal form, the VGPO jamming signal is a chirp function with a chirp rate of V f . This is the inspiration for the jamming identification in this paper.
As the specific form of s(t) is not assumed, we may have to do some preprocessing to prevent the influence of the potential time-varying frequency components in signal form on chirp rate estimation. The preprocessing may be different for various kinds of signals. Take the widely used linear frequency modulation (LFM) signal as an example, it has a quadratic term in terms of t, i.e., the inherent chirp rate. Generally speaking, the acceleration V f of the VGPO jamming (three times the acceleration of gravity in Xband is about 2 kHz/s) is much smaller than the chirp rate of LFM signal (typical value 10 10 − 10 14 Hz/s, [1]). As a result, the de-chirping processing should be performed first in the recognition procedure for the radar that emits LFM signal.
After being sampled with discrete-time t = nT s where T s denotes the sampling interval, and performing some appropriate preprocessing, the jamming signal s J (t) in mth PRI can be expressed as where where f d0 and ϕ 0 separately denote the time-invariant frequency and constant phase term in s J (n).
The phase of the discretized jamming signal in (4) contains three components: a quadratic item with respect to nT s with a coefficient V f , a linear item with respect to nT s with a coefficient f dj , and a time-independent item ϕ j . Ignoring the time-independent phase item, the coefficient f dj of the linear item is different in various PRIs, while the coefficient of the quadratic phase items is the same. However, the following analysis indicates that the quantization in the frequency domain will influence the estimated result in the chirp rate domain when performing joint estimation to the jamming signal.

Discrete chirp-Fourier transform and its scaling property
The discrete chirp-Fourier transform algorithm is used to jointly estimate the chirp rate and the frequency of a chirp-type signal. For an analog chirp signal , its discretized form with sample rate 1/N is where W N = exp(−j2π/N). The corresponding N-point DCFT of signal x I (n) is defined as [29] where k and l represent the frequency and chirp rate variables, respectively. In (7), the duration of the signal is normalized as 1, which makes both of k and l less equal than the length of signal N. In radar signal processing, the Doppler frequency of the signal is less than the sampling frequency of f s , and the chirp rate is less than f 2 s correspondingly. After sampling with discrete-time t = nT s , the discretized form of x(t) is rewritten as To maintain the consistency between the estimated values of frequency and chirp rate and the actual ones, the corresponding form of DCFT should be modified as Each bin in frequency domain and chirp rate domain represents f s /N and f 2 s /N, respectively. When |X mc (k, l)| maximizes, the estimated values of k and l are respectively, and the corresponding estimated values of frequency and chirp rate are k 0 and l 0 , respectively. Essentially, the time delay of a signal is equivalent to zero padding in front of the signal. It changes the phase of the signal's DFT result, while the estimated frequency and the envelope of DFT results are unchanged. In a very short timestamp, the positions of target and jammer are usually considered to be located in the same range cell. Under this assumption, both of the time delay t r and t J are constant, which means the number of zero paddings in each PRI is the same during this period. In other words, the time delay in these PRIs would not influence the estimated value of frequency through DFT. As the DCFT algorithm is a generalized form of the DFT algorithm, this conclusion is applicable for jointly estimating the values of chirp rate and frequency with the DCFT algorithm.
From (10), one can see that the chirp rate variable only appears in the factor W ln 2 N . It means that the chirp rate domain of DCFT can be broadened or compressed after dividing the variable l by a scaling factor C, which equivalently increases or decreases the minimum interval in the chirp rate axis. According to the scalability of the DCFT algorithm in the chirp rate domain, we further define the modified DCFT as follows: At this time, the maximum value in the chirp rate axis changes from f 2 s to f 2 s /C, and the corresponding minimum interval decreases to f 2 s /NC. Similarly, for the chirp signal in (9), the peak position in the modified DCFT domain is located at and the corresponding values of frequency and chirp rate are k 0 and l 0 , respectively. It indicates that the estimated values of chirp rate and frequency are unchanged in the modified DCFT compared with the DCFT in (10). As mentioned above, the value of V f is in a reasonable range. What is more, the sampling frequency f s is always known as prior information. As a result, the scaling factor C has a rough estimation for the radar side. The bigger the scaling factor C, the smaller the minimum interval in the chirp rate domain for the modified DCFT algorithm. In fact, if C is big enough or much bigger than N, the DCFT result for different l is almost the same. In this case, the DCFT algorithm is considered as degenerating into the DFT algorithm and could cause a huge estimation error to the chirp rate at low SNR.

The analysis of quantitative effect
In this part, we analyze the influence of a small frequency offset on the DCFT estimation results. In the DFT algorithm, we estimate the constant frequency f n0 of a signal by matching it with a quantitative value in frequency domain, i.e., f n0 = n 0 f s /N, where N is the total number of DFT and n 0 (0 ≤ n 0 < N − 1) is a certain nonnegative integer. When the signal's frequency f n1 has an offset f based on f n0 , i.e., the peak still exhibits at the position f = f n0 only if | f | < f s /2N. It means the frequency offset f is ignored due to the quantization in the frequency domain. Similarly, the quantization error must occur in the DCFT algorithm as well. From [23], it is pointed out that, for a discrete chirp signal its DCFT result has the similar peak property as the signal wherek 0 andl 0 are close to integers l 0 and k 0 , respectively, i.e.
where and η are two arbitrary positive values. However, this conclusion is drawn based on the condition thatl 0 andk 0 are close enough to l 0 and k 0 , i.e., , η ≈ 0. As the process of pulling off a velocity gate is continuous, the pseudo-Doppler frequency of the jamming signal changes in each PRI. Moreover, the minimum interval in the chirp rate axis decreases after adopting the modified DCFT in (12). As a result, it becomes very necessary to analyze the quantitative effect of the DCFT algorithm in this paper. For a discrete VGPO jamming in (4), its modified DCFT after ignoring the timeindependent phase item ϕ j can be calculated by Let the Doppler frequency wheref dj is an integer multiple of f s /NC, and f m denotes the remainder Doppler frequency besidesf dj in mth PRI. Here, we assume the total number of pulses is M when the estimated frequency through the modified DCFT is fixed asf dj . Then, we have (5) and (18). Then, (17) can be rewritten as follows: where In the complex coordinate system, each complex exponential can be regarded as a unit vector. Consequently, P can be considered as the modulus of a vector summed by N unit vectors, and it reaches to maximum of N when all of the vectors point in the same direction. According to [23], P maximizes when the frequency variation of DCFT matches to the quantized frequencyf dj , i.e., φ 2 = 0. Then, we have On the other side, φ 1 (n) differs by an integer to various n, where n 1 and n 2 are two arbitrary different integers between 0 and N − 1 and ϕ 0n is an integer. Substituting the expression of φ 1 into (22), we can yield Let α 1 = ϕ 0n /(n 1 − n 2 ), and α 2 = n 1 + n 2 . Then, if P reaches to maximum N, it implies that the variable l to arbitrary n 1 and n 2 in (23) is solvable, i.e., The effective value of φ 1 only exists in the fractional part from 0 to 1 due to the periodicity of the exponential function. Correspondingly, the true value of l is the remainder after dividing by N in (24). In other words, for a fixed m, the true value of l should be fixed, too. From (24), one can see that l and m are linearly related. In general, not all vectors in (20) point in the same direction, and P cannot reach to maximum N. ϕ 0n in (22) is not an integer anymore to a fixed pair of integers n 1 and n 2 , which changes the estimated chirp rate of jamming signal in mth PRI. Among various PRIs, m is the unique determining factor that influences the value P, and the estimated chirp rate l. From this point, ϕ 0n is not an integer either for various pairs of n 1 and n 2 . However, it can be found from the formulation of φ 1 that ϕ 0n is linearly related to m for fixed n 1 and n 2 . It means l is still proportional to m. As the influence of m to φ 1 (n) and l is small. Once V f is fixed, all integer pairs n 1 and n 2 that make P reach to the maximum value (not N) are basically determined. Therefore, the slope between l and m can be considered to be unchanged at this time.
More generally, f includes an extra frequency independent of m (part off dj ) which results in ϕ 0n is not an integer either. Then, the estimated chirp rate of VGPO jamming in mth PRI may differ from the one in (24), while the linear relationship between l and m still exists.
It can be seen that, when a joint estimation of the VGPO jamming with the DCFT algorithm is performed, the quantization in frequency domain will affect the estimated value of the chirp rate, and the estimated chirp rate is linearly related to the residue of the quantized frequency item.
Unlike in (3), we next discuss a more general scenario that the false target is maneuvering, which means f dj has a non-linear correlation with time t. Taking the quadratic relation as an example, f dj changes as follows over time t: where A f is the cubic coefficient of time t. Substituting the f dj in (2) with f dqj , and following the same derivation from (4)∼6 and (17)∼(19), we can calculate the modified DCFT of the jamming signal in non-linear scenario. Similar to (5), the frequency item of the preprocessed jamming signal f dqj can be expressed as follows: Still, we let φ 1 (n) be the function that includes high order terms of n and f m , i.e., Obviously, besides a cubic item of nT s , there is a quadratic item of n with respect to m, and a linear item of n with respect to m 2 in φ 1 (n). With the same condition as in (22), we can find that there is a quadratic relation between l and m in this situation. More generally, the relationship between l and m is the same with that between f dj and t. Here, we term this phenomenon as a quantitative effect. Therefore, it can be concluded that a series of small offsets in frequency would increase the values of the chirp rate in joint estimation. The quantitative effect can be interpreted as, within the timestamp when the estimated quantized frequency is unchanged, the estimated chirp rate of the VGPO jamming through the modified DCFT algorithm increases with a tendency that is the same as the relationship between f dj and t.

Feature extraction
In this part, we use an identification parameter to benchmark the quantitative effect. As the pull-off stage of VGPO jamming continues for several seconds [1], the quantitative effect will appear periodically. For the real target echo, the estimated chirp rate through the modified DCFT algorithm is essentially the estimation error caused by thermal noise. Its value equals 0 in an ideal case and randomly changes under noisy environment. By utilizing the distinctive quantitative effect, the VGPO jamming can be distinguished from the real target echo.
Here, we define the feature, the mean square to variance ratio (MSVR) of the normalized estimated chirp rate, and utilize it to distinguish the VGPO jamming and the real target echo. Let C e (m) represent the estimated chirp rate in mth PRI, and its normalized form is yield by where C max e and C min e represent the maximum and minimum of C e , respectively. Let C e =[ C e (1) C e (2) · · · C e (m)], and the MSVR is defined as where μ c and σ c stand for the mean and the variance of C e , respectively. From the aforementioned content, it is concluded that C e (m) of VGPO jamming linearly (or non-linearly) increases in M PRIs. Although the estimated chirp rate is normalized from 0 to 1, its variance becomes smaller and smaller as the JNR increases. On the contrary, all of the elements in C e of real target echo can be regarded as random, which means its MSVR is steady. Consequently, the VGPO jamming signal and the real target echo can be recognized by setting a proper threshold.
The flowchart of the identification procedure to VGPO jamming based on the DCFT algorithm is provided in Fig. 1. After being sampled, the received signal is preprocessed. Then, the chirp rate of the signal is estimated by the modified DCFT in each pulse. Finally, by calculating the value of MSVR with the estimated results in M PRIs, the VGPO jamming can be identified from the target signal.

Results and discussion
In this section, we conduct some simulation experiments to illustrate the performance of the proposed method in this article. We assume the waveform of the transmitted signal is LFM. It is worth noting that the pulse numbers are set relatively large in some simulations which are similar to [18]. This is because our method is performed from the time when the interference begins to implement, and the shifted frequency is small at the early pull-off stage and even within a single frequency interval.

Scaling property
Firstly, we demonstrate the effectiveness of the scaling property of the modified DCFT algorithm. With the acceleration V f 50 kHz/s for instance, we consider the scenario that f dj linearly changes with t as shown in (3). The carrier frequency is 10 GHz. The bandwidth and time width of the signal is 5 MHz and 100 μs, respectively, and the pulse repetition period is 500 μs. The velocity of the target is fixed as 200 m/s. The SNR and JNR are 5 dB and 10 dB, respectively, and the sampling frequency is 10 MHz.
Both of the DCFT results of jamming signal and target echo and their corresponding profiles at the frequency where the |X mc | 2 maximizes are shown in Fig. 2. In the first row, the scaling factor C is 1 (without scaling). In the second and the third rows, C is 4e3 and 2e7, respectively. In the profiles, both of the modules of the DCFT results of jamming and target are normalized to 1 to compare the estimated chirp rate more visually, and the peaks of DCFT results of jamming and target are respectively signed by a red circle and a blue star.
By comparing the results shown in various rows in Fig. 2, it is easy to find that the DCFT result magnifies in the chirp rate domain as the scaling factor C increases. The estimated chirp rates of target and jamming are distinguishable in Fig. 2i, but not in Fig. 2c and f. Although the estimated chirp rate is quite different from the parameter we set, it is not surprising to us, due to the existence of noise and the quantitative effect.

Quantitative effect
In the second simulation, 500 pulses of jamming signal and target echo are respectively considered to verify the effectiveness of the quantitative effect. The radar signal's parameter is the same as above. Both of the linear and non-linear (quadratic) scenarios in (3) and (26) are considered. The V f is 20 kHz/s, which is equivalent with the acceleration 200 m/s 2 in [28]. The A f is set to 100 kHz/s 2 which is relatively large to obtain an intuitional result. The scaling factor of C is 2e5. The JNR and SNR are 15 dB and 10 dB, respectively. The estimated values of frequency and chirp rate of the target and linear VGPO in various pulses are shown in Fig. 3a, and the ones of the target and quadratic VGPO are shown in Fig. 3b. For comparison, we also plot the real shifted Doppler frequency of the jamming signal in each pulse.
The estimated frequency of the jamming signal stepwise increases in both linear and quadratic scenarios. In Fig. 3a, the estimated chirp rate of the jamming signal increases linearly from 0 to 0.08 s and increases quadratically in Fig. 3b, which is coincident with our theoretical analysis. It should be noticed that there always appears a flat area where (2020) 2020:37 Page 11 of 18 Fig. 2 The DCFT results of the target echo and VGPO jamming signal with various scaling factor C. a Jamming (without scaling). b Target (without scaling). c Profiles of a and b after normalization. d Jamming (C = 4e3). e Target (C = 4e3). f Profiles of a and b after normalization. g Jamming (C = 2e7). h Target (C = 2e7). i Profiles of g and h after normalization the estimated chirp rate is 0 (e.g., from 0.062 to 0.074 s in Fig. 3b). This is because the estimated frequency increases at this point, that is to say,f dj in (18) is quantized into the next frequency bin. Correspondingly, the remaining Doppler frequency f m has to subtract f s /N, which causes the estimated chirp rate l to become negative. However, the value range of l is from 0 to N − 1. As a result, the estimated chirp rate l becomes 0 which is the closest one to the theoretical value.

Influence of different parameters to DCFT
We then explore the influence of different signal parameters, pulse numbers, and SNR (or JNR) on identifying the VGPO jamming by utilizing the proposed feature. We assume the bandwidth of the signal is fixed as 5 MHz, and the sampling frequency is 10 MHz, while the pulse duration changes from 20 to 100 μs with duty cycle 20%. The carrier frequency is 15 GHz, and the scaling factor C is 2e5.
To inquire about the question of how long it would take for the proposed method to identify the jamming signal successfully, we jointly consider the influence of the signal's time width and pulse numbers. The linear pull-off scenario is considered. V f is still set to 20 kHz/s. The MSVR of the target echo and the VGPO jamming under various signal parameters and pulse numbers are shown in Fig. 4a and b, respectively. The MSVR of the target echo is steady around 0.45 after 100 pulses, even under different signal time width from Fig. 4a. On the contrary, the MSVR of the VGPO jamming is various under different time width and different pulse numbers in Fig. 4b, and its maximum is about 2.96. By comparing Fig. 4a and b, we can find that the MSVR of jamming signal is always bigger than that of the target echo.
Overall, the MSVR of VGPO jamming increases as the number of pulses increases until the estimated frequency changes. The wider the time width, the faster the value of MSVR increases. It should be noticed that the ladder-shaped portion is caused due to a change in the estimated frequency.
We plot the profiles in Fig. 4b when time width are 50 μs, 80 μs, and 100 μs, respectively, and show them in Fig. 5. We also draw the lines corresponding to the number of pulses We then study the influence of SNR on the performance of the proposed feature. We consider the scenario when the pulse duration is 20 μs according to the result in Fig. 4b, the total time of interference implementation is 0.08 s (800 pulses).
Besides a target in uniform speed, an accelerated target with an acceleration 2.7 m/s 2 , a linear and non-linear pull-off VGPO jamming are also considered. The V f of the linear pull-off jamming is 2 kHz/s. A non-linear pull-off case is a simple form of the maneuvering scenario, and its A f is 20 kHz/s 2 . The JSR is fixed as 3 dB, and the SNR changes from 0 to 15 dB. The relationships between MSVR and SNR are shown in Fig. 6. It can be seen that two kinds of VGPO jamming almost have the same MSVR. It means that the quadratic term of the frequency pull-off function in (26) has less effect on MSVR in a very short time (0.1 s) if the coefficient of the quadratic term is not too large (e.g., 20 kHz/s 2 ).
Different from the results of the linear and non-linear VGPO jamming, the MSVRs of the target are almost steady which are similar to the result in Fig. 4a. It means that the MSVR of the target signal is unchanged even with various parameters. Moreover, the MSVR of the target and the accelerated target are almost the same. Combined with the MSVR of the two jamming signals, we can conclude that a tiny acceleration or small maneuvering on target hardly changes the value of MSVR.

Discrimination performance
Before we conduct the discrimination performance of the proposed method, we consider the influence of the speed of frequency offset, i.e., the value of V f , on MSVR. The sampling frequency is 20 MHz, and the SNR and JSR are 10 dB and 5 dB, respectively. The other parameters are the same as the previous subsection. The values of V f are separately assumed as 20 kHz/s and 50 kHz/s. The MSVR of target echo and jamming signal versus pulse numbers with various V f are shown in Fig. 7 where VGPO1 and VGPO2 separately denote the linear VGPO jamming with V f of 20 kHz/s and 50 kHz/s. The MSVR of the jamming signal with V f of 50 kHz/s is larger than that of the jamming with V f of 20 kHz/s. Furthermore, the more pulses involved in the calculation, the greater the MSVR of the jamming signal, which is consistent with the result in Fig. 4b.
To compare the performance with the results in [18] and [19], 800 and 850 pulses are separately considered, with the same parameters in Fig. 6. The value of V f is set as 20 kHz/s, and 300 Monte Carlo trials are considered. As the value of MSVR of the target signal is stable under different radar parameters, the threshold can be set as the average between the minimum MSVR of the jamming signal, which is usually the corresponding value when the SNR is 0 dB, and the stable value of MSVR of the target. In this simulation, the threshold is set as 0.65 as the MSVR of jamming signal is around 0.8 when the SNR is 0 dB in Fig. 6. Figure 8 shows the identification probability of jamming signal versus SNR based on the proposed method. From Fig. 8, when the SNR is 0 dB, the identification probability of VGPO jamming could reach 0.975 with 850 pulses and 0.81 with 800 pulses. When the SNR is 2 dB, the jamming signal can always be correctly recognized.
From the aspect of the total data needed to correctly identify the jamming signal, our method only uses the signal data within 0.085 s and is less than the method in [18] which needs the signal data within 0.096 s. From the aspect of the identification probability, our method can correctly recognize the jamming signal when the SNR is 2 dB, which is smaller than that in [18] (the blue dotted line in Fig. 8). It should be noted that the value of V f is 50 kHz/s in [18], which also means our method would be better if the same value of V f is used based on the result in Fig. 7.

Discussion
Although this method can recognize the VGPO jamming with less data and a higher recognition rate compared to the open literature, it is not suitable for the jamming of range gate pull-off (RGPO) with uniform acceleration. This is because its chirp rate after de-chirping is pretty small (about 0.2 for the false target with the same motion state with the same simulation parameters in this paper). However, as a more general form of VGPO, the range-velocity gate pull-off (R-VGPO) jamming can also be identified by this method. Even in mixed-signal including target echo and jamming, this method is still effective.

Conclusions
In this paper, the discrimination issue between the target echo and the VGPO jamming is studied. To this purpose, we resort the discrete chirp-Fourier transform to jointly estimate the chirp rate and frequency of a signal, and the scaling characteristic of the DCFT algorithm is studied and explored to increase the estimated accuracy in chirp rate domain. Based on the actual joint estimation results, the quantitative effect of the VGPO signal is further analyzed; it guarantees that the relationship between the pull-off acceleration and the time can be judged from the relationship between the estimated chirp rate and the pulse numbers. To distinguish VGPO jamming and target echo, the MSVR is defined. Simulation results indicate that for a timestamp when the estimated frequency is unchanged, the MSVR of VGPO jamming increases with the pulse numbers and time width (with the same duty cycle) increases and is always larger than that of a target which is stable. Compared to other methods, the proposed method can recognize the jamming signal with less data and have a higher recognition; it means a great potential in practical applications.