Next Article in Journal
Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter
Next Article in Special Issue
Performance Analysis of the Direct Position Determination Method in the Presence of Array Model Errors
Previous Article in Journal
Cell Selection Technique for Millimeter-Wave Cellular Systems with Hybrid Beamforming
Previous Article in Special Issue
Dual-Channel Cosine Function Based ITD Estimation for Robust Speech Separation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals

College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1460; https://doi.org/10.3390/s17061460
Submission received: 10 April 2017 / Revised: 13 June 2017 / Accepted: 14 June 2017 / Published: 21 June 2017

Abstract

:
For targets with complex motion, such as ships fluctuating with oceanic waves and high maneuvering airplanes, azimuth echo signals can be modeled as multicomponent quadratic frequency modulation (QFM) signals after migration compensation and phase adjustment. For the QFM signal model, the chirp rate (CR) and the quadratic chirp rate (QCR) are two important physical quantities, which need to be estimated. For multicomponent QFM signals, the cross terms create a challenge for detection, which needs to be addressed. In this paper, by employing a novel multi-scale parametric symmetric self-correlation function (PSSF) and modified scaled Fourier transform (mSFT), an effective parameter estimation algorithm is proposed—referred to as the Two-Dimensional product modified Lv’s distribution (2D-PMLVD)—for QFM signals. The 2D-PMLVD is simple and can be easily implemented by using fast Fourier transform (FFT) and complex multiplication. These measures are analyzed in the paper, including the principle, the cross term, anti-noise performance, and computational complexity. Compared to the other three representative methods, the 2D-PMLVD can achieve better anti-noise performance. The 2D-PMLVD, which is free of searching and has no identifiability problems, is more suitable for multicomponent situations. Through several simulations and analyses, the effectiveness of the proposed estimation algorithm is verified.

1. Introduction

In the past decades, radar signal processing has been a very important research field [1,2,3], especially estimating radar parameters [4,5,6]. The parameters of radar signals mainly include frequency [7], phase [8], and direction-of-arrival [9,10,11]. In this paper, we only focus on research into the frequency of the quadratic frequency modulation (QFM) signal [2,4]. Recently, many researchers have paid much attention to presenting effective inverse synthetic aperture radar (ISAR) imaging algorithms, which are important for high-resolution ISAR imaging and recognition of moving targets [12,13]. In general, the primary steps for ISAR imaging are range alignment [14,15] and phase adjustment [16,17], and translational-induced phase error correction. Then, the conventional range-Doppler (RD) method is utilized to reconstruct a well-focused ISAR image. However, in practical application, as the target’s motions are not smooth, the conventional RD method is not valid. For slow maneuvering targets, the azimuth echoes can be modeled as a linear frequency modulated (LFM) signal [18,19,20,21], and as the Doppler frequency shifts, we cannot obtain high-resolution imaging using the conventional RD. The chirp rate is identified as the cause of the target image defocus. Recently, because of this situation, many methods have been proposed, such as the fractional Fourier Transform (FrFT) [18], the stretch keystone-Wigner transform [19], the LPP-Hough transform [20], and the Lv’s distribution (LVD) [21]. However, for targets with complex motion, such as high maneuvering airplanes and ships fluctuating with oceanic waves, the azimuth echoes can be modeled as QFM signals, and methods for LFM signals are not suited for QFM signals. For ISAR imaging based on QFM signals, the chirp rate and the quadratic chirp rate induces a Doppler frequency shift and deteriorates the ISAR image.
Many algorithms about the parameter estimation of QFM signals for targets with complex motion have been proposed, and they generally fall into two categories: noncorrelation algorithms and correlation algorithms. Noncorrelation algorithms include the maximum likelihood algorithm [22], the quantification-based method [23], and the discrete chirp Fourier transform [24]. Correlation algorithms include the higher order ambiguity function (HAF) [2], the product higher order ambiguity function [25], the cubic phase function (CPF) [26,27], the product generalized cubic phase function [28], the product high-order matched phase transform (PHMT) [29], the HAF-integrated CPF (HAF-ICPF) [30], the scaled Fourier transform (SCFT)-based algorithm [31], the chirp rate-quadratic chirp rate distribution (CRQCRD) [32], the Modified Lv’s Distribution (MLVD) [33], and the new scaled Fourier transform (SCFT)-based algorithm [34]. Although noncorrelation algorithms can get high anti-noise performance with a lower signal-to-noise ratio (SNR) and have no cross-influence, they are not suitable for the QFM signal due to the high computation cost ( O ( M 2 N log 2 N ) ) [22,23,24], where N is the effective length of the slow time in ISAR, and M is the number of searching points, which is always greater than N . Thus, noncorrelation algorithms are worse in practical application. Compared to noncorrelation algorithms, although correlation algorithms are influenced by cross terms and have lower anti-noise performances, they have lower computational complexity and have advantages in radar system complexity. Thus, correlation algorithms are more practical than noncorrelation algorithms.
Through analyses of implementation, cross terms, error propagation, anti-noise performance and computational cost, the HAF-ICPF [30], the MLVD [33] and the algorithm presented in [34] will outperform other correlation algorithms [2,22,23,24,25,26,27,28,29,30,31,32,33,34]. Therefore, these algorithms are chosen to compare with the proposed estimation algorithm. However, some disadvantages still exist in these algorithms, as follows: (1) Due to requiring one-dimensional brute-force searching, the computation cost of the HAF-ICPF algorithm is still high. In addition, error propagation exists in HAF-ICPF, because it estimates the chirp rate and the quadratic chirp rate sequentially; (2) The algorithm presented in [34] overcomes brute-force searching and error propagation, but the algorithm does not have a high anti-noise performance, as redundancy information is not used; (3) Although the MLVD improves anti-noise performance by using redundancy information, it is still poor. In addition, the MLVD suffers from the identifiability problem when dealing with multiple components’ QFM signals with the same cubic order phase coefficients or the same cubic and quadratic order phase coefficients.
In this paper, according to the analyses above and the MLVD for the QFM signal, a novel algorithm, which is called 2D-PMLVD, is presented. The proposed estimation algorithm is based on a novel multi-scale parametric symmetric self-correlation function (PSSF) and the idea of keystone transform. By means of this processing, the QFM signal is transformed into different 2-D frequency domains, but at the same coordinate position. By multiplying all 2-D frequency domains, the proposed algorithm not only suppresses cross terms, but also improves anti-noise performance. Through analyses of the implementation, the 2D-PMLVD inherits advantages of the MLVD. Due to there being no brute-force searching, and its being implemented only using the FFT, the IFFT, and complex multiplication, the proposed estimation algorithm is easily implemented with low computation cost. Therefore, the 2D-PMLVD is more suitable in practical application.
The remainder of the paper is organized as follows. Section 2 introduces the model of QFM signals, and the principle of 2D-PMLVD for mono-QFM signals and multi-QFM signals. In addition, the performance of the 2D-PMLVD for multicomponent QFM signals is analyzed. Section 3 presents the details of the implementation in the discrete-time form. The anti-noise performance and the computational cost of the 2D-PMLVD are analyzed in Section 4. A conclusion is given in Section 5.

2. 2D-PMLVD with QFM Signal

In this section, the model of the QFM signals and the principle of the 2D-PMLVD are introduced, respectively.
For ISAR imaging of targets with complex motions, the geometry used here is based on the model in [31,32]. In this paper, range alignment and translational induced phase error correction, which can be referred to in [14,15,16,17], will not be discussed in detail. We will only focus on the processing of the Doppler frequency shift. Suppose that the radar transmits the LFM signal. According to the analyses in [31,32], after the range alignment and the translational induced phase error correction, azimuth echoes of the l th range cell take the form of multicomponent QFM signals
s l ( t ) = p = 1 P A p exp [ j 2 π ( ϕ 1 , p t + 1 2 ϕ 2 , p t 2 + 1 6 ϕ 3 , p t 3 ) ] + z ( t )
where P is the phase of the p th component, t denotes the slow time and its sampling frequency is the same as the pulse repetition frequency (PRF), z ( t ) denotes the additive complex white Gaussian noise with a variance of δ 2 , ϕ 1 , p , ϕ 2 , p , ϕ 3 , p denote the centroid frequency, the chirp rate and the quadratic chirp rate, respectively.

2.1. 2D-PMLVD with Mono-QFM Signal

In order to facilitate the understanding of the algorithm, we first consider a noise-free QFM signal. It is represented as
s ( t ) = A exp [ j ϕ ( t ) ] = A exp [ j 2 π ( ϕ 1 t + 1 2 ϕ 2 t 2 + 1 6 ϕ 3 t 3 ) ]
where ϕ ( t ) is the phase function ϕ 1 , ϕ 2 and ϕ 3 denote the centroid frequency, the chirp rate and the quadratic chirp rate, respectively.
According to analyses in [21,33], a novel multi-scale PSSF is defined as
R s ( t , τ , D ) = s ( t + τ + a 2 + D ) s ( t τ + a 2 D ) s * ( t + τ + a 2 D ) s * ( t τ + a 2 + D )
where a denotes a constant time-delay, τ denotes the lag time variable, and * denotes the complex conjugation. D = { d i } i = 1 L ( L > 1 ) is the set of scale factors, and the selection criteria will be presented in Section 2.3.1. According to [33,34], the scale factor d i and the corresponding lag-time τ are utilized to complete the order reduction. Besides that, d i can reduce the number of self-correlations, which benefits the anti-noise performance.
Substituting (2) into (3) yields
R s ( t , τ , D ) = A 4 exp { j 2 π [ 2 ϕ 2 D ( τ + a ) + 2 ϕ 3 D t ( τ + a ) ] }
It is easily seen from Equation (4) that the time variable t and lag variable τ couple with each other in exponential phase terms. For every scale factor, Equation (4) suffers the same influence of the coupling; we only need to analyze Equation (4) with one scale factor to explain the influence. Thus, Equation (4) is represented as
R s ( t , τ , d i ) = A 4 exp { j 2 π [ 2 ϕ 2 d i ( τ + a ) + 2 ϕ 3 d i t ( τ + a ) ] }
The form of Equation (5) is the same as the self-correlation function defined for the MLVD [33] and LVD [21]. According to [21,33], we have known that the coupling influence energy accumulation. In order to analyze the influence of the coupling for the energy accumulation of the proposed algorithm, we use the energy accumulation method presented in [21,33], referred to as EA. The EA method is also the basic energy accumulation method used for the proposed algorithm. The EA method is divided into two steps. Firstly, perform FFT on Equation (5) with respect to t . Secondly, perform FFT on the first step with respect to τ . After the first step, the signal energy will peak along the inclined line f t = 2 ϕ 3 d i ( τ + a ) . It is easily seen that the signal energy accumulation cannot be accomplished by the second step. Thus, in order to effectively accumulate the energy of the signal, the coupling needs to be removed. Here, the idea of keystone transformation is borrowed to rescale the time axis, which is defined as
Κ ( Φ ( t , τ ) ) Φ ( t n h ( τ + a ) , τ )
where a is a constant time-delay and h is a scaling factor different from d i . The choice of parameters a and h will be presented in Section 2.3.1. Φ ( t , τ ) is the phase function with respect to ( t , τ ) . Κ is the scaling operator. t n is denotes scaled time.
Applying the scaling operator Κ to (5), we have
K [ R x ( t , τ , d i ) ] = A 4 exp { j 2 π [ 2 ϕ 2 d i ( τ + a ) + 2 ϕ 3 d i h t n ] } = A 4 exp { j 4 π ϕ 2 d i a } exp [ j 2 π ( 2 ϕ 2 d i τ + 2 ϕ 3 d i h t n ) ]
It is easily seen from (7) that the coupling has been removed by introducing the new time variable t n . Then, by applying the first step of EA to Equation (7), the signal energy becomes a beeline f t n = 2 ϕ 3 d i / h . Thus, the signal energy can be accumulated by the second step. The process can be expressed as
L x ( f t n , i , f τ , i , d i ) = F F T τ ( F F T t n ( Κ [ R s ( t , τ , d i ) ] ) ) = A 4 exp ( j 4 π ϕ 2 a d i ) δ ( f τ , i 2 ϕ 2 d i ) δ ( f t n , i 2 ϕ 3 d i / h )
where δ ( · ) denotes the Dirac delta function, F F T ( · ) is the Fast Fourier transform operator. Equation (8) is similar to MLVD [33]. It can be easily seen from Equation (8) that a sole peak appears on the 2-D frequency domain f τ f t n . The coordinate of the peak is ( f τ , i = 2 ϕ 2 d i , f t n , i = 2 ϕ 3 d i / h ) . According to the coordinate of the peak, the parameters ϕ 1 , ϕ 2 , and ϕ 3 can be estimated.
However when the EA method is applied to the multi-scale PSSF, the process is represented as
L x ( f t n , f τ , D ) = τ ( t n ( Κ [ R s ( t , τ , D ) ] ) ) = A 4 exp ( j 4 π ϕ 2 a D ) δ ( f τ 2 ϕ 2 D ) δ ( f t n 2 ϕ 3 D h )
It is easily seen from Equation (9) that the mono-component signal will be translated into multi-peaks on f τ f t n . The number of peaks is the same as the number of scale factors. Suppose that d 1 is the reference factor, and the peak which is correspond to d 1 is the reference peak. The relationships between reference peak and other peaks are represented as
P 1 , i = f τ , 1 f τ , i = f t n , 1 f t n , i = d 1 d i ( 1 i L )
where P 1 , i is zoom factor. f τ , 1 and f t n , 1 are the coordinates of the reference peak on the 2-D frequency domain. f τ , i and f t n , i are the coordinates of the peak which corresponds to d i . It can be easily seen from Equation (10) that scale factor influences the location of the peak. Due to the values of scale factors have been selected by the selection criteria, which will be presented in Section 2.3.2, the relationships in Equation (10) can be easily obtained. According to these relationships, an improved energy accumulation method is presented, which is called m S F T t , i F F T τ , i . The m S F T t , i F F T τ , i not only inherits advantages of decoupling but also improves the performance of energy accumulation. By using the m S F T t , i F F T τ , i , the peaks of the QFM signal are transformed into different 2-D frequency domains, but at the same position under different scale factors. After multiplying all the 2-D frequency domains, a higher peak will appear. The m S F T t , i F F T τ , i will be presented in detail in Section 3. The processing is named 2D-PMLVD, which is represented as
Q 2 D P M L V D = i = 1 L F F T τ , i ( m S F T t , i [ R s ( t , τ , d i ) ] )
2D-PMLVD not only suppresses cross terms, but also improves the anti-noise ability, which will be proved in Section 2.2 and Section 4.1.
2D-PMLVD is divided into four steps. Firstly, select one of the scale factors as the reference one and calculate the relationships between reference factor and other factors by Equation (10). Secondly, transform the QFM signal into different 2-D frequency domains based on m S F T t , i F F T τ , i . Thirdly, multiply all the 2-D frequency domain to obtain the higher peak. Fourthly, detect the peak value to estimate the signal parameters.
For example, suppose that the set of scales contains two factors which are represented as D = { d 1 , d 2 } . According to Equation (10), we get
P 1 , 2 = d 1 / d 2
where d 1 is the reference scale. After applying m S F T t , i F F T τ , i to the QFM signal, two peaks will appear on different 2-D frequency domains. The reference peak will appear on f τ , 1 f t n , 1 , and the other peak will appear on f τ , 2 f t n , 2 . However, due to the m S F T t , i F F T τ , i operator, the relation between the different frequency domain is represented as
f τ , 2 f t n , 2 = ( f τ , 1 f t n , 1 ) / P 1 , 2
In order words, the new frequency interval becomes [ 0 , 1 / P 1 , 2 ] . Thus, because of the transformation, the two peaks appear at the same position. Then, a higher peak is obtained by multiplying the different frequency domains. According to the analysis in [25], the proposed algorithm can be improved when the number of scale factors is increased. Then, the computational cost will increase simultaneously. Hence, the number of scale factors will not be too large. In this paper, the number is selected as 2. By this selection, the cross terms can be effectively suppressed which are shown in Figure 1 and Figure 2.

2.2. 2D-PMLVD with Multi-QFM Signals

In this part, 2D-PMLVD with multi-QFM signals is analyzed. The proposed algorithm not only improves the ability to suppress cross terms, but also improves the performance of the EA method. The EA method, which is presented in MLVD [33] and LVD [21], has been proved to be effective for suppressing cross terms. Now, model the noise-free multi-QFM signals as
s m u l t i ( t ) = p = 1 P A p exp [ j ϕ ( t ) ] = p = 1 P A p exp [ j 2 π ( ϕ 1 , p t + 1 2 ϕ 2 , p t 2 + 1 6 ϕ 3 , p t 3 ) ]
To formulate the cross term problem arising from multi-QFM signals, we first consider two signals ( P = 2 ) in Equation (14), which is represented as
s m u l t i ( t ) = p = 1 2 A p exp [ j ϕ ( t ) ] = A 1 exp [ j 2 π ( ϕ 1 , 1 t + 1 2 ϕ 2 , 1 t 2 + 1 6 ϕ 3 , 1 t 3 ) ] + A 2 exp [ j 2 π ( ϕ 1 , 2 t + 1 2 ϕ 2 , 2 t 2 + 1 6 ϕ 3 , 2 t 3 ) ]
For ease of understanding, the multi-scale PSSF can be written as
R s ( t , τ , D ) = R ( t + τ + a 2 , D ) R * ( t τ + a 2 , D )
where
R ( t , D ) = s ( t + D ) s * ( t D )
Submitting (14) into (16b), we obtain
R ( t , D ) = R a u t o ( t , D ) + R c r o s s ( t , D )
where R a u t o ( t , D ) and R c r o s s ( t , D ) denote auto-terms and cross terms with different scale factors, respectively. In order to describe the shortcoming of the EA method and explain the principle of the energy accumulation method in the proposed algorithm, we firstly consider one scale factor in Equation (17). It is represented as
R ( t , d i ) = R a u t o ( t , d i ) + R c r o s s ( t , d i )
Submitting (14) into (18a),
R a u t o ( t , d i ) = A 1 2 exp [ j 2 π ( 2 ϕ 1 , 1 d i + 1 3 ϕ 3 , 1 d i 3 ) ] × exp [ j 2 π ( 2 ϕ 2 , 1 d i t + ϕ 3 , 1 d i t 2 ) ] + A 2 2 exp [ j 2 π ( 2 ϕ 1 , 2 d i + 1 3 ϕ 3 , 2 d i 3 ) ] × exp [ j 2 π ( 2 ϕ 2 , 2 d i t + ϕ 3 , 2 d i t 2 ) ] ,
R c r o s s ( t , d i ) = A 1 A 2 exp [ j 2 π ( ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 2 ϕ 2 , 1 d i 2 1 2 ϕ 2 , 2 d i 2 + 1 6 ϕ 3 , 1 d i 3 + 1 6 ϕ 3 , 2 d i 3 ) ] × exp { j 2 π [ ( ϕ 1 , 1 ϕ 1 , 2 + ϕ 2 , 1 d i + ϕ 2 , 2 d i + 1 2 ϕ 3 , 1 d i 2 1 2 ϕ 3 , 2 d i 2 ) t + 1 2 ( ϕ 2 , 1 ϕ 2 , 2 + ϕ 3 , 1 d i + ϕ 3 , 2 d i ) t 2 + 1 6 ( ϕ 3 , 1 ϕ 3 , 2 ) t 3 ] } + A 1 A 2 exp { j 2 π [ ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 2 d i 2 ( ϕ 2 , 2 ϕ 2 , 1 ) + 1 6 d i 3 ( ϕ 3 , 1 + ϕ 3 , 2 ) ] } × exp { j 2 π [ ( ϕ 1 , 2 ϕ 1 , 1 + ϕ 2 , 1 d i + ϕ 2 , 2 d i + 1 2 ϕ 3 , 2 d i 2 1 2 ϕ 3 , 1 d i 2 ) t + 1 2 ( ϕ 2 , 2 ϕ 2 , 1 + ϕ 3 , 1 d i + ϕ 3 , 2 d i ) t 2 + 1 6 ( ϕ 3 , 2 ϕ 3 , 1 ) t 3 ] } .
According to [21], Equation (16a) is similar to the form of the self-correlation function defined for the LVD. Thus, under multi-QFM signals, if Equation (18c) takes the form of the LFM signal, the cross terms will accumulate as an auto-term. Let Δ ϕ 1 = ϕ 1 , 1 ϕ 1 , 2 , Δ ϕ 2 = ϕ 2 , 1 ϕ 2 , 2 , and Δ ϕ 3 = ϕ 3 , 1 ϕ 3 , 2 . LFM signals will exist in the cross terms as in the following two cases [32,35]:
  • Δ ϕ 1 0 , Δ ϕ 2 0 , and Δ ϕ 3 = 0 .
  • Δ ϕ 1 0 , Δ ϕ 2 = 0 , and Δ ϕ 3 = 0 .
In these two cases, the spurious peaks will appear on f τ f t n , which greatly interferes with the performance of the detection. However, 2D-PMLVD overcomes the problem by introducing the multi-scale PSSF and m S F T t , i F F T τ , i . This will be analyzed in detail as follows.

2.2.1. Case One

When Δ ϕ 1 0 , Δ ϕ 2 0 , and Δ ϕ 3 = 0 , the R c r o s s ( t , d i ) is represented as
R c r o s s ( t , d i ) = A 1 A 2 exp { j 2 π [ ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 2 d i 2 ( ϕ 2 , 1 ϕ 2 , 2 ) + 1 3 ϕ 3 , : d i 3 ] } × exp { j 2 π [ ( ϕ 1 , 1 ϕ 1 , 2 + ϕ 2 , 1 d i + ϕ 2 , 2 d i ) t + 1 2 ( ϕ 2 , 1 ϕ 2 , 2 + 2 ϕ 3 , : d i ) t 2 ] } + A 1 A 2 exp { j 2 π [ ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 2 d i 2 ( ϕ 2 , 2 ϕ 2 , 1 ) + 1 3 ϕ 3 , : d i 3 ] } × exp { j 2 π [ ( ϕ 1 , 2 ϕ 1 , 1 + ϕ 2 , 1 d i + ϕ 2 , 2 d i ) t + 1 2 ( ϕ 2 , 2 ϕ 2 , 1 + 2 ϕ 3 , : d i ) t 2 ] }
where ϕ 3 , : = ϕ 3 , 1 = ϕ 3 , 2 . It can be easily seen from Equation (19) that R c r o s s ( t , d i ) takes the form of two LFM signals. This means that the cross terms will accumulate two spurious peaks. The coordinates of the spurious peaks are
( f τ , i C r 1 = ϕ 1 , 1 ϕ 1 , 2 + ϕ 2 , 1 d i + ϕ 2 , 2 d i , f t n , i C r 1 = ϕ 2 , 1 ϕ 2 , 2 + 2 ϕ 3 , : d i )   and   ( f τ , i C r 2 = ϕ 1 , 2 ϕ 1 , 1 + ϕ 2 , 1 d i + ϕ 2 , 2 d i , f t n , i C r 2 = ϕ 2 , 2 ϕ 2 , 1 + 2 ϕ 3 , : d i )
It can be easily seen from Equation (18b) that R a u t o ( t , d i ) takes the form of two LFM signals. The coordinates of the true peaks are
( f τ , i A u 1 = 2 ϕ 2 , 1 d i , f t n , i A u 1 = 2 ϕ 3 , 1 d i )   and   ( f τ , i A u 2 = 2 ϕ 2 , 2 d i , f t n , i A u 2 = 2 ϕ 3 , 2 d i )
For different scale factors, the coordinates of the true peaks meet Equation (10), and can be represented as
P i , j = f τ , i A u 1 f τ , j A u 1 = f t n , i A u 1 f t n , j A u 1 = f τ , i A u 2 f τ , j A u 2 = f t n , i A u 2 f t n , j A u 2
However, through the analysis of Equation (20), the coordinates of the spurious peaks do not meet Equation (10). In other words, the coordinates are not proportional to the scale factor. Thus, the spurious peaks cannot be transformed into the same position by the operator m S F T t , i F F T τ , i . After multiplying all the frequency domain, the spurious peaks will not be enhanced. Compared with the true peaks, the spurious peaks are “suppressed”. The result is based on the premise that the scale factors are different. This is because if the scale factors are equivalent, the spurious peaks are in the same position, which means that the cross terms cannot be suppressed. This is proved by Example 1. The selection criterion for the scale factor is introduced in Section 2.3.2.
Example 1.
We consider two QFM signals P 1 and P 2 . Signal parameters are set as follows: A 1 = 1 , ϕ 1 , 1 = 89 , ϕ 2 , 1 = 42 , ϕ 3 , 1 = 4 for P 1 ; and A 2 = 1 , ϕ 1 , 2 = 10 , ϕ 2 , 2 = 20 , ϕ 3 , 2 = 4 for P 2 . The sampling frequency F s is set to 256 Hz. The effective signal length N is equal to 256. Both a and h are set to one. The quadratic chirp rate of P 1 and P 2 are equal to each other. Figure 1a shows the results of MLVD with a scale factor equal to 56. A u 1 and A u 2 denote the true peaks of the auto-terms. Because ϕ 3 , 2 = ϕ 3 , 1 , two spurious peaks, denoted by C r 1 and C r 2 , exist in the 2-D frequency domain alongside the true peaks. It can be easily seen that the spurious peaks are even higher than the true ones in Figure 1a. This creates a serious challenge for signal detection. Figure 1b shows the results of 2D-PMLVD with two scale factors equal to 56 and 64. Because the multi-scale PSSF and the operator m S F T t n , i F F T τ , i were applied, the spurious peaks of the cross terms are suppressed, and the true peaks of the auto-terms are greatly strengthened. Figure 1c,d gives the contour of Figure 1a,b, respectively.

2.2.2. Case Two

When Δ ϕ 1 0 , Δ ϕ 2 = 0 , and Δ ϕ 3 = 0 , R a u t o ( t , d i ) is represented as
R a u t o ( t , d i ) = { A 1 2 exp [ j 2 π ( 2 ϕ 1 , 1 d i + 1 3 ϕ 3 , : d i 3 ) ] + A 2 2 exp [ j 2 π ( 2 ϕ 1 , 2 d i + 1 3 ϕ 3 , : d i 3 ) ] } × exp [ j 2 π ( 2 ϕ 2 , : d i t + ϕ 3 , : d i t 2 ) ]
where ϕ 2 , : = ϕ 2 , 1 = ϕ 2 , 2 . And R c r o s s ( t , d i ) is represented as
R c r o s s ( t , d i ) = A 1 A 2 exp [ j 2 π ( ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 3 ϕ 3 , : d i 3 ) ] × exp { j 2 π [ ( ϕ 1 , 1 ϕ 1 , 2 + 2 ϕ 2 , : d i ) t + ϕ 3 , : d i t 2 ] } + A 1 A 2 exp [ j 2 π ( ϕ 1 , 1 d i + ϕ 1 , 2 d i + 1 3 d i 3 ϕ 3 , : ) ] × exp { j 2 π [ ( ϕ 1 , 2 ϕ 1 , 1 + 2 ϕ 2 , : d i ) t + ϕ 3 , : d i t 2 ] }
Unlike case one, R a u t o ( t , d i ) only takes the form of one LFM signal in case two. That is because the two LFM signals are the same. It means that, by the proposed algorithm, only a higher peak of the auto-terms is obtained. The coordinates of the true peak are ( 2 ϕ 2 , : d i , 2 ϕ 3 , : d i ) . It can be easily seen that the coordinates are proportional to the scale factor.
It can be easily seen from Equation (24) that R c r o s s ( t , d i ) also takes the form of two LFM signals. The coordinates of the two peaks are ( ϕ 1 , 2 ϕ 1 , 1 + 2 ϕ 2 , : d i , 2 ϕ 3 , : d i ) and ( ϕ 1 , 1 ϕ 1 , 2 + 2 ϕ 2 , : d i , 2 ϕ 3 , : d i ) , respectively. Thus, the coordinates are not proportional to the scale factor. This is the same as in case one.
The process for suppressing the cross terms is the same as in case one. After applying 2D-PMLVD to case two, the spurious peaks can effectively be suppressed and the identifiability is improved. This is shown in Example 2.
Example 2.
We consider two QFM signals P 3 and P 4 . Signals parameters are set as follows: A 3 = 1 , ϕ 1 , 3 = 10 , ϕ 2 , 3 = 20 , ϕ 3 , 3 = 10 for P 3 ; and A 4 = 1 , ϕ 1 , 4 = 30 , ϕ 2 , 4 = 20 , ϕ 3 , 4 = 10 for P 4 . The sampling frequency F s is set to 256 Hz. The effective signal length N is equal to 256. Both a and h are set to one. The quadratic chirp rate and chirp rate of P 3 and P 4 are the same. Figure 2a shows the results of MLVD with a scale factor equal to 56. Because the quadratic chirp rate and chirp rate are the same, the auto-term is accumulated into one peak, which is represented as A u 3 ( A u 4 ) . C r 3 and C r 4 denote the spurious peaks. Figure 2b shows the results of 2D-PMLVD with two scale factors. By means of 2D-PMLVD, the spurious peaks of the cross terms are suppressed, and the true peak is greatly strengthened. Figure 2c,d gives the contour of Figure 2a,b, respectively.

2.3. Parameter Selection and Resolution

In this part, the selection criteria of the parameters a , h and scale factors are presented. In practical applications, the QFM signal is of finite-length and discrete. This means that the selection criteria only need to be discussed in a finite-length discrete setting.

2.3.1. Selection of a and h

We know that the optimal value of a h is equal to one [21,36]. Define the discrete parameter pair ( a , h ) as ( m / f s , h ) . In order to directly read the quadratic chirp rate from the 2-D frequency domain, the optimal value of h is set to one. According to the length of redundancy, the selection of the parameters follows two criteria.
Criterion 1: According to the analysis in [21], as long as the signal is long enough to guarantee at least 1 s redundancy, that is to say N a l l / f s N 2 D P M L V D / f s 1 , the parameters are
( m = f s ,   h = 1 )
N a l l and N 2 D P M L V D denote the length of the signal and the effective length of the signal, respectively. In this case, the parameters of a and h are equal to one, and the previously processed data can be used as redundancy information to perform the proposed algorithm. Because 2D-PMLVD includes a scale factor set that is different from the LVD, the effective length of 2D-PMLVD is shorter than that of LVD. It means that 2D-PMLVD is much more suited to fulfilling the condition of redundancy.
Criterion 2: However, in practical applications, the total sampling time may be very short, even less than 1 s, or the frequencies of the chirp signal may vary quickly with time even if the length of signal is long enough. The parameter selection criterion is not satisfied by the above conditions. The new condition is N a l l / f s N 2 D P M L V D / f s < 1 . To guarantee the condition that a is equal to one, the parameters are represented as
( m = N a l l N 2 D P M L V D ,   h = f s / ( N a l l N 2 D P M L V D ) )
The derivation of Equation (26) is in accordance with [21], which aims at making h as close as possible to one under the condition of a h = 1 .

2.3.2. Selection of d i

According to the analysis in [33], we know that the scale factors not only reduce order, but also benefit anti-noise performance. Thus, their selection is very important for 2D-PMLVD. According to the analysis in [37], we know that the optimal di value is obtained by minimizing the values for the MSE of chirp rate and quadratic chirp rate. This is represented as
d i 0.089 N a l l
where denotes the round up operator. We have seen in Section 2.2 that, in order to suppress the cross terms, the scale factors in D will be different. When we move away from the optimal choice, we will certainly obtain a sub-optimal set. As far as accuracy is concerned, it would be better to have the effective signal length be as close to the optimal length as possible, which suggests that the selection of the other scale factors in the set D should be close to d i o p t .
For a mono-QFM signal, the values of scale factors may be equivalent, as there is no cross term. However, in actual processing, we do not know the QFM signal is mono-component or multi-component. Thus, the scale factors in D are different.

2.3.3. Resolution

In the proposed algorithm, which is similar to the frequency resolution of Fourier transform, the resolution for representing the parameters of CR and QCR relate to the signal length, which makes sense, since the 1-D QFM signals are converted into the 2-D single frequency signals in the CR-QCR domain. In this paper, due to the sampling scheme of Claasen and Mecklenbräuker (CM) and the finite-length effect, the resolution for representing the parameters of CR and QCR are different. Moreover, the resolutions also relate to the scale factors. We have seen that the multi-scale PSSF can be divided into Equation (16a,b), and that Equation (16a) is the same as PSIAF, which is presented in [21]. Thus, some theorems in [21] can also be applied to this paper. Applying Equation (16b) to a QFM signal, the auto-terms can be represented as different LFM signals with different scale factors and the centroid frequency and the chirp rate of the i th LFM signal can be expressed as
f i = 2 d i φ 2 , p   and   γ i = 2 d i φ 3 , p
Suppose Δ denotes sampling interval. According to the Appendix B in reference [21], we have known that the resolutions of centroid frequency and the chirp rate are expressed as
δ f i = 1 / N i Δ   and   δ γ i = 2 / ( h N i Δ )
where N i denotes the length of the i th LFM signal; and N i is represented as N i = N a l l 2 d i , where N a l l denotes the length of the QFM signal. According to Equations (28) and (29), the corresponding resolution for φ 2 , p and φ 3 , p is represented as
δ φ 2 , p = 1 / ( 2 d i N i Δ )   and   δ φ 3 , p = 1 / ( h d i N i Δ )
It can be easily seen from Equation (30) that different scale factors correspond to different resolutions; and the resolutions are influenced by the worst resolutions in the proposed algorithm. Thus, the resolutions are influenced not only by the signal length but also by the scale factors.

3. Implementation

In this section, the discrete implementation of 2D-PMLVD will be discussed by using the sampling scheme of Claasen and Mecklenbräuker (CM) [38]. The discrete form of (4) expressed as
R s ( t , τ , D ) = A 4 exp { j 2 π [ 2 ϕ 2 D ( 2 m T s + a ) + 2 ϕ 3 D n T s ( 2 m T s + a ) ] }
The key step of the implementation is the scaling transform. Similar to the keystone transform, the operator Κ can be fulfilled by using DFT-IFFT or SFT-IFFT, which are analyzed in detail in [19]. According to the analysis in [19], the SFT-IFFT transformation only uses complex multiplications and FFT based on the scaling principle, and is more suitable for implementation than DFT-IFFT. If the Fourier transform is represented as
S ( f ) = ( r ( t ) ) = r ( t ) exp ( j 2 π f t ) d t
then the representation of SFT is expressed as
S ( λ , f ) = ¯ ( r ( t ) ) = r ( t ) exp ( j 2 π λ f t ) d t
where λ = ( 2 m T s + a ) h / N .
According to the analysis in Section 2.1 and Section 2.2, the energy of the auto-terms is expected to be accumulated into the same position under different scale factors. However, the EA method, which is presented in [21,33], is only valid for mono-scale factors, and is invalid for multi-scale factors. According to the analysis of SFT, an improved energy accumulation method is proposed, called modified SFT-FFT ( m S F T t , i F F T τ , i ) . The operator is divided into two parts: a performing m S F T t , i [ ] operator along the t axis and a performing F F T τ , i [ ] operator along the τ axis.
The m S F T t , i [ ] is different from Equation (33), which is represented as
m S F T t , i [ R ( t , τ , d i ) ] = R ( t , τ , d i ) exp ( j 2 π λ f i Δ t ) d t
where R ( t , τ , d i ) denotes the PSSF when scale factor is d i . f i Δ = f 1 / P 1 , i ( 1 i L ) , where f 1 denotes the reference frequency domain and f i Δ denotes the i th frequency domain, which corresponds to the i th scale factor. Compared with the interval of f 1 , the interval of f i Δ becomes [ 0 , 1 / P 1 , i ] by introducing P 1 , i . The variable λ in m S F T t , i [ ] is the same as in Equation (33). Thus, the coupling between t and τ can also be removed by the m S F T t , i [ ] operator. Based on Equation (10), we know that P 1 , i denotes the relationships between reference energy peak position and the i th energy peak position. By using the zoom factor P 1 , i , the m S F T t , i [ ] operator changes the frequency interval. Then, under different scale factors, the signal energy of the auto-term is accumulated at the same position on the f t n axis. In other words, under different scale factors, the signal energy of the auto-term is accumulated into the same beeline which is represented as f t n = X ( τ ) , where X ( τ ) denotes the function of τ . Then, applying the F F T τ , i [ ] operator along the τ axis, the signal energy of the auto-term is accumulated as a sharp peak in the 2-D frequency domain. For the i th scale factor, the number of FFT point is N u m 1 / P 1 , i , where N u m 1 denotes the number of the FFT point for the reference scale factor. Thus, m S F T t , i F F T τ , i is represented as
m S F T t , i F F T τ , i [ R s ( t , τ , d i ) ] = F F T τ , i ( m S F T t , i [ R s ( t , τ , d i ) ] )
Through analyses of the implementation, m S F T t , i [ ] not only inherits advantages of the SFT with regard to the transformation of the time-domain, but also improves the performance in the frequency-domain. Moreover, the computation cost does not increase too much. The flowchart of m S F T t , i F F T τ , i is shown in Figure 3.
The main steps of estimation parameters by 2D-PMLVD are illustrated as follows.
Input: QFM signals, D . The factors in D are obtained by the selection criterion introduced in Section 2.3.1.
Step1): Substitute QFM signals and D into Equation (3), and R s ( t , τ , D ) is obtained.
Step2):
for i = 1 : L do
Apply the m S F T t , i operator to R s ( t , τ , d i ) along the t axis by Equation (34).
Apply the F F T τ , i [ ] operator to Equation (34) along the τ axis, and we can obtain the data, represented as β i .
end for
Step3): Multiply all β i , we can obtain different peaks, which correspond to different QFM components.
Step4): Detect the peak value to estimate the signal parameter chirp rate and quadratic chirp rate.
Output: the QFM signal parameter chirp rate and quadratic chirp rate.

4. Analyses of Anti-Noise Performance and Computational Cost

In this section, the anti-noise performance and computational cost, which play important roles in parameters estimation [32,33,34], will be analyzed for 2D-PMLVD. HAF-ICPF [30], MLVD [33], and the algorithm presented in [34] are chosen as comparisons.

4.1. Anti-Noise Performance Analysis

In this simulation, we consider a noisy QFM signal. According to analysis methods employed in other estimation algorithms [19,21,32,33,34,36,39], we utilize input-output SNR [21,32] and mean square error (MSE) [37] to evaluate the anti-noise performance of the proposed algorithm with Example 3.
Example 3.
We consider one QFM signal denoted by P 5 . The parameters of P 5 are the same as P 1 . The sampling frequency F s is set to 256 Hz. The effective signal length N is equal to 256. Both a and h are set to one. The scale factors are equal to 56 and 64. The input SNRs tested in Figure 4 are S N R i n = [ 8 : 1 : 0 ] , and 100 trials are performed for each S N R i n value.
The input-output SNR comparison of the four algorithms is shown in Figure 4a, and the solid line denotes the result of the matched filter for P 5 . HAF-ICPF contains HAF and ICPF. HAF is used to reduce the order of the QFM signal, and then ICPF is used to accumulate the energy and estimate the parameter. We know that the threshold SNR of HAF is low [2]. However, according to [40], ICPF will improve the threshold SNR. Therefore, the threshold SNR of HAF-ICPF is −2dB. The evidence for this result can be obtained as described in [33,37]. The algorithm presented in [34], which employs a fourth-order autocorrelation function, is based on coherent accumulation [34]. The algorithm presented in [34] defines a novel parametric autocorrelation function to complete the order reduction and energy accumulation, which benefits anti-noise performance, although redundancy information is not utilized. Thus, in Figure 4a, the threshold SNR of the algorithm presented in [34] is −3 dB. In contrast to the algorithm presented in [34], MLVD utilizes redundancy information, and therefore the anti-noise performance is based on a better anti-noise algorithm. However, MLVD only utilizes one scale factor, which is not favorable for anti-noise performance. Thus, its threshold SNR is −4 dB, which conforms to the result in [33]. In Figure 4a, the threshold SNR of 2D-PMLVD is −5 dB, which is better than the other three algorithms. The advantages include: (1) the energy accumulation of 2D-PMLVD is coherent; (2) in the defined multi-scale PSSF, the redundancy information and multi-scale factors not only reduce the number of self-correlations, but also benefit anti-noise performance [33,37]; and (3) by the m S F T t n , i F F T τ , i and production, the cross terms are effectively suppressed. In addition, it is worthwhile noting that the error propagations exist in HAF-ICPF as a result of the sequential estimation of chirp rate and the quadratic chirp rate. However, this problem does not exist in the other three algorithms.
For 2D-PMLVD, the observed MSEs for the chirp rate and the quadratic chirp rate estimation are plotted in Figure 4b,c as functions of SNR. The corresponding Cramer-Rao bounds (CRBs), which can be obtained as described in [41,42], are also shown by the solid line in the Figure 4. Obviously, in Figure 4b,c, both the MSEs of the chirp rate and quadratic chirp rate estimations are close to CRB when S N R 5 d B . Therefore, the threshold SNR of the 2D-PMVLD is −5 dB for MSEs, and these results conform to those in Figure 4a.

4.2. Computational Cost

The anti-noise performance of PMVLD has been analyzed in the Section 4.1. In this section, the computational cost will be analyzed for PMVLD.
For HAF-ICPF, the main implementation procedures include HAF ( O ( N ) ) [2] and ICPF ( O ( N 3 ) ) [40], where N denotes signal length. Thus, the HAF-ICPF computational cost is ( O ( N 3 ) ) [30,33,34]. For the algorithm presented in [34], its main implementation procedures include a parametric autocorrelation function ( O ( N 2 ) ) , the GSCFT ( O ( N 2 log 2 N ) ) [34] and an FFT operation ( O ( N 2 log 2 N ) ) . Thus, its computational cost is ( O ( N 2 log 2 N ) ) [34]. For MLVD, the main implementation procedures include the parametric symmetric self-correlation function ( O ( N 2 ) ) , the keystone transform ( O ( N 2 log 2 N ) ) and the FFT operation ( O ( N 2 log 2 N ) ) . Thus, the MLVD computational cost is also ( O ( N 2 log 2 N ) ) [33]. For 2D-PMLVD, the main implementation procedures include the defined parametric symmetric self-correlation function ( O ( N 2 ) ) , the m S F T t , i F F T τ , i ( O ( N 2 log 2 N ) ) , and ( L 1 ) products of frequency domain. However, L is usually a small number (e.g., 2), meaning that the additional cost is not excessive. Thus, the computational cost of 2D-PMLVD is ( O ( N 2 log 2 N ) ) . Table 1 lists the computational costs of these four estimation algorithms.

5. Conclusions

In this paper, we propose 2D-PMLVD for QFM signals by employing multi-scale PSSF and m S F T t , i F F T τ , i . Due to the characteristics of multi-scale PSSF, 2D-PMLVD can utilize m S F T t , i F F T τ , i to remove coupling and accumulate energy to the same position under different scale factors. Based on the products of different 2-D frequency domains, cross terms are suppressed and auto-terms are enhanced. The main contributions of the proposed method include the following: (1) it eliminates the procedure of parameter searching; (2) it improves identifiability for multi-QFM signals; (3) it improves anti-noise performance; and (4) it achieves a balance between low computational cost and high performance with low SNR.

Acknowledgments

This work was supported by the China NSF Grant (61671168), the Fundamental Research Funds for the Central Universities (HEUCF160807) and the Natural Science Foundation of Heilongjiang Province of China (QC2016085). The authors wish to thank the anonymous reviewers for their valuable comments on improving this paper.

Author Contributions

The main idea was proposed by Fulong Jing, Shuhong Jiao and Changbo Hou. Weijian Si and Yu Wang performed the experiments and analyzed the simulation results. Fulong Jing wrote the manuscript and all the authors participated in amending the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2017, PP. [Google Scholar] [CrossRef]
  2. Barbarossa, S.; Petrone, V. Analysis of polynomial-phase signals by the integrated generalized ambiguity function. IEEE Trans. Signal Process. 1997, 45, 316–327. [Google Scholar] [CrossRef]
  3. Gu, Y.; Leshem, A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation. IEEE Trans. Signal Process. 2012, 60, 3881–3885. [Google Scholar]
  4. Song, Y.E.; Wang, C.G.; Shi, P. Algorithm based on the linear canonical transform for qfm signal parameters estimation. IET Signal Process. 2016, 10, 318–324. [Google Scholar] [CrossRef]
  5. Gu, Y.; Goodman, N.A. Information-theoretic compressive sensing kernel optimization and bayesian Cramér-Rao bound for time delay estimation. IEEE Trans. Signal Process. 2017, PP. [Google Scholar] [CrossRef]
  6. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2017, 17, 755–765. [Google Scholar] [CrossRef]
  7. Wang, P.; Li, H.; Djurovic, I.; Himed, B. Performance of instantaneous frequency rate estimation using high-order phase function. IEEE Trans. Signal Process. 2010, 58, 2415–2421. [Google Scholar] [CrossRef]
  8. Mo, S.; Wang, Y.; Liu, C. An estimation algorithm for phase errors in synthetic aperture radar imagery. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1818–1822. [Google Scholar]
  9. Wang, B.; Zhang, Y.D.; Wang, W. Robust doa estimation in the presence of mis-calibrated sensors. IEEE Signal Process. Lett. 2017, PP. [Google Scholar] [CrossRef]
  10. Guo, M.; Tao, C.; Wang, B. An improved doa estimation approach using coarray interpolation and matrix denoising. Sensors 2017, 17, 1140. [Google Scholar] [CrossRef] [PubMed]
  11. Zhou, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive Sensing Based Coprime Array Direction-of-Arrival Estimation. IET Commun. 2017, in press. [Google Scholar] [CrossRef]
  12. Berizzi, F.; Mese, E.D.; Diani, M.; Martorella, M. High-resolution isar imaging of maneuvering targets by means of the range instantaneous doppler technique: Modeling and performance analysis. IEEE Trans. Image Process. 2001, 10, 1880–1890. [Google Scholar] [CrossRef] [PubMed]
  13. Li, Y.; Xing, M.; Su, J.; Quan, Y. A new algorithm of isar imaging for maneuvering targets with low snr. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 543–557. [Google Scholar] [CrossRef]
  14. Zhu, D.; Li, Y.; Zhu, Z. A keystone transform without interpolation for sar ground moving-target imaging. IEEE Geosci. Remote Sens. Lett. 2007, 4, 18–22. [Google Scholar] [CrossRef]
  15. Xing, M.; Wu, R.; Bao, Z. High resolutio isar imaging of high speed moving targets. IEE Proc. Radar Sonar Navig. 2005, 152, 58–67. [Google Scholar] [CrossRef]
  16. Wahl, D.E.; Eichel, P.H.; Ghiglia, D.C.; Jakowatz, C.V., Jr. Phase gradient autofocus—A robust tool for high resolution sar phase correction. IEEE Trans. Aerosp. Electron. Syst. 1994, 30, 827–835. [Google Scholar] [CrossRef]
  17. Xi, L.; Liu, G.; Ni, J. Autofocusing of isar images based on entropy minimization. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 1240–1252. [Google Scholar] [CrossRef]
  18. Almeida, L.B. The fractional fourier transform and time-frequency representations. IEEE Trans. Signal Process. 1994, 42, 3084–3091. [Google Scholar] [CrossRef]
  19. Lv, X.; Xing, M.; Wan, C.; Zhang, S. Isar imaging of maneuvering targets based on the range centroid doppler technique. IEEE Trans. Image Process. 2010, 19, 141–153. [Google Scholar] [PubMed]
  20. Bi, G.; Li, X.; Samson See, C.M. Lfm signal detection using lpp-hough transform. Signal Process. 2011, 91, 1432–1443. [Google Scholar] [CrossRef]
  21. Lv, X.; Bi, G.; Wan, C.; Xing, M. Lv’s distribution: Principle, implementation, properties, and performance. IEEE Trans. Signal Process. 2011, 59, 3576–3591. [Google Scholar] [CrossRef]
  22. Abatzoglou, T.J. Fast maximum likelihood joint estimation of frequency and frequency rate. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP ’86, Tokyo, Japan, 7–11 April 1986; pp. 1409–1412. [Google Scholar]
  23. Zheng, J. Fast parameter estimation algorithm for cubic phase signal based on quantifying effects of doppler frequency shift. Prog. Electromagn. Res. 2013, 142, 57–74. [Google Scholar] [CrossRef]
  24. Wu, L.; Wei, X.; Yang, D.; Wang, H. Isar imaging of targets with complex motion based on discrete chirp fourier transform for cubic chirps. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4201–4212. [Google Scholar] [CrossRef]
  25. Barbarossa, S.; Scaglione, A.; Giannakis, G.B. Product high-order ambiguity function for multicomponent polynomial-phase signal modeling. IEEE Trans. Signal Process. 1998, 46, 691–708. [Google Scholar] [CrossRef]
  26. O’Shea, P. A new technique for instantaneous frequency rate estimation. IEEE Signal Process. Lett. 2002, 9, 251–252. [Google Scholar] [CrossRef]
  27. O’Shea, P. A fast algorithm for estimating the parameters of a quadratic fm signal. IEEE Trans. Signal Process. 2004, 52, 385–393. [Google Scholar] [CrossRef]
  28. Wang, Y.; Jiang, Y. Inverse synthetic aperture radar imaging of maneuvering target based on the product generalized cubic phase function. IEEE Geosci. Remote Sens. Lett. 2011, 8, 958–962. [Google Scholar] [CrossRef]
  29. Wang, Y.; Jiang, Y. Isar imaging of a ship target using product high-order matched-phase transform. IEEE Geosci. Remote Sens. Lett. 2009, 6, 658–661. [Google Scholar] [CrossRef]
  30. Wang, Y. Inverse synthetic aperture radar imaging of manoeuvring target based on range-instantaneous- doppler and range-instantaneous-chirp-rate algorithms. IET Radar Sonar Navig. 2012, 6, 921–928. [Google Scholar] [CrossRef]
  31. 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]
  32. Zheng, J.; Su, T.; Zhang, L.; Zhu, W. 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]
  33. Li, Y.; Su, T.; Zheng, J.; He, X. ISAR imaging of targets with complex motions based on modified lv’s distribution for cubic phase signal. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4775–4784. [Google Scholar] [CrossRef]
  34. Zheng, J.; Liu, H.; Liao, G.; Su, T. Isar imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis. IEEE Sens. J. 2016, 16, 2509–2518. [Google Scholar] [CrossRef]
  35. Zheng, J.; Su, T.; Zhu, W.; Zhang, L. Isar imaging of nonuniformly rotating target based on a fast parameter estimation algorithm of cubic phase signal. IEEE Trans. Geosci. Remote Sens. 2015, 53, 1–14. [Google Scholar]
  36. Lv, X.; Xing, M.; Zhang, S.; Bao, Z. Keystone transformation of the wigner-ville distribution for analysis of multicomponent lfm signals. Signal Process. 2009, 89, 791–806. [Google Scholar] [CrossRef]
  37. 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]
  38. Claasen, T.A.C.M.; Mecklenbräuker, W.F.G. The wigner distribution—A tool for time-frequency signal analysis—Part II: Discrete time signals. Philips J. Res. 1980, 35, 276–300. [Google Scholar]
  39. Luo, S.; Bi, G.; Lv, X.; Hu, F. Performance analysis on lv distribution and its applications. Digit. Signal Process. 2013, 23, 797–807. [Google Scholar] [CrossRef]
  40. Wang, P.; Li, H.; Djurović, I.; Himed, B. Integrated cubic phase function for linear fm signal analysis. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 963–977. [Google Scholar] [CrossRef]
  41. Peleg, S.; Porat, B. The cramer-rao lower bound for signals with constant amplitude and polynomial phase. IEEE Trans. Signal Process. 1991, 39, 749–752. [Google Scholar] [CrossRef]
  42. Peleg, S.; Porat, B. Linear fm signal parameter estimation from discrete-time observations. IEEE Trans. Aerosp. Electron. Syst. 1991, 27, 607–616. [Google Scholar] [CrossRef]
Figure 1. Simulation of Example 1. (a) Stereogram with MLVD; (b) Stereogram with 2D-PMLVD; (c) Contour with MLVD; (d) Contour with 2D-PMLVD.
Figure 1. Simulation of Example 1. (a) Stereogram with MLVD; (b) Stereogram with 2D-PMLVD; (c) Contour with MLVD; (d) Contour with 2D-PMLVD.
Sensors 17 01460 g001
Figure 2. Simulation of Example 2. (a) Stereogram with MLVD; (b) Stereogram with 2D-PMLVD; (c) Contour with MLVD; (d) Contour with 2D-PMLVD.
Figure 2. Simulation of Example 2. (a) Stereogram with MLVD; (b) Stereogram with 2D-PMLVD; (c) Contour with MLVD; (d) Contour with 2D-PMLVD.
Sensors 17 01460 g002
Figure 3. Flowchart of m S F T t , i F F T τ , i .
Figure 3. Flowchart of m S F T t , i F F T τ , i .
Sensors 17 01460 g003
Figure 4. Anti-noise performance analysis. (a) Input-output SNR comparison; (b) MSE of the chirp rate estimation; (c) MSE of the quadratic chirp rate estimation.
Figure 4. Anti-noise performance analysis. (a) Input-output SNR comparison; (b) MSE of the chirp rate estimation; (c) MSE of the quadratic chirp rate estimation.
Sensors 17 01460 g004
Table 1. Computational Cost.
Table 1. Computational Cost.
Estimation AlgorithmComputational Cost
HAF-ICPF ( O ( N 3 ) )
The algorithm presented in [34] ( O ( N 2 log 2 N ) )
MLVD ( O ( N 2 log 2 N ) )
2D-PMLVD ( O ( N 2 log 2 N ) )

Share and Cite

MDPI and ACS Style

Jing, F.; Jiao, S.; Hou, C.; Si, W.; Wang, Y. Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals. Sensors 2017, 17, 1460. https://doi.org/10.3390/s17061460

AMA Style

Jing F, Jiao S, Hou C, Si W, Wang Y. Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals. Sensors. 2017; 17(6):1460. https://doi.org/10.3390/s17061460

Chicago/Turabian Style

Jing, Fulong, Shuhong Jiao, Changbo Hou, Weijian Si, and Yu Wang. 2017. "Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals" Sensors 17, no. 6: 1460. https://doi.org/10.3390/s17061460

APA Style

Jing, F., Jiao, S., Hou, C., Si, W., & Wang, Y. (2017). Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals. Sensors, 17(6), 1460. https://doi.org/10.3390/s17061460

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