Next Article in Journal
A General Paradigm for Retrieving Soil Moisture and Surface Temperature from Passive Microwave Remote Sensing Data Based on Artificial Intelligence
Next Article in Special Issue
An Optimization Method for Collaborative Radar Antijamming Based on Multi-Agent Reinforcement Learning
Previous Article in Journal
Can the Perception Data of Autonomous Vehicles Be Used to Replace Mobile Mapping Surveys?—A Case Study Surveying Roadside City Trees
Previous Article in Special Issue
DRFM-Based Repeater Jamming Reconstruction and Cancellation Method with Accurate Edge Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radar Maneuvering Target Detection Based on Product Scale Zoom Discrete Chirp Fourier Transform

1
Electronic Information School, Wuhan University, Wuhan 430072, China
2
Mingfei Weiye Technology Co., Ltd., Wuhan 430061, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(7), 1792; https://doi.org/10.3390/rs15071792
Submission received: 8 February 2023 / Revised: 19 March 2023 / Accepted: 22 March 2023 / Published: 27 March 2023
(This article belongs to the Special Issue Advanced Radar Signal Processing and Applications)

Abstract

:
Long-time coherent integration works to significantly increase the detection probability for maneuvering targets. However, during the observation time, the problems that often tend to occur are range cell migration (RCM) and Doppler frequency cell migration (DFCM), due to the high velocity and acceleration of the maneuvering target, which can reduce the detection of the maneuvering targets. In this regard, we propose a new coherent integration approach, based on the product scale zoom discrete chirp Fourier transform (PSZDCFT). Specifically, by introducing the zoom operation into the modified discrete chirp Fourier transform (MDCFT), the zoom discrete chirp Fourier transform (ZDCFT) can correctly estimate the centroid frequency and chirp rate of the linear frequency-modulated signal (LFM), regardless of whether the parameters of the LFM signal are outside the estimation scopes. Then, the scale operation, combined with ZDCFT, is performed on the radar echo signal in the range frequency slow time domain, to remove the coupling. Thereafter, a product operation is executed along the range frequency to inhibit spurious peaks and reinforce the true peak. Finally, the velocity and acceleration of the target estimated from the true peak position, are used to construct a phase compensation function to eliminate the RCM and DFCM, thus achieving coherent integration. The method is a linear transform without energy loss, and is suitable for low signal-to-noise (SNR) environments. Moreover, the method can be effectively fulfilled based on the chirp-z transform (CZT), which prevents the brute-force search. Thus, the method reaches a favorable tradeoff between anti-noise performance and computational load. Intensive simulations demonstrate the effectiveness of the proposed method.

Graphical Abstract

1. Introduction

In recent years, radar maneuvering target detection has attracted extensive research interest in spaceborne, airborne, and ground-based radar [1,2,3], due to the growing demand for practical applications, such as for high-resolution remote sensing imaging and high precision tracking [4,5,6,7,8]. However, it is often difficult for radar to detect high maneuvering targets, due to the weak radar echoes caused by their low radar cross section (RCS) [9,10,11,12,13]. Long-time coherent integration can enhance the signal-to-noise ratio (SNR) of the echo signal and increase the detection probability of maneuvering targets [14,15,16]. However, during coherent processing, the high velocity and acceleration of the maneuvering target can induce linear range cell migration (LRCM), quadratic range cell migration (QRCM), and Doppler frequency cell migration (DFCM). Not only that, when the radar pulse repetition frequency (PRF) is low, the excess motion parameters of the maneuvering target may not be estimated correctly. These adverse effects can seriously degrade the detection performance of traditional methods (e.g., moving target detection, MTD). Therefore, in order to improve the detection performance, it is imperative to propose an effective method to solve these problems.
For the elimination of LRCM, many successful methods have been proposed. Representative methods involve the keystone transform (KT) [17,18], Radon Fourier transform (RFT) [19,20], axis rotation transform (ART) [21], scaled inverse Fourier transform (SIFT) [22], and frequency-domain deramp-keystone transform (FDDKT) [23], where the KT is a linear transform for eliminating LRCM of multiple targets without prior information. However, the KT may suffer from Doppler ambiguity, due to the high velocity and low PRF. The RFT and ART eliminate the LRCM through the two-dimensional search procedure, which may require a heavy computational burden. The SIFT and FDDKT can achieve the cancellation of LRCM with the help of the symmetric autocorrelation function, which can maintain a low computational burden at the cost of anti-noise performance. Nevertheless, the coherent gain of the above method can decline sharply, due to overlooking the QRCM and DFCM arising from the target’s acceleration.
To address this issue, various methods have been introduced. The generalized RFT (GRFT) [24], is a classical method which extracts the target trajectory through a three-dimensional search, and constructs a corresponding filter to accomplish coherent integration. However, the consequent problem is that the GRFT has severe blind speed side lobes (BSSLs) interference, due to the nature of its filter construction. Based on this, the Radon-Lv’s distribution (RLVD) [25] and Radon-fractional Fourier transform (RFRFT) [26] were developed, and they can successfully avert the effect of BSSLs due to their adoption of LVD and FRFT as filters, respectively. Unfortunately, the aforementioned methods are prohibitively computationally intensive and are not conducive to real-time processing. In this respect, improved axis rotation fractional Fourier transform (IAR-FRFT) [27], modified axis rotation transform and Lv’s transform (MART-LVT) [28], and modified location rotation transform and improved discrete chirp Fourier transform (MLRT-IDCFT) [29], can slightly diminish the computational burden, by reducing the search dimension. However, these three methods come at the cost of not being able to eliminate the QRCM, which may degrade some of the accumulation gain. Moreover, the computational burden of these methods remain considerable. To further drastically reduce the computational cost for real-time processing, some rank-reduction based methods have been widely investigated. Typical methods include, the adjacent cross-correlation function (ACCF) [30], three-dimensional scaled transform (TDST) [31], second-order Wigner–Ville distribution (SoWVD) [32], and frequency autocorrelation function and Lv’s distribution (FAF-LVD) [33]. These methods can reduce the rank of the echo signal by carrying out correlation operations, thus avoiding the brute-force search and reducing the computational load. Although these methods have low computational complexity and facilitate real-time processing, they are only suitable for high SNR environments due to the loss of signal energy and cross terms caused by correlation operations. Besides, time-reversal transform (TRT)-based methods, such as improved axis rotation and time reversal transform (IAR-TRT) [34], keystone transform and time reversal transform (KT-TRT) [35], and phase compensation and time-reversal transform (PC-TRT) [36], can also obtain efficient calculations, by implementing TRT operations. However, similar to the rank-reduction-based approaches, the anti-noise performance of TRT-based methods can also drop tremendously in a low SNR environment, due to the fact that TRT operations are nonlinear.
With the aim of striking a better balance between computational cost and anti-noise performance, this paper proposes a novel coherent integration approach, based on the product scale zoom discrete chirp Fourier transform (PSZDCFT). The basic idea of the method, is that the velocity and acceleration of the maneuvering target are estimated via PSZDCFT, then the phase compensation function is constructed to eliminate the RCM and DFCM, and finally, the slow-time Fourier transform is applied, to accomplish coherent integration. The zoom operation is used to extend the estimation scopes of the modified discrete chirp Fourier transform (MDCFT) [37] when the velocity and acceleration of the target are superabundant, which also results in the occurrence of other false peaks in the two-dimensional spectrum. The scale operation is used to eliminate the coupling between the true peak location and the range frequency, while making the false peak locations couple with the range frequency. The product operation is used to identify and reinforce the true peak, while suppressing the false peaks. The most important features of the proposed method are, extending the estimation ranges of MDCFT by zoom operation, and improving the parameter estimation performance by product operation. The former can be applied to the uniform motion target with velocity ambiguity, and the latter can be developed into other parameter estimation methods, such as dechirp estimation. Moreover, in contrast to the rank-reduction-based and time-reverse-based methods, the proposed method is a linear transform without loss of target energy and no cross terms for multiple targets. Compared with the traditional methods of brute-force search, such as GRFT and RLVD, this method is search-free, thanks to its efficient implementation by fast Fourier transform (FFT), inverse fast Fourier transform (IFFT), and complex multiplication, which greatly reduces the computational burden, without significant degradation in detection performance. Thus overall, the presented approach reaches a favorable tradeoff between anti-noise performance and computational complexity.
The rest of this paper is arranged as follows. In Section 2, the signal model is formulated and the problem to be addressed is presented. Section 3 details the development and rationale of the proposed method. In Section 4, the relevant analysis of the proposed approach is given. Section 5 provides the simulation results. Finally, the conclusions are drawn in Section 6.

2. Signal Model and Problem Formulation

Assume that a pulsed radar transmits a narrow-band linear frequency-modulated signal (LFM), as follows
s ( t ^ , η m ) = rect t ^ T p exp j π γ t ^ 2 exp j 2 π f c t ^ + η m ,
where T p and rect ( · ) represent the pulse width and rectangular function, respectively. γ and t ^ are the frequency modulated rate and fast time, respectively. f c is the carrier frequency. η m = m / f r ( m = M / 2 , , 1 , 0 , 1 , , M / 2 1 ) is the slow time. f r = 1 / T r and M indicate the PRF and pulse number, respectively. T r is the pulse repetition interval (PRI).
We consider a second-order motion model and assume that the instantaneous slant range R ( η m ) , between the target and the radar is
R η m = R 0 + v 0 η m + 1 2 a 0 η m 2 ,
where R 0 , v 0 , and a 0 are the initial range, radial velocity, and acceleration of the target, respectively.
Then, after the signal demodulation, the received baseband signal is stated as
s r t ^ , η m = σ 0 rect t ^ 2 R η m / c T p × exp j 4 π f c c R η m × exp j π γ t ^ 2 R η m c 2 ,
where σ 0 and c indicate the echo amplitude and velocity of light, respectively. After pulse compression (PC) using matched filters, the received signal can be noted as
s c t ^ , η m = σ 1 sinc B t ^ 2 R η m c exp j 4 π f c R η m c ,
where σ 1 and B are the signal amplitude after PC and the signal bandwidth, respectively.
Then, performing FFT on Equation (4) along the fast time t ^ , yields
S c f , η m = σ 2 rect f B exp j 4 π f + f c R η m c ,
where f is the range frequency corresponding to the fast time t ^ , and σ 2 indicates the signal amplitude after performing FFT.
Substituting Equation (2) into Equation (5), one obtains
S c f , η m = σ 3 exp j 4 π f c R 0 + v 0 η m + 1 2 a 0 η m 2 × exp j 4 π λ R 0 + v 0 η m + 1 2 a 0 η m 2 ,
where σ 3 = σ 2 rect f B . λ = c / f c is the wavelength.
Since the range frequency is coupled to slow time, the first exponential term of Equation (6) could bring about the LRCM due to velocity and QRCM due to acceleration. The second exponential term may result in the DFCM induced by acceleration. Therefore, it is critical to accurately estimate the velocity and acceleration of the maneuvering target, to eliminate the adverse effects of the RCM and DFCM and improve the detection performance.

3. Description of the Proposed Method

3.1. MDCFT and Its Limitations

In [38], Xia proposes the DCFT to estimate the centroid frequency and chirp rate of the LFM signal. However, this method has the requirement that the total sampling number of the signal is prime and the parameters of the signal are integers. With this in mind, the MDCFT is proposed, which eliminates two constraints of DCFT. The principle of MDCFT is as follows.
The analog form of the LFM signal can be expressed as
z η m = σ exp j 2 π f 0 η m + μ 0 η m 2 ,
where σ and η m denote the signal amplitude and slow time, respectively. f 0 and u 0 are the centroid frequency and chirp rate, respectively. After discretization, Equation (7) can be written as
z m = σ exp j 2 π f 0 T r m + μ 0 T r 2 m 2 = σ W M k 0 m + l 0 M m 2 ,
where k 0 = f 0 T r M and l 0 = μ 0 T r 2 M 2 represent the digital centroid frequency and chirp rate of the signal, respectively. W M = exp ( j 2 π / M ) .
The MDCFT of z ( m ) is defined as
Z M D C F T ( k , l ) = 1 M m = M / 2 M / 2 1 z ( m ) W M k m + l M m 2 ,
where k [ M / 2 , M / 2 1 ] and l [ M / 2 , M / 2 1 ] denote the centroid frequency index and chirp frequency index, respectively. It is obvious that the estimation scopes of both k 0 and l 0 are [ M / 2 , M / 2 ) . With the substitution of Equation (8) into Equation (9), we can obtain
Z M D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M 1 M l l 0 m 2 W M k k 0 m .
Depending on whether k 0 and l 0 exceed the estimation scopes, the peak position of Z M D C F T ( k , l ) can be divided into the following cases.
Case 1: k 0 < M / 2 and l 0 < M / 2 . This means that both k 0 and l 0 are within their estimation scopes, and as shown in Equation (10), the peak of the two-dimensional spectrum can be observed at ( k 0 , l 0 ) .
Case 2: k 0 M / 2 and l 0 < M / 2 . In this case, k 0 is over its valid estimation range, which implies that Doppler ambiguity would occur. k 0 can be written as
k 0 = k 0 + p M ,
where p = r o u n d ( k 0 / M ) is the number of the Doppler ambiguity. k 0 ( M / 2 , M / 2 ) . Inserting Equation (11) into Equation (10), one obtains
Z M D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M 1 M l l 0 m 2 W M k k 0 m .
Here, we use the constant equation W M p M m = 1 . As shown in Equation (12), we can see that the peak is located at ( k 0 , l 0 ) in the two-dimensional spectrum, which is the folded peak. From the coordinates of the folded peak, we cannot get the correct digital centroid frequency.
Case 3: M / 2 l 0 < M 2 / 2 . At this point, since l 0 is out of the estimation scope, k 0 and l 0 cannot be accurately estimated, regardless of the value of k 0 . In the two-dimensional spectrum, Z M D C F T ( k , l ) has no distinctive sharp peak.
Case 4: M 2 / 2 l 0 . In this case, l 0 is still out of the estimation scope, but similar to case 2, chirp rate ambiguity can occur. l 0 can be written as
l 0 = l 0 + h M 2 ,
where h = r o u n d ( l 0 / M 2 ) is the number of the chirp rate ambiguity. l 0 ( M 2 / 2 , M 2 / 2 ) . By taking Equation (13) into Equation (10), we have
Z M D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M 1 M l l 0 m 2 W M k k 0 m .
W M 1 M h M 2 m 2 = 1 is used here. The analysis of Equation (14) is the same as the previous three cases, and it is worth noting that the correct value of l 0 cannot be obtained in this case. In practice, chirp rate ambiguity is almost impossible, that is, l 0 M 2 / 2 . Therefore, we will not analyze this case later.

3.2. Zoom DCFT

In the previous subsection, we expounded the principle and limitations of the MDCFT. In order to remove these limitations and make MDCFT more practical, we proposed the zoom DCFT method (ZDCFT).
The ZDCFT of z ( m ) is defined as follows
Z Z D C F T ( k , l ) = ZDCFT [ z ( m ) ] = 1 M m = M / 2 M / 2 1 z ( m ) W M α k m + β l M m 2 ,
where k [ M / 2 , M / 2 ) and l [ M / 2 , M / 2 ) represent the centroid frequency index and chirp frequency index, respectively. α N * and β N * denote the zoom factors of the centroid frequency and chirp rate, respectively.
In order to clarify the role of the ZDCFT, we analyze the first three cases in the previous subsection one by one.
Inserting Equation (8) into Equation (15), we obtain
Z Z D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M l l l 0 / β m 2 W M k k k 0 / α m ,
where M l = M 2 / β and M k = M / α are the new implied periods of Z Z D C F T ( k , l ) in the k and l domains, respectively.
For case 1, i.e., k 0 < M / 2 and l 0 < M / 2 , we should just take α = β = 1 . Then, the ZDCFT degenerates into the MDCFT. The peak of the theoretical two-dimensional spectrum is located at ( k 0 , l 0 ) .
For case 2, i.e., k 0 M / 2 and l 0 < M / 2 , we know that l 0 < M / 2 , so β is taken as 1. Since k 0 M / 2 exceeds the estimation scope, we should take the appropriate α ( α > 1 ) to extend the estimation scope of k 0 , i.e.,
[ M 2 α , M 2 α ) .
In other words, it is to make k 0 / α lie within the computational range of k, i.e., k 0 / α [ M / 2 , M / 2 ) . However, the ensuing problem is that Z M D C F T ( k , l ) would generate α peaks in the k domain. This is because the length of the computational interval in the k-domain is α times longer than its implied period by the zoom operation, i.e., M / M k = α . Therefore, with β = 1 , Equation (16) can be rewritten as follows
Z Z D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M l l l 0 m 2 W M k k k 0 / α F M k m ,
where F Z , M l = M 2 , and M k = M / α . It can be seen from Equation (18) that these α peaks are located at
k ^ = k 0 / α + F M k l ^ = l 0 k ^ [ M / 2 , M / 2 ) l ^ [ M / 2 , M / 2 ) F Z .
When F = 0 , we can obtain the correct digital centroid frequency and chirp rate, i.e., ( k ^ 0 , l ^ 0 ) = ( α k ^ , l ^ ) .
Example 1.
Figure 1 illustrates the role of the zoom factor of the centroid frequency α. The simulation parameters of the LFM signal are: μ 0 = 45 Hz/s, f 0 = 630 Hz, f r = 400 Hz, and M = 800 . Then, we can get k 0 = f 0 M / f r = 1260 and l 0 = μ 0 M 2 / f r 2 = 180 . It can be seen that k 0 is out of the estimation scope while l 0 is not. Therefore, we choose α = 4 and β = 1 , which means that four peaks would be observed in the two-dimensional spectrum, of which the true peak is located at (315, 180). With the help of the new implied period M k = M / α = 200 , we can also derive the theoretical x-coordinate of one of the false peaks according to Equation (19), i.e., 85 = 315 2 × 200 . The value of the y-coordinate is the same as the true peak. The result of MDCFT is given in Figure 1a. It can be seen that, although there is a peak in the two-dimensional spectrum, the x-coordinate of this peak is the folded digital centroid frequency, i.e., k 0 = 340 . Figure 1b shows the result of ZDCFT, which is in agreement with the theory.
For case 3, i.e., M / 2 l 0 < M 2 / 2 , as in case 2, we should select a suitable β ( β > 1 ) to make l 0 / β lie within M / 2 . To put it another way, the estimation scope of l 0 is extended by β , i.e.,
[ M 2 β , M 2 β ) .
Note that, it is sufficient that β is less than M according to Equation (20) and that l 0 is less than M 2 / 2 . Thus, unlike the analysis in case 2, the length of the computational interval in the l domain is still less than its new implied period, i.e., M < M l = M 2 / β , s . t . β < M . From another perspective, the computational range of l is a subset of [ M l / 2 , M l / 2 ) , i.e., [ M / 2 , M / 2 ) [ M l / 2 , M l / 2 ) , which implies that Z M D C F T ( k , l ) would generate only one peak in the l domain, independent of β . Considering k 0 , Equation (16) is rewritten as follows
Z Z D C F T ( k , l ) = σ 1 M m = M / 2 M / 2 1 W M l l l 0 / β m 2 W M k k k 0 / α F M k m ,
where M l = M 2 / β , and M k = M / α . From Equation (21), we can see that in the two-dimensional spectrum, there are still α peaks, which appear at
k ^ = k 0 / α + F M k l ^ = l 0 / β k ^ [ M / 2 , M / 2 ) l ^ [ M / 2 , M / 2 ) F Z .
When F = 0 , we can obtain the correct digital centroid frequency and chirp rate, i.e., ( k ^ 0 , l ^ 0 ) = ( α k ^ , β l ^ ) .
Example 2.
In Figure 2, we show the joint action of the zoom factors α and β. The simulation parameters for the LFM signal are set as: μ 0 = 145 Hz/s, f 0 = 630 Hz, f r = 400 Hz, and M = 800 . Then, both k 0 and l 0 are out of the estimation scopes, i.e., k 0 = 1260 > M / 2 = 400 and l 0 = 580 > M / 2 . At this point, we have that α equals 4 and β equals 2, which indicates that there would be four peaks along the k-dimension and one peak along the l-dimension, in the two-dimensional spectrum. The true peak should be formed at (315, 290) in the two-dimensional spectrum. Similar to Example 1, the location of one of the three false peaks should appear at (−85, 290), based on Equation (22). Figure 2a shows the correct result of ZDCFT, in accordance with the theoretical analysis. For comparison, Figure 2b,c display the wrong results of ZDCFT, when α = 1 and β = 2 and when α = 4 and β = 1 , respectively. It is evident that the former ( α = 1 , β = 2 ) is blurred in the Doppler domain, while the latter ( α = 4 , β = 1 ) is invalid because the chirp rate cannot be correctly estimated.
Remark 1.
The zoom factors for α and β should be as small as possible, with the guarantee that the estimation ranges contain k 0 and l 0 [39]. The reason for this is that the spatial resolutions in the k and l domains are proportional to α and β, respectively, i.e., ρ k = α / M and ρ l = β / M 2 .

3.3. Scale ZDCFT

The previous subsection shows that ZDCFT can obtain correct parameter estimations, regardless of whether the centroid frequency and chirp rate of the LFM signal exceed the estimation scopes of MDCFT. For the two-dimensional radar echo signal in the range frequency-slow time domain, it is an LFM signal in the slow-time dimension with a certain range frequency cell. Therefore, we can use the ZDCFT to estimate the velocity and acceleration of the maneuvering target. However, we may obtain different values of the velocity and acceleration for different range frequency cells via ZDCFT. Considering this problem, we propose the scale ZDCFT method (SZDCFT).
The discretized form of Equation (6), for the slow-time variable η m , is as follows
S c f , m = σ 3 exp j 2 π ψ ( f ) 2 R 0 λ + 2 v 0 λ T r m + a 0 λ T r 2 m 2 = σ 4 W M ψ ( f ) K 0 m + L 0 M m 2 ,
where σ 4 = σ 3 exp j 4 π λ ψ f R 0 denotes the complex amplitude of the signal. K 0 = 2 v 0 λ M T r and L 0 = a 0 λ M 2 T r 2 represent the digital centroid frequency and chirp rate of the target, respectively. ψ ( f ) = 1 + f / f c is the scale factor.
Then, the SZDCFT of Equation (23) is defined as
S S Z D C F T ( f , k , l ) = SZDCFT [ S c ( f , m ) ] = 1 M m = M / 2 M / 2 1 S c ( f , m ) W M ψ ( f ) α k m + β l M m 2 = σ 4 1 M m = M / 2 M / 2 1 W M l f ( l L 0 / β ) m 2 W M k f ( k K 0 / α F M k f ) m ,
where M l f = M 2 β ψ ( f ) and M k f = M α ψ ( f ) denote the implied periods of S S Z D C F T ( f , k , l ) in the l and k domains associated with the range frequency, respectively.
Similarly, in the two-dimensional spectrum, we can observe from Equation (24) that there are α peaks, which are situated at
k ^ = K 0 / α + F M k f l ^ = L 0 / β k ^ [ M / 2 , M / 2 ) l ^ [ M / 2 , M / 2 ) F Z .
In Equation (25), the true peaks will appear at the location of ( K 0 / α , L 0 / β ) when F = 0 , which are irrespective of the range frequency. Therefore, it can be seen that the SZDCFT can neutralize the effect of the range frequency.
Example 3.
Figure 3 simulates the results of SZDCFT for two different scale factors, ψ 1 = 1.01 and ψ 2 = 0.99 . The parameters of the LFM signal are the same as in Example 2. Correspondingly, the zoom factors of the centroid frequency α, and chirp rate β, are equal to 4 and 2, respectively. In the two-dimensional spectrum, the true peak would be formed at (315, 290). The results of SZDCFT for the scale factors ψ 1 and ψ 2 , are depicted in Figure 3a,b, respectively, which are consistent with the theoretical derivation. As a comparison, Figure 3c,d show the parameter estimations via ZDCFT, in which it can be seen that the coordinates of the peak are coupled to the scale factors.
Remark 2.
It should be noted that the selection criteria for the values of α and β, are the same as that of the ZDCFT. This is because, although the existence of the scale factor ψ ( f ) could make the implied periods of the SZDCFT different from those of the ZDCFT, fortunately, for the narrowband radar system, there are f f c and ψ ( f ) 1 . Therefore, when judging whether the parameters to be estimated exceed the estimation scopes of SZDCFT, the changes in the implied periods caused by the scale factor ψ ( f ) can be ignored.
Here, we give more specific selection criteria for the values of α and β , combined with the target motion parameters.
The velocity of the maneuvering target can be written as
v 0 = N a m b v a m b + v r e s ,
where N a m b and v a m b = λ f r 2 are the ambiguity integers and blind speed, respectively. v r e s < v a m b / 2 is the unambiguous velocity.
When the target’s velocity is not ambiguous, i.e., N a m b = 0 and v 0 = v r e s < v a m b / 2 , then we have K 0 = 2 v r e s λ M T r < M / 2 . Evidently, α is taken as 1.
When the target’s velocity is ambiguous, i.e., N a m b 0 , we have K 0 = 2 v 0 λ M T r = N a m b M K 0 , where K 0 = 2 v r e s λ M T r ( M / 2 , M / 2 ) . This is similar to Equation (11). Therefore, in practical applications, we usually determine α according to the velocity scope of interest. For example, we assume that the velocity scope of interest is [ v m a x , v m a x ] . Then, we can obtain the maximum value of K 0 as K 0 , m a x = 2 v m a x λ M T r . According to Equation (25), in order to ensure that K 0 , m a x can be estimated, α should satisfy the following inequality
α 4 v m a x λ f r α N * ,
where · is the ceiling operation. In order to have a better estimation accuracy of the velocity, we generally take α = 4 v m a x λ f r .
Similar to the analysis of velocity, for the acceleration of the maneuvering target a 0 , we assume that the acceleration scope of interest is [ a m a x , a m a x ] . Then, we can also obtain the maximum value of L 0 , m a x , i.e., L 0 , m a x = a m a x λ M 2 T r 2 . Similarly, based on Equation (25), we can also obtain the constraint of β as follows
2 a m a x M λ f r 2 β < M β N * .
From the analysis of β in the previous subsection, it can be known that the SZDCFT would only produce one peak in the l domain under the constraint of β < M . Here, we default to the acceleration being unambiguous, which is often the case in practice.

3.4. Product Operation and Coherent Integration

The SZDCFT uncouples the slow time and range frequency, and extends the parameter estimation scope. However, we are unable to identify the true peak when there are multiple peaks in the k domain. To solve this problem, we propose the product SZDCFT method (PSZDCFT).
By re-examining Equations (24) and (25), we can find that: (1) For different range frequency cells, S S Z D C F T ( f , k , l ) has different implied periods in the k domain ( M k f = M α ψ ( f ) ), which means that when there are multiple peaks, the interval between peaks in the two-dimensional spectrum varies with the change in the range frequency. (2) The true peak position is always constant in the two-dimensional spectrum, i.e., ( K 0 / α , L 0 / β ) , while the positions of the false peaks are changed with the range frequency. The comparison between Figure 3a,b supports the above assertions very well.
Therefore, it naturally occurs to us to perform a product operation on the SZDCFT results for all range frequency cells, which can curb the false peaks and heighten the true peak. Then, the PSZDCFT method is defined as
S P S Z D C F T k , l = Product [ S S Z D C F T ( f , k , l ) ] = f = B / 2 B / 2 S S Z D C F T f , k , l σ P S Z D C F T · δ k K 0 / α , l = L 0 / β ,
where σ P S Z D C F T is the amplitude of accumulation.
From Equation (29), it can be seen that there will be only true peak in the two-dimensional spectrum, while the false peaks are eliminated. Then, based on the peak location, we can obtain the velocity and acceleration of the maneuvering target as follows
k m a x , l m a x = arg max k , l Product SZDCFT S c ( f , m ) v ^ 0 = α λ f r k m a x 2 M a ^ 0 = β λ f r 2 l m a x M 2 .
After obtaining the estimated velocity and acceleration of the maneuvering target, we can construct the phase compensation function to eliminate the RCM and DFCM, i.e.,
H P C F f , η m = exp j 4 π c f v ^ 0 η m exp j 2 π c ( f + f c ) a ^ 0 η m 2 .
Finally, we can achieve the coherent integration by applying the IFFT and FFT with respect to f and η m .
S C I t ^ , f η m = FFT η m IFFT f S c f , η m · H P C F f , η m = σ C I sinc B t ^ 2 R 0 c × sinc C P I f η m + 2 v r e s λ ,
where σ C I and C P I denote the amplitude and coherent processing interval, respectively. f η m is the folded Doppler frequency with respect to η m .

4. Analysis of the Proposed Method

4.1. PSZDCFT for Multi-Targets

The SZDCFT is a linear transform, which can avert the interference of cross terms. Therefore, for multi-target scenarios, the proposed approach is still applicable. It should be noted that the product operation may annihilate the weak target when there are significant differences in the echo amplitudes between targets. Fortunately, the “CLEAN” technique [40] can be used, to remove the strong target effect.

4.2. Implementation of SZDCFT

The SZDCFT of S c ( f , m ) can be rewritten as follows
S S Z D C F T ( f , k , l ) = 1 M m = M / 2 M / 2 1 S c ( f , m ) W M ψ ( f ) α k m + β l M 2 m 2 = 1 M m = M / 2 M / 2 1 S c ( f , m ) W M ψ ( f ) β l M 2 m 2 × exp j 2 π M ψ ( f ) α k m .
Equation (33) shows that, for each fixed l, S S Z D C F T ( f , k , l ) is the scaled Fourier transform (SFT) of the signal S c ( f , m ) W M ψ ( f ) β l M 2 m 2 . As we know, the SFT can be efficiently implemented by FFT-based chirp-z transform (CZT).
Assuming a discrete signal x ( n ) , n = N / 2 , , N / 2 1 , the SFT of x ( n ) is defined as follows
X ξ k = n = N / 2 N / 2 1 x ( n ) exp j 2 π N ξ n k = n = N / 2 N / 2 1 x ( n ) W ˜ n k ,
where k = N / 2 , N / 2 1 and W ˜ = exp j 2 π N ξ . ξ is the scale factor. Then, the Bluestein equation is applied here, i.e.,
n k = 1 2 n 2 + k 2 ( n k ) 2 .
Inserting Equation (35) into Equation (34), one obtains
X ( ξ k ) = n = N / 2 N / 2 1 x ( n ) W ˜ n k = W ˜ k 2 2 n = N / 2 N / 2 1 x ( n ) W ˜ n 2 2 W ˜ ( k n ) 2 2 = W ˜ k 2 2 x ( n ) W ˜ n 2 2 W ˜ n 2 2 ,
where ⊛ represents convolution. Figure 4 illustrates the calculation process of Equation (36). Thus, the SFT of N points can be efficiently realized by complex multiplication, FFT and IFFT, which requires a computational load of O ( 3 N l o g 2 N ) .

4.3. Computational Complexity

The computational complexities of GRFT, RLVD, ACCF, PC-TRT, FAF-LVD, and the proposed method, are analyzed and compared in this subsection. Suppose M, N, N v , and N a represent the number of pulses, range cells, searching velocity, and searching acceleration, respectively.
The GRFT accomplishes coherent integration through a three-dimensional parameter search. Thus, its computational load is about O ( N N v N a M ) .
The RLVD firstly obtains the target trajectory through a three-dimensional search, which requires a computational cost of O ( N N v N a ) , and then accomplishes coherent integration via LVD, with a computational cost of O ( 3 M 2 l o g 2 M ) . Therefore, its total computational complexity is O ( 3 N N v N a M 2 l o g 2 M ) .
For the ACCF method, it uses an adjacent cross-correlation operation to correct the RCM and eliminate the DFCM, which requires a computational cost of O ( 2 M N l o g 2 N ) . Then, the computational cost required to achieve coherent integration is O ( M N ( l o g 2 M + l o g 2 N ) ) . Hence, its total computational cost is about O ( M N ( 3 l o g 2 N + l o g 2 M ) .
As for PC-TRT, it decouples the slow time and range frequency by constructing the phase compensation function with the searched velocity, and performing the time reversal transform, and finally realizes coherent integration using IFFT and FFT. Therefore, its computational cost is about O ( N v M N l o g 2 ( M N ) ) .
The main calculations of FAF-LVD contain second-order KT O ( 4 M N l o g 2 M ) , FAF O ( M N l o g 2 N ) , CZT-based SFT O ( 3 M N l o g 2 M ) , phase compensation O ( M N l o g 2 N ) , and LVD O ( 3 M 2 l o g 2 M ) . Therefore, its total computational burden is about O ( 3 M 2 l o g 2 M + M N ( 7 l o g 2 M + 2 l o g 2 N ) ) .
For each fixed range cell of the proposed method, the computational costs required by the SZDCFT operation and the product operation are O ( M 2 + 3 M 2 l o g 2 M ) and O ( M 2 ) , respectively. Therefore, the overall computational load of the PSZDCFT, is in the order of O ( N M 2 ( 2 + 3 l o g 2 M ) ) .
Table 1 shows the computational complexities of the above methods. We assume that N = N v = N a = M . Then, Figure 5 depicts the computational complexity curves. The RLVD and GRFT take up too much time and are not conducive to real-time processing. It can be seen that although the computational complexity of PC-TRT remains of the same order of magnitude as that of the present method, i.e., O ( M 3 l o g 2 M ) , the PC-TRT has an inferior detection performance under the condition of low SNR, which is due to the fact that the TRT operation loses the bulk of the signal energy. The FAF-LVD and ACCF obtain the lowest computational loads via the correlation operation, but there are also noticeable drops in the anti-noise performance. Hence, it can be inferred that the proposed method accomplishes a better compromise between computational burden and anti-noise performance.

4.4. Procedure of the Proposed Method

The flow chart of the proposed method is shown in Figure 6 and its main procedures are described in the following subsection.
Step 1: Perform PC on the received baseband signal s r ( t ^ , η m ) , then employ the FFT along the range dimension, to obtain S c ( f , η m ) .
Step 2: For a fixed range frequency f, apply the SZDCFT on S c ( f , η m ) , to obtain S S Z D C F T ( f , k , l ) .
Step 3: Determine the scope of the next range frequency, f. If f [ B / 2 , B / 2 ] , execute step 2. If f [ B / 2 , B / 2 ] , execute product operation, to obtain S P S Z D C F T ( k , l ) .
Step 4: Estimate the target’s velocity and acceleration according to Equation (30), and then construct the phase compensation function based on Equation (31).
Step 5: Perform the IFFT and FFT along the range dimension and slow-time dimension, respectively, to accomplish coherent integration.

5. Numerical Results

In this section, the validity and reliability of the presented method are verified by numerical simulations, with the use of the Matlab numerical software; the radar system parameters are listed in Table 2.

5.1. Coherent Integration for a Single Target

Figure 7 illustrates the coherent integration ability for a single target, where the motion parameters of the target are: initial range cell number n 0 = 451 st, radial velocity v 0 = 120 m/s, and radial acceleration a 0 = 45 m/s 2 . The input SNR after PC is 20 dB. Based on Equations (27) and (28), we take α = 4 and β = 2 , to extend the estimation scopes of velocity and acceleration to [−150, 150) m/s and [−75, 75) m/s 2 , respectively. Figure 7a displays the target motion trajectory after PC. It is apparent that serious RCM occurs to the target, due to the high velocity and acceleration. Figure 7b provides the result of MDCFT when f = 4 MHz. It can be seen that the MDCFT fails because the motion parameters exceed the estimation scopes of MDCFT. Figure 7c,d show the results of SZDCFT when f = ± 4 MHz. It can be observed that there are four peaks in each two-dimensional spectrum diagram. The true peak positions are the same, i.e., ( v ^ 0 , a ^ 0 ) = ( 120 , 45 ) , while the false peak positions are interlaced. The result of PSZDCFT is depicted in Figure 7e. As predicted, the false peaks are removed, while the true peak is reinforced and recognized. Then, based on the location of the true peak, we obtain the estimated velocity v ^ 0 and acceleration a ^ 0 as being 120 m/s and 45 m/s 2 , respectively. Finally, the estimated parameters are used to construct a phase compensation function, to eliminate the RCM and DFCM, and then the target energy can be well focused, in Figure 7f.
Figure 8 intuitively demonstrates the excellent ability of the product operation with respect to focusing energy and estimating parameters, in a low SNR environment. The input SNR is −8 dB after PC. Other parameters are identical to those in Figure 7. Figure 8a shows the result after PC, from which it can be observed that the target trajectory is completely annihilated in the noise. Figure 8b provides the result of SZDCFT at f = 0 MHz. It can be noticed that the peak is not visible, which is not favorable for parameter estimation. In this regard, we perform a product operation for the SZDCFT results along the range frequency cells, to accumulate the signal energy dispersed in different range frequency cells, which is illustrated in Figure 8c. Comparing Figure 7c,d and Figure 8b, it is easy to see that the product operation still has a remarkable parameter estimation ability in the low SNR environment, which is also confirmed in Section 5.4. On this basis, coherent integration can naturally also be accomplished well, in Figure 8d.

5.2. Coherent Integration for Multi-Targets

We also evaluate the coherent integration capability of the presented approach for multiple targets in Figure 9. The motion parameters of the three maneuvering targets are given in Table 3. Again, after PC, the input SNR is 20 dB. Similarly, depending on the motion parameters of the three maneuvering targets, we choose α = 4 and β = 2 . Three curved trajectories can be observed in Figure 9a, which indicates that the three targets undergo severe RCMs. Figure 9b,c depict the results of SZDCFT when f = ± 4 MHz. We can find that after performing the SZDCFT operation, the true peak positions are independent of the range frequency, while the false peak positions are coupled with the range frequency. Then, we perform a product operation on the results of SZDCFT for all range frequency cells, to obtain three spikes, shown in Figure 9d. Based on the three peak locations, we obtain the motion parameters of targets A, B, and C, i.e., v ^ A = 120 m/s, a ^ A = 45 m/s 2 , v ^ B = 108 m/s, a ^ B = 30 m/s 2 , v ^ C = 117 m/s, and a ^ C = 42 m/s 2 . Finally, the coherent integration results for targets A, B, and C are provided in Figure 9e,f,g, respectively.

5.3. Detection Performance

In this section, we evaluate the target detection capability of the presented method through Monte Carlo experiments. As a comparison, five other representative methods, i.e., RLVD, GRFT, FAF-LVD, PC-TRT, and ACCF, are also simulated. The constant false alarm rate (CFAR) technique, with false alarm rate P f a = 10 6 , is used. The input SNR after PC varies from −25 to 10 dB, with a step of 1 dB. For each SNR value, 500 independent simulations are run. Figure 10 depicts the detection probability curves.
It can be seen that the GRFT has the best detection performance, thanks to the multidimensional parameter search. The RLVD achieves a detection performance close to GRFT’s, again attributed to the parameter search and excellent LVD operation. Using GRFT as a baseline, the proposed method suffers an SNR loss of about 3 dB, when the detection probability p d = 0.9. This is because the proposed method estimates the parameters by the product operation, which introduces more random errors compared with GRFT, which is based on the maximum likelihood estimation (MLE), in an extremely low SNR environment. For a clearer explanation, Figure 11 offers the results of parameter estimation for the two methods at an extremely low SNR, of −12 dB. The simulation parameters are the same as those in Figure 7. Figure 11a gives the results of the proposed method by the product operation. It can be seen that at extremely low SNRs, multiple randomly distributed peaks appear in the plot after the product operation, due to the influence of noise, which deteriorates the parameter estimation performance. As a contrast, the results of GRFT, based on the MLE, is displayed in Figure 11b. Fortunately, due to the fact that the proposed approach is a linear transform with no energy loss, it can be observed from Figure 10 that the method has a superior anti-noise performance at low SNRs, compared with FAF-LVD, PC-TRT, and ACCF. Therefore, considering the processing time and detection performance, we can deduce that the presented approach is more applicable for maneuvering targets in low SNR environments.

5.4. Parameter Estimation Performance

In this section, we evaluate the parameter estimation performance of GRFT, FAF-LVD, PC-TRT, and the presented method, by Monte Carlo experiments. For the sake of comparison, the ranges of the estimated motion parameters by the four approaches are set to be identical. The input SNR after PC varies from −25 to 10 dB, in steps of 1 dB. For each SNR value, 500 independent experiments are carried out. The root mean square error (RMSE) is used as the evaluation criterion. The RMSE curves for velocity and acceleration are plotted in Figure 12.
The GRFT is demonstrated to have optimal parameter estimation performance, owing to the MLE. The presented method has a performance loss of about 3 dB compared to GRFT, due to the product operation. The PC-TRT has the worst RMSE, because it only uses one time slice for nonlinear correlation, which loses much signal energy. Although the FAF-LVD can reduce the energy loss through a variable latency, its RMSE is dramatically increased when the SNR is less than −5 dB, compared to the presented method. Thus, the presented approach acquires a better equilibrium between computational cost, detection performance, and parameter estimation.

6. Conclusions

In this paper, a new coherent integration method, based on PSZDCFT, is proposed for maneuvering targets, which effectively addresses the RCM and DFCM. For maneuvering targets with velocity ambiguity and excessive acceleration, the proposed method introduces a zoom operation into MDCFT, to extend its parameter estimation ranges, and it is interesting to note that this operation is also suitable for uniform motion targets with velocity ambiguity. Considering the coupling of range frequency and slow time, the scale operation is combined with ZDCFT to decouple, which not only obtains the correct parameter estimates, but also paves the way for the product operation, excluding false peaks. Naturally, the product operation is employed along the range frequency, which can inhibit false peaks and accumulate energy dispersed in different range cells. Finally, by constructing the phase compensation function with the estimated parameters, coherent integration is achieved. The computational complexity analysis shows that the computational cost of the proposed method is comparatively small, thanks to its fast implementation using the CZT-based SFT. Intensive numerical simulation results demonstrate that the proposed method is robust in a low SNR environment, due to the fact that the proposed method is a linear transform with no energy loss and no cross terms for multiple targets. Therefore, the proposed method obtains a superior equilibrium between computational complexity and detection performance, and is more suitable for maneuvering target detection in comparison with most current methods.

Author Contributions

Conceptualization, L.X. and H.G.; methodology, L.X. and H.G.; software, L.X. and T.L.; validation, L.X., T.L. and B.F.; formal analysis, L.X. and B.F.; resources, L.L.; writing—original draft preparation, L.X.; writing—review and editing, L.X. and H.G.; supervision, B.F. and L.L.; project administration, L.X., L.L. and T.L.; funding acquisition, H.G. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (under Grant 61671333), in part by the Natural Science Foundation of Hubei Province (under Grant 2014CFA093), in part by the Fundamental Research Funds for the Central Universities (under Grant 2042019K50264, Grant 2042019gf0013, and Grant 2042020gf0003), and in part by the Fundamental Research Funds for the Wuhan Maritime Communication Research Institute (under Grant 2017J-13 and Grant KCJJ2019-05).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the fact that it is currently privileged information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, M.; Huang, P.; Liu, X.; Liao, G.; Zhang, Z.; Wang, Z.; Hou, Q. Space Maneuvering Target Integration Detection and Parameter Estimation for a Spaceborne Radar System with Target Doppler Aliasing. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2020, 13, 3579–3594. [Google Scholar] [CrossRef]
  2. Wang, Z.; Shi, S.; Cheng, Z.; He, Z. A Modified Sequential Multiplexed Method for Detecting Airborne and Sea Targets with Over-the-Horizon Radar. IEEE Access 2020, 8, 84082–84092. [Google Scholar] [CrossRef]
  3. Huang, X.; Tang, S.; Zhang, L.; Li, S. Ground-Based Radar Detection for High-Speed Maneuvering Target via Fast Discrete Chirp-Fourier Transform. IEEE Access 2019, 7, 12097–12113. [Google Scholar] [CrossRef]
  4. Zhao, L.; Tao, H.; Chen, W.; Song, D. Maneuvering Target Detection Based on Subspace Subaperture Joint Coherent Integration. Remote Sens. 2021, 13, 1948. [Google Scholar] [CrossRef]
  5. Lin, L.; Sun, G.; Cheng, Z.; He, Z. Long-Time Coherent Integration for Maneuvering Target Detection Based on ITRT-MRFT. IEEE Sens. J. 2020, 20, 3718–3731. [Google Scholar] [CrossRef]
  6. Zheng, J.; Zhu, K.; Niu, Z.; Liu, H.; Liu, Q. Generalized Dechirp-Keystone Transform for Radar High-Speed Maneuvering Target Detection and Localization. Remote Sens. 2021, 13, 3367. [Google Scholar] [CrossRef]
  7. Jin, K.; Lai, T.; Wang, Y.; Li, G.; Liu, Q.; Zhao, Y. Coherent Integration for Radar High-Speed Maneuvering Target Based on Frequency-Domain Second-Order Phase Difference. Electronics 2019, 8, 287. [Google Scholar] [CrossRef] [Green Version]
  8. Suo, P.; Tao, S.; Tao, R.; Nan, Z. Detection of high-speed and accelerated target based on the linear frequency modulation radar. IET Radar Sonar Nav. 2014, 8, 37–47. [Google Scholar] [CrossRef]
  9. Zhang, Z.; Liu, N.; Hou, Y.; Zhang, S.; Zhang, L. A Coherent Integration Segment Searching Based GRT-GRFT Hybrid Integration Method for Arbitrary Fluctuating Target. Remote Sens. 2022, 14, 2695. [Google Scholar] [CrossRef]
  10. Jin, K.; Li, G.; Lai, T.; Jin, T.; Zhao, Y. A Novel Long-Time Coherent Integration Algorithm for Doppler-Ambiguous Radar Maneuvering Target Detection. IEEE Sens. J. 2020, 20, 9394–9407. [Google Scholar] [CrossRef]
  11. Tian, M.; Liao, G.; Zhu, S.; He, X.; Liu, Y.; Li, Y. An Efficient Method for Ground Maneuvering Target Refocusing and Motion Parameter Estimation Based on DPT–KT–MFP. Remote Sens. 2021, 13, 1092. [Google Scholar] [CrossRef]
  12. Huang, X.; Zhang, L.; Li, S.; Zhao, Y. Radar high-speed small target detection based on keystone transform and linear canonical transform. Digit. Signal Process. 2018, 82, 203–215. [Google Scholar] [CrossRef]
  13. Swerling, P. Radar probability of detection for some additional fluctuating target cases. IEEE Trans. Aerosp. Electron. Syst. 1997, 33, 698–709. [Google Scholar] [CrossRef]
  14. Tian, M.; Liao, G.; Zhu, S.; Liu, Y.; He, X.; Li, Y. Long-time coherent integration and motion parameters estimation of radar moving target with unknown entry/departure time based on SAF-WLVT. Digit. Signal Process. 2020, 107, 102854. [Google Scholar] [CrossRef]
  15. Kong, L.; Li, X.; Cui, G.; Yi, W.; Yang, Y. Coherent Integration Algorithm for a Maneuvering Target with High-Order Range Migration. IEEE Trans. Signal Process. 2015, 63, 4474–4486. [Google Scholar] [CrossRef]
  16. Zheng, J.; Zhang, J.; Xu, S.; Liu, H.; Liu, Q. Radar Detection and Motion Parameters Estimation of Maneuvering Target Based on the Extended Keystone Transform. IEEE Access 2018, 6, 76060–76074. [Google Scholar] [CrossRef]
  17. Perry, R.; DiPietro, R.; Fante, R. SAR imaging of moving targets. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 188–200. [Google Scholar] [CrossRef]
  18. Zhu, D.; Li, Y.; Zhu, Z. A keystone transform without interpolation for SAR ground moving-target imaging. IEEE Geosci. Remote Sens. Lett. 2007, 4, 18–22. [Google Scholar] [CrossRef]
  19. Xu, J.; Yu, J.; Peng, Y.; Xia, X. Radon-Fourier transform for radar target detection, I: Generalized Doppler filter bank. IEEE Trans. Aerospace Electr. Syst. 2011, 47, 1186–1202. [Google Scholar] [CrossRef]
  20. Xu, J.; Yu, J.; Peng, Y.; Xia, X. Radon-Fourier transform for radar target detection (II): Blind speed sidelobe suppression. IEEE Trans. Aerospace Electr. Syst. 2011, 47, 2473–2489. [Google Scholar] [CrossRef]
  21. Rao, X.; Tao, H.; Su, J.; Guo, X.; Zhang, J. Axis rotation MTD algorithm for weak target detection. Digit. Signal Process. 2014, 26, 81–86. [Google Scholar] [CrossRef]
  22. Zheng, J.; Su, T.; Zhu, W.; He, X.; Liu, Q. Radar high-speed target detection based on the scaled inverse Fourier transform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1108–1119. [Google Scholar] [CrossRef]
  23. Zheng, J.; Su, T.; Liu, H.; Liao, G.; Liu, Z.; Liu, Q. Radar high-speed target detection based on the frequency-domain deramp-keystone transform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 285–294. [Google Scholar] [CrossRef]
  24. Xu, J.; Xia, X.; Peng, S.; Yu, J.; Peng, Y.; Qian, L. Radar Maneuvering Target Motion Estimation Based on Generalized Radon-Fourier Transform. IEEE Trans. Signal Process. 2012, 60, 6190–6201. [Google Scholar]
  25. Li, X.; Cui, G.; Yi, W.; Kong, L. Coherent Integration for Maneuvering Target Detection Based on Radon-Lv’s Distribution. IEEE Signal Process. Lett. 2015, 22, 1467–1471. [Google Scholar] [CrossRef]
  26. Chen, X.; Guan, J.; Liu, N.; He, Y. Maneuvering Target Detection via Radon-Fractional Fourier Transform-Based Long-Time Coherent Integration. IEEE Trans. Signal Process. 2014, 62, 939–953. [Google Scholar] [CrossRef]
  27. Rao, X.; Tao, H.; Su, J.; Xie, J.; Zhang, X. Detection of constant radial acceleration weak target via IAR-FRFT. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 3242–3253. [Google Scholar] [CrossRef]
  28. Li, X.; Sun, Z.; Yi, W.; Cui, G.; Kong, L. Radar detection and parameter estimation of high-speed target based on MART-LVT. IEEE Sens. J. 2019, 19, 1478–1486. [Google Scholar] [CrossRef]
  29. Lin, L.; Cheng, Z.; He, Z. Coherent detection and parameter estimation for ground moving target based on MLRT-IDCFT. Digit. Signal Process. 2022, 120, 103259. [Google Scholar] [CrossRef]
  30. Li, X.; Cui, G.; Yi, W.; Kong, L. A fast maneuvering target motion parameters estimation algorithm based on ACCF. IEEE Signal Process Lett. 2015, 22, 270–274. [Google Scholar] [CrossRef]
  31. Zheng, J.; Liu, H.; Liu, J.; Du, X.; Liu, Q. Radar high-speed maneuvering target detection based on three-dimensional scaled transform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2821–2833. [Google Scholar] [CrossRef]
  32. Huang, P.; Liao, G.; Yang, Z.; Xia, X.; Ma, J.; Zhang, X. A fast SAR imaging method for ground moving target using a second-order WVD transform. IEEE Trans. Geosci. Remote Sens. 2016, 54, 1940–1956. [Google Scholar] [CrossRef]
  33. Jin, K.; Lai, T.; Zhu, S.; Li, G.; Jin, T.; Zhao, Y. Coherent Detection and Parameter Estimation for Radar High-Speed Maneuvering Target Based on FAF-LVD. Circuits Syst. Signal Process. 2020, 39, 2600–2622. [Google Scholar] [CrossRef]
  34. Huang, P.; Xia, X.; Liu, X.; Liao, G. Refocusing and Motion Parameter Estimation for Ground Moving Targets Based on Improved Axis Rotation-Time Reversal Transform. IEEE Trans. Comput. Imaging 2018, 4, 479–494. [Google Scholar] [CrossRef]
  35. Huang, P.; Liao, G.; Yang, Z.; Xia, X.; Ma, J.; Zhang, X. An approach for refocusing of ground moving target without target motion parameter estimation. IEEE Trans. Geosci. Remote Sens. 2017, 55, 336–350. [Google Scholar] [CrossRef]
  36. Li, X.; Huang, L.; Zhang, S.; Wang, W. Phase Compensation and Time-Reversal Transform for High-Order Maneuvering Target Detection. IEEE Geosci. Remote Sens. Lett. 2022, 19, 4025905. [Google Scholar] [CrossRef]
  37. Guo, X.; Sun, H.; Wang, S.; Liu, G. Comments on discrete chirp-Fourier transform and its application to chirp rate estimation. IEEE Trans. Signal Process. 2002, 50, 3115–3116. [Google Scholar]
  38. Xia, X. Discrete chirp-Fourier transform and its application to chirp rate estimation. IEEE Trans. Signal Process. 2000, 48, 3122–3133. [Google Scholar] [CrossRef] [Green Version]
  39. Wan, J.; Zhou, Y.; Zhang, L.; Chen, Z.; Yu, H. Efficient Algorithm for SAR Refocusing of Ground Fast-Maneuvering Targets. Remote Sens. 2019, 11, 2214. [Google Scholar] [CrossRef] [Green Version]
  40. Misiurewicz, J.; Kulpa, K.; Czekala, Z.; Filipek, T. Radar detection of helicopters with application of CNEAN method. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3525–3537. [Google Scholar] [CrossRef]
Figure 1. Simulation results of Example 1. (a) The result of MDCFT. (b) The result of ZDCFT ( α = 4 , β = 1 ).
Figure 1. Simulation results of Example 1. (a) The result of MDCFT. (b) The result of ZDCFT ( α = 4 , β = 1 ).
Remotesensing 15 01792 g001
Figure 2. Simulation results of Example 2. (a) The correct result of ZDCFT ( α = 4 , β = 2 ). (b) The blurred result of ZDCFT ( α = 1 , β = 2 ). (c) The invalid result of ZDCFT ( α = 4 , β = 1 ).
Figure 2. Simulation results of Example 2. (a) The correct result of ZDCFT ( α = 4 , β = 2 ). (b) The blurred result of ZDCFT ( α = 1 , β = 2 ). (c) The invalid result of ZDCFT ( α = 4 , β = 1 ).
Remotesensing 15 01792 g002
Figure 3. Simulation results of Example 3. (a) The result of SZDCFT ( ψ 1 = 1.01 ). (b) The result of SZDCFT ( ψ 2 = 0.99 ). (c) The result of ZDCFT ( ψ 1 = 1.01 ). (d) The result of ZDCFT ( ψ 2 = 0.99 ).
Figure 3. Simulation results of Example 3. (a) The result of SZDCFT ( ψ 1 = 1.01 ). (b) The result of SZDCFT ( ψ 2 = 0.99 ). (c) The result of ZDCFT ( ψ 1 = 1.01 ). (d) The result of ZDCFT ( ψ 2 = 0.99 ).
Remotesensing 15 01792 g003
Figure 4. The calculation process of CZT using FFT.
Figure 4. The calculation process of CZT using FFT.
Remotesensing 15 01792 g004
Figure 5. Computational complexity curves.
Figure 5. Computational complexity curves.
Remotesensing 15 01792 g005
Figure 6. The flowchart of the proposed method.
Figure 6. The flowchart of the proposed method.
Remotesensing 15 01792 g006
Figure 7. Coherent integration for a single target at SNR = 20 dB. (a) The result of PC. (b) The result of MDCFT at f = 4 MHz. (c) The result of SZDCFT at f = 4 MHz. (d) The result of SZDCFT at f = 4 MHz. (e) The result of PSZDCFT. (f) The result of coherent integration.
Figure 7. Coherent integration for a single target at SNR = 20 dB. (a) The result of PC. (b) The result of MDCFT at f = 4 MHz. (c) The result of SZDCFT at f = 4 MHz. (d) The result of SZDCFT at f = 4 MHz. (e) The result of PSZDCFT. (f) The result of coherent integration.
Remotesensing 15 01792 g007
Figure 8. The role of product operation in a low SNR environment (a) The result of PC. (b) The result of SZDCFT at f = 0 MHz. (c) The result of PSZDCFT. (d) The result of coherent integration.
Figure 8. The role of product operation in a low SNR environment (a) The result of PC. (b) The result of SZDCFT at f = 0 MHz. (c) The result of PSZDCFT. (d) The result of coherent integration.
Remotesensing 15 01792 g008
Figure 9. Coherent integration for multi-targets. (a) The result of PC. (b) The result of SZDCFT at f = 4 MHz. (c) The result of SZDCFT at f = 4 MHz. (d) The result of PSZDCFT. (e) The result of coherent integration for target A. (f) The result of coherent integration for target B. (g) The result of coherent integration for target C.
Figure 9. Coherent integration for multi-targets. (a) The result of PC. (b) The result of SZDCFT at f = 4 MHz. (c) The result of SZDCFT at f = 4 MHz. (d) The result of PSZDCFT. (e) The result of coherent integration for target A. (f) The result of coherent integration for target B. (g) The result of coherent integration for target C.
Remotesensing 15 01792 g009
Figure 10. Detection probability curves of different methods.
Figure 10. Detection probability curves of different methods.
Remotesensing 15 01792 g010
Figure 11. The results of parameter estimation for the proposed method and GRFT at SNR = −12 dB. (a) The result of of parameter estimation for the proposed method. (b) The result of of parameter estimation for GRFT.
Figure 11. The results of parameter estimation for the proposed method and GRFT at SNR = −12 dB. (a) The result of of parameter estimation for the proposed method. (b) The result of of parameter estimation for GRFT.
Remotesensing 15 01792 g011
Figure 12. Parameter estimation performance of different methods. (a) RMSE of velocity estimation. (b) RMSE of acceleration estimation.
Figure 12. Parameter estimation performance of different methods. (a) RMSE of velocity estimation. (b) RMSE of acceleration estimation.
Remotesensing 15 01792 g012
Table 1. Computational complexity comparison.
Table 1. Computational complexity comparison.
MethodsComputational Complexity
GRFT O ( N v N a M N )
RLVD O ( 3 N N v N a M 2 l o g 2 M )
ACCF O ( M N ( 2 l o g 2 N + l o g 2 M ) )
PC-TRT O ( N v M N l o g 2 ( M N ) )
FAF-LVD O ( 3 M 2 l o g 2 M + M N ( 7 l o g 2 M + 2 l o g 2 N ) )
Proposed O ( N M 2 ( 2 + 3 l o g 2 M ) )
Table 2. Simulation parameters of radar.
Table 2. Simulation parameters of radar.
Carrier frequency, f c 1 GHz
Bandwidth, B10 MHz
Sample frequency, f s 20 MHz
Pulse repetition frequency, P R F 500 Hz
Pulse duration, T p 20 μ s
Number of pulses, M1000
Table 3. Motion parameters of three maneuvering targets.
Table 3. Motion parameters of three maneuvering targets.
Motion ParametersTarget ATarget BTarget C
Initial range cell number451st476st501st
Radial velocity120 m/s−108 m/s117 m/s
Radial acceleration45 m/ s 2 −30 m/ s 2 42 m/ s 2
SNR (after PC)20 dB20 dB20 dB
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xia, L.; Gao, H.; Liang, L.; Lu, T.; Feng, B. Radar Maneuvering Target Detection Based on Product Scale Zoom Discrete Chirp Fourier Transform. Remote Sens. 2023, 15, 1792. https://doi.org/10.3390/rs15071792

AMA Style

Xia L, Gao H, Liang L, Lu T, Feng B. Radar Maneuvering Target Detection Based on Product Scale Zoom Discrete Chirp Fourier Transform. Remote Sensing. 2023; 15(7):1792. https://doi.org/10.3390/rs15071792

Chicago/Turabian Style

Xia, Lang, Huotao Gao, Lizheng Liang, Taoming Lu, and Boning Feng. 2023. "Radar Maneuvering Target Detection Based on Product Scale Zoom Discrete Chirp Fourier Transform" Remote Sensing 15, no. 7: 1792. https://doi.org/10.3390/rs15071792

APA Style

Xia, L., Gao, H., Liang, L., Lu, T., & Feng, B. (2023). Radar Maneuvering Target Detection Based on Product Scale Zoom Discrete Chirp Fourier Transform. Remote Sensing, 15(7), 1792. https://doi.org/10.3390/rs15071792

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop