## 1.

## Introduction

The inverse synthetic aperture radar (ISAR) imaging technique for moving targets has attracted the attention of many radar researchers in the past three decades due to its significance in both civil and military applications.^{1}2.3.4.5.^{–}^{6} Generally, in order to obtain a well-focused ISAR image, the first procedure is to implement motion compensation, which includes range alignment (the translational and rotational migrations) and phase adjustment (the Doppler phase caused by the translation).^{6} Then the conventional range-Doppler (RD) algorithm can be used to generate a focused ISAR image, which indicates that, during the coherent integration interval, the Doppler frequency shift must remain constant. However, instances of targets that are usually uncooperative and engaged in complicated maneuvers often occurs.^{7}8.9.10.11.12.13.^{–}^{14} In this case, the RD algorithm does not work well for the time-varying characteristic of Doppler frequency. To deal with this problem, many articles have developed different solution for different degree of maneuverability, and they can fall into two categories: linear frequency modulation (LFM) signal model^{6}7.8.^{–}^{9} for targets with smooth motion and cubic phase signal (CPS) model^{2}^{,}^{4}^{,}^{5}^{,}^{10}11.12.13.^{–}^{14} for targets with complex motion. For the slow-maneuvering target, the assumption that the azimuth echo is characterized as LFM signal is appropriate and practical, and has been demonstrated by Ref. 6 with experimental ship data. Many representative algorithms based on LFM signal model have been proposed, such as Radon Wigner transform,^{3} the modified Wigner-Ville distribution (WVD),^{6} the stretch keystone-Wigner transform,^{7} the Lv’s distribution,^{8} and the fractional Fourier transform.^{9} For targets with complex motion, in which rotational motion may cause higher-order phase terms, the traditional LFM signal model is not appropriate. For a more reasonable and accurate approximation, many articles have characterized the azimuth echo in a range bin as multicomponent CPSs, which is identified to be much more appropriate than the LFM model.^{10}11.12.13.^{–}^{14}

With regard to the parameter estimation of CPS, numerous algorithms, including the cubic phase function (CPF),^{15}^{,}^{16} the higher-order ambiguity function,^{17} the product generalized CPF (PGCPF),^{18} the product high-order matched-phase transform,^{19} and the modified version of CPF (MCPF),^{14} have been addressed, but all of them involve multilinear transformation and discrete Fourier transform for the nonuniformly spaced signal sample, which result in extensive cross-terms under multi-CPSs and undesirable computational complexity.^{10}11.12.13.^{–}^{14}^{,}^{18}^{,}^{19}20.21.22.^{–}^{23} In Ref. 20, the parameter estimation algorithm based on the local polynomial Wigner distribution (LPWD) is proposed and has been successfully applied to ISAR imaging.^{21}^{,}^{22} However, after compensating the third-order term with the estimated parameter, the LPWD algorithm estimates the second-order coefficient using CPF method, which does lead to the heavy computational burden due to the Fourier transform with respect to the nonuniformly spaced data. In general, the performance and computational cost are conflicting requirements, e.g., the modified discrete chirp Fourier transform for CPS can obtain robust performance against noise and handle multi-CPSs, but it requires a multidimensional search.

Inspired by the proposed LPWD in Refs. 21 and 22 and considering the tradeoff between the computational cost and the performance of parameter estimation, an effective parameter estimation algorithm for multi-CPSs is introduced based on the local polynomial ambiguity function (LPAF) in this paper, which does not suffer from the considerable cross-term interference and heavy computational burden due to its moderate order of nonlinearity, coherent integration via fast Fourier transform (FFT). Then, the high-quality ISAR image can be achieved by the corresponding ISAR imaging algorithm based on LPAF.

This paper is organized as follows. ISAR imaging model with complex motion is given in Sec. 2. In Sec. 3, the parameter estimation algorithm for CPS is proposed based on LPAF. In Sec. 4, the performance of LPAF is discussed in detail, which mainly includes the computational cost and the robustness in the presence of noise. The LPAF-based ISAR imaging algorithm is presented in Sec. 5. Finally, the application of the LPAF-based ISAR imaging algorithm and the conclusion are given in Secs. 6 and 7, respectively.

## 2.

## Inverse Synthetic Aperture Radar Imaging Model with Complex Motion

ISAR imaging model on the imaging project plane is shown in Fig. 1, where the radar is supposed to locate at the same plane with the target, and to orient along the $X$ axis. ${P}_{i}({x}_{i},{y}_{i})$ is the location on the $i$’th point scatterer in the Cartesian coordinates with its origin at the geometric center $O$ of the target, and the $i$’th point scatterer can also be described as ${P}_{i}({r}_{i},{\theta}_{i})$, where ${x}_{i}={r}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\theta}_{i}$, ${y}_{i}={r}_{i}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\theta}_{i}$. ${R}_{t}$ and $\theta ({t}_{m})$ are the translational distance and rotational displacement, respectively.

The instant distance $R({t}_{m})$ between the radar and the $i$’th point scatterer at the instant slow time ${t}_{m}$ can be approximated as

## (1)

$$R({t}_{m})={\{{R}_{t}^{2}+{r}_{i}^{2}+2{R}_{t}{r}_{i}\text{\hspace{0.17em}}\mathrm{cos}[\theta ({t}_{m})+{\theta}_{i}]\}}^{1/2}\phantom{\rule{0ex}{0ex}}\approx {R}_{t}+\frac{{r}_{i}^{2}+2{R}_{t}{r}_{i}\text{\hspace{0.17em}}\mathrm{cos}[\theta ({t}_{m})+{\theta}_{i}]}{2{R}_{t}}\phantom{\rule{0ex}{0ex}}\approx {R}_{t}+{r}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\theta}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta ({t}_{m})-{r}_{i}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\theta}_{i}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta ({t}_{m})\phantom{\rule{0ex}{0ex}}={R}_{t}+{x}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta ({t}_{m})-{y}_{i}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta ({t}_{m}).$$For targets with complex motion, e.g., the highly maneuvering airplane and ship fluctuation with oceanic waves, ${R}_{t}$ and $\theta ({t}_{m})$ are time-varying and can be expressed in Taylor series as follows:

## (2)

$${R}_{t}={R}_{0}+{v}_{t}{t}_{m}+\frac{1}{2}{a}_{t}{t}_{m}^{2}+\frac{1}{6}{\kappa}_{t}{t}_{m}^{3}\cdots ,$$## (3)

$$\theta ({t}_{m})=\alpha {t}_{m}+\frac{1}{2}\beta {t}_{m}^{2}+\frac{1}{6}\gamma {t}_{m}^{3}\cdots ,$$## (4)

$$R({t}_{m})\approx {R}_{0}+{x}_{i}+({v}_{t}-{y}_{i}\alpha ){t}_{m}+\frac{1}{2}({a}_{t}-{y}_{i}\beta ){t}_{m}^{2}+\frac{1}{6}({\kappa}_{t}-{y}_{i}\gamma ){t}_{m}^{3}.$$Assume the transmitted LFM signal takes the following form:

## (5)

$${s}_{t}(\widehat{t})=\mathrm{rect}\left(\frac{\widehat{t}}{T}\right)\mathrm{exp}\left[j2\pi \right({f}_{c}\widehat{t}+\frac{1}{2}\mu {\widehat{t}}^{2}\left)\right],$$## (6)

$${s}_{r}(\widehat{t},{t}_{m})={\sigma}_{i}\text{\hspace{0.17em}}\mathrm{rect}\left[\frac{\widehat{t}-2R({t}_{m})/c}{T}\right]\mathrm{exp}\left(j2\pi \right\{-{f}_{c}\frac{2R({t}_{m})}{c}+\frac{1}{2}\mu {[\widehat{t}-\frac{2R({t}_{m})}{c}]}^{2}\left\}\right),$$After the pulse compression by using the matched filter (MF) $H(\widehat{t})=\mathrm{rect}(\widehat{t}/T)\mathrm{exp}(j\pi \mu {\widehat{t}}^{2})$, we have

## (7)

$${s}_{\mathrm{pc}}(\widehat{t},{t}_{m})={\sigma}_{i}B\text{\hspace{0.17em}}\mathrm{sinc}\left\{B\right[\widehat{t}-\frac{2R({t}_{m})}{c}\left]\right\}\mathrm{exp}[-j\frac{4\pi}{\lambda}R({t}_{m})],$$In Eq. (7), the complex motion will result in range migration and Doppler frequency shift. The translational range migration has the same feature for all scatterers and has no contribution for ISAR imaging. Hence, in this paper, the standard range alignment method^{24} and the phase gradient autofocus method^{25} are utilized to eliminate the effect caused by translation, and keystone transform to correct the rotational range migration. Here, we only concentrate on the Doppler frequency shift caused by rotation.

Consequently, after the motion compensation (the range migration compensation and the translational induced phase error correction), the azimuth echo of the $i$’th scatter in the range-azimuth domain can be rewritten as

## (8)

$${s}_{\mathrm{pc}}(\widehat{t},{t}_{m})={\sigma}_{i}B\text{\hspace{0.17em}}\mathrm{sinc}\left\{B\right[\widehat{t}-\frac{2({R}_{o}+{x}_{i})}{c}\left]\right\}\mathrm{exp}\left\{j2\pi \right[\left(\frac{2{y}_{i}\alpha}{\lambda}\right){t}_{m}+\frac{1}{2}\left(\frac{2{y}_{i}\beta}{\lambda}\right){t}_{m}^{2}+\frac{1}{6}\left(\frac{2{y}_{i}\gamma}{\lambda}\right){t}_{m}^{3}\left]\right\}.$$In realistic ISAR imaging applications, there will be multiple scatterers in a range bin simultaneously, which is more accurate and has been confirmed in many literatures. Thus, the azimuth echo signals in a range bin can be rewritten as

## (9)

$$s({t}_{m})=\sum _{i=1}^{Q}{A}_{i}\text{\hspace{0.17em}}\mathrm{exp}\left[j2\pi \right({\varphi}_{i,1}{t}_{m}+\frac{1}{2}{\varphi}_{i,2}{t}_{m}^{2}+\frac{1}{6}{\varphi}_{i,3}{t}_{m}^{3}\left)\right],$$## 3.

## Parameters Estimation for Cubic Phase Signal

## 3.1.

### Introduction of Local Polynomial Wigner Distribution

LPWD,^{20} a special kind of time-frequency distribution proposed by Stankovic, is widely used for the parameters estimation of multi-CPSs. Then the ISAR imaging algorithm^{21}^{,}^{22} is also proposed based on LPWD. LPWD is defined in Ref. 20 as follows:

## (10)

$$W(\overrightarrow{\omega},t)={\int}_{-\infty}^{+\infty}s(t+\frac{\tau}{2}){s}^{*}(t-\frac{\tau}{2})\mathrm{exp}[-j2\pi \theta (\overrightarrow{\omega},\tau )]\mathrm{d}\tau ,$$## (11)

$$\theta (\overrightarrow{\omega},\tau )={\omega}_{1}\tau +{\omega}_{2}\frac{{\tau}^{2}}{2!}+{\omega}_{3}\frac{{\tau}^{3}}{3!}+\dots +{\omega}_{m}\frac{{\tau}^{m}}{m!},$$Thereby, the instantaneous frequency estimation of polynomial phase signal (PPS) can be achieved with a significantly reduced distribution dimension.

## 3.2.

### Definition of Local Polynomial Ambiguity Function

Considering an analytic signal $s(t)$, the ambiguity function (AF) is defined as

## (12)

$$A(\tau ,\omega )={\int}_{-\infty}^{+\infty}s(t+\frac{\tau}{2}){s}^{*}(t+\frac{\tau}{2})\mathrm{exp}(-j2\pi \omega t)\mathrm{d}t.$$Borrowing the idea from the definitions of LPWD and AF, the LPAF is defined as follows:

## (13)

$$\mathrm{LPAF}(\overrightarrow{\omega},\tau )={\int}_{-\infty}^{+\infty}s(t+\frac{\tau}{2}){s}^{*}(t-\frac{\tau}{2})\mathrm{exp}[-j2\pi \theta (\overrightarrow{\omega},t)]\mathrm{d}t.$$Similar to the definition of LPWD, we have

## (14)

$$\theta (\overrightarrow{\omega},t)={\omega}_{1}t+{\omega}_{2}\frac{{t}^{2}}{2!}+{\omega}_{3}\frac{{t}^{3}}{3!}+\dots +{\omega}_{m}\frac{{t}^{m}}{m!},$$Consider mono-CPS with the following form:

## (15)

$$s(t)=A\text{\hspace{0.17em}}\mathrm{exp}\left[j2\pi \right({\varphi}_{1}t+\frac{1}{2}{\varphi}_{2}{t}^{2}+\frac{1}{6}{\varphi}_{3}{t}^{3}\left)\right],$$## (16)

$$\mathrm{LPAF}({\omega}_{1},{\omega}_{2},\tau )={\int}_{-\infty}^{+\infty}s(t+\frac{\tau}{2}){s}^{*}(t-\frac{\tau}{2})\mathrm{exp}[-j2\pi ({\omega}_{1}t+{\omega}_{2}\frac{{t}^{2}}{2}\left)\right]\mathrm{d}t.$$Substituting Eq. (15) in Eq. (16), we obtain

## (17)

$$\mathrm{LPAF}({\omega}_{1},{\omega}_{2},\tau )={A}^{2}\text{\hspace{0.17em}}\mathrm{exp}[j\mathrm{\Phi}(\tau )]{\int}_{-\infty}^{+\infty}\mathrm{exp}[-j2\pi ({\omega}_{1}-{\varphi}_{2}\tau )t]\mathrm{exp}[-j2\pi (\frac{1}{2}{\omega}_{2}-\frac{1}{2}{\varphi}_{3}\tau \left){t}^{2}\right]\mathrm{d}t,$$## (18)

$$\mathrm{LPAF}({\omega}_{1},{\omega}_{2})={D}^{2}{\int}_{-\infty}^{+\infty}\mathrm{exp}[-j2\pi ({\omega}_{1}-{\varphi}_{2}{\tau}_{0})t]\mathrm{exp}[-j2\pi (\frac{1}{2}{\omega}_{2}-\frac{1}{2}{\varphi}_{3}{\tau}_{0}\left){t}^{2}\right]\mathrm{d}t,$$## (19)

$$\{\begin{array}{l}{\omega}_{1}={\varphi}_{2}{\tau}_{0}\\ {\omega}_{2}={\varphi}_{3}{\tau}_{0}\end{array}.$$From Eq. (19), it can be seen that the different lag-time instance ${\tau}_{0}$ indicates the different zoom factor, which can enlarge the estimation ranges of CR and QCR to avoid the spectrum aliasing and estimation error^{26} in some extreme ISAR imaging case, e.g., higher maneuverability. On the other hand, the symmetric instantaneous self-correlation operation will shorten the effective signal in Eq. (16) when ${\tau}_{0}\ne 0$. Based on the recommendation from Ref. 27, the optimal discrete lag is ${N}_{{\tau}_{0}/2}^{\mathrm{opt}}=(N+1)/2\text{\hspace{0.17em}\hspace{0.17em}}H$, where $H$ is the PPS order. For convenience to analyze, we choose ${\tau}_{0}=1$ so that the values of CR and QCR can be read directly from CRQCRD.

Thereby, we can estimate the parameters ${\varphi}_{2}$ and ${\varphi}_{3}$ simultaneously by searching the position of the maximum of $|\mathrm{LPAF}({\omega}_{1},{\omega}_{2})|$.

## (20)

$$({\widehat{\varphi}}_{2},{\widehat{\varphi}}_{3})=\underset{{\omega}_{1},{\omega}_{2}}{\mathrm{arg}\mathrm{max}}|\mathrm{LPAF}({\omega}_{1},{\omega}_{2})|.$$With the estimated parameters $({\widehat{\varphi}}_{2},{\widehat{\varphi}}_{3})$ and dechirp technique, the other parameters can be estimated via FFT after compensating the high-order phase terms.

According to the aforementioned analyses on LPWD and LPAF for CPS, we know that LPAF is similar to LPWD except for the definition of $\theta (\overrightarrow{\omega},\mathrm{\Delta})$ in Eqs. (11) and (14), where $\mathrm{\Delta}=\tau $ for LPWD and $\mathrm{\Delta}=t$ for LPAF. However, instead of the special time instance as $t=0$,^{21} LPAF can choose the different lag-time instance ${\tau}_{0}$ to introduce a zoom factor adjusting the estimation ranges of CR and QCR in Eq. (19). It must be emphasized that such a difference circumvents the problems of LPWD including spectrum aliasing, error propagation, and computational burden, which will be discussed in Sec. 4.

## 3.3.

### Local Polynomial Ambiguity Function for Multicomponent Cubic Phase Signals

In this section, the performance of LPAF for multi-CPSs is shown. As we know, correlation operation will result in the cross-term effects among multi-CPSs, which determines whether the proposed algorithm is suitable for multi-CPSs or not. For simplicity, we choose special case as ${\tau}_{0}=1$ and consider two-component CPSs in the following structure:

## (21)

$$s(t)={s}_{p}(t)+{s}_{q}(t)\phantom{\rule{0ex}{0ex}}={A}_{p}\text{\hspace{0.17em}}\mathrm{exp}[j2\pi ({\varphi}_{p,1}t+{\varphi}_{p,2}{t}^{2}+{\varphi}_{p,3}{t}^{3})]+{A}_{q}\text{\hspace{0.17em}}\mathrm{exp}[j2\pi ({\varphi}_{q,1}t+{\varphi}_{q,2}{t}^{2}+{\varphi}_{q,3}{t}^{3})].$$Substituting Eq. (21) in Eq. (16), we obtain the LPAF for $s(t)$ in Eq. (21) as follows:

## (22)

$$\mathrm{LPAF}({\omega}_{1},{\omega}_{2})={D}_{p}^{2}{\int}_{-\infty}^{+\infty}\mathrm{exp}[-j2\pi ({\omega}_{1}-{\varphi}_{p,2})t]\times \mathrm{exp}[-j2\pi (\frac{1}{2}{\omega}_{2}-\frac{1}{2}{\varphi}_{p,3}\left){t}^{2}\right]\mathrm{d}t\phantom{\rule{0ex}{0ex}}+{D}_{q}^{2}{\int}_{-\infty}^{+\infty}\mathrm{exp}[-j2\pi ({\omega}_{1}-{\varphi}_{q,2})t]\times \mathrm{exp}[-j2\pi (\frac{1}{2}{\omega}_{2}-\frac{1}{2}{\varphi}_{q,3}\left){t}^{2}\right]\mathrm{d}t\phantom{\rule{0ex}{0ex}}+{\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})+{\mathrm{LPAF}}_{qp,\text{cross}}({\omega}_{1},{\omega}_{2}).$$Note that we have two cross-terms ${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})$ and ${\mathrm{LPAF}}_{qp,\text{cross}}({\omega}_{1},{\omega}_{2})$, and we only take ${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})$ as an example to analyze. Obviously, ${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})$ is with the form

## (23)

$${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})={D}_{pq,\text{cross}}{\int}_{-\infty}^{+\infty}\mathrm{exp}\{-j2\pi [{\omega}_{1}-({\varphi}_{p,1}-{\varphi}_{q,1})-\frac{1}{2}({\varphi}_{p,2}+{\varphi}_{q,2})-\frac{1}{8}({\varphi}_{p,3}-{\varphi}_{q,3})\left]t\right\}\phantom{\rule{0ex}{0ex}}\times \mathrm{exp}\{-j2\pi [{\omega}_{2}-\frac{1}{2}({\varphi}_{p,2}-{\varphi}_{q,2})-\frac{1}{4}({\varphi}_{p,3}+{\varphi}_{q,3})\left]{t}^{2}\right\}\phantom{\rule{0ex}{0ex}}\times \mathrm{exp}\left\{j2\pi \right[\frac{1}{6}({\varphi}_{p,3}-{\varphi}_{q,3})\left]{t}^{3}\right\}\mathrm{d}t,$$As is obvious from Eq. (23), when ${\varphi}_{p,3}\ne {\varphi}_{q,3}$, the auto-terms can be detected correctly by peak detection, and the phase coefficients can be estimated easily. Here, assume that the received signal contains the same QCR, i.e., ${\varphi}_{p,3}={\varphi}_{q,3}$. Then, ${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})$ can be further simplified as

## (24)

$${\mathrm{LPAF}}_{pq,\text{cross}}({\omega}_{1},{\omega}_{2})={D}_{pq,\text{cross}}{\int}_{-\infty}^{+\infty}\mathrm{exp}\{-j2\pi [{\omega}_{1}-({\varphi}_{p,1}-{\varphi}_{q,1})-\frac{1}{2}({\varphi}_{p,2}+{\varphi}_{q,2})\left]t\right\}\phantom{\rule{0ex}{0ex}}\times \mathrm{exp}\{-j2\pi [{\omega}_{2}-\frac{1}{2}({\varphi}_{p,2}-{\varphi}_{q,2})-\frac{1}{2}({\varphi}_{p,3})\left]{t}^{2}\right\}\mathrm{d}t.$$From Eq. (24), we can find that the spurious peak will be located at

## (25)

$$\{\begin{array}{l}{\omega}_{1}=({\varphi}_{p,1}-{\varphi}_{q,1})+\frac{1}{2}({\varphi}_{p,2}+{\varphi}_{q,2})\\ {\omega}_{2}=\frac{1}{2}({\varphi}_{p,2}-{\varphi}_{q,2})+\frac{1}{2}({\varphi}_{p,3})\end{array}.$$It is worthwhile to mention that only when ${\varphi}_{p,3}={\varphi}_{q,3}$, the LPAF for $s(t)$ causes spurious peak, which affects the detection and estimation for multi-CPSs. Fortunately, the phase coefficients ${\varphi}_{j}(j=\mathrm{0,1},2\dots )$ are not only related to target rotational parameter (e.g., CF corresponding to the initial rotational velocity, CR corresponding to the rotational acceleration, and QCR corresponding to the rotational acceleration rate), but also proximately proportional to the cross-range coordinates. Thus, for different scatterers with different coordinates in a range bin, we have ${\varphi}_{p3}\ne {\varphi}_{q3}(p\ne q)$.^{14} Accordingly, the cross-terms of the LPAF induced by multi-CPSs can be avoided effectively.

*Example 1.* Consider a mono-component CPS and two-component CPSs, respectively. The sampling frequency and the effective length of signal are 256 Hz and 256, respectively. The signal parameters are shown in Table 1.

## Table 1

Parameters for CPS.

Ai | ϕi,1 | ϕi,2 | ϕi,3 | |||
---|---|---|---|---|---|---|

i=1 | 1 | 40 Hz | 30 Hz/s | −60 Hz/s2 | ||

i=2 | 1.05 | −30 Hz | −80 Hz/s | 40 Hz/s2 |

The WVD of CPS-A1 is shown in Fig. 2(a), in which the curves appear due to QCR of CPS. With the proposed parameter estimation algorithm, a sole peak appears in the CRQCRD in Figs. 2(b) and 2(c), where CR and QCR can be obtained by peak detection. Moreover, we obtain the CRQCRD of the two-component CPSs in Fig. 2(d). For easy visualization, the stereogram, the right view ($XZ$ plane), and the left view ($XZ$ plane) of CRQCRD are depicted in Figs. 2(d)–2(f), where two different peaks appear in CRQCRD. Therefore, the corresponding parameters of each component can be easily estimated via searching the peaks and CLEAN technique,^{7}^{,}^{8} sequentially.

## 4.

## Performance of Local Polynomial Ambiguity Function

In this section, we will discuss the performance of LPAF, mainly including the computational cost and the robustness in the presence of noise, which determine whether the proposed algorithm can be utilized for the parameter estimation of multi-CPSs in realistic ISAR imaging applications or not. In order to verify the advantages of the proposed estimation algorithm, other two representative algorithms, MCPF and LPWD, are chosen as reference to compare with the proposal.

## 4.1.

### Analysis of the Computational Cost

For MCPF, the estimation procedures of the parameter ${\varphi}_{3}$ require $O(3{N}_{{t}_{m}}^{2})$ complex multiply operations (CMs) for definition of multilinear function and $O({N}_{{t}_{m}}^{3})$ operations for discrete Fourier transform (DFT) along the nonuniformly spaced lag-time axis. For estimation of the parameter ${\varphi}_{2}$, the implementation procedures require $O({N}_{{t}_{m}}^{2})$ operations for CPF and $O({N}_{{t}_{m}}^{3})$ operations for DFT due to the two-dimensional (2-D) kernel. Thus, the computational complexity with MCPF is in the order of $O({N}_{{t}_{m}}^{3})$. According to Refs. 11 and 13, the nonuniform fast Fourier transform (NUFFT)^{11}^{,}^{13} can be utilized to speed up the implementation of the Fourier transform along the nonuniformly spaced lag-time axis. Unfortunately, the interpolation necessary for NUFFT is not preferable in some realistic applications.

The LPAF algorithm, which can directly obtain CRQCRD, estimates CR and QCR with 2-D maximizations and can reduce error propagation effect to some extent. According to the aforementioned analyses and discussion on the modified discrete chirp-Fourier transform,^{28}^{,}^{29} which retains capability of efficient FFT implementation, the implementation procedures of parameters estimation based on LPAF include (1) the self-correlation operation in Eq. (16), after which the signal has the form of the linear frequency modulated signal, (2) FFT operation after matching the QCR ${\varphi}_{3}$ in Eq. (18), and (3) the parameters estimation by peak detection technique. In this paper, the searching range of the parameter ${\varphi}_{3}$ is $[-\mathrm{ROUND}(F{t}_{m}/2),\mathrm{ROUND}(F{t}_{m}/2)]$, where $\mathrm{ROUND}(\xb7)$ and $F{t}_{m}$ denote the round-up operator and pulse repetition frequency (PRF), respectively. The resolution of the parameter ${\varphi}_{3}$ is $\mathrm{PRF}/{N}_{{t}_{m}}$, where ${N}_{{t}_{m}}$ is the effective length of ${t}_{m}$ (the same as the number of effective echoes); thus, the searching point of the parameter ${\varphi}_{3}$ is exactly ${N}_{{t}_{m}}$ in ISAR imaging. In conclusion, the computational complexity with this approach is about in the order of $O({N}_{{t}_{m}}^{2}\text{\hspace{0.17em}}{\mathrm{log}}_{2}\text{\hspace{0.17em}}{N}_{{t}_{m}})$.

For LPWD, the implementation procedures for the estimation of parameters ${\varphi}_{1}$ and ${\varphi}_{3}$ are analogous to the LPAF algorithm. However, the parameter ${\varphi}_{2}$ is estimated by the integrated form of cubic phase function, which includes the defined bilinear function [$O({N}_{{t}_{m}}^{2})$] and Fourier transform [$O({N}_{{t}_{m}}^{3})$] along the nonuniformly spaced lag-time axis. Therefore, the computational complexity with the LPWD-based algorithm is also in the order of $O({N}_{{t}_{m}}^{3})$, which is obviously greater than the LPAF algorithm.

It is obvious that, compared with MCPF and LPWD, the proposed algorithm is computationally more efficient than the other two algorithms since the latter involves nonuniformly spaced data computation. Furthermore, as mentioned in Ref. 21, MCPF estimates phase coefficients one by one via dechirp technique, which will reduce the estimation precision due to error propagation. On the contrary, the LPAF algorithm can obtain CR and QCR simultaneously. Table 2 gives the computational costs of the three algorithms.

## Table 2

Computational cost.

Computation cost | LPWD | MCPF | LPAF |
---|---|---|---|

Number of CMs | O(Ntm3) | O(Ntm3) | O(Ntm2 log2 Ntm) |

## 4.2.

### Analysis of Performance in Noise

In order to further illustrate the practicability of the proposed parameter estimation algorithm, the robustness of performance against noise, which is evaluated by input-output signal-to-noise ratio (SNR)^{9} and mean square error (MSE),^{30}^{,}^{31} will be analyzed in this section.

*Example 2.* Consider a mono-CPS as shown in Eq. (15); the sampling frequency and the effective length of signal are 256 Hz and 256, respectively. The signal parameters are set as follows: $A=1$, ${\varphi}_{1}=40\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}$, ${\varphi}_{2}=24\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/\mathrm{s}$, ${\varphi}_{3}=60\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/{\mathrm{s}}^{2}$. Hereon, the CPS is embedded in complex white Gaussian noise, and the value of ${\mathrm{SNR}}_{\mathrm{in}}$ varies from $-5$ to 2 dB in steps of 1 dB, and 100 Monte Carlo simulations are performed for each ${\mathrm{SNR}}_{\mathrm{in}}$ value. Meanwhile, the MSEs for the parameters ${\varphi}_{2}$ and ${\varphi}_{3}$ with the corresponding Cramer-Rao bounds (CRBs)^{30} are also obtained and shown in Figs. 3(b) and 3(c), respectively.

In Fig. 3(a), MF, LPWD, and MCPF are chosen to compare with the proposed LPAF. For MCPF, its implementation is similar to PGCPF,^{18} which defines a multilinear function of fourth order and estimates unknown parameters one by one with two one-dimensional maximizations. Thus, the performance of MCPF closely matches that of MF until the input SNR falls below 0 dB, and the evidence for MCPF can be found in Ref. 10. For LPWD, its performance also closely matches that of MF until the input SNR falls below 0 dB, and we can find the corresponding evidence in Ref. 21. According to the analyses in Refs. 8 and 11, compared to the self-correlation function [$s(\tau /2){s}^{*}(-\tau /2)$, for $t=0$] of LPWD, the self-correlation function [$s(t+{\tau}_{0}/2){s}^{*}(t-{\tau}_{0}/2)$, for $\tau ={\tau}_{0}$] of LPAF can reduce the influence of nonlinearity and benefit the performance in the presence of noise. Thereby, the proposed LPAF can obtain robust performance against noise and its threshold SNR is $-3\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$. In addition, the MSEs for the estimations ${\varphi}_{2}$ and ${\varphi}_{3}$ are plotted in Figs. 3(b) and 3(c) as a function of SNR, respectively. Obviously, the MSEs of CR and QCR are close to the corresponding CRB when ${\mathrm{SNR}}_{\mathrm{in}}\ge -3\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$, which conforms to the result in Fig. 3(a).

## 5.

## Inverse Synthetic Aperture Radar Imaging Algorithm Based on Local Polynomial Ambiguity Function

On the above analyses and confirmation, a new ISAR imaging algorithm for target with complex motion is proposed based on LPAF. Detail implementation procedures are given as follows:

Step 1: Complete the range compression and motion compensation, and obtain data ${s}^{k}({t}_{m})$ of the $k$’th (where $1\le k\le K$, and $K$ is the total number of range bins) range bin.

Step 2: Substitute ${s}^{k}({t}_{m})$ in Eq. (16).

where ${\mathrm{LPAF}}_{\text{cross}}^{k}({\omega}_{1},{\omega}_{2})$ and ${\mathrm{LPAF}}_{\text{noise}}^{k}({\omega}_{1},{\omega}_{2})$ denote the cross-terms and noise, respectively.## (26)

$${\mathrm{LPAF}}^{k}({\omega}_{1},{\omega}_{2})=\sum _{i=1}^{Q}{D}_{i}^{2}{\int}_{-\infty}^{+\infty}\mathrm{exp}[-j2\pi ({\omega}_{1}-{\varphi}_{i,2})t]\mathrm{exp}[-j2\pi (\frac{1}{2}{\omega}_{2}-\frac{1}{2}{\varphi}_{i,3}\left){t}^{2}\right]\mathrm{d}t\phantom{\rule{0ex}{0ex}}+{\mathrm{LPAF}}_{\text{cross}}^{k}({\omega}_{1},{\omega}_{2})+{\mathrm{LPAF}}_{\text{noise}}^{k}({\omega}_{1},{\omega}_{2}),$$Step 3: Complete estimation of the parameters ${\varphi}_{i,2}$ and ${\varphi}_{i,3}$ of the $i$’th CPS by finding the location of peak.

Step 4: Dechirp ${s}^{k}({t}_{m})$ with the estimated parameters ${\varphi}_{i,2}$ and ${\varphi}_{i,3}$ to estimate ${\varphi}_{i,1}$ and ${A}_{i}$.

## (28)

$$({\widehat{A}}_{i}=\frac{G}{N},{\widehat{\varphi}}_{i,1}=f)=\underset{G,f}{\mathrm{arg}\mathrm{max}}\left|\mathrm{FFT}\right\{{s}^{k}({t}_{m})\times \mathrm{exp}[-j2\pi (\frac{1}{2}{\widehat{\varphi}}_{i,2}{t}_{m}^{2}+\frac{1}{6}{\widehat{\varphi}}_{i,3}{t}_{m}^{3}\left)\right]\left\}\right|.$$Step 5: Subtract the estimated $i$’th CPS from the original signal ${s}^{k}({t}_{m})$ in frequency domain.

where $\mathrm{WIN}({\widehat{\varphi}}_{i,1})=\{\begin{array}{ll}0,& f-\mathrm{\Delta}f/2<{\widehat{\varphi}}_{i,1}<f+\mathrm{\Delta}f/2\\ 1,& \text{otherwise}\end{array}$ is a notch filter centred at the CF ${\widehat{\varphi}}_{i,1}$, and $\mathrm{\Delta}f$ is the corresponding bandwidth.## (29)

$${s}^{k}({t}_{m})=\mathrm{IFFT}(\mathrm{WIN}({\widehat{\varphi}}_{i,1})\times \mathrm{FFT}\{{s}^{k}({t}_{m})\mathrm{exp}[-j2\pi (\frac{1}{2}{\widehat{\varphi}}_{i,2}{t}_{m}^{2}+\frac{1}{6}{\widehat{\varphi}}_{i,3}{t}_{m}^{3}\left)\right]\left\}\right)\phantom{\rule{0ex}{0ex}}\times \mathrm{exp}\left[j2\pi \right(\frac{1}{2}{\widehat{\varphi}}_{i,2}{t}_{m}^{2}+\frac{1}{6}{\widehat{\varphi}}_{i,3}{t}_{m}^{3}\left)\right],$$Step 6: Renew $i=i+1$ and repeat steps 2 to 5 until the residual energy of the signal is less than the threshold $T$.

Step 7: Renew $k=k+1$ and repeat steps 1 to 6 until $k=K$.

## 6.

## Inverse Synthetic Aperture Radar Imaging Results Based on Local Polynomial Ambiguity Function

In this section, we will utilize the synthetic ship target and aircraft target with the model shown in Figs. 4(a) and 6(a) to verify the proposed ISAR imaging algorithm successively. Motion parameters and corresponding radar parameters are listed in Tables 3 and 4, respectively.

## Table 3

Simulation parameters of ship.

Radar | Ship target | |||
---|---|---|---|---|

Carrier frequency | 15 GHz | Target distance | 24 km | |

Bandwidth | 200 MHz | Motions | Translation | Rotation |

Pulse repetition frequency (PRF) | 256 Hz | Velocity | 30 m/s | 0.012 rad/s |

Effective echo | 256 | Acceleration | 3 m/s2 | 0.02 rad/s2 |

Sample frequency | 200 MHz | Acceleration rate | 2 m/s3 | 0.04 rad/s3 |

## Table 4

Simulation parameters of aircraft.

Radar | Aircraft target | |||
---|---|---|---|---|

Carrier frequency | 10 GHz | Target distance | 40 km | |

Bandwidth | 150 MHz | Motions | Translation | Rotation |

PRF | 256 Hz | Velocity | 20 m/s | 0.03 rad/s |

Echo number | 512 | Acceleration | 1 m/s2 | 0.04 rad/s2 |

Sample frequency | 300 MHz | Acceleration rate | 2 m/s3 | 0.08 rad/s3 |

## 6.1.

### Ship Target

Similar to Refs. 10 and 13, a fluctuating ship target modeled as a set of ideal scatterers is provided to demonstrate the effectiveness of the proposed algorithm in this paper. After pulse compression and motion compensation, the received signal can be simplified as the 2-D reflectivity function shown in Eq. (8). Figure 4(b) gives the conventional RD imaging results with FFT operation along the azimuth slow time axis under the situation of ${\mathrm{SNR}}_{\mathrm{in}}=-3\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$, which is blurred severely by CR and QCR. In this section, we choose the special lag-time instance ${\tau}_{0}=1$ to achieve the LPAF so that CR and QCR can be read directly.

To demonstrate the accuracy of the parameters estimation by the proposed LPAF, the data of the 153rd range bin in Fig. 4(b) are extracted to analyze, and the corresponding CRs and QCRs of these scatterers are estimated as shown in Table 5. Then, with the estimated CRs and QCRs, we compensate the Doppler frequency shift and obtain Figs. 4(d)–4(f), where the ideal matched results are also shown to make a comparison with each mentioned algorithm. Compared with the results in Figs. 4(c), 4(e), and 4(f), the result obtained by LPAF in Fig. 4(d) is very close to the ideal results, which demonstrate that the accuracy of the parameters estimation by the proposed LPAF is high. Although some spurious peaks appear in Fig. 4(d), the same as in Table 5, the LPAF is also better than other algorithms in the estimation accuracy and suppression on cross-terms.

## Table 5

Results of parameters estimation.

Signal parameters | A1 | A2 | A3 | A4 | A5 | A6 | A7 | A8 | Spurious peak | |
---|---|---|---|---|---|---|---|---|---|---|

CR (Hz/s) | 60 | 40 | 20 | 0 | −20 | −40 | −60 | −80 | ||

Estimated CR (Hz/s) | 60 | 40 | 20 | 0 | −20 | −40 | −61 | −75 | −54 | −42 |

QCR (Hz/s2) | 120 | 80 | 40 | 0 | −40 | −80 | −120 | −160 | ||

Estimated QCR (Hz/s2) | 120 | 80 | 40 | 0 | −38 | −80 | −124 | −84 | −38 | 82 |

In addition, the comparison of the LPAF-based algorithm with the MCPF-based algorithm and the LPWD-based algorithm for ISAR imaging is discussed. Generally, to quantitatively evaluate the performance of the proposed algorithm, the MSEs of 2-D reconstructed coordinates are defined as $\mathrm{MSE}={\Vert {\mathfrak{R}}_{\mathrm{est}}-\mathfrak{R}\Vert}_{2}$, where $\mathfrak{R}$ and ${\mathfrak{R}}_{\mathrm{est}}$ represent the original data and the reconstructed data obtained by the proposed algorithm, that is, the corresponding relation between $\mathfrak{R}$ and ${\mathfrak{R}}_{\mathrm{est}}$ is determinate and known ahead of time. However, the number and the position of the scatterers on target are usually unknown and cannot be obtained in advance. Entropy can be used to measure the smoothness of a distribution function. The smoother a distribution function, the larger is its entropy.^{10}^{,}^{13} Owing to this property, we use the entropy of Eq. (30) as a criterion to measure the quality of the image in Table 6. Better focus results in a sharper image and thus smaller entropy. Meanwhile, simulation times of Fig. 5 (only including effective range bins) are also listed to certify the high efficiency of LPAF-based ISAR imaging algorithm.

## Table 6

Entropies and time of ISAR images in Fig. 5.

Fig. 5(a) | Fig. 5(b) | Fig. 5(c) | |
---|---|---|---|

Entropies | 3.4421 | 4.8379 | 4.2688 |

Simulation time (s) | 1.2318 | 110.0910 | 202.9427 |

From Fig. 5, we see that LPAF-based ISAR imaging algorithm has better suppression on cross-terms compared with the LPWD-based ISAR imaging algorithm and MCPF-based ISAR imaging algorithm. Evidently, most scatterers are relocated correctly in Fig. 5(a). Due to the four order of nonlinearity of MCPF and error propagation effect of LPWD (the estimation precision for CR is influenced by QCR), the spurious scatterers in Figs. 5(b) and 5(c) are more than those in Fig. 5(a). Furthermore, the entropy of Fig. 5(a) is also smaller than those of Figs. 5(b) and 5(c) in Table 6, which further shows that the LPAF-based ISAR imaging algorithm can obtain a well-focused image.

## 6.2.

### Aircraft Target

In this section, the synthetic aircraft target as shown in Fig. 6(a), which consists of 141 scatterers, is chosen to demonstrate the validity of the proposed ISAR imaging algorithm in this paper. We know that the CR and QCR of the azimuth echo in ISAR imaging are determined by the target’s size and the effective rotating vector. When the target’s size is a little larger or its maneuverability is a little higher, we have to adjust the parameter estimation range to meet the practical requirements. Here, we choose the special instance ${\tau}_{0}=1/2$ to enlarge the parameter estimation ranges and avoid spectrum aliasing and estimation error.

Obviously, the image obtained by the traditional RD algorithm is blurred severely in Fig. 6(b). In order to confirm the validity of the proposed parameter estimation algorithm, when all scatterers have been corrected into right slant range bins after motion compensation, we extract the data of the 211st range bin to analyze in detail. Under the situation of ${\mathrm{SNR}}_{\mathrm{in}}=4\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$, the WVD as shown in Fig. 6(c) is computed for the radar echo signal of the 211st range bin. It is not difficult to find that the Doppler frequency is time-varying and nonlinear, which also demonstrates that the CPS model is very appropriate. Here we apply the LPAF-based parameter estimation algorithm to the extracted data. In Figs. 6(d) and 6(e), two peaks appear and can be easily picked out. Then the corresponding CR and QCR of two peaks can be correctly estimated via peak detection technique and CLEAN technique.

With the parameter estimation algorithm, we can estimate $({\varphi}_{p,2},{\varphi}_{p,3})$ and $({\varphi}_{q,2},{\varphi}_{q,3})$ as $(4\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/\mathrm{s},8\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/{\mathrm{s}}^{2})$ and $(-4\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/\mathrm{s},-8\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{Hz}/{\mathrm{s}}^{2})$, respectively. After the Doppler spread compensation with the estimated parameters, we obtain the results in Fig. 6(f) by performing FFT, where the results of the conventional RD algorithm are also provided. Simulation shows that signal energy cannot be focused efficiently for Doppler frequency shift induced by CR and QCR. Consequently, the effective parameter estimation and compensation are quite essential.

Generally speaking, the better the focused quality of image, the smaller is the entropy.^{5}^{,}^{10}^{,}^{11}^{,}^{13}^{,}^{14} As before, the entropy is applied to measure the quality of the image. Compared with the LPWD-based imaging algorithm and the MCPF-based imaging algorithm, this proposed LPAF-based imaging algorithm can obtain robust performance against noise. Therefore, most scatterers have been reconstructed correctly in Fig. 7(a) and spurious scatterers in Fig. 7(a) are fewer than those in Figs. 7(b) and 7(c). Also, the entropy of image in Fig. 7(a) is smaller than those in other figures in Table 7, which indicates that the image quality for the LPAF-based ISAR imaging algorithm is better than those for other algorithms.

## Table 7

Entropies of ISAR images in Fig. 7.

Fig. 7(a) | Fig. 7(b) | Fig. 7(c) | |
---|---|---|---|

Entropies | 5.6388 | 5.9089 | 5.8016 |

In order to further validate the effectiveness of the LPAF-based ISAR imaging algorithm for targets with complex motion, we provide the results of the LPAF-based ISAR imaging algorithm under a range of different SNRs in Fig. 8, where the scope of ${\mathrm{SNR}}_{\mathrm{in}}$ is from 2 to $-8\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$ with an increment of $-2\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$. It is worthwhile to mention that the less spurious scatterers appear in Figs. 8(a)–8(c), and the entropies of Figs. 8(a)–8(c) are small than those of Figs. 8(d)–8(f) in Table 8. Especially, the images in Figs. 8(e) and 8(f) are blurred severely by spurious scatterers. Thus, we conclude that the proposed ISAR imaging algorithm fails when the input SNR falls below $-4\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{dB}$.

## Table 8

Entropies of ISAR images in Fig. 8.

Fig. 8(a) | Fig. 8(b) | Fig. 8(c) | Fig. 8(d) | Fig. 8(e) | Fig. 8(f) | |
---|---|---|---|---|---|---|

Entropies | 5.6527 | 5.7258 | 5.8516 | 6.5244 | 7.9125 | 8.9486 |

## 7.

## Conclusion

For ISAR imaging of targets with complex motion, the azimuth echo in a range bin can be modeled as multicomponent CPSs, which have been proved more appropriate than the LFM model. This paper presents a new ISAR imaging algorithm for targets with complex motion based on LPAF. LPAF can directly obtain CRQCRD only by FFT operation, where two parameters can be estimated simultaneously via peak detection. With the analyses and the simulations, we demonstrate that the LPAF outperforms MCPF and LPWD for the parameter estimation of CPS in the robustness of performance against noise, the computational cost, and the cross-term suppression. Finally, the effectiveness of the proposed ISAR imaging algorithm is validated by the experimental results of simulated data.

## Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant Nos. 61271024 and 61201292 and in part by the Science and Technology Foundation of Shaanxi Province under Grant 2012JM8015. Qian Lv conceived the work in this paper that led to the submission, experiment design, and accomplishment of the writing of the manuscript. Jibin Zheng played an important role in revising the manuscript and providing English language support. Tao Su is mainly responsible for approving the final version.

## References

F. Berizzi et al., “High-resolution ISAR imaging of maneuvering targets by means of the range instantaneous Doppler technique: modeling and performance analysis,” IEEE Trans. Image Process. 10(12), 1880–1890 (2001).IIPRE41057-7149http://dx.doi.org/10.1109/83.974573Google Scholar

Y. Li et al., “Inverse synthetic aperture radar imaging of ship target with complex motion,” IET Radar Sonar Navig. 2(6), 395–403 (2008).http://dx.doi.org/10.1049/iet-rsn:20070101Google Scholar

J. Wood and D. Barry, “Radon transformation of time-frequency distributions for analysis of multicomponent signals,” IEEE Trans. Signal Process. 42(11), 3166–3177 (1994).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.330375Google Scholar

Y. Wang, “Inverse synthetic aperture radar imaging of manoeuvring target based on range-instantaneous-Doppler and range-instantaneous-chirp-rate algorithms,” IET Radar Sonar Navig. 6(9), 921–928 (2012).http://dx.doi.org/10.1049/iet-rsn.2012.0091Google Scholar

J. Zheng et al., “ISAR imaging of targets with complex motions based on the keystone time-chirp rate distribution,” IEEE Geosci. Remote Sens. Lett. 11(7), 1275–1279 (2014).http://dx.doi.org/10.1109/LGRS.2013.2291992Google Scholar

M. Xing et al., “New ISAR imaging algorithm based on modified Wigner-Ville distribution,” IET Radar Sonar Navig. 3(1), 70–80 (2009).http://dx.doi.org/10.1049/iet-rsn:20080003Google Scholar

X. Lv et al., “ISAR imaging of maneuvering targets based on the range centroid Doppler technique,” IEEE Trans. Image Process. 19(1), 141–153 (2010).IIPRE41057-7149http://dx.doi.org/10.1109/TIP.2009.2032892Google Scholar

X. Lv et al., “Lv’s distribution: principle, implementation, properties, and performance,” IEEE Trans. Signal Process. 59(8), 3576–3591 (2011).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2011.2155651Google Scholar

R. Tao, N. Zhang and Y. Wang, “Analysing and compensating the effects of range and Doppler frequency migrations in linear frequency modulation pulse compression radar,” IET Radar, Sonar Navig. 5(1), 12–22 (2011).http://dx.doi.org/10.1049/iet-rsn.2009.0265Google Scholar

X. Bai et al., “ISAR imaging of a ship target based on parameter estimation of multicomponent quadratic frequency-modulated signals,” IEEE Trans. Geosci. Remote Sens. 52(2), 1418–1429 (2014).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2013.2251348Google Scholar

J. Zheng et al., “ISAR imaging of targets with complex motion based on the chirp rate-quadratic chirp rate distribution,” IEEE Trans. Geosci. Remote Sens. 52(11), 7276–7289 (2014).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2014.2310474Google Scholar

Y. Li et al., “A new algorithm of ISAR imaging for maneuvering targets with low SNR,” IEEE Trans. Aerosp. Electron. Syst. 49(1), 543–557 (2013).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2013.6404119Google Scholar

J. Zheng et al., “ISAR imaging of nonuniformly rotating target based on a fast parameter estimation algorithm of cubic phase signal,” IEEE Trans. Geosci. Remote Sens. 53(9), 4727–4740 (2015).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2015.2408350Google Scholar

Y. Wang and B. Zhao, “Inverse synthetic aperture radar imaging of nonuniformly rotating target based on the parameters estimation of multicomponent quadratic frequency-modulated signals,” IEEE Sensors J. 15(7), 4053–4061 (2015).ISJEAZ1530-437Xhttp://dx.doi.org/10.1109/JSEN.2015.2409884Google Scholar

P. O’shea, “A fast algorithm for estimating the parameters of a quadratic FM signal,” IEEE Trans. Signal Process. 52(2), 385–393 (2004).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2003.821097Google Scholar

P. O’shea, “A new technique for instantaneous frequency rate estimation,” IEEE Signal Process. Lett. 9(8), 251–252 (2002).http://dx.doi.org/10.1109/LSP.2002.803003Google Scholar

S. Peleg and B. Porat, “Linear FM signal parameter estimation from discrete-time observations,” IEEE Trans. Aerosp. Electron. Syst. 27(4), 607–616 (1991).IEARAX0018-9251http://dx.doi.org/10.1109/7.85033Google Scholar

Y. Wang and Y. Jiang, “Inverse synthetic aperture radar imaging of maneuvering target based on the product generalized cubic phase function,” IEEE Geosci. Remote Sens. Lett. 8(5), 958–962 (2011).http://dx.doi.org/10.1109/LGRS.2011.2143387Google Scholar

Y. Wang and Y. Jiang, “ISAR imaging of a ship target using product high-order matched-phase transform,” IEEE Geosci. Remote Sens. Lett. 6(4), 658–661 (2009).http://dx.doi.org/10.1109/LGRS.2009.2013876Google Scholar

L. Stankovic, “Local polynomial Wigner distribution,” Signal Process. 59(1), 123–128 (1997).http://dx.doi.org/10.1016/S0165-1684(97)00086-8Google Scholar

Y. Wang, J. Kang and Y. Jiang, “ISAR imaging of maneuvering target based on the local polynomial Wigner distribution and integrated high-order ambiguity function for cubic phase signal model,” IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 7(7), 2971–2991 (2014).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/JSTARS.2014.2301158Google Scholar

Y. Wang, B. Zhao and J. Kang, “Asymptotic statistical performance of local polynomial Wigner distribution for the parameters estimation of cubic-phase signal with application in ISAR imaging of ship target,” IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 8(3), 1087–1098 (2015).http://dx.doi.org/10.1109/JSTARS.2014.2355219Google Scholar

L. Wu et al., “ISAR imaging of targets with complex motion based on discrete chirp Fourier transform for cubic chirps,” IEEE Trans. Geosci. Remote Sens. 50(10), 4201–4212 (2012).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2012.2189220Google Scholar

M. Xing, R. Wu and Z. Bao, “High resolution ISAR imaging of high speed moving targets,” IEE Proc. Radar Sonar Navig. 152(2), 58–67 (2005).http://dx.doi.org/10.1049/ip-rsn:20045084Google Scholar

D. Wahl et al., “Phase gradient autofocus—a robust tool for high resolution SAR phase correction,” IEEE Trans. Aerosp. Electron. Syst. 30(3), 827–835 (1994).IEARAX0018-9251http://dx.doi.org/10.1109/7.303752Google Scholar

Y. Li et al., “Inverse synthetic aperture radar imaging of targets with nonsevere maneuverability based on the centroid frequency chirp rate distribution,” J. Appl. Remote Sens. 9, 095065 (2015).http://dx.doi.org/10.1117/1.JRS.9.095065Google Scholar

S. Barbarossa, A. Scaglione and G. B. Giannakis, “Product high-order ambiguity function for multicomponent polynomial-phase signal modeling,” IEEE Trans. Signal Process. 46(3), 691–708 (1998).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.661336Google Scholar

X. Xia, “Discrete chirp-Fourier transform and its application to chirp rate estimation,” IEEE Trans. Signal Process. 48(11) 3122–3133 (2000).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.875469Google Scholar

X. Guo et al., “Comments on discrete chirp-Fourier transform and its application to chirp rate estimation,” IEEE Trans. Signal Process. 50(12), 3115–3116 (2002).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2002.805492Google Scholar

S. Peleg and B. Porat, “The Cramer-Rao lower bound for signals with constant amplitude and polynomial phase,” IEEE Trans. Signal Process. 39(3), 749–752 (1991).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.80864Google Scholar

P. Wang et al., “Integrated cubic phase function for linear FM signal analysis,” IEEE Trans. Aerosp. Electron. Syst. 46(3), 963–977 (2010).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2010.5545167Google Scholar

## Biography

**Qian Lv** received her BS degree in measuring and control techniques and instruments from Xi’an Shiyou University, Shaanxi, China, in 2013. She is currently working toward her PhD with the National Laboratory of Radar Signal Processing, Xidian University, Xi’an, China. Her research interests include synthetic aperture radar (SAR) and inverse SAR signal processing, time-frequency analysis, and interferometric inverse synthetic aperture radar (InISAR) image.

**Tao Su** received his BS degree in information theory, his MS degree in mobile communication, and his PhD in signal and information processing from Xidian University, Xi’an, China, in 1990, 1993, and 1999, respectively. Since 1993, he has been a professor with the National Laboratory of Radar Signal Processing, School of Electronic Engineering. His research interests include high-speed real-time signal processing on radar, sonar and telecommunications, digital signal processing, parallel processing system design, and FPGA IP design.

**Jibin Zheng** received his BS degree in electronic information science and technology from Shandong Normal University, Shandong, China, in 2009 and his PhD in signal and information processing from Xidian University, Xi’an, China, in 2015. From September 2012 to September 2014, he worked as a visiting PhD student at the Department of Electrical Engineering, Duke University, Durham, North Carolina. His research interests include SAR and inverse SAR signal processing, cognitive radar.