Next Article in Journal
Application of Lagrange Relaxation to Decentralized Optimization of Dispatching a Charging Station for Electric Vehicles
Next Article in Special Issue
Remote Sensing Image Fusion Based on Sparse Representation and Guided Filtering
Previous Article in Journal
An Efficient 3-D Ray Tracing Method: Prediction of Indoor Radio Propagation at 28 GHz in 5G Network
Previous Article in Special Issue
Retrieval of Three-Dimensional Surface Deformation Using an Improved Differential SAR Tomography System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coherent Integration for Radar High-Speed Maneuvering Target Based on Frequency-Domain Second-Order Phase Difference

1
National Digital Switching System Engineering and Technological Research Center, Zhengzhou 450000, China
2
Air Traffic Control and Navigation College, Air Force Engineering University, Xi’an 710000, China
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(3), 287; https://doi.org/10.3390/electronics8030287
Submission received: 18 January 2019 / Revised: 24 February 2019 / Accepted: 27 February 2019 / Published: 4 March 2019

Abstract

:
In recent years, target detection has drawn increasing attention in the field of radar signal processing. In this paper, we address the problem of coherent integration for detecting high-speed maneuvering targets, involving range migration (RM), quadratic RM (QRM), and Doppler frequency migration (DFM) within the coherent processing interval. We propose a novel coherent integration algorithm based on the frequency-domain second-order phase difference (FD-SoPD) approach. First, we use the FD-SoPD operation to reduce the signal from three to two dimensions and simultaneously eliminate the effects of QRM and DFM, which leads to signal-to-noise ratio improvement in the velocity-acceleration domain. Next, we estimate the target motion parameters from the peak position without the need for a search procedure. We show that this algorithm can be easily implemented by using complex multiplications combined with fast Fourier transform (FFT) and inverse FFT (IFFT) operations. We perform comparisons with several representative algorithms and show that the proposed technique can be used to achieve a good trade-off between computational complexity and detection performance. We present both simulated and experimental data to demonstrate the effectiveness of the proposed method.

1. Introduction

With the increasing requirements for space target detection and high-resolution imaging, radar high-speed maneuvering target detection has drawn growing attention [1,2,3,4,5,6,7,8,9,10,11]. Normally, a low-speed target is located in the same range cell during the short observation time, and the conventional moving target detection (MTD) algorithm [12] can complete coherent integration by using fast Fourier transform (FFT). It is well known that in order to improve the detection ability in far-range and low radar cross section (RCS) targets, a long-term coherent integration is always required [13]. In this case, for high-speed maneuvering targets, the linear range migration (LRM), quadratic range migration (QRM), and Doppler frequency migration (DFM) effects will inevitably occur, thereby deteriorating integration performance.
As for radar coherent integration, many successful detection algorithms have been proposed, such as the keystone transform (KT) [14,15], scaled inverse Fourier transform (SCIFT) [16,17], frequency-domain deramp-keystone transform (FDDKT) [18], modified location rotation transform (MLRT) [19], and Radon Fourier transform (RFT) [20]. For a moving target with linear range migration, these algorithms achieve satisfactory antinoise performance, parameter estimation accuracy, and detection ability with reasonable computational cost. Nevertheless, they may suffer from integration performance degradation due to ignoring the effects of QRM and DFM caused by the target’s acceleration.
To address these issues, many advanced methods have been recently proposed. They can be roughly divided into three categories.
(a)
Radon transform-based algorithms, such as generalized Radon Fourier transform (GRFT) [21], Radon-fractional Fourier transform (RFRFT) [22], and Radon-Lv’s Distribution (RLVD) [23,24,25,26]. These kinds of algorithms implement phase compensation and parameter estimation by searching the maneuvering target motion trajectory. Although they can obtain coherent integration under a low signal-to-noise ratio (SNR), the huge computational load seriously limits their practical application.
(b)
KT based algorithms, such as second-order keystone transform (SoKT) [27], Doppler keystone transform (DKT) [28], keystone-Lv’s distribution (KT-LVD) [29], and so on. The KT is used to correct the QRM blindly, which reduces the calculation cost to a certain extent, but it still needs to use parameter searching to eliminate the Doppler ambiguity.
(c)
Correlation-based algorithms: The representative adjacent cross-correlation function and Lv’s distribution (ACCF-LVD) algorithm proposed in References [30,31,32] reduces the migration order by ACCF and quickly estimates the motion parameters without any searching procedure, which greatly reduces the computational burden and benefits practical applications. Unfortunately, this method is only effective when the input SNR is high [33]. The three-dimensional scaled transform (TDST) method was then presented to realize coherent integration and motion parameters estimations for maneuvering targets under a low SNR background [34]. This method eliminates the coupling effectively among spatial frequency, slow time, and time delay. However, the complex three-dimensional transform is usually less suitable for realistic applications.
Aiming to realize the coherent integration of radar high-speed maneuvering targets with low computational complexity, we propose a novel frequency-domain second-order phase difference (FD-SoPD) algorithm in this paper. First, the SoPD is used in the spatial frequency domain to eliminate the impact of acceleration. Then, we can simultaneously estimate the velocity and acceleration from the peak position, followed by phase compensation and coherent integration. The proposed technique has the following contributions: (a) It reduces the signal from three to two dimensions, thus avoiding the complex operation of TDST; (b) the target motion parameters can be easily estimated by FFT without any searching process; (c) the phase difference eliminates the Doppler ambiguity, thus the high speed of target can be accurately estimated; (d) it achieves a good balance between the computational cost and detection ability. Finally, we present both simulated and experimental data to demonstrate the effectiveness of the proposed method.
The remainder of this paper is organized as follows. In Section 2, the signal model for the maneuvering target is established. In Section 3, we deduce the principle of FD-SoPD in detail and discuss the situations of single target and multi-targets, respectively. Section 4 analyses the computational burden. In Section 5, we evaluate the performance via several numerical experiments. Finally, conclusions are drawn in Section 6.

2. Signal Model and Problem Formulation

Suppose the radar transmits a linear frequency modulated (LFM) signal, which can be expressed as:
s t ( t ^ ) = rect ( t ^ T p ) exp ( j 2 π f c t ^ + j π γ t ^ 2 )
where,
rect ( t ^ T p ) = { 1 ,    | t ^ | < T p / 2 0 ,    | t ^ | > T p / 2
is the rectangular window function, and t ^ is the fast time. T p , f c and γ indicate the pulse width, carrier frequency, and frequency modulation rate, respectively. Assume that there are K targets in the scene of radar observation. During the accumulation time, the distance between the maneuvering target and radar can be approximated as second order, i.e.,
R i ( t m ) = r i + v i t m + a i t m 2 / 2
where t m = m / P R F is the slow time, m and P R F denote the transmitted pulse number index and the pulse repetition frequency (PRF). r i , v i , and a i are respectively the initial slant range, radial velocity, and acceleration of the i th target.
Ignoring the influence of noise, the received signal of K targets after down conversion can be expressed as [35]:
s r ( t ^ , t m ) = i = 1 K A r , i rect ( t ^ 2 R i ( t m ) / c T p ) exp { j 4 π f c R i ( t m ) c } exp { j π γ ( t ^ 2 R i ( t m ) c ) 2 }
where A r , i is the target reflectivity, and c is the light speed.
After pulse compression, the radar echoes are written as:
s c ( t ^ , t m ) = i = 1 K A c , i sinc [ B ( t ^ 2 R i ( t m ) c ) ] exp { j 4 π f c R i ( t m ) c }
where A c , i denotes the amplitude after compression and B = γ T p is the bandwidth of transmitted signal.
Substituting Equation (3) into Equation (5), we obtain:
s c ( t ^ , t m ) = i = 1 K A c , i sinc [ B ( t ^ 2 ( r i + v i t m + a i t m 2 / 2 ) c ) ] exp { j 4 π f c ( r i + v i t m + a i t m 2 / 2 ) c }
As can be seen from Equation (6), the signal envelope indicates the target range, which changes nonlinearly with the slow time. When the offset exceeds the range resolution Δ r = c / 2 B , the LRM will occur. If the target has a large acceleration, the QRM can be seen in the envelope. In this case, the conventional MTD is invalid. In addition, the phase in Equation (6) indicates a linear change in the Doppler frequency of the i th target, i.e.,
f d , i = 2 λ d ( r i + v i t m + a i t m 2 / 2 ) d t m = 2 v i + 2 a i t m λ
where λ = c / f c is the wave length. Similarly, when the offset exceeds a Doppler resolution, DFM would occur and defocus the target energy in the Doppler domain. Moreover, for high-speed targets, we often have f d , i P R F , which induces the Doppler ambiguity and makes it hard to estimate the target’s velocity. Therefore, the coherent accumulation of high-speed maneuvering targets can only be achieved by effectively eliminating the LRM, QRM, and DFM.

3. The Principle of the FD-SoPD

3.1. FD-SoPD with Mono-Target

According to Equation (6), the compressed signal of the i th target is
s c ( t ^ , t m ) = A c , i sinc [ B ( t ^ 2 ( r i + v i t m + a i t m 2 / 2 ) c ) ] exp { j 4 π f c ( r i + v i t m + a i t m 2 / 2 ) c }
Performing the Fourier transform (FT) along the t ^ axis, we obtain the spatial spectrum of the signal, i.e.,
S ( f r , t m ) = A f r , i r e c t ( f r B ) exp ( j 4 π ( f r + f c ) R i ( t m ) c )                          = A f r , i r e c t ( f r B ) exp { j 4 π ( f r + f c ) c ( r i + v i t m + a i t m 2 / 2 ) }
where f r is the frequency of the spatial harmonic from the spatial spectrum of the fast time signal record, and A f r , i is the amplitude of the spatial harmonic obtained by FT of the fast time signal record.
As shown in Equation (9), the coupling between f r and t m (or t m 2 ) is the fundamental cause of LRM or QRM. Moreover, the existence of t m 2 broadens the Doppler spectrum and makes the signal energy defocused. If the velocity v i and acceleration a i are accurately estimated, it is easy to perform phase compensation and coherent integration.
The proposed FD-SoPD is defined as:
R S o P D ( t m , τ , τ 1 ; f r ) = S ( f r ; t m + τ 2 ) S * ( f r ; t m τ 2 ) [ S ( f r ; t m + τ 1 2 ) S * ( f r ; t m τ 1 2 ) ] *
where τ and τ 1 are two lag time variables. Substituting Equation (9) into Equation (10) yields:
R S o P D ( f r , t m , τ ; τ 0 ) = | A f r , i | 4 r e c t ( f r B ) exp [ j 4 π f r + f c c v i ( τ 1 τ ) ]                                                     × exp [ j 4 π f r + f c c a i t m ( τ 1 τ ) ]
When τ and τ 1 have fixed nonzero lag time difference, the coupling between τ and τ 1 will be eliminated, i.e.,
2 τ 0 = τ 1 τ
where τ 0 is the fixed lag time. Equation (11) can be further expressed as:
R S o P D ( f r , t m , τ ; τ 0 ) = | A f r , i | 4 r e c t ( f r B ) exp ( j 8 π f r + f c c v i τ 0 )                                                     × exp ( j 8 π f r + f c c a i τ 0 t m )
As shown in Equation (13), three axes, f r , t m and τ , exist in R S o P D ( f r , t m , τ ; τ 0 ) . However, the signal energy is constant along the τ axis, and thus can be accumulated coherently by the addition operation as follows:
R A ( f r , t m ; τ 0 ) = ADD τ [ R S o P D ( f r , t m , τ ; τ 0 ) ]                                        = G m | A f r , i | 4 r e c t ( f r B ) exp ( j 8 π f r + f c c v i τ 0 ) exp ( j 8 π f r + f c c a i τ 0 t m )
where ADD τ ( · ) is the addition function along the τ axis and G m denotes the corresponding integration gain.
Remark 1.
From Equation (14), we may find three features of the FD-SoPD. (a) The signal is reduced from three to two dimensions, which avoids the multidimensional scaled transform in TDST algorithm. (b) The QRM and DRM are simultaneously eliminated. (c) Equation (14) is equivalent to a uniform motion model, whose velocity is 2 a i τ 0 . Thus, the Doppler ambiguity is eliminated. If the envelope migration caused by the velocity 2 a i τ 0 exceeds a range cell, the KT is needed, i.e.,
( f r + f c ) t m = f c τ m
where τ m denotes the scaled slow-time variable.
After performing the KT on Equation (14), we have:
R A ( f r , τ m ; τ 0 ) = G m | A f r , i | 4 r e c t ( f r B ) exp ( j 8 π f c c v i τ 0 )                                      × exp ( j 8 π f r c v i τ 0 ) exp ( j 8 π f c c a i τ 0 τ m )
Applying the FT with respect to f r and τ m , we get:
S F ( t ^ , f τ m ) = FT τ m { FT f r [ R A ( f r , τ m ; τ 0 ) ] } = A F , i exp ( j 8 π f c v i τ 0 c ) sinc [ B ( t ^ + 4 v i τ 0 c ) ] sinc [ C I ( f τ m 4 a i τ 0 λ ) ]
where A F , i is amplitude after two-dimensional FT, f τ m is the frequency with respect to τ m , and C I denotes the coherent integration time.
From Equation (17), we can simultaneously estimate the velocity v i and acceleration a i of the i th target, i.e.,
( v ^ i = c t ^ max 4 τ 0 , a ^ i = λ f τ m , max 4 τ 0 )
Notice that the fixed lag time constant τ 0 is important in the implementation of the SoPD. A large fixed lag time τ 0 will improve the parameter estimation accuracy, whereas spectrum aliasing may occur. Therefore, the compromise consideration usually chooses τ 0 < T a / 4 , where T a is the observation time.
Utilizing the estimated parameters, we can construct the phase compensation function to compensate the LRM, QRM, and DFM in Equation (9)
H c o m ( f r , t m ) = exp ( j 4 π f r v i t m + a i t m 2 / 2 c ) exp ( j 4 π f c a i t m 2 / 2 c )
Finally, the signal energy will be integrated by the IFT and the FT operations,
E ( t ^ , f d ) = FT t m { IFT f r [ S ( f r , t m ) H c o m ( f r , t m ) ] }                            = A E , i sinc [ B ( t ^ 2 r i c ) ] δ ( f d + f d 0 , i )
where f d is the Doppler frequency with respect to t m , and f d 0 , i is the Doppler frequency of the target.
In Equation (20), the signal energy of a high-speed maneuvering target is integrated into a single peak in the range-Doppler domain. The peak position is ( 2 r i / c , f d 0 , i ) and the peak value is | E ( 2 r i / c , f d 0 , i ) | . Here, the constant false alarm rate (CFAR) [36] technique can be used for the target detection, i.e.,
| E ( 2 r i / c , f d 0 , i ) |    H 1 H 0    η
where η is the CFAR threshold. If | E ( 2 r i / c , f d 0 , i ) | is larger than the threshold, there will be a moving target. Otherwise, no target is detected.
Remark 2.
Different from the Radon transform based algorithms in [13,21,22,23], the proposed FD-SoPD method avoids the brute-force searching procedure of unknown motion parameters. In addition, it can be easily implemented by FFT and IFFT, which significantly reduces the computational complexity.
In the following, we will give an example to demonstrate how the FD-SoPD works to accomplish target motion parameter estimation and coherent integration.
Example 1.
We use an ideal maneuvering point target in this example. The parameters of frequency-modulated continuous-wave (FMCW) radar are set as: The carrier frequency f c = 1 GHz , the bandwidth B = 100   MHz , the sampling frequency f s = 2   MHz , pulse repetition frequency PRF = 128   Hz , and the number of integration pulses M = 256 and τ 0 = T a / 5 . The motion parameters of the maneuvering point target are: A r , i = 1 , r i = 3 km , v i = 15 m / s , a i = 1 m / s 2 . Simulation results are shown in Figure 1.
Figure 1a shows the result of pulse compression, where serious LRM occurs due to the target’s high speed and the radar’s high resolution. Figure 1b shows the result of FD-SoPD transform. It is obvious that, after the FD-SoPD with respect to slow time, the envelope migration is effectively eliminated, and the target energy is located in the same range cell. Performing the FT with respect to τ m , we can get the velocity-acceleration distribution, as shown in Figure 1c. The target energy is well accumulated as one peak, and we can estimate v ^ i = 15.01 m / s and a ^ i = 1.031 m / s 2 from its position. Finally, coherent integration is easily accomplished with FFT, and the result is given in Figure 1d.

3.2. FD-SoPD with Multi-Targets

In this subsection, we will analyze the performance of FD-SoPD under multi-targets in detail. Assume that there are K maneuvering targets in the scene of radar observation. The compressed signal can be expressed as Equation (6). Accordingly, the signal in the spatial frequency domain is:
S ( f r , t m ) = i = 1 K A f r , i r e c t ( f r B ) exp ( j 4 π ( f r + f c ) R i ( t m ) c )                            = i = 1 K A f r , i r e c t ( f r B ) exp { j 4 π ( f r + f c ) c ( r i + v i t m + a i t m 2 / 2 ) }
Substituting Equation (22) into Equation (10) along the slow time, we obtain the FD-SoPD of S ( f r , t m )
R S o P D ( f r , t m , τ ; τ 0 ) = R a u t o ( f r , t m , τ ; τ 0 ) + R c r o s s ( f r , t m , τ ; τ 0 )
where R a u t o ( f r , t m , τ ; τ 0 ) and R c r o s s ( f r , t m , τ ; τ 0 ) denote the auto-terms and cross terms, and can be written as
R a u t o ( f r , t m , τ ; τ 0 ) = i = 1 K | A f r , i | 4 r e c t ( f r B ) exp ( j 8 π f r + f c c v i τ 0 )                                                  × exp ( j 8 π f r + f c c a i τ 0 t m )
R c r o s s ( f r , t m , τ ; τ 0 ) = R 2 + R 3 + R 4 + R 5
The summation R 5 can be further expanded as the following three parts:
R 5 = R 6 + R 7 + R 8
The detailed expressions of R i ( i = 2 , 3 , 8 ) are given in the Appendix A.
After the addition, KT and two-dimensional FT, the velocity and acceleration of target will be estimated simultaneously, i.e.,
S F ( t ^ , f τ m ) = FT τ m { FT f r { KT { ADD τ [ R S o P D ( f r , t m , τ ; τ 0 ) ] } } }
According to the specific motion of the maneuvering target, we consider the cross-terms resulting from the following two cases.
Case 1.
The acceleration of any two targets is different, i.e., a i a j , i , j = 1 , 2 , , K , i j . In this case, R a u t o has a similar form with Equation (13), which can be integrated after FFT. R 2 has the linear term of τ and the coupling term between t m and τ , which cannot be accumulated in the addition operation and FFT of Equation (27). R 3 , R 4 , R 7 and R 8 have the quadratic term τ 2 and a coupling term between t m and τ , and thus cannot be accumulated as well. It is known from Equation (34) that R 6 has a symmetric property about τ , which will become a sinusoidal oscillation term after the addition in Equation (27). Therefore, the energy of R 6 will be smeared after performing two-dimensional FFT. In summary, the cross terms can be ignored compared to the auto-terms. Here, we give an example to illustrate the discussion of Case 1.
Example 2.
In this example, we use two maneuvering targets designated as Tr1 and Tr2, respectively. Radar parameters are the same as those in Example 1. Target motion parameters are set as: A r , 1 = 1 , r 1 = 3 km , v 1 = 15 m / s , a 1 = 1 m / s 2 for target Tr1; A r , 2 = 1 , r 2 = 3.2 km , v 2 = 12 m / s , a 2 = 0.6 m / s 2 for target Tr2.
Figure 2a is the result of pulse compression. Figure 2b gives the result of FD-SoPD. It is obvious that the auto-terms are corrected into beelines, while the cross-terms cannot be corrected. Thus, after the addition and two-dimensional FT, only the auto-terms are accumulated into two peaks, as shown in Figure 2c. We can estimate the velocity and acceleration of targets as v ^ 1 = 15.01 m / s , a ^ 1 = 1.067 m / s 2 , v ^ 2 = 12 m / s , a ^ 2 = 0.5822 m / s 2 . After compensating the RM, QRM, and DFM with the estimated motion parameters, these two targets are coherently integrated, as shown in Figure 2d,e.
Case 2.
Some of the accelerations coincide, i.e., a i = a j or c i 2 = c j 2 , i , j = 1 , 2 , , K , i j . In this case, the coupling term between t m and τ in R 2 is eliminated, but there are still linear terms of τ . The energy of R 2 will be accumulated only when c i 1 = c j 1 , which means the two targets have the same velocity and acceleration. i.e., R 1 = R 2 . R 3 , R 4 and R 8 have the quadratic term τ 2 and coupling term between t m and τ , and thus cannot be accumulated. R 7 has linear terms of τ and a random initial phase regarding target reflectivity and the initial range, which defocuses the target energy.
As for R 6 , when c i 2 = c j 2 , it can be simplified as:
R 6 = i = 1 K j = 1 K i j | A i | 2 | A j | 2 exp { j 4 π ε [ ( c i 1 + c j 1 ) τ 0 + 4 c i 2 τ 0 t m ] }
Substituting Equation (28) into Equation (27), we can see that R 6 can achieve energy accumulation, and the peak position is in the middle of the auto-terms, that is, the acceleration is the same as the real value, while the velocity is estimated as the average of the two targets.
Example 3.
In this example, two maneuvering targets designated as Tr1 and Tr2 have the same acceleration. Radar parameters are the same as those in Example 1. Target motion parameters are set as: A r , 1 = 1 , r 1 = 2.9 km , v 1 = 15 m / s , a 1 = 1 m / s 2 for Target Tr1; A r , 2 = 1 , r 2 = 3.1 km , v 2 = 15 m / s , a 2 = 1 m / s 2 for target Tr2.
Figure 3a shows the target trajectories after pulse compression. Figure 3b is the result of FD-SoPD. It is obvious that, in addition to the auto-terms, the cross term R 6 is also corrected as a beeline, which locates in the middle of them. Thus, the energy of R 6 is accumulated into Peak 1 in Figure 3c, and the auto-terms form Peak 2 and Peak 3. Moreover, the motion parameters of Peak 1 also confirm the theoretical analysis in Case 2.
The next step is to determine whether all of these peaks are real maneuvering targets. After phase compensation with Peak 1, the coherent integration result is shown in Figure 3d, where no target will be detected by CFAR detection. Therefore, Peak 1 belongs to a cross-term peak. In contrast, integration with Peak 2 or Peak 3 can both produce a single sharp peak in the range-Doppler domain. Thus, Peak 2 and Peak 3 belong to the auto-term peaks. This also provides us with a method for estimating the target motion parameters combined with CFAR detection. The detailed flowchart of the FD-SoPD algorithm is given in Figure 4.

4. Computational Burden Analysis of the FD-SoPD Algorithm

In this section, we will analyze the computational burden of the algorithm. The SCIFT [16], TDST [34], and ACCF-LVD [30] are selected for comparisons.
We denote the number of range cells and pulses by N r and M . For SCIFT, its main procedures include the symmetric autocorrelation function ( O ( 2 N r M log 2 M ) ), the chirp-z based SCIFT ( O ( 3 M N r log N r ) ), and FFT along the lag time axis ( O ( N r M log 2 M ) ). Therefore, the computational complexity is about O ( 3 M N r ( log 2 M + log N r ) ) .
For TDST, to complete the two steps of scaled Fourier transform (SCFT), the computational complexities are O ( 3 N r M 2 log 2 M ) and O ( 3 M 2 N r log 2 N r ) , respectively. Thus, the total computational complexity is in the order of O ( 3 N r M 2 log 2 N r M ) .
For ACCF-LVD, its main procedures include ACCF operation ( O ( 2 M N r log 2 N r ) ) and chirp-z based LVD algorithm ( O ( 3 M 2 log 2 M ) ). Therefore, its computational cost is about O ( 2 M N r log 2 N r + 3 M 2 log 2 M ) .
The implementation of the proposed algorithm needs the calculation of R S o P D ( f r , t m , τ ; τ 0 ) ( O ( 2 M 2 N r ) ), chirp-z based KT ( O ( 3 N r M log 2 M ) ), and two dimensional FFT ( O ( M N r ( log 2 N r + log 2 M ) ) ). Therefore, the overall computational cost of the proposed method is in the order of O ( 2 M 2 N r ) .
The computational complexities are listed in Table 1. Under the assumption of N r = M , Figure 5 shows the computational complexities of the above four methods. Obviously, the TDST takes too much time and is not suitable for real-time processing. Table 1 also gives the detailed values of computational resources. The TDST takes up much more memory to store the three-dimensional matrix [34]. In comparison, the SCIFT, ACCF-LVD, and FD-SoPD show great advantages in this aspect. Therefore, we could conclude that the proposed FD-SoPD cost moderates computational time and resources, which helps practical applications.

5. Numerical Results

In the section, we will give several numerical experiments to demonstrate the effectiveness of the proposed algorithm. The simulation parameters are given in Table 2.

5.1. Coherent Integration Performance

In this part, the coherent integration performance of the proposed method for a maneuvering target is evaluated. We choose the representative MTD, SCIFT, TDST, and ACCF-LVD algorithms as references. Complex zero-mean white Gaussian noise is added to radar echoes, and the SNR is set to be 5dB after compression. Figure 6a shows the target trajectory, and the result of MTD is given in Figure 6b. It is obvious that the MTD cannot integrate the target energy due to ignoring the LRM, QRM, and DFM. The velocity estimation and coherent integration results of the SCIFT are shown in Figure 6c,d, respectively. Unfortunately, the SCIFT is also invalid due to ignoring the target’s acceleration.
Figure 7a shows the integration result of LVD, where no significant peak can be found. Thus, the ACCF-LVD cannot integrate target energy in such a low SNR, as shown in Figure 7b. Figure 7c–f give respectively the parameter estimation and integration results of the FD-SoPD and TDST. Although both algorithms can accurately estimate the target motion parameters and perform coherent accumulation, the proposed algorithm has much more advantages in computational efficiency and resources.
Detailed results of parameter estimation and target detection are given in Table 3.

5.2. Detection Performance

The detection ability of the above five algorithms is evaluated combined with the CFAR detector. Assume the radar data is contaminated by the zero-mean white Gaussian noise and input SNRs after pulse compression are set as [-20:1:20] dB. 200 trials are done for each SNR value. The false alarm rate is set as P f a = 10 6 . Figure 8 shows the simulation result, where one can see that the MTD and SCIFT have the poorest detection probability due to ignoring the QRM or DFM. The adjacent cross-correlation function suffers more energy loss than the SoPD in the slow time domain [7]. Thus, the required SNR of FD-SoPD is about 4 dB less than ACCF-LVD. However, compared with TDST, the proposed algorithm suffers about 7 dB loss due to two-order bilinear transformation in SoPD. Considering the advantages of FD-SoPD, we can conclude that the proposed coherent detection algorithm achieves a good balance between the computational burden and detection ability.

5.3. Parameter Estimation Performance

We also evaluate the motion parameters estimation performance of FD-SoPD. The SNR after range compression varies from −15dB–20dB. The parameters for the radar and target are given in Table 2. two-hundred Monte Carlo simulations are performed for each SNR value. The root mean square error (RMSE) is utilized as a benchmark. The ACCF-LVD and TDST, which can estimate the velocity and acceleration of target, are selected for comparisons. Figure 9a,b show the RMSEs of the estimated velocity and acceleration. It can be seen that the TDST has the best estimation performance at the cost of huge computational burden. The performance of the proposed method is about 4dB better than those of ACCF-LVD on the input SNR threshold. However, compared with TDST, the FD-SoPD suffers from about 8dB SNR loss due to the constant delay in Equation (12). Overall, the proposed technique strikes a better balance between parameter estimation performance and computational cost.

5.4. Experimental Data Processing

In this subsection, we adopt the measured data of a DJI Phantom 3 commercial UAV to demonstrate the proposed FD-SoPD method. The data was collected in March 2017 by the National University of Defense Technology, Hunan, China. Figure 10a–c show the experimental scene, FMCW radar system, and radar antennas, respectively. Radar parameters are given in Table 4. Figure 10d shows the target trajectory after pulse compression, where the UAV moves across 7 range cells during the observation time. Figure 10e gives the parameter estimation result, where we could read the velocity and acceleration of the UAV, i.e., v ^ = 1.217 m / s and a ^ = 0.2145 m / s 2 . Finally, coherent integration of FD-SoPD can be obtained with the estimated velocity and acceleration, as shown in Figure 10f. At the same time, the integration results of SCIFT and MTD are also given in Figure 10g,h. Due to ignoring the LRM, QRM, or DFM, the target energy is distributed in the range-Doppler domain. However, the proposed method can estimate the acceleration of the target accurately. Thus, a well-focused peak is obtained, which is beneficial to target detection. Detailed results of parameter estimation and coherent integration are given in Table 5.

6. Conclusions

A frequency-domain second-order phase difference method is proposed to achieve coherent integration and parameter estimation. First, the FD-SoPD is performed to eliminate the QRM, DFM, and Doppler ambiguity simultaneously. After that, parameter estimation and coherent integration are accomplished. Compared with ACCF-LVD, the FD-SoPD could obtain better detection performance with moderate computation complexity. Simulations and experimental data processing results demonstrate the effectiveness of the proposed algorithm.

Author Contributions

K.J., T.L. and Y.Z. conceived and designed the experiments; G.L. performed the experiments; Y.W. analyzed the data and contributed analysis tools; K.J. and Y.W. wrote the paper.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Tian Jin for his useful suggestions as well as the help of the real radar data processing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this appendix, we will give the expressions of cross terms R 2 - R 8 . In order to simplify the expression form, we define
ε = c / ( f r + f c ) c i 1 = v i c i 2 = a i / 2 A i = A f r , i r e c t ( f r B ) exp ( j 4 π ( f r + f c ) c r i )
Then the cross terms R 2 - R 8 can be written as follows.
R 2 = i = 1 K j = 1 K i j | A i | 2 | A j | 2 exp { j 4 π ε [ ( c i 1 c j 1 ) τ + 2 c i 1 τ 0 + 2 ( c i 2 c j 2 ) t m τ + 4 c i 2 τ 0 t m ] }
R 3 = p = 1 K | A p | 2 exp [ j 4 π ε ( c p 1 τ + 2 c p 2 t m τ ) ] × k = 1 K l = 1 K k l A k * A l × exp { j 4 π ε [ c k 1 ( t m + τ 2 + τ 0 ) + c k 2 ( t m + τ 2 + τ 0 ) 2 c l 1 ( t m τ 2 τ 0 ) + c l 2 ( t m τ 2 τ 0 ) 2 ] }
R 4 = q = 1 K | A q | 2 exp [ j 4 π ε ( c q 1 τ + 2 c q 2 t m τ + 4 c q 2 t m τ 0 ) ] × i = 1 K j = 1 K i j A i A j * exp { j 4 π ε [ c i 1 ( t m + τ 2 ) + c i 2 ( t m + τ 2 ) 2 c j 1 ( t m τ 2 ) c j 2 ( t m τ 2 ) 2 ] }
R 5 = i = 1 K j = 1 K i j A i A j * exp { j 4 π ε [ c i 1 ( t m + τ 2 ) + c i 2 ( t m + ) 2 ] }        × exp { j 4 π ε [ c j 1 ( t m τ 2 ) + c j 2 ( t m τ 2 ) 2 ] }         × k = 1 K l = 1 K k l A k * A l exp { j 4 π ε [ c k 1 ( t m + τ 2 + τ 0 ) + c k 2 ( t m + τ 2 + τ 0 ) 2 ] }         × exp { j 4 π ε [ c l 1 ( t m τ 2 τ 0 ) + c l 2 ( t m τ 2 τ 0 ) 2 ] }
R 6 = i = 1 K j = 1 K i j | A i | 2 | A j | 2 exp { j 4 π ε [ ( c i 1 + c j 1 ) τ 0 + ( c i 2 c j 2 ) τ 0 2 ] }         × exp { j 4 π ε [ 2 ( c i 2 + c j 2 ) τ 0 t m + ( c i 2 c j 2 ) τ 0 τ ] }
R 7 = i = 1 K j = 1 K i j A i A j * exp { j 4 π ε [ c i 1 ( t m + τ 2 ) + c i 2 ( t m + τ 2 ) 2 ] }         × exp { j 4 π ε [ c j 1 ( t m τ 2 ) + c j 2 ( t m τ 2 ) 2 ] }         × k = 1 K l = 1 K k i , l j , k l A k * A l exp { j 4 π ε [ c k 1 ( t m + τ 2 + τ 0 ) + c k 2 ( t m + τ 2 + τ 0 ) 2 ] }        × exp { j 4 π ε [ c l 1 ( t m τ 2 τ 0 ) + c l 2 ( t m τ 2 τ 0 ) 2 ] }
R 8 = i = 1 K j = 1 K i j A i A j * exp { j 4 π ε [ c i 1 ( t m + τ 2 ) + c i 2 ( t m + τ 2 ) 2 c j 1 ( t m τ 2 ) c j 2 ( t m τ 2 ) 2 ] } × k = 1 K l = 1 K k = i , l j , k l A k * A l exp { j 4 π ε [ c k 1 ( t m + τ 2 + τ 0 ) + c k 2 ( t m + τ 2 + τ 0 ) 2 c l 1 ( t m τ 2 τ 0 ) c l 2 ( t m τ 2 τ 0 ) 2 ] } + i = 1 K j = 1 K i j A i A j * exp { j 4 π ε [ c i 1 ( t m + τ 2 ) + c i 2 ( t m + τ 2 ) 2 c j 1 ( t m τ 2 ) c j 2 ( t m τ 2 ) 2 ] } × k = 1 K l = 1 K k i , l = j , k l A k * A l exp { j 4 π ε [ c k 1 ( t m + τ 2 + τ 0 ) + c k 2 ( t m + τ 2 + τ 0 ) 2 c l 1 ( t m τ 2 τ 0 ) c l 2 ( t m τ 2 τ 0 ) 2 ] }

References

  1. Huang, P.H.; Liao, G.S.; Yang, Z.Z.; Xia, X.G.; Ma, J.T.; Ma, J.T. Long-Time Coherent Integration for Weak Maneuvering Target Detection and High-Order Motion Parameter Estimation Based on Keystone Transform. IEEE Trans. Signal Process. 2016, 64, 4013–4026. [Google Scholar] [CrossRef]
  2. Suo, P.C.; Tao, S.; Tao, R.; Nan, Z. Detection of high-speed and accelerated target based on the linear frequency modulation radar. IET Radar Sonar Navig. 2014, 8, 37–47. [Google Scholar] [CrossRef]
  3. Zhu, S.Q.; Liao, G.S.; Yang, D.; Tao, H.H. A New Method for Radar High-Speed Maneuvering Weak Target Detection and Imaging. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1175–1179. [Google Scholar]
  4. Li, X.L.; Sun, Z.; Yi, W.; Cui, G.L.; Kong, L.J.; Yang, X.B. Computationally efficient coherent detection and parameter estimation algorithm for maneuvering target. Signal Process. 2018, 155, 130–142. [Google Scholar] [CrossRef]
  5. Li, X.L.; Cui, G.L.; Yi, W.; Kong, L.J. Fast coherent integration for maneuvering target with high-order range migration via TRT-SKT-LVD. IEEE Trans. Aerospace Electr. Syst. 2016, 52, 2803–2814. [Google Scholar] [CrossRef]
  6. Wu, W.; Wang, G.H.; Sun, J.P. Polynomial Radon-Polynomial Fourier Transform for Near Space Hypersonic Maneuvering Target Detection. IEEE Trans. Aerospace Electr. Syst. 2018, 54, 1306–1322. [Google Scholar] [CrossRef]
  7. Zhang, J.C.; Su, T.; Zheng, J.B.; He, X.H. 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]
  8. Lao, G.C.; Yin, C.B.; Ye, W.; Sun, Y.; Li, G.J. An SAR-ISAR Hybrid Imaging Method for Ship Targets Based on FDE-AJTF Decomposition. Electronics 2018, 7, 46. [Google Scholar] [CrossRef]
  9. Lv, Y.K.; Wu, Y.H.; Wang, H.Y.; Qiu, L.; Jiang, J.W.; Sun, Y. An Inverse Synthetic Aperture Ladar Imaging Algorithm of Maneuvering Target Based on Integral Cubic Phase Function-Fractional Fourier Transform. Electronics 2018, 7, 148. [Google Scholar] [CrossRef]
  10. Lazarov, A.; Minchev, C. ISAR Geometry, Signal Model, and Image Processing Algorithms. IET Radar Sonar Navig. 2017, 11, 1425–1434. [Google Scholar] [CrossRef]
  11. Li, X.L.; Kong, L.J.; Cui, G.L.; Yi, W. A low complexity coherent integration method for maneuvering target detection. Dig. Signal Process. 2016, 49, 137–147. [Google Scholar] [CrossRef]
  12. Rao, X.; Tao, H.H.; Su, J.; Xie, J.; Zhang, X.Y. Detection of Constant Radial Acceleration Weak Target via IAR-FRFT. IEEE Trans. Aerospace Electr. Syst. 2015, 51, 3242–3253. [Google Scholar] [CrossRef]
  13. Xu, J.; Yu, J.; Peng, Y.N.; Xia, X.G. Radon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank. IEEE Trans. Aerospace Electr. Syst. 2011, 47, 1186–1202. [Google Scholar] [CrossRef]
  14. 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]
  15. Sun, G.C.; Xing, M.D.; Xia, X.G.; Wu, Y.R.; Bao, Z. Robust Ground Moving-Target Imaging Using Deramp–Keystone Processing. IEEE Trans. Geosci. Remote Sens. 2013, 51, 966–982. [Google Scholar] [CrossRef]
  16. Zheng, J.B.; Su, T.; Zhu, W.T.; He, X.H.; 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]
  17. Niu, Z.Y.; Zheng, J.B.; Su, T.; Zhang, J.C. Fast implementation of scaled inverse Fourier transform for high-speed radar target detection. Electr. Lett. 2017, 53, 1142–1144. [Google Scholar] [CrossRef]
  18. Zheng, J.B.; Su, T.; Liu, H.W.; Liao, G.S.; Liu, Z.; Liu, Q.H. 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]
  19. Sun, Z.; Li, X.L.; Yi, W.; Cui, G.L.; Kong, L.J. A Coherent Detection and Velocity Estimation Algorithm for the High-Speed Target Based on the Modified Location Rotation Transform. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2018, 11, 2346–2361. [Google Scholar] [CrossRef]
  20. Xu, J.; Yu, J.; Peng, Y.N.; Xia, X.G. 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. Xu, J.; Xia, X.G.; Peng, S.B.; Yu, J.; 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]
  22. Chen, X.L.; Guan, J.; Liu, N.B.; 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]
  23. Li, X.L.; Cui, G.L.; Yi, W.; Kong, L.J. Coherent Integration for Maneuvering Target Detection Based on Radon-Lv’s Distribution. IEEE Signal Process. Lett. 2015, 22, 1467–1471. [Google Scholar] [CrossRef]
  24. Lv, X.L.; Bi, G.A.; Wan, C.R.; Xing, M.D. Lv’s Distribution: Principle, Implementation, Properties, and Performance. IEEE Trans. Signal Process. 2011, 59, 3576–3591. [Google Scholar] [CrossRef]
  25. Lv, X.L.; Xing, M.D.; Zhang, S.H.; Bao, Z. Keystone transformation of the Wigner-Ville distribution for analysis of multicomponent LFM signals. Signal Process. 2009, 59, 791–806. [Google Scholar] [CrossRef]
  26. Zheng, J.B.; Liu, H.W.; Liu, Q.H. Parameterized Centroid Frequency-Chirp Rate Distribution for LFM Signal Analysis and Mechanisms of Constant Delay Introduction. IEEE Trans. Signal Process. 2017, 65, 6435–6447. [Google Scholar] [CrossRef]
  27. Kirkland, D. Imaging moving targets using the second-order keystone transform. IET Radar Sonar Navig. 2011, 5, 902–910. [Google Scholar] [CrossRef]
  28. Li, G.; Xia, X.G.; Peng, Y.N. Doppler Keystone Transform: An Approach Suitable for Parallel Implementation of SAR Moving Target Imaging. IEEE Geosci. Remote Sens. Lett. 2008, 5, 573–577. [Google Scholar] [CrossRef]
  29. Li, X.L.; Cui, G.L.; Yi, W.; Kong, L.J. Manoeuvring target detection based on keystone transform and Lv’s distribution. IET Radar Sonar Navig. 2016, 10, 1234–1242. [Google Scholar] [CrossRef]
  30. Li, X.L.; Cui, G.L.; Yi, W.; Kong, L.J. A Fast Maneuvering Target Motion Parameters Estimation Algorithm Based on ACCF. IEEE Signal Process. Lett. 2015, 22, 270–274. [Google Scholar] [CrossRef]
  31. Li, X.L.; Cui, G.L.; Kong, L.J.; 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]
  32. Li, X.L.; Cui, G.L.; Yi, W.; Kong, L.J. Radar Maneuvering Target Detection and Motion Parameter Estimation Based on TRT-SGRFT. Signal Process. 2017, 133, 107–116. [Google Scholar] [CrossRef]
  33. He, X.P.; Liao, G.S.; Zhu, S.Q.; Xu, J.W.; Guo, Y.F.; Wei, J.Q. Fast Non-Searching Method for Ground Moving Target Refocusing and Motion Parameters Estimation. Digital Signal Process. 2018, 79, 152–163. [Google Scholar] [CrossRef]
  34. Zheng, J.B.; Liu, H.W.; Liu, J.; Du, X.L.; 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]
  35. Zheng, J.B.; Su, T.; Zhang, L.; Zhu, W.T.; Liu, Q.H. ISAR Imaging of Targets with Complex Motion Based on the Chirp Rate–Quadratic Chirp Rate Distribution. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7276–7289. [Google Scholar] [CrossRef]
  36. Guida, M.; Longo, M.; Lops, M. Biparametric CFAR procedures for lognormal clutter. IEEE Trans. Aerospace Electr. Syst. 1993, 29, 798–809. [Google Scholar] [CrossRef]
Figure 1. Simulation results of Example 1. (a) The result of pulse compression; (b) the result of the frequency-domain second-order phase difference (FD-SoPD); (c) velocity–acceleration distribution; (d) the result of coherent integration.
Figure 1. Simulation results of Example 1. (a) The result of pulse compression; (b) the result of the frequency-domain second-order phase difference (FD-SoPD); (c) velocity–acceleration distribution; (d) the result of coherent integration.
Electronics 08 00287 g001aElectronics 08 00287 g001b
Figure 2. Simulation results of Example 2. (a) The result of pulse compression; (b) the result of FD-SoPD; (c) the velocity-Acceleration distribution; (d) coherent integration result of Tr1; (e) coherent integration result of Tr2.
Figure 2. Simulation results of Example 2. (a) The result of pulse compression; (b) the result of FD-SoPD; (c) the velocity-Acceleration distribution; (d) coherent integration result of Tr1; (e) coherent integration result of Tr2.
Electronics 08 00287 g002aElectronics 08 00287 g002b
Figure 3. Simulation results of Example 3; (a) the result of pulse compression; (b) the result of the FD-SoPD; (c) the velocity-Acceleration distribution; (d) the coherent integration result with peak 1; (e) Coherent integration result with peak 2; (f) the coherent integration result with peak 3.
Figure 3. Simulation results of Example 3; (a) the result of pulse compression; (b) the result of the FD-SoPD; (c) the velocity-Acceleration distribution; (d) the coherent integration result with peak 1; (e) Coherent integration result with peak 2; (f) the coherent integration result with peak 3.
Electronics 08 00287 g003
Figure 4. Flowchart of the proposed coherent detection algorithm.
Figure 4. Flowchart of the proposed coherent detection algorithm.
Electronics 08 00287 g004
Figure 5. Computational complexity comparison.
Figure 5. Computational complexity comparison.
Electronics 08 00287 g005
Figure 6. Coherent integration for a maneuvering target. (a) The result after pulse compression; (b) the integration result of moving target detection (MTD); (c) the velocity estimation result of the scaled inverse Fourier transform (SCIFT); (d) the integration result of the SCIFT.
Figure 6. Coherent integration for a maneuvering target. (a) The result after pulse compression; (b) the integration result of moving target detection (MTD); (c) the velocity estimation result of the scaled inverse Fourier transform (SCIFT); (d) the integration result of the SCIFT.
Electronics 08 00287 g006
Figure 7. Integration performance comparison. (a) Parameter estimation result of LVD; (b) integration result of the ACCF-LVD; (c) parameter estimation result of the proposed algorithm; (d) integration result of the proposed algorithm; (e) parameter estimation result of the TDST; (f) integration result of the TDST.
Figure 7. Integration performance comparison. (a) Parameter estimation result of LVD; (b) integration result of the ACCF-LVD; (c) parameter estimation result of the proposed algorithm; (d) integration result of the proposed algorithm; (e) parameter estimation result of the TDST; (f) integration result of the TDST.
Electronics 08 00287 g007aElectronics 08 00287 g007b
Figure 8. Detection probability of five algorithms.
Figure 8. Detection probability of five algorithms.
Electronics 08 00287 g008
Figure 9. Motion parameters estimation performance of the three methods. (a) Estimation of the root mean square error (RMSE) of velocity; (b) estimation RMSE of acceleration.
Figure 9. Motion parameters estimation performance of the three methods. (a) Estimation of the root mean square error (RMSE) of velocity; (b) estimation RMSE of acceleration.
Electronics 08 00287 g009
Figure 10. Experimental data processing results. (a) The experimental scene; (b) the FMCW radar system; (c) the radar antennas; (d) moving trajectory of the UAV; (e) parameter estimation result of FD-SoPD; (f) the coherent integration result of FD-SoPD; (g) the coherent integration result of the MTD; (h) the coherent integration result of the SCIFT.
Figure 10. Experimental data processing results. (a) The experimental scene; (b) the FMCW radar system; (c) the radar antennas; (d) moving trajectory of the UAV; (e) parameter estimation result of FD-SoPD; (f) the coherent integration result of FD-SoPD; (g) the coherent integration result of the MTD; (h) the coherent integration result of the SCIFT.
Electronics 08 00287 g010aElectronics 08 00287 g010b
Table 1. The computational burden comparisons of different algorithms.
Table 1. The computational burden comparisons of different algorithms.
MethodComputational ComplexityTime Cost (s) 1)Computational Resources
SCIFT O ( 3 M N r ( log 2 M + log N r ) ) 6.35 O ( 2 M N r )
TDST O ( 3 N r M 2 log 2 N r M ) 335.43 O ( 2 M 2 N r )
ACCF-LVD O ( 2 M N r log 2 N r + 3 M 2 log 2 M ) 6.06 O ( M N r )
FD-SoPD O ( 2 M 2 N r ) 12.24 O ( M N r )
1) The main configuration of the computer. CPU: Intel Core i7-6700HQ 2.60 GHz; RAM: 16.00G; Operating System: Windows 7; Software: Matlab 2015a.
Table 2. Simulation parameters for the radar and target.
Table 2. Simulation parameters for the radar and target.
ParametersValueParametersValue
Carrier frequency1 GHzBandwidth100 MHz
Sample frequency2 MHzPRF128 Hz
Pulse duration 2   ms Pulse number256
Initial slant range3 kmRadial velocity15 m/s
Radial acceleration1 m/s2--
Table 3. Comparisons of simulated parameter estimation and target detection results.
Table 3. Comparisons of simulated parameter estimation and target detection results.
Initial Range (km)Velocity (m/s)Acceleration (m/s2)Detection Result
MTD2.95−1.65-No target
SCIFT3.07655.72-No target
ACCF-LVD2.881043.8411.386No target
Proposed3.0015.011.031Detected
TDST3.0015.061.012Detected
Table 4. Frequency-modulated continuous-wave (FMCW) radar parameters.
Table 4. Frequency-modulated continuous-wave (FMCW) radar parameters.
Radar ParameterValueRadar ParameterValue
Carrier frequency9.5 GHzPRF50 Hz
Bandwidth1 GHzSampling frequency1 MHz
Pulse width0.0102 sCoherent time0.92 s
Transmit power25 dbmWeight7 kg
Radar length35 cmRadar width24 cm
Radar height20 cm--
Table 5. Comparisons of experimental parameter estimation results.
Table 5. Comparisons of experimental parameter estimation results.
Range CellVelocity (m/s)Acceleration (m/s2)Peak Value
MTD207−0.279- 2.573 × 10 7
SCIFT204−0.283- 4.013 × 10 7
FD-SoPD2041.2170.2145 8.089 × 10 7

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/electronics8030287

AMA Style

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(3):287. https://doi.org/10.3390/electronics8030287

Chicago/Turabian Style

Jin, Ke, Tao Lai, Yubing Wang, Gongquan Li, and Yongjun Zhao. 2019. "Coherent Integration for Radar High-Speed Maneuvering Target Based on Frequency-Domain Second-Order Phase Difference" Electronics 8, no. 3: 287. https://doi.org/10.3390/electronics8030287

APA Style

Jin, K., Lai, T., Wang, Y., Li, G., & Zhao, Y. (2019). Coherent Integration for Radar High-Speed Maneuvering Target Based on Frequency-Domain Second-Order Phase Difference. Electronics, 8(3), 287. https://doi.org/10.3390/electronics8030287

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