1. Introduction
With the development of stealth technology, modern ground and aerial targets are able to present the characteristics of high speed, strong maneuverability, long range, and low radar cross-section (RCS) [
1]. The traditional moving target detection (MTD) algorithms do not have enough capacity to detect these high-speed maneuvering weak targets. The long-time integration technique is an effective way to increase the signal-to-noise ratio (SNR) and improve radar detection performance [
2]. However, when the integration time increases, the high speed and acceleration of the target cause range migration (RM) and Doppler frequency migration (DFM) [
3], which limit the performance of classical integration algorithms [
4]. Therefore, new approaches to eliminate RM and DFM have been investigated [
5].
The typical long-time integration techniques are mainly divided into two categories: incoherent integration and coherent integration. The incoherent integration methods only use the amplitude of echoes to accumulate target signal, which causes poor detection performance in low SNR scenarios [
6]. Classical incoherent integration methods include the Hough transform, Radon transform, dynamic programming, and particle filtering methods [
7].
Coherent integration performs better than incoherent integration by compensating the phase fluctuation among different sampling pulses. The keystone transform (KT) method corrects the range walk by rescaling the slow time for each range frequency [
8]. KT and several improved versions of KT have been widely used, as they can correct the RM effectively without any prior knowledge about the target motion parameter. The Radon Fourier transform (RFT) method eliminates the RM via joint searching along the range and velocity directions of the moving target and integrates coherently via Doppler filtering [
9]. The axis rotation MTD (AR-MTD) method eliminates the linear RM by rotating the two-dimensional echoes data plane and realizes coherent integration via the classical MTD algorithm [
10]. A fast coherent integration method based on sequence reversing transform was presented in [
11], providing a good balance between computational cost and detection ability. By employing the symmetric autocorrelation function and the scaled inverse Fourier transform (SCIFT), a coherent detection algorithm was introduced in [
12]; this algorithm can detect high-speed targets without brute-force searching of unknown motion parameters. However, the above methods only consider RM correction, and suffer performance degradation when DFM appears.
KT and second-order KT (SOKT) [
13]-based methods have been proposed to regulate both RM and DFM. The KT-matched filtering processing (KT-MFP) method corrects the linear RM via KT and then jointly searches in the fold factor and acceleration domain to remove the residual RM and compensate the DFM. KT-MFP achieves coherent integration through the slow-time Fourier transform (FT) [
14]. The SKT-DLVD method uses the segmented keystone transform (SKT) to correct the range walks of targets and then applies the Doppler Lv’s transform (DLVD) to estimate the velocities of targets [
15]. The KT-LCT method employs KT to eliminate RM. After this, the linear canonical transform (LCT) is applied to compensate DFM and realize coherent integration [
16]. In the low-frequency ultra-wideband synthetic aperture radar (SAR), [
17] utilizes the first-order KT to correct the range walk and then uses the SOKT to compensate for the range curvature. However, when Doppler ambiguity occurs due to the high speed of the target or limited pulse repetition frequency (PRF), the Doppler ambiguity number has to be estimated before KT and SOKT processing, which increases the computational burden.
In order to integrate coherently under the condition of Doppler ambiguity, the SOKT-RFT method utilizes SOKT for range curvature correction and the improved de-chirping method for DFM compensation. Then, RFT is applied to correct the linear RM. Because SOKT-RFT eliminates RM and DFM in steps, the integration performance can deteriorate due to compensation errors in the previous steps [
18]. The SOKT-MFrRT method uses the SOKT to eliminate quadratic range cell migration and the modified fractional Radon transform (MFrRT) to estimate the ambiguity number of Doppler frequency [
19].
To eliminate the RM and DFM effects simultaneously, the Radon-fractional Fourier transform (RFRFT) removes the RM effect via three-dimensional searching within the parameter space and realizes integration using the fractional Fourier transform (FRFT) [
20]. Inspired by this, the Radon-Lv’s distribution (RLVD) method eliminates the RM via jointly searching in the target motion parameter space and achieves coherent integration via Lv’s distribution. The RLVD obtains better integration and detection performance than RFRFT [
21]. The computational complexity of RLVD is quite large due to the need for multi-dimensional joint searching [
22]. In this regard, IAR-FRFT eliminates the linear RM via the improved axis rotation (IAR) transform and realizes coherent integration by FRFT [
23]. An approach combining the modified axis rotation transform (MART) and Lv’s transform (LVT) is presented in [
24]. Compared with RFRFT and RLVD, the computational cost of IAR-FRFT and MART-LVT is decreased; nevertheless, these two methods suffer integration loss because the range curvature induced by the target acceleration is ignored.
For cases of increasing target maneuverability and long observation time, the short-time generalized radon Fourier transform (STGRFT) method was presented in [
25] to detect a maneuvering weak target with multiple motion models. The STGRFT was able to eliminate RM and DFM as well as to estimate the model changing-point time and accumulate the target energy distributed in different motion stages. In the wideband radar scenario, a coherent integration algorithm based on the sub-band keystone transform and extended Lv’s distribution (ELVD) was proposed to estimate the motion parameters and reconstruct the high-resolution range profile (HRRP) of a maneuvering weak target [
26]. In [
27], high-order motion parameter estimation was modeled as an under-estimated linear regression and the complex-field Bayesian compression sensing (BCS) algorithm was designed to resolve the sparse recovery.
In this paper, a low computational complexity coherent integration algorithm named SOKT-IAR-LVD is proposed to eliminate the RM and DFM caused by the constant radial acceleration of the target. First, the SOKT-IAR-LVD employs SOKT to eliminate the range curvature caused by target acceleration and alleviate linear range migration. Second, the IAR is applied to regulate linear range migration by rotating the fast time axis and the target envelope is aligned along the slow time axis with a quadratic phase characteristic. Third, based on the quadratic phase characteristic of the target signal, the LVD is adopted to accumulate the target signal into a well-focused peak in the centroid frequency–chirp rate (CFCR) domain. The motion parameters of the target are estimated by the rotation angle of IAR and the peak position of LVD outputs. The coherent integration gain and computational complexity of SOKT-IAR-LVD are analyzed.
Without needing to estimate the Doppler ambiguity number and target acceleration, SOKT-IAR-LVD possesses a much lower computational cost than RLVD and achieves better integration performance than the IAR-FRFT, SOKT-RFT, and AR-MTD methods.
The rest of this paper is organized as follows. In
Section 2, the echo of a constant radial acceleration target is modeled and the phase terms which cause the RM and DFM are analyzed. In
Section 3, the procedure of the proposed SOKT-IAR-LVD method is detailed, and the integration gain and computational complexity are analyzed.
Section 4 presents simulation results that demonstrate the efficiency of the SOKT-IAR-LVD method.
2. Signal Model
Assume that the radar transmitted waveform is a linear frequency modulated (LFM) signal, as follows:
where
is the carrier frequency,
is the pulse repetition interval (PRI),
is the pulse width,
B is the bandwidth of the LFM waveform,
is the frequency modulation rate,
and
denote the fast time and slow time, respectively,
M is the number of transmitted LFM pulses in a coherent processing interval (CPI), and
For a target moving towards the radar, the instantaneous range between radar and target is
where
is the initial range between the radar and target and
v and
a represent the radial velocity and acceleration of the target, respectively. For simplicity, we define
Because
, the baseband target echo can be written as
where
c is the speed of light,
is the Doppler frequency, and
is the carrier wavelength. The matched filtered output of the received target echo is
where sinc
denotes the
function. In (
5), the peak position of the target envelope in the
mth PRI is
which varies with the increase of slow time. Denote
as the range offset during a CPI, which has
When
is less than half of the range resolution, that is,
the peak position of the target envelope can be approximately regarded as in the same range cell during the integration time, meaning that the traditional MTD method can be employed for coherent integration. However, the condition of (
8) is frequently not met for high-speed maneuvering targets, and the integration performance deteriorates when using the MTD method.
The two range migration terms at the right side of (
7) should be eliminated before coherent integration. Applying FT on (
5) to the fast time
, the matched filtered output in the range frequency domain is
Because of the high target speed or low radar PRF, under-sampling induces ambiguity of the measured target velocity. The velocity of the target is written as
where
is the Doppler ambiguity number,
is the blind velocity,
is the PRF, and
is the measured velocity and satisfies
. Because
, we can ignore the amount of
in
and write (
9) as
where
. Thus, it has
There are six exponential terms in (
12):
is the initial range term;
is the Doppler term determined by the measured velocity
, and results in a linear range migration;
is the phase term induced by the blind velocity, and also causes a linear range migration;
is the frequency modulation term induced by the target acceleration, and leads to the range curvature and DFM;
indicates the characteristic of the high-speed motion and causes a fixed target envelope offset from its nominal position; and
is a constant term.
Here, and suffer from the first-order coupling between f and m, suffers from the second-order coupling between f and m, which causes both range curvature and DFM, and , , and should be decoupled prior to coherent integration.
3. SOKT-IAR-LVD Method
In this section, the SOKT-IAR-LVD method to improve coherent integration performance is detailed. After matched filtering in the range frequency domain, SOKT is employed to eliminate the second-order coupling in and correct the range curvature caused by the target acceleration. Then, the IAR is used to eliminate the coupling between range frequency and slow time in the terms and . At last, the LVD is applied to coherently integrate the target echo. The target motion parameters are estimated by the IAR and LVD results. The details of the SOKT-IAR-LVD method are as follows.
3.1. Range Curvature Correction via SOKT
SOKT is utilized to correct the range curvature caused by the target acceleration. SOKT is a process of rescaling the slow time axis for each range frequency. The scaling formula of SOKT is defined as
where
denotes the new slow-time variable.
Substituting (
13) into (
12), the SOKT output in the fast time frequency domain is
As shown in (
14), SOKT removes the second-order coupling in
for range curvature elimination after matched filtering. As
, the following approximations hold:
Therefore, (
14) is rewritten as
As shown in (
16), the second-order coupling between
f and
m in (
12) has been eliminated. After performing the inverse Fourier transform (IFT) on (
16) into the
domain, we have
The set
(
17) can now be rewritten as
where
r represents the range corresponding to the fast time
. In (
18), the range offset of the target envelope varies linearly with the slow time
. In addition, the linear RM caused by the measured velocity
is reduced to half of its value through SOKT processing.
With the fast time sampling frequency
, the discrete form of (
18) is
where
,
, and
and
represent the range cell numbers of
r and
, respectively.
3.2. Range Migration Correction via IAR
The AR-MTD method concentrates the target echoes in a range cell via the axis rotation transform. However, the Doppler resolution may vary with the axis rotation angle. In the SOKT-IAR-LVD method, the IAR regulates the linear range migration in (
19) by rotating the fast time axis. The slow time axis remains unchanged in order to maintain a constant Doppler resolution.
Figure 1 shows the diagram of the IAR transform. In (
19), the target signal envelope after SOKT is distributed along a straight line with a slope of
in the coordinate system
. In
Figure 1, the angle between the target signal envelope and slow time axis
is defined as
, and has
The IAR processing is
where
represents a new coordinate system after axis rotation,
denotes the axis rotation angle, and
is the angle rotation step. We denote
as the number of searching angles. Because the target velocity is frequently limited to a certain range, the angle rotation searching area
can be predetermined in order to reduce the computational burden when the target’s velocity is limited in the region
with
According to (
20), the angle rotation searching area
is computed as follows:
Substituting (
20) and (
21) into (
19) gives
When
, (
26) is written as
As shown in (
27), the peak position of the target envelope is concentrated at the range cell
. The set
(
27) is simplified as follows:
with
In (
28), the linear range migration has been corrected by IAR and the output of IAR is a chirp signal with a chirp rate of
. Inspired by the excellent performance of LVD in extracting the parameters of chirp signals [
28], the LVD method is adopted to integrate the target signal in (
28) and then estimate the target motion parameters.
3.3. Coherent Integration and Parameter Estimation with LVD
In LVD processing, the symmetric instantaneous autocorrelation function (SIAF) of (
28) is
where ∗ denotes the complex conjugation. Because
when
or
, the number of valid elements in the SIAF matrix is
.
The variables
and
l in (
30) are coupled with each other in the exponential phase term. We set
where
is the LVD slow time variable and
h is a scaling factor that determines the chirp rate estimation range of
. Substituting (
31) into (
30), we have
where the coupling between
and
l has been eliminated. Performing two-dimensional (2D) FT on (
32), the output of LVD is
with
where
represents the centroid frequency–chirp rate (CFCR) domain with
. The target echo is coherently accumulated at
,
. Therefore, the velocity and acceleration of target are estimated by the peak position of the LVD results:
According to (
37), the scaling factor
h should be no less than
, where
denotes the maximal target radial acceleration.
Finally, based on the IAR and LVD results, the Doppler ambiguity number
and target velocity
v are estimated as follows:
3.4. Procedure of the SOKT-IAR-LVD Method
Based on the above analysis, the flow chart of the SOKT-IAR-LVD method is shown in
Figure 2. The procedure of SOKT-IAR-LVD is divided into the following steps.
Step 1. Initialize the parameters of the SOKT-IAR-LVD method, including the rotation angle searching area , angle search step , and scaling factor h.
Step 2. Perform matched filtering on the radar echoes in the frequency domain.
Step 3. Apply SOKT to the matched filtered outputs to eliminate the range curvature, then perform range IFT on the SOKT outputs.
Step 4. For each , apply IAR to the SOKT outputs to remove the residual linear RM, then perform LVD transform on the IAR results along the slow time.
Step 5. After computing the LVD results for all rotation angles in , search for the peak of the LVD results at each range cell and the corresponding values of , , and .
Step 6. Apply constant false alarm rate (CFAR) detection to the peak of the LVD results. If a target is detected, the range, velocity, and acceleration of the target are estimated using the values of
,
, and
using (
35)–(
39), respectively.
3.5. Integration Gain Analysis
The coherent integration gain of SOKT-IAR-LVD mainly depends on the LVD transform. We denote the output SNR of matched filtering as
, and the output SNR of the IAR transform is the same as
. Adding the complex Gaussian noise to the IAR output in (
28) with
, we have
where
represents the noise with power
and the amplitude
is
. The SIAF of
is
In (
41), the amplitude of the target signal changes to
, while the noise power is
. Because the number of valid elements in the SIAF matrix is
, the output SNR of LVD is
Therefore, the integration gain of LVD is
The integration gain of the RFT- and FRFT-based methods only depends on the pulse number
M, which is
when the output SNR of matched filtering satisfies
The SOKT-IAR-LVD method achieves a higher integration gain than the RFT and FRFT based methods. Condition (
45) frequently holds for long coherent integration scenarios.
3.6. Computation Complexity Analysis
Next, the computational load of SOKT-IAR-LVD is analyzed. Because SOKT and LVD can be realized by Chirp-z transform–inverse fast FT (CZT-IFFT) [
29] and scaled FT–IFFT (SFT-IFFT) [
30], respectively, the computational complexities of SOKT and IAR are
and
, where
is the number of range cells in a PRI. The computation complexity of LVD at a rotation angle
is
. Therefore, the overall computational cost of SOKT-IAR-LVD is approximately
.