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
where
and
represent the pulse width and rectangular function, respectively.
and
are the frequency modulated rate and fast time, respectively.
is the carrier frequency.
is the slow time.
and
M indicate the PRF and pulse number, respectively.
is the pulse repetition interval (PRI).
We consider a second-order motion model and assume that the instantaneous slant range
, between the target and the radar is
where
,
, and
are the initial range, radial velocity, and acceleration of the target, respectively.
Then, after the signal demodulation, the received baseband signal is stated as
where
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
where
and
B are the signal amplitude after PC and the signal bandwidth, respectively.
Then, performing FFT on Equation (
4) along the fast time
, yields
where
f is the range frequency corresponding to the fast time
, and
indicates the signal amplitude after performing FFT.
Substituting Equation (
2) into Equation (
5), one obtains
where
.
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
where
and
denote the signal amplitude and slow time, respectively.
and
are the centroid frequency and chirp rate, respectively. After discretization, Equation (
7) can be written as
where
and
represent the digital centroid frequency and chirp rate of the signal, respectively.
.
The MDCFT of
is defined as
where
and
denote the centroid frequency index and chirp frequency index, respectively. It is obvious that the estimation scopes of both
and
are
. With the substitution of Equation (
8) into Equation (
9), we can obtain
Depending on whether and exceed the estimation scopes, the peak position of can be divided into the following cases.
Case 1:
and
. This means that both
and
are within their estimation scopes, and as shown in Equation (
10), the peak of the two-dimensional spectrum can be observed at
.
Case 2:
and
. In this case,
is over its valid estimation range, which implies that Doppler ambiguity would occur.
can be written as
where
is the number of the Doppler ambiguity.
. Inserting Equation (
11) into Equation (
10), one obtains
Here, we use the constant equation
. As shown in Equation (
12), we can see that the peak is located at
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: . At this point, since is out of the estimation scope, and cannot be accurately estimated, regardless of the value of . In the two-dimensional spectrum, has no distinctive sharp peak.
Case 4:
. In this case,
is still out of the estimation scope, but similar to case 2, chirp rate ambiguity can occur.
can be written as
where
is the number of the chirp rate ambiguity.
. By taking Equation (
13) into Equation (
10), we have
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
cannot be obtained in this case. In practice, chirp rate ambiguity is almost impossible, that is,
. 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
is defined as follows
where
and
represent the centroid frequency index and chirp frequency index, respectively.
and
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
where
and
are the new implied periods of
in the
k and
l domains, respectively.
For case 1, i.e., and , we should just take . Then, the ZDCFT degenerates into the MDCFT. The peak of the theoretical two-dimensional spectrum is located at .
For case 2, i.e.,
and
, we know that
, so
is taken as 1. Since
exceeds the estimation scope, we should take the appropriate
to extend the estimation scope of
, i.e.,
In other words, it is to make
lie within the computational range of
k, i.e.,
. However, the ensuing problem is that
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.,
. Therefore, with
, Equation (
16) can be rewritten as follows
where
,
, and
. It can be seen from Equation (
18) that these
peaks are located at
When , we can obtain the correct digital centroid frequency and chirp rate, i.e., .
Example 1. Figure 1 illustrates the role of the zoom factor of the centroid frequency α. The simulation parameters of the LFM signal are: Hz/s, Hz, Hz, and . Then, we can get and . It can be seen that is out of the estimation scope while is not. Therefore, we choose and , 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 , we can also derive the theoretical x-coordinate of one of the false peaks according to Equation (19), i.e., . 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., . Figure 1b shows the result of ZDCFT, which is in agreement with the theory. For case 3, i.e.,
, as in case 2, we should select a suitable
(
) to make
lie within
. To put it another way, the estimation scope of
is extended by
, i.e.,
Note that, it is sufficient that
is less than
M according to Equation (
20) and that
is less than
. 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.,
. From another perspective, the computational range of
l is a subset of
, i.e.,
, which implies that
would generate only one peak in the
l domain, independent of
. Considering
, Equation (
16) is rewritten as follows
where
, and
. From Equation (
21), we can see that in the two-dimensional spectrum, there are still
peaks, which appear at
When , we can obtain the correct digital centroid frequency and chirp rate, i.e., .
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: Hz/s, Hz, Hz, and . Then, both and are out of the estimation scopes, i.e., and . 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 and and when and , respectively. It is evident that the former () is blurred in the Doppler domain, while the latter () 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 and [39]. The reason for this is that the spatial resolutions in the k and l domains are proportional to α and β, respectively, i.e., and . 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
, is as follows
where
denotes the complex amplitude of the signal.
and
represent the digital centroid frequency and chirp rate of the target, respectively.
is the scale factor.
Then, the SZDCFT of Equation (
23) is defined as
where
and
denote the implied periods of
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
In Equation (
25), the true peaks will appear at the location of
when
, 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, and . 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 and , 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 could make the implied periods of the SZDCFT different from those of the ZDCFT, fortunately, for the narrowband radar system, there are and . 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 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
where
and
are the ambiguity integers and blind speed, respectively.
is the unambiguous velocity.
When the target’s velocity is not ambiguous, i.e., and , then we have . Evidently, is taken as 1.
When the target’s velocity is ambiguous, i.e.,
, we have
, where
. 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
. Then, we can obtain the maximum value of
as
. According to Equation (
25), in order to ensure that
can be estimated,
should satisfy the following inequality
where
is the ceiling operation. In order to have a better estimation accuracy of the velocity, we generally take
.
Similar to the analysis of velocity, for the acceleration of the maneuvering target
, we assume that the acceleration scope of interest is
. Then, we can also obtain the maximum value of
, i.e.,
. Similarly, based on Equation (
25), we can also obtain the constraint of
as follows
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 . 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,
has different implied periods in the
k domain (
), 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.,
, 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
where
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
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.,
Finally, we can achieve the coherent integration by applying the IFFT and FFT with respect to
f and
.
where
and
denote the amplitude and coherent processing interval, respectively.
is the folded Doppler frequency with respect to
.
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
can be rewritten as follows
Equation (
33) shows that, for each fixed
l,
is the scaled Fourier transform (SFT) of the signal
. As we know, the SFT can be efficiently implemented by FFT-based chirp-z transform (CZT).
Assuming a discrete signal
,
, the SFT of
is defined as follows
where
and
.
is the scale factor. Then, the Bluestein equation is applied here, i.e.,
Inserting Equation (
35) into Equation (
34), one obtains
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
.
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, , and 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 .
The RLVD firstly obtains the target trajectory through a three-dimensional search, which requires a computational cost of , and then accomplishes coherent integration via LVD, with a computational cost of . Therefore, its total computational complexity is .
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 . Then, the computational cost required to achieve coherent integration is . Hence, its total computational cost is about .
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 .
The main calculations of FAF-LVD contain second-order KT , FAF , CZT-based SFT , phase compensation , and LVD . Therefore, its total computational burden is about .
For each fixed range cell of the proposed method, the computational costs required by the SZDCFT operation and the product operation are and , respectively. Therefore, the overall computational load of the PSZDCFT, is in the order of .
Table 1 shows the computational complexities of the above methods. We assume that
. 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.,
, 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 , then employ the FFT along the range dimension, to obtain .
Step 2: For a fixed range frequency f, apply the SZDCFT on , to obtain .
Step 3: Determine the scope of the next range frequency, f. If , execute step 2. If , execute product operation, to obtain .
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.
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.