Next Article in Journal
Three-Stage Operational Amplifier with Frequency Compensation Using Cascade Zero
Next Article in Special Issue
Domain Adaptive Few-Shot Learning for ISAR Aircraft Recognition with Transferred Attention and Weighting Importance
Previous Article in Journal
On Dynamic Node Cooperation Strategy Design for Energy Efficiency in Hierarchical Federated Learning
Previous Article in Special Issue
Vision-Based Support for the Detection and Recognition of Drones with Small Radar Cross Sections
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Coherent Integration Algorithm for Maneuvering Target Detection Based on Symmetric Instantaneous Autocorrelation Function

1
CAS Key Laboratory of Microwave Remote Sensing, National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(11), 2363; https://doi.org/10.3390/electronics12112363
Submission received: 1 April 2023 / Revised: 20 May 2023 / Accepted: 22 May 2023 / Published: 23 May 2023

Abstract

:
Detection and parameter estimation of maneuvering targets having a jerking motion are some of the challenging problems for modern radar systems. Such targets usually introduce range migration (RM) and Doppler frequency migration (DFM) problems leading to serious performance degradation in detection. To address these problems, a novel coherent integration (CI) algorithm is proposed based on a new symmetric instantaneous autocorrelation function (NSIAF), which can be utilized to reduce the order on the slow time and to eliminate the linear range migration (LRM) first. Then, the jerk and acceleration of the target are estimated after applying the keystone transform (KT) and the scaled Fourier transform (SFT); both of these are then used to construct the reference function for matched filtering. Finally, CI and target detection can be accomplished by the scaled inverse Fourier transform (SCIFT) after matched filtering. Both simulation data (this work) and practical radar experiment data (data set of others) were processed to validate the proposed algorithm. Compared with other representative algorithms, our algorithm can achieve a good balance between computational complexity and detection performance.

1. Introduction

The method and algorithm for coherent integration (CI) of echoes from highly maneuvering targets (MaTs), along with motion parameter estimation, have been attracting growing attention, and significant efforts have made in research as part of the rapid development of modern radars [1,2,3,4,5]. Due to the feasible low-radar cross-section of maneuvering targets, the returned radar signal can be quite weak, which can make the detection more difficult. With the help of long-time CI techniques, the MaT detection performance of radar can be significantly improved [6,7,8,9]. However, target maneuvering can cause high-order range migration (RM) and Doppler frequency migration (DFM), which creates considerable challenges for CI.
In order to address the linear range migration (LRM) caused by the target velocity, several algorithms have been developed, e.g., the keystone transform (KT) can correct the LRM without requiring knowledge about the target’s velocity beforehand [10,11]. However, its performance might be limited due to velocity ambiguity resulting from the resampling process. Additionally, the interpolation operation can also affect the detection performance while resulting in a significant computational burden. The Radon–Fourier transform (RFT) can remove the LRM and realize CI by jointly searching in the target’s motion parameter space. The axis-rotation moving-target detection (AR-MTD) [12] corrects the range migration (RM) by rotating the fast-time and slow-time axes at the same time and completes the CI by moving target detection (MTD). Both the RFT and the AR-MTD suffer from heavy computational burdens because of motion parameter searching. The scaled inverse Fourier transform (SCIFT) [13] can eliminate the LRM and accomplish the CI without any searching procedure. However, the aforementioned algorithms cannot deal with a target experiencing acceleration or a jerking motion.
The acceleration of MaTs will induce quadratic range migration (QRM) and linear Doppler frequency migration (LDFM). The Radon-fractional Fourier transform (RFRFT) [14] was proposed to deal with MaTs experiencing acceleration. The RFRFT can extract the target trajectory and eliminate the LDFM by applying the Radon transform and the fractional Fourier transform together. Nonetheless, the RFRFT is hindered by high computational cost resulting from the three-dimensional searching. The second-order KT (SOKT) [15], and the KT-Lv’s distribution (KT-LVD) [16] algorithms were developed to reduce the dimension of parameter searching. However, these algorithms still require one-dimensional searching. The discrete polynomial-phase transform (DPT) [17], although it can achieve fast CI without any searching, is not suitable for low signal-to-noise-ratio (SNR) scenarios due to the cross-correlation applied. Additionally, all the above algorithms fail to deal with a maneuvering target experiencing a jerking motion because cubic range migration (CRM) and quadratic Doppler frequency migration (QDFM) are induced.
Although the generalized Radon–Fourier transform (GRFT) [18] is an effective algorithm to accomplish CI regarding acceleration, the brute-force searching of motion parameters significantly increases the computation complexity. In addition, the blind speed side lobe (BSSL) effect may appear due to discrete pulse sampling, finite range resolution, and limited integration time, leading to serious false alarms and CI degradation. The adjacent cross-correlation function (ACCF) applied iteratively [19] is an efficient CI algorithm for MaTs, but the detection performance is sacrificed since the multiple nonlinear operations are conducted with a fixed lag time. The multiple nonlinear operations using variable lag time were studied in [20], demonstrating higher anti-noise performance than the ACCF-based algorithms. Unfortunately, this algorithm is not suitable for maneuvering targets with a jerking motion.
In this paper, we proposed a novel fast CI algorithm for the detection of MaTs with a jerking motion. We first define a new symmetric instantaneous autocorrelation function (NSIAF) with variable lag time, which can be applied to reduce the order on the slow time and eliminate the LRM. Afterwards, the KT and the scaled Fourier transform (SFT) [21] are applied to decouple the range frequency, lag time, and slow time, and obtain the acceleration and jerk estimations, according to which the matched function is constructed and applied. The CI can then be accomplished using SCIFT [13], and in the final step, the target is detected through peak searching. The proposed algorithm can be implemented with a low computational cost by circumventing parameter searching. Comparisons with other representative algorithms indicate that this proposed algorithm can achieve a favorable balance between the computational complexity and the detection performance. Finally, both simulated data and practical radar data are processed to demonstrate the effectiveness of the proposed algorithm.
The remainder of this paper is organized as follows. Section 2 provides a brief introduction to the signal model. In Section 3, the proposed algorithm is described in detail. In Section 4, both simulated and practical radar data are processed to validate the proposed algorithm. Finally, the conclusion is drawn in Section 5.

2. Signal Model

Suppose that a pulse-Doppler radar transmits a linear frequency modulation (LFM) signal as
s ( t , t n ) = rect t T p exp ( j 2 π f c t + j π K t 2 )
where  rect  denotes a rectangle function with pulse duration of  T p K = B / T p  is the frequency modulation rate,  B  is the signal bandwidth,  f c  is the carrier frequency, and  t  denotes the fast time.
In this work, we focus on a maneuvering target with a jerking motion, so the instantaneous distance of the target relative to the radar can be expressed as [6,7,22]
R ( t n ) = R 0 + v t n + a t n 2 2 + b t n 3 6
where  t n = n T r   n = 1 , 2 ,   , N  is the slow time,  N  and  T r  are respectively the pulse number and the pulse repetition interval,  R 0 v a , and  b  denote the initial radial range, radial velocity, radial acceleration, and radial jerk.
The received echo signal after down conversion can be expressed as
s r ( t , t n ) = A 0 rect t 2 R ( t n ) / c T p exp j π K t 2 R ( t n ) c 2 exp j 4 π R ( t n ) λ
where  A 0  is the complex amplitude,  c  is the light speed,  f c  is the carrier frequency of the transmitted signal, and  λ = c / f c  is the wavelength.
The signal after pulse compression can be expressed as
s c ( t , t n ) = A 1 sin c B t 2 ( R 0 + v t n + a t n 2 / 2 + b t n 3 / 6 ) c × exp j 4 π R 0 + v t n + a t n 2 / 2 + b t n 3 / 6 λ
where  A 1  denotes the signal amplitude after pulse compression and  sinc x = sin π x / π x  denotes the SINC function.
As one can see from (4), the peak position and the phase of the range compressed signal are time varied with slow time  t n . Once the position variation exceeds the range resolution  ρ r = c / 2 B  in the SINC function, the compressed echo will be situated at cells having different ranges, i.e., the terms  t n t n 2 , and  t n 3  induce LRM, QRM, and CRM, respectively.
From the exponential term in (4), the instantaneous Doppler frequency can be derived as
f d ( t n ) = 2 λ d ( R 0 + v t n + a t n 2 / 2 + b t n 3 / 6 ) d t n = 2 v λ + 2 a t n λ + b t n 2 λ
As clearly shown in (5), the Doppler frequency varies with the slow time  t n  due to the radial acceleration  a  and jerk  b , which respectively induce the LDFM and QDFM.
The LRM, QRM, CRM, LDFM, and QDFM all can lead to serious CI loss and should be corrected to achieve good performance in MaT detection.

3. Novel Coherent Integration Algorithm for Maneuvering Target Detection

3.1. Principle of the Proposed Algorithm

As we know, the result of (4) is obtained by performing inverse fast Fourier transform (IFFT) on the following signal with respect to the range frequency  f r :
S c ( f r , t n ) = A 2 rect f r B exp j 4 π ( f r + f c ) R 0 c × exp j 4 π ( f r + f c ) v t n + a t n 2 / 2 + b t n 3 / 6 c
where  A 2  denotes the signal amplitude after FFT.
In the following, the NSIAF is defined with variable lag time regarding  S c ( f r , t n )  in the  ( f r , t n )  domain as
R ( f r , t n , τ ) = S f r , t n + τ 2 + τ 0 S f r , t n τ 2 τ 0 × S f r , t n + τ 2 τ 0 S f r , t n τ 2 + τ 0
where “ * ” denotes the complex conjugation operation, and  τ  and  τ 0  are the variable lag time and a constant time delay, respectively. Here  τ 0  is set as  0.089 T c  [23,24],  T c = N T r  is the total integration time.
By substituting (6) into (7), we obtain
R f r , t n , τ = A 3 exp j 8 π f r + f c c τ 0 τ a + b t n
where  A 3 = A 2 4 rect f r / B  is the signal amplitude of NSIAF.
Equation (8) shows the LRM has been eliminated and the highest order of the term  t n  has been reduced from three to one. However, the coupling between  f r τ , and  t n  still exists, which can be decoupled by the KT [10,11], via variable substitution as follows:
f c + f r τ = f c τ n
where  τ n  is the new lag-time variable.
Then, the signal becomes
R K T f r , t n , τ n = A 3 exp j 8 π f c c τ n τ 0 a + b t n exp j 2 π f r 0
As (10) shows,  f r  has been decoupled from  τ  and  t n , and the echo energies of pulses are all distributed along the  f r  axis [15], which can be accumulated directly via summation:
Q t n , τ n = SUM f r R K T f r , t n , τ n = A 4 exp j 8 π f c c a τ n τ 0 exp j 8 π f c c b τ n τ 0 t n
where  SUM f r  denotes the addition operation along the  f r  axis, and  A 4  is the signal amplitude after the addition operation.
In the following, we decouple  τ n  and  t n  by conducting SFT on (11) as follows:
Q f n , τ n = SFT t n Q t n , τ n = t n Q t n , τ n exp j 2 π τ n f n t n d t n = A 5 exp j 8 π f c c τ n τ 0 a δ f n + 4 λ τ 0 b
where  A 5  denotes the signal amplitude after SFT,  f n  is the scaled frequency variable corresponding to  t n , and  SFT t n  denotes the SFT operation with respect to  t n .
One can see from (12) that  t n  has been decoupled from  τ n . By applying the FFT on (12) with respect to  τ n , we have
Q f n , f τ n = A 6 δ f n + 4 λ τ 0 b δ f τ n + 4 λ τ 0 a
where  f τ n  is the frequency variable corresponding to  τ n , and  A 6  denotes the signal amplitude after FFT.
According to (13), the acceleration  a  and the jerk  b  can be estimated by searching for the peak position of  4 a τ 0 / λ , 4 b τ 0 / λ  in the  f n f τ n  domain, as
a ^ = λ f ^ τ n 4 τ 0 b ^ = λ f ^ n 4 τ 0
where  f ^ n  and  f ^ τ n  are the peak values along the  f n  and  f τ n  axes, respectively, and  a ^  and  b ^  are the estimated  a  and  b , respectively.
Based on (14), we can construct the corresponding matched function used to compensate for the acceleration and jerk as follows:
H c f r , t n = exp j 4 π f r + f c a ^ t n 2 / 2 + b ^ t n 3 / 6 c
By multiplying (6) with (15), we get
S 1 f r , t n = S c f r , t n H c f r , t n = A 2 rect f r B exp j 4 π ( f r + f c ) v t n c
As shown in (16) the QRM, LDFM, and QDFM have all been eliminated; LRM has not been eliminated, but can be easily eliminated by SCIFT [13]. So, the final CI result can be achieved:
S 2 f r , f d = SCIFT S 1 f r , t n = A 7 sin c B t 2 R 0 c δ f d + 2 v λ
where  SCIFT  denotes the SCIFT operation.
The whole flowchart of the proposed algorithm is given in Figure 1, which is explained in detail as follows:
Step 1: Perform the pulse compression and range FFT to obtain  S c ( f r , t n )  in (6).
Step 2: Calculate the NSIAF based on (7), and denote the result as  R f r , t n , τ .
Step 3: Conduct KT on  R f r , t n , τ , followed by the addition operation along the  f r  axis to get  Q t n , τ n .
Step 4: Apply the SFT and FFT on  Q t n , τ n  and estimate the acceleration  a ^  and the jerk  b ^  through searching the peak position.
Step 5: Construct the matched function based on (15) and multiply it with  S c ( f r , t n )  to get  S 1 f r , t n .
Step 6: Perform the SCIFT on  S 1 f r , t n  to realize coherent integration.
Step 7: Carry out the constant false alarm (CFAR) detection to accomplish target detection.

3.2. Cross-Term Analysis

The above analysis is valid only for a single target scenario. Cross-terms can appear when multiple targets are considered, and these shall be analyzed in the following.
Similar to (6), the received signal backscattered from  K  targets after pulse compression in the range frequency domain can be expressed as
S c ( f r , t n ) = k = 1 K A 2 k rect f r B exp j 4 π ( f r + f c ) R 0 c × exp j 4 π ( f r + f c ) v k t n + a k t n 2 / 2 + b k t n 3 / 6 c
By performing NSIAF on (18), we obtain
R f r , t n , τ = k = 1 K A 3 k exp j 8 π f r + f c c τ 0 τ ( a k + b k t n ) self - term + i = 1 K 1 j = i + 1 K R cross , i , j f r , t n , τ
where  R cross , i , j f r , t n , τ  denotes the cross term generated by the  i th  and  j th  target.
After the KT of (9) and addition operation of (11) applied to (19), we obtain
Q t n , τ n = k = 1 K A 4 k exp j 8 π f c c a k τ n τ 0 exp j 8 π f c c b k τ n τ 0 t n self - term + i = 1 K 1 j = i + 1 K Q cross , i , j t n , τ n
where  Q cross , i , j t n , τ n  is the outcome of  R cross , i , j f r , t n , τ  after KT and addition operations.
After performing SFT and FFT on (20) with respect to  t n  and  τ n , respectively, we have
Q f n , f τ n = k = 1 K A 6 k δ f n + 4 λ τ 0 b k δ f τ n + 4 λ τ 0 a k + i = 1 K 1 j = i + 1 K Q cross , i , j f n , f τ n
where  Q cross , i , j f n , f τ n  denotes the cross-term of  Q f n , f τ n . The detailed expressions of  R cross , i , j f r , t n , τ Q cross , i , j t n , τ n , and  Q cross , i , j f n , f τ n  can be found in Appendix A.
As (21) shows, all the  K  maneuvering targets can be well focused in the  f n f τ n  domain, whose peak locations correspond to  4 a k τ 0 / λ , 4 b k τ 0 / λ , while, as shown in Appendix A, the cross-terms are dispersed in the  f n f τ n  domain, so the cross-terms’ energies are much smaller than those of the self-terms.
Additionally, we should point out that the proposed algorithm may suffer from performance degradation in scenarios when both stronger and weaker maneuvering targets exist since the cross-terms associated with the strong targets may submerge the self-terms of the weak targets. In this case, the “CLEAN” technique [25] can be adopted to adaptively remove the influence of strong target.

3.3. SNR Loss of the NSIAF

It was pointed out in [22,26,27,28] that the use of autocorrelation functions will result in SNR loss due to the emergence of additional noise terms. In this subsection, we quantitatively analyze the SNR loss of the NSIAF defined by (7). For the convenience of analysis, we assume the input signal of the NSIAF is included with independent white Gaussian noise. Let
S 1 = S f r , t n + τ 2 + τ 0 , S 2 = S f r , t n τ 2 τ 0 S 3 = S f r , t n + τ 2 τ 0 , S 4 = S f r , t n τ 2 + τ 0 n 1 = n t n + τ 2 + τ 0 , n 2 = n t n τ 2 τ 0 n 3 = n t n + τ 2 τ 0 , n 4 = n t n τ 2 + τ 0
where  n 1 n 2 n 3 , and  n 4  are the Gaussian white noise with variance of  σ 2 . We define the input SNR as  SNR in = A 2 2 / σ 2 . Substituting (22) into (7) yields
R n f r , t n , τ = S 1 + n 1 S 2 + n 2 S 3 + n 3 S 4 + n 4 = R n s f r , t n , τ + R n c f r , t n , τ
where
R n s f r , t n , τ = A 2 4 exp j 8 π f r + f c c τ 0 τ a + b t n
R n c f r , t n , τ = n 1 n 2 n 3 n 4 + n 1 n 2 n 3 S 4 + n 1 n 2 S 3 n 4 + n 1 S 2 n 3 n 4 + S 1 n 2 n 3 n 4 + n 1 n 2 S 3 S 4 + n 1 S 2 n 3 S 4 + n 1 S 2 S 3 n 4 + S 1 n 2 n 3 S 4 + S 1 n 2 S 3 n 4 + S 1 S 2 n 3 n 4 + S 1 S 2 S 3 n 4 + n 1 S 2 S 3 S 4 + S 1 n 2 S 3 S 4 + S 1 S 2 n 3 S 4
According to (24) and (25), we can derive
E R n s f r , t n , τ 2 = A 2 8
E R n c f r , t n , τ 2 = σ 8 + 4 A 2 2 σ 6 + 6 A 2 4 σ 4 + 4 A 2 6 σ 2
Thus, the output SNR of the NSIAF is
SNR out = E R n s f r , t n , τ 2 E R n c f r , t n , τ 2 = A 2 8 σ 8 + 4 A 2 2 σ 6 + 6 A 2 4 σ 4 + 4 A 2 6 σ 2
Then the SNR loss can be calculated by
SNR out SNR in = 1 4 + 1 SNR in 3 + 4 1 SNR in 2 + 6 1 SNR in
As shown in (29), when  SNR in 1 SNR out / SNR in 1 / 4 , the maximum loss of the output SNR is about 6 dB [17].

4. Simulations and Experiments

In this subsection, we evaluate the performance of the proposed algorithm via simulations for both single target and multiple targets. The simulated radar parameters are listed in Table 1.

4.1. Simulation Results for Single MaT

In this simulation, let us consider a MaT whose motion parameters are set as follows: the initial range  R 0 = 150   km , the radial velocity  v = 480   m / s , the radial acceleration  a = 35   m / s 2 , and the jerk  b = 16   m / s 3 . This experiment is conducted without considering any noise for the purpose of better revealing the properties of the proposed algorithm. The simulation results are presented in Figure 2.
Figure 2a shows the pulse compression result, where serious RM is shown. Figure 2b presents the  t n τ  distribution at  f r = 30   MHz , i.e.,  R 30   MHz , t n , τ . After performing SFT on (11) with respect to  t n , we obtain the result of  Q f n , τ n , as shown in Figure 2c. After performing FFT on  Q f n , τ n  with respect to  τ n , one peak appears, as shown in Figure 2d, according to which the radial acceleration and jerk can be estimated as  a ^ = 35.02   m / s 2  and  b ^ = 16.04   m / s 3 , respectively. The final CI by SCIFT is presented in Figure 2e. Obviously, both  a  and  b  have been very accurately estimated.

4.2. Simulation Results for Multiple MaTs

We now evaluate the performance of the proposed algorithm on detection of four MaTs with noise considered, whose motion parameters are all different and are listed in Table 2. The echoes are added with complex additive Gaussian white noise, and the different SNRs of the four targets’ raw echoes are also listed in Table 2. Figure 3a shows the pulse compression result; as can be seen, four targets are presented with different trajectories. Figure 3b,c present the  r e s u l t s   o f   R 30   MHz , t n , τ  and  Q f n , τ n , respectively. Figure 3d shows the result of  Q f n , f τ n , from which one can see four apparent peaks. All the targets can be identified simultaneously even if their amplitudes are different due to different input SNRs. According to (14), the radial accelerations and jerks of the four targets can be estimated, and they are listed in Table 3. According to (15)–(17), the final CI result of the proposed algorithm is given in Figure 3e. The result of iterative ACCF is given in Figure 3f, where the echo energies are all defocused. Figure 3g shows the result of GRFT, where although four peaks corresponding to these four targets can still be located, the BSSLs indicated by the red circle can also result in false targets, besides lowering the detection probability of the truth target. The MTD [29] result is presented in Figure 3h, where the echo energies are distributed in several range cells and in several Doppler frequency cells, as well due to serious RM and DFM.

4.3. Computational Complexity and Detection Performance

In this subsection, the computation complexity and the detection performance of the proposed are compared with that of GRFT and iterative ACCF.
Suppose the numbers of range bins and pulses are  M  and  N , respectively. The computational complexity for the main procedures of the proposed algorithm is composed of four parts, i.e., calculation of  R f r , t n , τ   [ O ( M N 2 ) ] , the operation of Chirp-z-based KT  [ O ( M N 2 log N ) ] , the FFT and SFT operations  [ O ( 2 M N log N ) ] , and the SCIFT  [ O ( M N 2 ) ] . Thus, the overall computation complexity of the proposed algorithm is  O ( 2 M N 2 + M N 2 log N + 2 M N log N ) . For the computational complexities of GRFT and iterative ACCF, one can refer to [18,19], and these are not repeated here. All complexities are listed in Table 4, where  N v N a  and  N b  represent the searching numbers of radial velocity, acceleration, and jerk, respectively, in GRFT. Obviously, if  N v N a , and  N b  are greater than  N  (in fact, this is usually the case, especially when accurate estimates are required), the computational burden of our algorithm can be significantly lower than that of GRFT. Although the iterative ACCF algorithm has the lowest computational complexity, its anti-noise performance is inferior to ours, as shall be shown later.
In the following, we evaluate the detection performances of the proposed algorithm via Monte Carlo trials. The radar system parameters are the same as those in Table 1, and the target in Section 4.1 is selected. Complex additive white Gaussian noises are added to make the SNRs before pulse compression as [−35:1:0] dB, and 100 Monte Carlo trials are conducted for each SNR value. The detection performances of our algorithm along with the GRFT, iterative ACCF, and MTD are presented in Figure 4, all using the cell averaging constant false alarm (CA-CFAR) by setting the false alarm probability as  P f a = 10 5 . As shown, the GRFT performs the best, whereas our algorithm performs better than the MTD and iterative ACCF.
According to Figure 4 and the above computational complexity analysis, we can conclude that our proposed algorithm can strike a balance between the computational complexity and the detection ability.

4.4. Practical Radar Experiment

In this subsection, the practical Ka-band radar data [30] are processed to validate the proposed algorithm. The basic parameters of the Ka-band radar are provided as follows: the carrier frequency is 35 GHz, the PRF is 32 kHz, and the range cell is about 1.875 m.
The target trajectory is shown in Figure 5a; as can be seen, serious RM exists. Figure 5b shows the MTD result, which shows that the CI performance is degraded due to RM and DFM. The  Q f n , f τ n  result is presented in Figure 5c, from which the acceleration and jerk can be estimated as  a ^ = 0.31   m / s 2  and  b ^ = 0.99   m / s 3 , respectively, according to the peak position. From the true instant velocities provided by the GPS, the true acceleration and jerk corresponding to the selected data segment can be fitted as  a = 0.26   m / s 2  and  b = 1.02   m / s 3 , respectively. The true velocity and the fitted curve are presented in Figure 5d, which shows that the estimated  a  and  b  are very close to the truth values. Figure 5e shows the final CI result. The Doppler frequency profiles are plotted in Figure 5f; as can be seen, our result is about 7.52 dB better than that of MTD.

5. Conclusions

This paper proposes a novel CI algorithm for MaTs having a jerking motion without the need to search for motion parameters. The NSIAF with variable lag time is defined and used to reduce the order on the slow time and to eliminate the LRM. Then the KT and the SFT are applied to obtain the acceleration and jerk estimations, according to which the matched function is constructed and applied, and the CI and target detection can be finally accomplished using the SCIFT. The searching for motion parameters can be avoided in the whole processes. The performance of the proposed algorithm is compared with that of other representative algorithms, and shows that our algorithm can achieve a good balance between the computational complexity and the detection ability. Processing results of both simulated and practical radar data validate the proposed algorithm very well.
Furthermore, we should point out that the proposed algorithm is not specifically designed for the detection of high- or low-velocity targets. For targets with much lower velocities, the applicability depends on whether the LRM, QRM, CRM, LDFM, and QDFM occur. For some radar systems working at much higher carrier frequencies, the low-speed maneuvering targets can also lead to these problems, in which case our algorithm is still applicable. If this is not the case, i.e., the above migrations do not exist simultaneously, then some other algorithms with comparatively lower computational complexity can be used to accomplish the coherent integration.

Author Contributions

Conceptualization, Y.M., J.Y. and Y.Z.; methodology, Y.M. and J.Y.; software, Y.M.; validation, Y.M., J.Y. and Y.Z.; formal analysis, Y.M.; investigation, Y.M. and J.Y.; resources, Y.Z.; writing—original draft preparation, Y.M.; writing—revising and editing, Y.Z.; visualization, Y.M.; supervision, Y.Z.; project administration, Y.Z.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant 61971402.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data of practical radar experiment are from Science Data Bank (https://cstr.cn/31253.11.sciencedb.908).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In Appendix A, we discuss the behaviors of the cross-terms in NSIAF processing. The cross-term  R cross , i , j f r , t n , τ  in (19) can be expressed as
R cross , i , j f r , t n , τ = p = 1 6 R cross , i , j , p f r , t n , τ
where
R cross , i , j , 1 f r , t n , τ = 2 A 2 i 3 A 2 j exp j 4 π f r + f c c 3 R 0 i + R 0 j + 3 v i + v j t n + 1 2 a i a j τ τ 0 + 1 2 b i b j τ τ 0 t n + 3 a i + a j 1 2 t n 2 + 1 8 τ 2 + 1 2 τ 0 2 + 3 b i + b j 1 6 t n 3 + 1 8 τ 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c ( v i v j ) 1 2 τ τ 0 + a i a j 1 2 τ t n τ 0 t n + b i b j 1 48 τ 3 1 6 τ 0 3 + 1 4 τ τ 0 2 1 8 τ 2 τ 0 + 1 4 τ t n 2 1 2 τ 0 t n 2
R cross , i , j , 2 f r , t n , τ = 2 A 2 i 3 A 2 j exp j 4 π f r + f c c R 0 i R 0 j + v i v j t n 1 2 3 a i + a j τ τ 0 1 2 3 b i + b j τ τ 0 t n + a i a j 1 2 t n 2 + 1 8 τ 2 1 2 τ 0 2 + b i b j 1 6 t n 3 + 1 8 τ 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 τ + τ 0 + a i a j 1 2 τ t n + τ 0 t n + b i b j 1 48 τ 3 + 1 6 τ 0 3 + 1 4 τ τ 0 2 + 1 8 τ 2 τ 0 + 1 4 τ t n 2 + 1 2 τ 0 t n 2
R cross , i , j , 3 f r , t n , τ = 2 A 2 i 2 A 2 j 2 exp j 4 π f r + f c c τ τ 0 t n b i + b j + τ τ 0 a i + a j × cos 4 π f r + f c c 2 R 0 i R 0 j + a i a j t n 2 + 1 4 τ 2 + τ 0 2 + 2 v i v j t n + b i b j 1 3 t n 3 + 1 4 τ 2 t n + τ 0 2 t n
R cross , i , j , 4 f r , t n , τ = 2 A 2 i 2 A 2 j 2 exp j 4 π f r + f c c 2 τ τ 0 a i + a j + 2 τ τ 0 t n b i + b j × cos 4 π f r + f c c v i v j 1 2 τ τ 0 + a i a j 1 2 τ t n τ 0 t n + b i b j 1 48 τ 3 1 6 τ 0 3 + 1 4 τ τ 0 2 1 8 τ 2 τ 0 + 1 4 τ t n 2 1 2 τ 0 t n 2
R cross , i , j , 5 f r , t n , τ = 2 A 2 i A 2 j 3 exp j 4 π f r + f c c R 0 i R 0 j v i v j t n + 1 2 a i + 3 a j τ τ 0 + 1 2 b i + 3 b j τ τ 0 t n + a i a j 1 2 t n 2 + 1 8 τ 2 + 1 2 τ 0 2 + b i b j 1 6 t n 3 + 1 8 τ 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 τ + τ 0 + a i a j 1 2 τ t n + τ 0 t n + b i b j 1 48 τ 3 + 1 6 τ 0 3 + 1 4 τ τ 0 2 + 1 8 τ 2 τ 0 + 1 4 τ t n 2 + 1 2 τ 0 t n 2
R cross , i , j , 6 f r , t n , τ = 2 A 2 i A 2 j 3 exp j 4 π f r + f c c R 0 i R 0 j v i v j t n + 1 2 a i + 3 a j τ τ 0 + 1 2 b i + 3 b j τ τ 0 t n a i a j 1 2 t n 2 + 1 8 τ 2 + 1 2 τ 0 2 b i b j 1 6 t n 3 + 1 8 τ 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 τ τ 0 + a i a j 1 2 τ t n τ 0 t n + b i b j 1 48 τ 3 1 6 τ 0 3 + 1 4 τ τ 0 2 1 8 τ 2 τ 0 + 1 4 τ t n 2 1 2 τ 0 t n 2
The cross-term  Q cross , i , j t n , τ n  in (20) can be expressed as
Q cross , i , j ( t n , τ n ) = p = 1 6 SUM f r Q cross , i , j , p f r , t n , τ n
where
Q cross , i , j , 1 f r , t n , τ n = 2 A 2 i 3 A 2 j exp j 4 π f r + f c c 3 R 0 i + R 0 j + 3 v i + v j t n + 1 2 a i a j f c τ n f r + f c τ 0 + 1 2 b i b j f c τ n f r + f c τ 0 t n + 3 a i + a j 1 2 t n 2 + 1 8 f c τ n f r + f c 2 + 1 2 τ 0 2 + 3 b i + b j 1 6 t n 3 + 1 8 f c τ n f r + f c 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c ( v i v j ) 1 2 f c τ n f r + f c τ 0 + a i a j 1 2 f c τ n f r + f c t n τ 0 t n + b i b j 1 48 f c τ n f r + f c 3 1 6 τ 0 3 + 1 4 f c τ n f r + f c τ 0 2 1 8 f c τ n f r + f c 2 τ 0 + 1 4 f c τ n f r + f c t n 2 1 2 τ 0 t n 2
Q cross , i , j , 2 f r , t n , τ n = 2 A 2 i 3 A 2 j exp j 4 π f r + f c c R 0 i R 0 j + v i v j t n 1 2 3 a i + a j f c τ n f r + f c τ 0 1 2 3 b i + b j f c τ n f r + f c τ 0 t n + a i a j 1 2 t n 2 + 1 8 f c τ n f r + f c 2 1 2 τ 0 2 + b i b j 1 6 t n 3 + 1 8 f c τ n f r + f c 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 f c τ n f r + f c + τ 0 + a i a j 1 2 f c τ n f r + f c t n + τ 0 t n + b i b j 1 48 f c τ n f r + f c 3 + 1 6 τ 0 3 + 1 4 f c τ n f r + f c τ 0 2 + 1 8 f c τ n f r + f c 2 τ 0 + 1 4 f c τ n f r + f c t n 2 + 1 2 τ 0 t n 2
Q cross , i , j , 3 f r , t n , τ n = 2 A 2 i 2 A 2 j 2 exp j 4 π f c c τ n τ 0 t n b i + b j + τ n τ 0 a i + a j × cos 4 π f r + f c c 2 R 0 i R 0 j + a i a j t n 2 + 1 4 f c τ n f r + f c 2 + τ 0 2 + 2 v i v j t n + b i b j 1 3 t n 3 + 1 4 f c τ n f r + f c 2 t n + τ 0 2 t n
Q cross , i , j , 4 f r , t n , τ n = 2 A 2 i 2 A 2 j 2 exp j 4 π f c c 2 τ n τ 0 a i + a j + 2 τ n τ 0 t n b i + b j × cos 4 π f r + f c c v i v j 1 2 f c τ n f r + f c τ 0 + a i a j 1 2 f c τ n f r + f c t n τ 0 t n + b i b j 1 48 f c τ n f r + f c 3 1 6 τ 0 3 + 1 4 f c τ n f r + f c τ 0 2 1 8 f c τ n f r + f c 2 τ 0 + 1 4 f c τ n f r + f c t n 2 1 2 τ 0 t n 2
Q cross , i , j , 5 f r , t n , τ n = 2 A 2 i A 2 j 3 exp j 4 π f r + f c c R 0 i R 0 j v i v j t n + 1 2 a i + 3 a j f c τ n f r + f c τ 0 + 1 2 b i + 3 b j f c τ n f r + f c τ 0 t n + a i a j 1 2 t n 2 + 1 8 f c τ n f r + f c 2 + 1 2 τ 0 2 + b i b j 1 6 t n 3 + 1 8 f c τ n f r + f c 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 f c τ n f r + f c + τ 0 + a i a j 1 2 f c τ n f r + f c t n + τ 0 t n + b i b j 1 48 f c τ n f r + f c 3 + 1 6 τ 0 3 + 1 4 f c τ n f r + f c τ 0 2 + 1 8 f c τ n f r + f c 2 τ 0 + 1 4 f c τ n f r + f c t n 2 + 1 2 τ 0 t n 2
Q cross , i , j , 6 f r , t n , τ n = 2 A 2 i A 2 j 3 exp j 4 π f r + f c c R 0 i R 0 j v i v j t n + 1 2 a i + 3 a j f c τ n f r + f c τ 0 + 1 2 b i + 3 b j f c τ n f r + f c τ 0 t n a i a j 1 2 t n 2 + 1 8 f c τ n f r + f c 2 + 1 2 τ 0 2 b i b j 1 6 t n 3 + 1 8 f c τ n f r + f c 2 t n + 1 2 τ 0 2 t n × cos 4 π f r + f c c v i v j 1 2 f c τ n f r + f c τ 0 + a i a j 1 2 f c τ n f r + f c t n τ 0 t n + b i b j 1 48 f c τ n f r + f c 3 1 6 τ 0 3 + 1 4 f c τ n f r + f c τ 0 2 1 8 f c τ n f r + f c 2 τ 0 + 1 4 f c τ n f r + f c t n 2 1 2 τ 0 t n 2
Compared with the self-term in (20), the following three aspects can be observed:
(a) The cosine function terms will degrade the accumulation of the cross-terms. In fact, only when  v i = v j a i = a j  and  b i = b j , the cosine function terms equal to 1 [31,32,33,34].
(b) After KT, the new coupling between  f r  and  τ n  is generated, which means the KT does not work for the cross-term.
(c) In contrast to (20), where  f r  has been decoupled from  τ  and  t n , the echo energies will be dispersed by direct summation along  f r .
The cross-term  Q cross , i , j f n , f τ n  in (21) can be expressed as
Q cross , i , j f n , f τ n = FFT τ n SFT t n Q cross , i , j ( t n , τ n )
It is obvious that the high-order terms related to the time variable  τ n  and  t n  exist in the cosine terms and the exponential terms in  Q cross , i , j t n , τ n , such as  t n 2 t n 3 τ n 2 , and  τ n 3 . After performing FFT with respect to  τ n  and SFT with respect to  t n , all the cross-terms cannot be integrated along straight lines parallel to the  t n  axis in the  t n f τ n  domain or to the  τ n  axis in the  f n τ n  domain due to the complex coupling between  t n  and  τ n  [24]. Therefore, the cross-term  Q cross , i , j f n , f τ n  cannot form any peaks in the  f n f τ n  domain. Compared with the self-term in (21), the cross terms will be smeared in this transformed domain.
The above analyses show that the cross-terms’ contribution cannot be accumulated well, so they are much smaller than those of the self-terms and thus can be ignored.

References

  1. Zheng, J.; Liu, H.; Liu, J.; Du, X.; Liu, Q.H. Radar High-Speed Maneuvering Target Detection Based on Three-Dimensional Scaled Transform. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2018, 11, 2821–2833. [Google Scholar] [CrossRef]
  2. Zhang, J.; Su, T.; Zheng, J.; He, X. Novel Fast Coherent Detection Algorithm for Radar Maneuvering Target with Jerk Motion. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2017, 10, 1792–1803. [Google Scholar] [CrossRef]
  3. 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, 1–5. [Google Scholar] [CrossRef]
  4. You, P.; Ding, Z.; Liu, S.; Zhang, G. Dechirp-receiving Radar Target Detection Based on Generalized Radon-Fourier Transform. IET Radar Sonar Navig. 2021, 15, 1096–1111. [Google Scholar] [CrossRef]
  5. Jin, K.; Lai, T.; Wang, Y.; Li, G.; 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]
  6. Gao, X.; Zhang, H.; Dang, T. Radar Detection and Fast Motion Parameter Estimation for Complex Manoeuvering Targets at High Speed and Acceleration. IET Radar Sonar Navig. 2022, 16, 1977–1996. [Google Scholar] [CrossRef]
  7. Wan, J.; He, Z.; Tan, X.; Li, D.; Liu, H.; Shu, Y.; Chen, Z. Coherent Integration for Maneuvering Target Detection via Fast Nonparametric Estimation Method. Signal Process. 2023, 203, 108820. [Google Scholar] [CrossRef]
  8. 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]
  9. Yao, D.; Zhang, X.; Sun, Z. Long-Time Coherent Integration for Maneuvering Target Based on Second-Order Keystone Transform and Lv’s Distribution. Electronics 2022, 11, 1961. [Google Scholar] [CrossRef]
  10. Zhu, D.Y.; Li, Y.; Zhu, Z.D. A Keystone Transform Without Interpolation for SAR Ground Moving-Target Imaging. IEEE Geosci. Remote Sens. Lett. 2007, 4, 18–22. [Google Scholar] [CrossRef]
  11. Su, J.; Xing, M.; Wang, G.; Bao, Z. High-Speed Multi-Target Detection with Narrowband Radar. IET Radar Sonar Navig. 2010, 4, 595. [Google Scholar] [CrossRef]
  12. Rao, X.; Tao, H.; Su, J.; Guo, X.; Zhang, J. Axis Rotation MTD Algorithm for Weak Target Detection. Digit. Signal Prog. 2014, 26, 81–86. [Google Scholar] [CrossRef]
  13. Zheng, J.; Su, T.; Zhu, W.; He, X.; Liu, Q.H. Radar High-Speed Target Detection Based on the Scaled Inverse Fourier Transform. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2015, 8, 1108–1119. [Google Scholar] [CrossRef]
  14. 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]
  15. Li, X.; Sun, Z.; Yi, W.; Cui, G.; Kong, L.; Yang, X.B. Computationally Efficient Coherent Detection and Parameter Estimation Algorithm for Maneuvering Target. Signal Process. 2019, 155, 130–142. [Google Scholar] [CrossRef]
  16. Li, X.; Cui, G.; Yi, W.; Kong, L. Manoeuvring Target Detection Based on Keystone Transform and Lv’s Distribution. IET Radar Sonar Navig. 2016, 10, 1234–1242. [Google Scholar] [CrossRef]
  17. Pang, C.; Liu, S.; Han, Y. Coherent Detection Algorithm for Radar Maneuvering Targets Based on Discrete Polynomial-Phase Transform. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2019, 12, 3412–3422. [Google Scholar] [CrossRef]
  18. Xu, J.; Xia, X.G.; Peng, S.B.; Ji, Y.; Peng, Y.N.; Qian, L.C. Radar Maneuvering Target Motion Estimation Based on Generalized Radon-Fourier Transform. IEEE Trans. Signal Process. 2012, 60, 6190–6201. [Google Scholar] [CrossRef]
  19. Li, X.; Cui, G.; Kong, L.; Yi, W. Fast Non-Searching Method for Maneuvering Target Detection and Motion Parameters Estimation. IEEE Trans. Signal Process. 2016, 64, 2232–2244. [Google Scholar] [CrossRef]
  20. 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 Observ. Remote Sens. 2016, 9, 285–294. [Google Scholar] [CrossRef]
  21. Bai, X.; Tao, R.; Wang, Z.; Wang, Y. ISAR Imaging of a Ship Target Based on Parameter Estimation of Multicomponent Quadratic Frequency-Modulated Signals. IEEE Trans. Geosci. Remote Sens. 2014, 52, 1418–1429. [Google Scholar] [CrossRef]
  22. Cao, R.; Li, M.; Zuo, L.; Sun, H. Fast Parameter Estimation Method for Maneuvering Target by Using Non-Uniformly Resampling Reducing Order Technique. Signal Process. 2019, 160, 299–315. [Google Scholar] [CrossRef]
  23. Djurovic, I.; Simeunovic, M.; Djukanovic, S.; Wang, P. A Hybrid CPF-HAF Estimation of Polynomial-Phase Signals: Detailed Statistical Analysis. IEEE Trans. Signal Process. 2012, 60, 5010–5023. [Google Scholar] [CrossRef]
  24. Huang, P.; Xia, X.; Liao, G.; Yang, Z.; Zhou, J.; Liu, X. Ground Moving Target Refocusing in SAR Imagery Using Scaled GHAF. IEEE Trans. Geosci. Remote Sens. 2018, 56, 1030–1045. [Google Scholar] [CrossRef]
  25. Li, X.; Kong, L.; Cui, G.; Yi, W. CLEAN-based Coherent Integration Method for High-speed Multi-targets Detection. IET Radar Sonar Navig. 2016, 10, 1671–1682. [Google Scholar] [CrossRef]
  26. Li, Y.; Zhang, J.; Niu, J.; Zhou, Y.; Wang, L. Computational Implementation and Asymptotic Statistical Performance Analysis of Range Frequency Autocorrelation Function for Radar High-Speed Target Detection. IEEE Trans. Comput. Imaging 2020, 6, 1297–1308. [Google Scholar] [CrossRef]
  27. Barbarossa, S. Analysis of Multicomponent LFM Signals by a Combined Wigner-Hough Transform. IEEE Trans. Signal Process. 1995, 43, 1511–1515. [Google Scholar] [CrossRef]
  28. Li, Y.; Zhang, J.; Zheng, J.; Niu, J.; Zhou, Y.; Meng, N. ISAR Imaging of Nonuniformly Rotating Targets With Low SNR Based on Third Order Autocorrelation Function. IEEE Access. 2020, 8, 74707–74719. [Google Scholar] [CrossRef]
  29. Barton, D.K. Radar System Analysis and Modeling; Artech House: Norfolk County, MA, USA, 2004; ISBN 1-58053-681-6. [Google Scholar]
  30. Song, Z.; Hui, B.; Fan, H.; Zhou, J.; Zhu, Y.; Da, K.; Zhang, X.; Su, H.; Jin, W.; Zhang, Y. A dataset for dim target detection and tracking of aircraft in radar echo sequences. Sci. Data Bank 2019. [Google Scholar] [CrossRef]
  31. Zhang, J.; Li, Y.; Su, T.; He, X. Quadratic FM Signal Detection and Parameter Estimation Using Coherently Integrated Trilinear Autocorrelation Function. IEEE Trans. Signal Process. 2020, 68, 621–633. [Google Scholar] [CrossRef]
  32. Li, D.; Zhan, M.; Su, J.; Liu, H.; Zhang, X.; Liao, G. Performances Analysis of Coherently Integrated CPF for LFM Signal Under Low SNR and Its Application to Ground Moving Target Imaging. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6402–6419. [Google Scholar] [CrossRef]
  33. Gao, X.; Zhang, H. Radar Coherent Integration Algorithm for Detection of Complexly Maneuvering Target with Extended Velocity and Acceleration Scopes. Digit. Signal Process. 2022, 129, 103689. [Google Scholar] [CrossRef]
  34. Huang, P.; Liao, G.; Yang, Z.; Xia, X.-G.; Ma, J.; Zheng, J. Ground Maneuvering Target Imaging and High-Order Motion Parameter Estimation Based on Second-Order Keystone and Generalized Hough-HAF Transform. IEEE Trans. Geosci. Remote Sens. 2017, 55, 320–335. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed algorithm.
Figure 1. Flowchart of the proposed algorithm.
Electronics 12 02363 g001
Figure 2. Simulation results for a single target without noise considered. (a) Pulse compression. (b) NSIAF at  f r = 30   MHz . (c Q f n , τ n . (d Q f n , f τ n . (e) CI.
Figure 2. Simulation results for a single target without noise considered. (a) Pulse compression. (b) NSIAF at  f r = 30   MHz . (c Q f n , τ n . (d Q f n , f τ n . (e) CI.
Electronics 12 02363 g002
Figure 3. Simulation results for four targets with noise considered. (a) Pulse compression. (b) NSIAF at  f r = 30   MHz . (c Q f n , τ n . (d Q f n , f τ n . (e) CI by our algorithm. (f) CI by ACCF iteratively. (g) CI by GRFT. (h) MTD.
Figure 3. Simulation results for four targets with noise considered. (a) Pulse compression. (b) NSIAF at  f r = 30   MHz . (c Q f n , τ n . (d Q f n , f τ n . (e) CI by our algorithm. (f) CI by ACCF iteratively. (g) CI by GRFT. (h) MTD.
Electronics 12 02363 g003aElectronics 12 02363 g003b
Figure 4. Comparison of detection performances of our algorithm along with the GRFT, iterative ACCF, and MTD.
Figure 4. Comparison of detection performances of our algorithm along with the GRFT, iterative ACCF, and MTD.
Electronics 12 02363 g004
Figure 5. Experimental results on practical radar data. (a) Pulse compression. (b) MTD. (c Q f n , f τ n  result. (d) True value and fitted curve. (e) Integrated result. (f) Doppler frequency profiles.
Figure 5. Experimental results on practical radar data. (a) Pulse compression. (b) MTD. (c Q f n , f τ n  result. (d) True value and fitted curve. (e) Integrated result. (f) Doppler frequency profiles.
Electronics 12 02363 g005
Table 1. Radar parameters.
Table 1. Radar parameters.
ParametersValue
Carrier frequency1 GHz
Bandwidth30 MHz
Sampling frequency60 MHz
Pulse duration5 μs
Pulse repetition frequency400 Hz
Pulse number500
Table 2. Simulated targets’ parameters.
Table 2. Simulated targets’ parameters.
ParametersTarget ATarget BTarget CTarget D
SNR (dB)−11.61−12.43−13.85−14.94
Initial range (km)150150.1150.2150.05
Radial velocity (m/s)−500420300−250
Radial acceleration (m/s2)40−2732−10
Radial jerk (m/s3)−181216−5
Table 3. Estimated targets’ accelerations and jerks.
Table 3. Estimated targets’ accelerations and jerks.
Estimated ParametersTarget ATarget BTarget CTarget D
Estimated radial acceleration (m/s2)39.99−27.1032.06−9.92
Estimated radial jerk (m/s3)−18.1811.9016.20−4.63
Table 4. Computational complexity.
Table 4. Computational complexity.
AlgorithmsComputational Complexity
GRFT O ( N v N a N b M N )
Iterative ACCF O ( M N log N + M N log M )
The proposed algorithm O ( 2 M N 2 + M N 2 log N + 2 M N log N )
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

Mi, Y.; Zhang, Y.; Yang, J. A Novel Coherent Integration Algorithm for Maneuvering Target Detection Based on Symmetric Instantaneous Autocorrelation Function. Electronics 2023, 12, 2363. https://doi.org/10.3390/electronics12112363

AMA Style

Mi Y, Zhang Y, Yang J. A Novel Coherent Integration Algorithm for Maneuvering Target Detection Based on Symmetric Instantaneous Autocorrelation Function. Electronics. 2023; 12(11):2363. https://doi.org/10.3390/electronics12112363

Chicago/Turabian Style

Mi, Yunpeng, Yunhua Zhang, and Jiefang Yang. 2023. "A Novel Coherent Integration Algorithm for Maneuvering Target Detection Based on Symmetric Instantaneous Autocorrelation Function" Electronics 12, no. 11: 2363. https://doi.org/10.3390/electronics12112363

APA Style

Mi, Y., Zhang, Y., & Yang, J. (2023). A Novel Coherent Integration Algorithm for Maneuvering Target Detection Based on Symmetric Instantaneous Autocorrelation Function. Electronics, 12(11), 2363. https://doi.org/10.3390/electronics12112363

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