27 June 2014 Velocity measurement and radar imaging of high-speed target by using stepped-chirp signal
Abstract
Stepped-chirp signal (SCS) is widely used by wideband high-resolution radars; however, its bandwidth synthesizing suffers from motion-induced phase errors resulting from the radial velocity of the target, especially for high-speed targets. So motion compensation is very crucial in signal processing for this kind of radar using SCS. Based on the cross-correlation inner frame method for velocity measurement of a single-scattering-center target, the multiple cross-correlation method (MCCM) is proposed for measuring the velocity of a complex target with multiple scattering centers (MSC), which is called SCS-MCCM. By this algorithm, not only the radial velocity, but also the real velocity of a target can be measured under the assumption that it is straight moving. After obtaining the velocity, motion compensation is conducted so as to realize high-resolution imaging. Theoretical analysis and simulations show that the proposed method is feasible and effective for an MSC target.

## CCIFM

The transmitted waveform consists of a series of chirp pulses (subchirps or subpulses), whose carrier frequencies increase (or decrease) pulse by pulse in a step of $△f$. If the initial distance between radar and the target is $R$, the radial target velocity is $V$, and the pulse repeat interval is ${T}_{r}$, then the time delay of the $i$’th subpulse is

## (1)

${\tau }_{i}=\frac{2R}{c}-\frac{2iV{T}_{r}}{c}.$

The echo of the $i$’th subpulse after being coherently received is

## (2)

$s\left(t,i\right)=u\left(t,i\right)\mathrm{exp}\left[-j2\pi \left({f}_{0}+i\mathrm{\Delta }f\right)\left(\frac{2R}{c}-\frac{2iV{T}_{r}}{c}\right)\right],$
where $u\left(t,i\right)=\mathrm{rect}\left[\left(t-{\tau }_{i}\right)/{T}_{p}\right]\mathrm{exp}\left[jK\pi {\left(t-{\tau }_{i}\right)}^{2}\right]$ is the baseband chirp pulse, ${f}_{0}$ is the carrier frequency of the first pulse, $△f$ is the frequency step, $K$ is the slope of the subpulse, ${T}_{p}$ is the subpulse time duration, $i=0,1\dots N-1$, and $N$ is the number of subpulses in a burst. In the following, we ignore $u\left(t,i\right)$ because it has no significant role in the derivation, as will be shown in the following. After cross-correlation processing on two subpulses spaced by $p$ subpulses, we can get

## (3)

$X\left(i\right)=s\left(i\right)·{s}^{*}\left(i+p\right)\phantom{\rule{0ex}{0ex}}=\mathrm{exp}\left(-j2\pi {f}_{0}p{T}_{r}\frac{2V}{c}\right)\mathrm{exp}\left(-j2\pi {p}^{2}\mathrm{\Delta }f{T}_{r}\frac{2V}{c}\right)\mathrm{exp}\left(j2\pi p\mathrm{\Delta }f\frac{2R}{c}\right)\mathrm{exp}\left(-j2\pi 2p\mathrm{\Delta }f\frac{2V}{c}i{T}_{r}\right),\phantom{\rule{0ex}{0ex}}i=0,1\dots N-p, 1\le p\le N-1.$

Equation (3) shows that $X\left(i\right)$ contains four terms: the first three terms are constant phases and the last term is a linear item of $i$, which indicates a frequency shift. According to Eq. (3), the Doppler frequency of a moving target can be formulated as19

## (4)

${f}_{d}=-2p{T}_{r}\mathrm{\Delta }f\frac{2V}{c},$
which is proportional to $V$. So the target’s velocity can be measured by fast Fourier transform (FFT) of $X\left(i\right)$. The velocity measurement accuracy is related to the resolution of Doppler frequency. If we increase the FFT dots to ${N}_{0}$ by zero-padding, the minimum measurable velocity is then

## (5)

$\mathrm{\Delta }V=\frac{c}{4{N}_{0}p{T}_{r}\mathrm{\Delta }f},$
while the velocity measurement scope is

## (6)

$V\in \left(-\frac{c}{8p{T}_{r}\mathrm{\Delta }f},\frac{c}{8p{T}_{r}\mathrm{\Delta }f}\right).$

Now, we shall show why $u\left(t,i\right)$ can be ignored in Eq. (3) in detail as follows. We can easily obtain the following expression:

## (7)

$u\left(t,i\right)·{u}^{*}\left(t,i+p\right)\phantom{\rule{0ex}{0ex}}=\mathrm{exp}\left\{jK\pi {\left[t-\left(\frac{2R}{c}-\frac{2iV{T}_{r}}{c}\right)\right]}^{2}\right\}\mathrm{exp}\left(-jK\pi {\left\{t-\left[\frac{2R}{c}-\frac{2\left(i+p\right)V{T}_{r}}{c}\right]\right\}}^{2}\right)\phantom{\rule{0ex}{0ex}}=\mathrm{exp}\left(-j\pi K\frac{4{p}^{2}{V}^{2}{T}_{r}^{2}}{{c}^{2}}\right)\mathrm{exp}\left(j\pi K\frac{8pV{T}_{r}}{{c}^{2}}R\right)\mathrm{exp}\left(-j2\pi t\frac{2pKV{T}_{r}}{c}\right)\mathrm{exp}\left(-j2\pi \frac{4pK{V}^{2}{T}_{r}}{{c}^{2}}i{T}_{r}\right).\phantom{\rule{0ex}{0ex}}$

For argument $i$, the first three terms in Eq. (7) are constant phases and the last term is a linear item, which indicates a frequency, i.e.,

## (8)

${f}_{d}^{\prime }=-\frac{4pK{V}^{2}{T}_{r}^{2}}{{c}^{2}}.$

Generally speaking, ${f}_{d}^{\prime }$ is much smaller than ${f}_{d}$, and the ratio $\left({f}_{d}^{\prime }/{f}_{d}\right)=\left({\mathrm{KVT}}_{r}/c\mathrm{\Delta }f\right)$ is usually ${<10}^{-3}$, so it can be ignored.

## SCS-MCCM

When the target is complex with MSCs, the correlation signals become much more complicated. Let us assume the number of scattering centers is $M$; the echo of the $i$’th subpulse is thus

## (9)

$s\left(t,i\right)=\sum _{m=1}^{M}\mathrm{exp}\left\{jK\pi {\left[t-\left(\frac{2{R}_{m}}{c}-\frac{2iV{T}_{r}}{c}\right)\right]}^{2}\right\}\mathrm{exp}\left[-j2\pi \left({f}_{0}+i\mathrm{\Delta }f\right)\left(\frac{2{R}_{m}}{c}-\frac{2iV{T}_{r}}{c}\right)\right].$

Because $s\left(t,i\right)$ is the sum of echoes from $M$ scatterers, the correlation signals will have self-terms (self-correlation of the same scatterer’s echo) and cross-terms (cross-correlation of two different scatterers’ echoes), which lead to different frequency shifts. Therefore, we cannot estimate the velocity just by FFT as is done in Sec. 2.

The cross-correlation between $s\left(t,i\right)$ and $s\left(t,i+p\right)$ is formulated as

## (10)

$X\left(t,i\right)=s\left(t,i\right)·{s}^{*}\left(t,i+p\right)\phantom{\rule{0ex}{0ex}}=\sum _{m,n=1}^{M}{\varphi }_{1}·{\varphi }_{2}·{\varphi }_{3}·{\varphi }_{4}·{\varphi }_{5}·{\varphi }_{6}·{\varphi }_{7}·{\varphi }_{8}·{\varphi }_{9}·{\varphi }_{10}·{\varphi }_{11}·{\varphi }_{12}·{\varphi }_{13},$
where
${\varphi }_{1}=\mathrm{exp}\left(-j2\pi {f}_{0}p{T}_{r}\frac{2V}{c}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{2}=\mathrm{exp}\left(-j2\pi {p}^{2}\mathrm{\Delta }f{T}_{r}\frac{2V}{c}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{3}=\mathrm{exp}\left(-jK\pi \frac{4{p}^{2}{V}^{2}{T}_{r}^{2}}{{c}^{2}}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{4}=\mathrm{exp}\left(-j2\pi \frac{2pKV{T}_{r}}{c}t\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{5}=\mathrm{exp}\left(j2\pi \frac{2K\mathrm{\Delta }{R}_{m,n}}{c}t\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{6}=\mathrm{exp}\left(-j2\pi 2p\mathrm{\Delta }f\frac{2V}{c}i{T}_{r}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{7}=\mathrm{exp}\left(j2\pi \mathrm{\Delta }f\frac{2\mathrm{\Delta }{R}_{m,n}}{c{T}_{r}}i{T}_{r}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{8}=\mathrm{exp}\left(-j2\pi \frac{4pK{V}^{2}{T}_{r}}{{c}^{2}}i{T}_{r}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{9}=\mathrm{exp}\left(j2\pi \frac{4KV\mathrm{\Delta }{R}_{m,n}}{{c}^{2}}i{T}_{r}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{10}=\mathrm{exp}\left(-j\pi K\frac{4{R}_{m}^{2}-4{R}_{n}^{2}}{{c}^{2}}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{11}=\mathrm{exp}\left(j\pi K\frac{8pV{T}_{r}{R}_{m}}{{c}^{2}}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{12}=\mathrm{exp}\left(j2\pi p\mathrm{\Delta }f\frac{2{R}_{m}}{c}\right)\phantom{\rule{0ex}{0ex}}{\varphi }_{13}=\mathrm{exp}\left(j2\pi {f}_{0}\frac{2\mathrm{\Delta }{R}_{m,n}}{c}\right)\phantom{\rule{0ex}{0ex}}\mathrm{\Delta }{R}_{m,n}={R}_{m}-{R}_{n},1$

Among the above phase terms, ${\varphi }_{1}$, ${\varphi }_{2}$, and ${\varphi }_{3}$ are constant phases, ${\varphi }_{10}$, ${\varphi }_{11}$, ${\varphi }_{12}$, and ${\varphi }_{13}$ are constant phases related to $R$, ${\varphi }_{4}$ and ${\varphi }_{5}$ are linear items about $t$, and ${\varphi }_{6}$, ${\varphi }_{7}$, ${\varphi }_{8}$, and ${\varphi }_{9}$ are linear items about $i$. As is shown, there are no coupling items about $t$ and $i$.

So, $X\left(t,i\right)$ has ${M}^{2}$ terms, including $M$ self-terms(when $m=n$) and $M\left(M-1\right)$ cross-terms(when $m\ne n$), which have different linear coefficients. That means the two-dimensional (2-D) spectrum of $X\left(t,i\right)$ has many frequency components, and each of them corresponds to a self-term or a cross-term of $X\left(t,i\right)$.

For self-terms, ${\varphi }_{5}$, ${\varphi }_{7}$, ${\varphi }_{9}$, ${\varphi }_{10}$, and ${\varphi }_{13}$ do not exist. As ${\varphi }_{8}$ is very small, as has been discussed previously, it can be ignored. With the constant phases ${\varphi }_{1}$, ${\varphi }_{2}$, and ${\varphi }_{3}$ further ignored, we get

## (11)

$X{\left(t,i\right)}_{\mathrm{self}}={\varphi }_{4}·{\varphi }_{6}\sum _{m=1}^{M}{\varphi }_{11}·{\varphi }_{12}.$

## (12)

$X{\left(t,i\right)}_{\mathrm{cross}}={\varphi }_{4}·{\varphi }_{6}\sum _{m,n=1,m\ne n}^{M}{\varphi }_{5}·\left({\varphi }_{7}·{\varphi }_{9}\right)·\left({\varphi }_{10}·{\varphi }_{11}·{\varphi }_{12}·{\varphi }_{13}\right).$

Let ${f}_{t}$ represent the frequency along the $t$ direction and ${f}_{d}$ represent the frequency along the $i$ direction. Because ${\varphi }_{4}$ (the linear item about $t$) and ${\varphi }_{6}$ (the linear item about $i$) are not related to the scatterer’s position ${R}_{m}$, ${f}_{t}$ and ${f}_{d}$ of the $M$ self-terms should be the same, and they are denoted as ${f}_{t}^{0}$ and ${f}_{d}^{0}$, respectively.

## (13)

${f}_{t}^{0}=-\frac{2pKV{T}_{r}}{c},$

## (14)

${f}_{d}^{0}=-2p\mathrm{\Delta }f\frac{2V}{c}$

The $M\left(M-1\right)$ cross-terms have different ${f}_{t}$ and different ${f}_{d}$, which are related to the relative positions of scattering centers $\mathrm{\Delta }{R}_{m,n}$. As shown in Eq. (12), ${f}_{t}$ is decided by ${\varphi }_{4}$ and ${\varphi }_{5}$, and ${f}_{d}$ is decided by ${\varphi }_{6}$, ${\varphi }_{7}$, and ${\varphi }_{9}$, i.e.,

## (15)

${f}_{t}={f}_{t}^{0}+\frac{2K\mathrm{\Delta }{R}_{m,n}}{c}\triangleq {f}_{t}^{0}+{f}_{t}^{1},$

## (16)

${f}_{d}={f}_{d}^{0}+\mathrm{\Delta }f\frac{2\mathrm{\Delta }{R}_{m,n}}{c{T}_{r}}+\frac{4KV\mathrm{\Delta }{R}_{m,n}}{{c}^{2}}\triangleq {f}_{d}^{0}+{f}_{d}^{1}.$

From Eqs. (15) and (16), we can see that ${f}_{t}$ and ${f}_{d}$ are symmetrical pairs about ${f}_{t}^{0}$ and ${f}_{d}^{0}$, respectively.

When ${f}_{d}^{1}$ is less than a frequency resolution cell, i.e., $|{f}_{d}^{1}|<\left[1/\left({T}_{r}·{N}_{0}\right)\right]$, we can have ${f}_{d}\approx {f}_{d}^{0}$. In this situation, the target can be regarded as a point target, so the velocity can be measured using the algorithm mentioned in Sec. 2. It means that the size of target should meet the requirement of $|\mathrm{\Delta }{R}_{m,n}|<\left[c/\left(2\mathrm{\Delta }f{N}_{0}\right)\right]$, which is usually not satisfied in practical situations. That is to say the performance of the algorithm of Sec. 2 cannot be as good as we expect for bigger targets.

In fact, we can also see from Eqs. (15) and (16) that the frequencies of cross-terms ${f}_{t}$ and ${f}_{d}$ are determined not only by $V$, but also by $\mathrm{\Delta }{R}_{m,n}$. However, from Eqs. (13) and (14), we can see that the frequencies of self-terms ${f}_{t}^{0}$ and ${f}_{d}^{0}$ are determined only by $V$. Therefore, if we can differentiate them in the 2-D spectrum of $X\left(t,i\right)$, the velocity can then be estimated by ${f}_{t}^{0}$ and ${f}_{d}^{0}$.

If a signal is sampled at rate ${f}_{s}$ and the sampling length is $L$, then the resolution of the spectrum is ${f}_{s}/L$. When $|{f}_{t}^{1}|>{f}_{s}/L$, i.e., the distances between scattering centers of the target satisfy the following condition:

## (17)

$|\mathrm{\Delta }{R}_{m,n}|>\frac{c{f}_{s}}{2KL},\phantom{\rule[-0.0ex]{2em}{0.0ex}}1

After performing FFT on $X\left(t,i\right)$ along the $t$ direction, ${f}_{t}$ and ${f}_{t}^{0}$ will be at least one frequency cell apart, so they can be easily distinguished. In fact, Eq. (17) is easy to meet via increasing $L$ by interpolation. Since ${f}_{t}$ is symmetrically distributed around ${f}_{t}^{0}$, we can just choose the central frequency of the spectrum to get rid of all other frequency components of the cross-terms. At the last step, we get ${f}_{d}^{0}$ by conducting FFT along the $i$ direction and, finally, obtain the velocity by

## (18)

$\stackrel{^}{V}=-\frac{c}{4p\mathrm{\Delta }f}{f}_{d}^{0},$
where $p\in \left[1,N-1\right]$. It is clear from Eq. (18) that the larger the $p$, the higher the measurement accuracy but the less the unambiguous measurement range of velocity. We must choose an appropriate $p$ to balance the measurement range and accuracy. In fact, we can use a smaller $p$ to conduct a coarse measurement with a large unambiguous range first and then use a larger $p$ to conduct an accurate measurement.

Figure 1 presents the flow chart of the proposed method, including the motion compensation step. The first step is to calculate a series of cross-correlations between two paired subpulses spaced by $p$ subpulses ($i=1$, $N-p$; the second step is to perform FFT on the correlation signals along the $t$ direction for $i=1$, $N-p$; the third step is to take the central spectral line of each pair to form a vector; the fourth step is to perform FFT on the above vector; the fifth step is to find the frequency corresponding to the peak value of the spectrum; the sixth step is to calculate the radial velocity according to Eq. (18); the seventh step is to conduct motion compensation according to the obtained velocity; the last step is to coherently synthesize the spectrums of all subpulses using the algorithm of Ref. 20 and, finally, get the high-resolution range profile through Inverse FFT (IFFT).

## Fig. 1

The flow chart of stepped-chirp signal multiple cross-correlation method (MCCM).

## Real Velocity Measurement

In this section, we shall show that besides the radial velocity, SCS-MCCM can also be used to measure the real velocity of the target. SCS-MCCM measures the radial velocity from the echoes of subpulses within a burst, i.e., it measures the instantaneous radial velocity at each burst. Then the real velocity can be estimated by combining the radial velocity measurement with the range measurement. Here, the target is supposed to be well tracked by radar.

We assume that the target moves straight during the observation time, as shown in Fig. 2. The real target velocity is ${V}_{0}$. When the radar is transmitting the $j$’th burst ($j=1,2,3\dots$), the distance between radar and target is ${R}_{j}$, and the angle between the target moving direction and radar boresight is ${\theta }_{j}$, so the radial velocity is ${V}_{j}={V}_{0}\mathrm{cos}{\theta }_{j}$ (${V}_{j}$ is positive when ${\theta }_{j}<90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$ and negative when ${\theta }_{j}>90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$).

## Fig. 2

The observation geometry.

If the burst repeat interval is ${T}_{b}$, then the distance between target position 1 and target position $j$ is ${V}_{0}\left(j-1\right){T}_{b}$. According to the cosine law, we can get Eqs. (19) and (20), respectively, as follows:

## (19)

$\mathrm{cos}\text{\hspace{0.17em}}{\theta }_{1}=\frac{{R}_{1}^{2}+{\left[{V}_{0}\left(j-1\right){T}_{b}\right]}^{2}-{R}_{j}^{2}}{2{R}_{1}{V}_{0}\left(j-1\right){T}_{b}},$

## (20)

$\mathrm{cos}\left(\pi -{\theta }_{j}\right)=\frac{{R}_{j}^{2}+{\left[{V}_{0}\left(j-1\right){T}_{b}\right]}^{2}-{R}_{1}^{2}}{2{R}_{j}{V}_{0}\left(j-1\right){T}_{b}}.$

By substituting $\mathrm{cos}\text{\hspace{0.17em}}{\theta }_{1}=\left({V}_{1}/{V}_{0}\right)$ into Eq. (19), we can get

## (21)

${R}_{j}^{2}-{R}_{1}^{2}={\left[{V}_{0}\left(j-1\right){T}_{b}\right]}^{2}-2{R}_{1}{V}_{1}\left(j-1\right){T}_{b}.$

By substituting $\mathrm{cos}\left(\pi -{\theta }_{j}\right)=-\left({V}_{j}/{V}_{0}\right)$ into Eq. (20), we can get

## (22)

${R}_{j}^{2}-{R}_{1}^{2}=-2{R}_{j}{V}_{j}\left(j-1\right){T}_{b}-{\left[{V}_{0}\left(j-1\right){T}_{b}\right]}^{2}.$

Finally, ${V}_{0}$ can be derived simultaneously from Eqs. (21) and (22).

## (23)

${V}_{0}=\sqrt{\frac{{R}_{1}{V}_{1}-{R}_{j}{V}_{j}}{\left(j-1\right){T}_{b}}}.$

Equation (23) indicates that the real velocity measurement can be easily obtained by using the results of the range measurement as well as the radial velocity measurement, where the range measurement can be done according to the time delay of the radar echoes.

In the following, let us analyze the impact of measurement errors of ${R}_{j}$ and ${V}_{j}$ on ${V}_{0}$ estimation. The measurement error (or accuracy) of ${V}_{0}$, i.e., $\mathrm{\Delta }{V}_{0}$, can be formulated as follows:

## (24)

$\mathrm{\Delta }{V}_{0}=\sqrt{{\left(\frac{\partial {V}_{0}}{\partial {R}_{j}}\mathrm{\Delta }{R}_{j}\right)}^{2}+{\left(\frac{\partial {V}_{0}}{\partial {V}_{j}}\mathrm{\Delta }{V}_{j}\right)}^{2}},$
where $\mathrm{\Delta }{R}_{j}$ and $\mathrm{\Delta }{V}_{j}$ are the measurement errors of ${R}_{j}$ and ${V}_{j}$, respectively. Equation (24) shows that when $\mathrm{\Delta }{R}_{j}$ and $\mathrm{\Delta }{V}_{j}$ are fixed (as discussed above), the smaller the partial derivatives with respect to ${R}_{j}$ and ${V}_{j}$, the smaller the measurement error of ${V}_{0}$, i.e., the higher the measurement accuracy of ${V}_{0}$.

The partial derivatives of Eq. (23) with respect to ${R}_{j}$ and ${V}_{j}$, respectively, are as follows:

## (25)

$\frac{\partial {V}_{0}}{\partial {R}_{j}}=-\frac{{V}_{j}}{2\sqrt{\left({R}_{1}{V}_{1}-{R}_{j}{V}_{j}\right)\left(j-1\right){T}_{b}}},$

## (26)

$\frac{\partial {V}_{0}}{\partial {V}_{j}}=-\frac{{R}_{j}}{2\sqrt{\left({R}_{1}{V}_{1}-{R}_{j}{V}_{j}\right)\left(j-1\right){T}_{b}}}.$

In Eqs. (25) and (26), when ${\theta }_{j}<90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, ${V}_{j}$ is positive. Along with the increase in $j$, $\left(j-1\right){T}_{b}$ becomes larger, while ${V}_{j}$ and ${R}_{j}$ become smaller, so $\left({R}_{1}{V}_{1}-{R}_{j}{V}_{j}\right)$ is positive and becomes larger, too; thus, the absolute values of $\partial {V}_{0}/\partial {R}_{j}$ and $\partial {V}_{0}/\partial {V}_{j}$ both become smaller. When ${\theta }_{j}>90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, ${V}_{j}$ is negative, and its absolute value becomes larger as $j$ increases. Along with increase in $j$, $\left(j-1\right){T}_{b}$ and ${R}_{j}$ also become larger, so $\left({R}_{1}{V}_{1}-{R}_{j}{V}_{j}\right)$ is still positive and becomes larger, too; once again, the absolute values of $\partial {V}_{0}/\partial {R}_{j}$ and $\partial {V}_{0}/\partial {V}_{j}$ both become smaller.

From the above analysis, we observed that the larger the $j$, the smaller $\partial {V}_{0}/\partial {R}_{j}$ and $\partial {V}_{0}/\partial {V}_{j}$. Therefore, in order to achieve high measurement accuracy of ${V}_{0}$, we should choose a $j$ as large as possible.

## Simulations

The simulation parameters are as follows: ${f}_{0}=13\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GHz}$, $△f=100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MHz}$, $B=120\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MHz}$, ${T}_{r}=16\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{us}$, $N=20$, $p=6$. Monte-Carlo simulations are conducted 100 times when the target velocities are produced with a uniform distribution in ($5000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $5500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$). For the SSC target, CCIFM and SCS-MCCM have the same performance for the velocity measurement. Figure 3 presents the error curves of the velocity measurements for the SCC target with different SNRs. Results show that in the case of $\mathrm{SNR}=20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$, the maximum measurement error is $<4\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, while in the case of $\mathrm{SNR}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$, the maximum measurement error is $<5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. That is to say the measurement error does not change significantly as the SNR decreases, i.e., MCCM exhibits a satisfied antinoise performance for the SSC target.

## Fig. 3

The velocity measurement errors for single-scattering-center (SSC) target.

Figures 4(a) and 4(b) present the high-resolution range profiles without and with motion compensation, respectively. When the target is moving at a high speed, the motion-induced phase errors will lead to image displacement and the image will be out of focus; therefore, high resolution cannot be achieved, as clearly shown in Fig. 4(a). Figure 4(b) shows that after velocity measurement and motion compensation, the burst of stepped-chirp can be synthesized correctly and a high-resolution range profile is, thus, obtained. The resolution in Fig. 4(b) is $\sim 0.075\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$, which is close to the theoretical value of a 2.04-GHz bandwidth signal. In the following, we shall demonstrate the good performance of the SCS-MCCM for an MSC target. Figure 5 gives the MSC target model used for the velocity measurement and imaging simulations, which is composed of 41 ideal point targets.

## Fig. 4

The synthesized range profiles for SSC target (a) without motion compensation and (b) with motion compensation.

## Fig. 5

The multiple scattering centers (MSC) target model.

We conduct Monte-Carlo simulation 100 times as before for the MSC target of Fig. 5 by using CCIFM with an $\mathrm{SNR}=20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$ and plot the obtained velocity measurement errors in Fig. 6, from which one can clearly see that they are very large and fluctuate dramatically. It means that CCIFM breaks down for the MSC target velocity measurement. Figure 7 presents the error curves of the velocity measurement by SCS-MCCM for the same MSC target with different SNRs. The target velocities in Figs. 7(a), 7(b), and 7(c) are varied in the range of ($2000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $2500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), ($5000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $5500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), and ($7000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $7500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), respectively. We can see from Fig. 7 that the measurement errors of the SCS-MCCM do not change as the target velocity changes, and they are all within $6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ when $\mathrm{SNR}=20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$, while the maximum measurement errors are $\sim 17\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ when $\mathrm{SNR}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$. Compared to the SSC target case, the measurement accuracy for the MSC target degrades remarkably as the SNR decreases. It means the performance is more sensitive to SNR. However, the accuracy still satisfies the requirement for motion compensation processing. Berizzi and Martorella derived the velocity accuracy requirement for SCS signal motion compensation12 as $\mathrm{\Delta }V<\left[c/\left(4N{f}_{0}{T}_{r}\right)\right]$, i.e., it is $\sim 18\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ with the parameters of this paper.

## Fig. 6

The velocity measurement error of cross-correlation inner frame method (CCIFM) for MSC target (SNR = 20 dB).

## Fig. 7

The velocity measurement errors of MCCM for MSC target: (a) $V\in$ ($2000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $2500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), (b) $V\in$ ($5000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $5500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), and (c) $V\in$ ($7000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, $7500\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$).

Figures 8 and 9 present the stepped frequency radar imaging results for the target of Fig. 5 using the system parameters as listed in Table 1 under the geometry of Fig. 2, where the initial distance ${R}_{1}$ between the radar and the target is 10 km, the initial view angle ${\theta }_{1}=60\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, and the real target velocity ${V}_{0}=7000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$.

## Fig. 8

The synthesized range profiles for MSC target (a) without motion compensation and (b) with motion compensation using CCIFM, and (c) with motion compensation using MCCM.

## Fig. 9

The imaging results (a) without motion compensation and (b) with motion compensation using CCIFM, and (c) with motion compensation using MCCM.

## Table 1

Simulation parameters.

 Initial carrier frequency 13 GHz Frequency step 100 MHz Bandwidth 120 MHz Number of pulses 20 Pulse width 8 us Subpulse repetition interval 16 us Burst repetition interval 2 ms Number of burst 128 SNR 20 dB

Figure 8(a) presents the synthesized range profile without motion compensation, from which one can see that the high-speed motion-induced phase errors make the range profile seriously smear. Figure 8(b) shows the result with motion compensation conducted using CCIFM. Because CCIFM cannot correctly measure the velocity of the MSC target, Fig. 8(b) scarcely improves compared to Fig. 8(a). Figure 8(c) gives the synthesized range profile with motion compensation conducted using SCS-MCCM; it is very well profiled as we expected. Figure 9(a) shows the 2-D imaging result without motion compensation. Due to the high speed of the target, the motion-induced phase errors make the image seriously blurred. Figures 9(b) and 9(c) are the images with motion compensation conducted using CCIFM and SCS-MCCM, respectively. As before, one can see that CCIFM does not work in the MSC target case. However, SCS-MCCM can handle this case and the radar image is very well focused, as shown in Fig. 9(c). The results of Fig. 9 also indicate that the high-speed motion-induced errors affect not only the range profiles, but also azimuth focusing. The above simulation experiments demonstrate that the proposed velocity measurement method for a high-speed target based on SCS is valid and effective, and the resulting measurement error can meet the requirement for high-resolution imaging.

Figures 10 and 11 present the real velocity measurement results using the same parameters and geometry as in the previous imaging simulation. The solid curves represent the estimated values versus the used number of received pulses. The dotted lines represent the real velocity set in the simulation, which is $7000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. From these two figures, one can see that the curves fluctuate severely and the measurement errors are quite large for both SSC and MSC targets when only several pulses are received and used. With an increase in the number of received pulses, the fluctuations of the curves quickly reduce, and the estimated values gradually approach the true value. The variance of the curve in Fig. 11 is larger than that in Fig. 10 because the velocity measurement accuracy for the SSC target is better than that for the MSC target. If we just take into account the parts of curves with a pulse number $>64$, we get the mean values of $6981\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ in Fig. 10 and $7147\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ in Fig. 11, respectively, i.e., the relative measurement errors are $\sim 0.27$ and 2.1%, respectively.

## Fig. 10

The real velocity measurement for SSC target.

## Fig. 11

The real velocity measurement for MSC target.

## Conclusion

High-resolution radar imaging of a high-speed target using SCS must properly solve the problem of motion estimation and compensation. We propose the SCS-MCCM to accurately measure the velocity of a high-speed target with MSCs, which is the first step for high-resolution imaging. Detailed analysis and derivations of the proposed method are presented. Although it is not the most accurate method of velocity measurement with SCS, the advantages of SCS-MCCM make it practically applicable and easy to adopt. It can measure the target radial velocity by using only the echoes of subpulses within a burst, so it does not require a high repetition frequency of burst and, thus, has a large unambiguous measurement range. If the target moves straight, this method can also be used to measure the real velocity of the target. To show the validation and effectiveness of the method, simulations for an MSC target are conducted, which clearly show that the radar image is very well focused after performing motion compensation using the proposed method, but is not focused by using CCIFM; at the same time, the real target velocity is also well estimated. We should emphasize that there is still a room for improving the velocity measurement accuracy for an MSC target with the proposed method, and this is our next effort.

## References

1. D. R. Wehner, High-Resolution Radar, 2nd ed., Chapter 5, Artech House, Boston/London (1995). Google Scholar

2. N. Levanon, “Stepped-frequency pulse-train radar signal,” IEE Proc. Radar, Sonar Navig. 149(6), 297–309 (2002).IRSNE21350-2395 http://dx.doi.org/10.1049/ip-rsn:20020432 Google Scholar

3. Q. ZhangY.-Q. Jin, “Aspects of radar imaging using frequency-stepped chirp signals,” EURASIP J. Appl. Signal Process. 2006, 1–8 (2006).1110-8657 http://dx.doi.org/10.1155/ASP/2006/85823 Google Scholar

4. N. LevanonE. Mozeson, Radar Signals, Chapter 7, John Wiley & Sons, Hoboken, New Jersey (2003). Google Scholar

5. J. S. Son, Range-Doppler Radar Imaging and Motion Compensation, Artech House, Boston/London (2000). Google Scholar

6. T. LongL. Ren, “HPRF pulse Doppler stepped frequency radar,” Sci. China 52(5), 883–893 (2009).1009-2757 http://dx.doi.org/10.1007/s11432-009-0087-8 Google Scholar

7. B. D. RiglingG. S. Goley, “On the feasibility of ISAR processing of stepped frequency waveforms,” in 2013 IEEE Radar Conf., pp. 1–6, IEEE & AES, Ottawa, Canada (2013). Google Scholar

8. T. J. AbatzoglouG. O. Gheen, “Range radical velocity and acceleration MLE using radar LFM pulse train,” IEEE Trans. Aerosp. Electron Syst. 34(4), 1–4 (1998).IEARAX0018-9251 http://dx.doi.org/10.1109/7.722676 Google Scholar

9. Y. Liuet al., “Motion compensation of moving targets for high range resolution stepped-frequency radar,” Sensors 8(5), 3429–3437 (2008).SNSRES0746-9462 http://dx.doi.org/10.3390/s8053429 Google Scholar

10. P. FrancescoT. Enrico, “Motion compensation for a frequency stepped radar,” in Proc. of Waveform Diversity and Design Conf., pp. 255–259, IEEE, Pisa, Italy (2007). Google Scholar

11. B. LiuW. Chang, “Range alignment and motion compensation for missile-borne frequency stepped chirp radar,” Prog. Electromagn. Res. 136, 523–542 (2013).PELREX1043-626X http://dx.doi.org/10.2528/PIER12110809 Google Scholar

12. F. BerizziM. Martorella, “A Capria. A contrast-based algorithm for synthetic range-profile motion compensation,” IEEE Trans. Geosci. Remote Sens. 46(10), 3053–3061 (2008).IGRSD20196-2892 http://dx.doi.org/10.1109/TGRS.2008.2002576 Google Scholar

13. S. HuixiaL. Zheng, “Compound velocity measurement based on modulated frequency stepped radar,” Syst. Eng. Electron. 33(3), 539–543 (2011). (In Chinese).1004-4132 Google Scholar

14. H. YuanS. WenZ. Cheng, “Simultaneous radar imaging and velocity measuring with step frequency waveforms,” J. Syst. Eng. Electron. 20(4), 741–747 (2009).1004-4132 Google Scholar

15. W. FeiL. Teng, “A new method of velocity estimation for inverse V-shape stepped frequency signal,” in CIE'06 Int. Conf. on Radar, pp. 1–3, CIE & IEEE, Shanghai, China (2006). Google Scholar

16. G. Liet al., “Range and velocity estimation of moving targets using multiple stepped-frequency pulse trains,” Sensors 8(2), 1343–1350 (2008).SNSRES0746-9462 http://dx.doi.org/10.3390/s8021343 Google Scholar

17. B. Yun-xiaet al., “Velocity measurement and compensation method based on range profile cross-correlation in stepped-frequency radar,” Syst. Eng. Electron. 30(11), 2112–2115 (2008). (In Chinese).1004-4132 Google Scholar

18. L. Zhaoet al., “A novel method of ISAR imaging for high speed targets using stepped-frequency chirp waveform,” in 2011 IEEE CIE Int. Conf. on Radar, Vol. 2, pp. 1604–1607, IEEE & CIE, Chengdu, China (2011). Google Scholar

19. G. F. XiaH. Y. SuP. K. Huang, “Velocity compensation methods for a LPRF modulated frequency stepped frequency (MFSF) radar,” J. Syst. Eng. Electron. 21(5), 746–751 (2010). http://dx.doi.org/10.3969/j.issn.1004-4132.2010.05.005 Google Scholar

20. Z. YunhuaW. JieL. Haibin, “Two simple and efficient approaches for compressing stepped chirp signals,” in 2005 Asia-Pacific Microwave Conf., Vol. 2, pp. 690–693, IEEE, Suzhou, China (2005). Google Scholar

## Biography

Wenshuai Zhai is an assistant professor with the Key Laboratory of Microwave Remote Sensing, Center for Space Science and Applied Research, Chinese Academy of Sciences. She received her BS degree in electronic engineering from Zhejiang University in 2004, and her MS and PhD degrees in electromagnetic field and microwave technology in 2007 and 2013, respectively. Her current research interests include SAR/ISAR imaging and high-resolution radar signal processing.

Yunhua Zhang is a professor with the Key Laboratory of Microwave Remote Sensing, Center for Space Science and Applied Research, Chinese Academy of Sciences. He received his BS degree from Xidian University, Xi’an, and his MS and PhD degrees from Zhejiang University, Hangzhou, China, in 1989, 1993, and 1995, respectively, all in electrical engineering. His research interests cover the system design and signal processing of microwave sensors and high-resolution radar, as well as computational electromagnetics.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Wenshuai Zhai, Wenshuai Zhai, Yunhua Zhang, Yunhua Zhang, } "Velocity measurement and radar imaging of high-speed target by using stepped-chirp signal," Journal of Applied Remote Sensing 8(1), 083606 (27 June 2014). https://doi.org/10.1117/1.JRS.8.083606 . Submission:
JOURNAL ARTICLE
13 PAGES

SHARE
KEYWORDS