A General Range-Velocity Processing Scheme for Discontinuous Spectrum FMCW Signal in HFSWR Applications

Discontinuous spectrum signal which has separate subbands distributed over a wide spectrum band is a solution to synthesize a wideband waveform in a highly congested spectrum environment. In this paper, we present a general range-velocity processing scheme for the discontinuous spectrum-frequency modulated continuous wave (DS-FMCW) signal specifically. In range domain, we propose a simple time rearrangement operation which converts the range transform problem of the DS-FMCW signal to a general spectral estimation problem of nonuniformly sampled data. Conventional periodogram results in a dirty range spectrum with high sidelobes which cannot be suppressed by traditional spectral weighting. In this paper, we introduce the iterative adaptive approach (IAA) in the estimation of the range spectrum. IAA is shown to have the ability to provide a clean range spectrum. On the other hand, the discontinuity of the signal spectrum has little impact on the velocity processing. However, with the range resolution improved, the influence of the target motion becomes nonnegligible. We present a velocity compensation strategy which includes the intersweep compensation and in-sweep compensation. Our processing scheme with the velocity compensation is shown to provide an accurate and clean range-velocity image which benefits the following detection process.


Introduction
High frequency surface wave radar (HFSWR) refers to a classification of radar that operates in the HF band  and utilizes the surface wave mode of propagation.These systems can be used both in surveillance to detect and track vessels and in oceanography to measure surface current and ocean wave parameters [1][2][3][4][5][6].
However, the HF band is a heavily congested part of the radio spectrum which makes it difficult to find a continuous silent frequency band to transmit radar signals.This limits the signal bandwidth of the HF radar system and results in a poor range resolution [7,8].In [7], a continuous measurement of noise and interference data in the frequency band of 3-6 MHz at Cape Race, Newfoundland, Canada, in the period between August 1, 1998, and May 10, 2000, shows that channels with a bandwidth of 20 kHz are readily available and a bandwidth of 100 kHz is not available.
To overcome the trouble of limited band, Green and Kutuzov proposed a waveform design idea which utilizes several discontinuous but clean subbands to synthesize a signal with the desired bandwidth [9][10][11].The discontinuity of the signal spectrum leads to high sidelobe levels of the matched filter output which cannot be suppressed by conventional spectral weighting.Therefore, most of the efforts at present focused on applying the sidelobe suppression technique in the matched filtering output.The first category of the methods is the mismatched filter (or instrumental variable filter) approach.It is a data-independent method which gets the filter coefficients by optimizing the autocorrelation function (ACF) and applies the filter to the matched filter output.The second category uses the CLEAN algorithm to reconstruct the target signal [12].The third category is based on spectrum reconstruction; this includes using AR model [9,13] or regularization interpolation [12] to fill the "gap" parts of the spectrum.And, on the other hand, [14][15][16] considered 2 International Journal of Antennas and Propagation suppressing the sidelobes from a perspective of waveform designing rather than signal processing.
While some papers use the model of frequency hopped pulse signal [16][17][18], some other papers consider a more abstract discontinuous spectrum signal model and concentrate on the analysis of the autocorrelation function and the power spectrum density (PSD) [14,15]; only a few papers use a signal model of FMCW [12,13].As the FMCW signal [19] and its variation FMICW signal [20] are widely used in HF radar applications such as the SeaSonde [3] and WERA [5] systems and they are quite different from the pulsed signal both in receiver structure and following signal processing method, we focus on the processing of discontinuous spectrum FMCW (DS-FMCW) signal specifically.Similar to the conventional FMCW signal, DS-FMCW uses stretch processing instead of correlation processing used in pulsed radar in the range domain.And we call this process range transform compared to the commonly used term "pulse compression" in pulsed radar.
Different from the conventional FMCW radar which outputs a constant-frequency complex sinusoid (beat signal) after mixing with the transmitted reference signal, the DS-FMCW radar outputs a constant-frequency signal but with phase hopped at specific locations.In this paper, we propose a kind of transformation applied on the time variable, by which the hopped phase is converted to a kind of "hopped" sampling instants.This process may be implicit in former works [12,13] and we choose to call it time rearrangement operation here.After the time rearrangement operation, the sampled postmixing signal becomes a nonuniformly sampled sequence of a complex sinusoid.Therefore, we set up a general spectral analysis model for the range processing of the DS-FMCW signal and the crux of the matter now becomes the spectral estimation of a nonuniformly sampled data sequence.
The spectral analysis of nonuniformly sampled data has drawn much attention since the nonuniformity in the data is common in various applications [21].And a large number of methods have emerged recently.The classical Fourier transform based periodogram suffers from the heavy leakage problem which leads to a low resolution and high sidelobe levels.Therefore, it is not applicable in the processing of DS-FMCW signal.In this paper, we apply the iterative adaptive approach (IAA) algorithm in the range spectrum estimation of the DS-FMCW signal.IAA was first proposed in [22] for target direction of arrival (DOA) estimation.It can be interpreted as an iteratively weighted least-square periodogram which eliminates almost completely the leakage problems of the conventional periodogram method in a fully data-adaptive manner [22,23].It is quite suitable to be applied in the range spectrum estimation of DS-FMCW signal.However, other spectral analysis methods suitable for the line spectrum estimation can also be taken into consideration.We use IAA as an example in this paper.
We consider the range-velocity processing of the DS-FMCW signal afterwards.With the range resolution improved, some subtle issues may occur due to target motion.Accordingly, we propose a velocity compensation strategy in Power Threshold this paper.Finally, we present a general range-velocity processing scheme for the DS-FMCW signal.This paper is organized as follows.Section 2 formulates the signal model and the range processing model of the DS-FMCW signal.The time rearrangement operation is introduced in this section.Section 3 extends the model to range-velocity processing and gives a velocity compensation strategy.A general range-velocity processing scheme for the DS-FMCW signal is presented at the end of this section.Section 4 gives a design example of a DS-FMCW signal and shows the effectiveness of our proposed processing scheme by several experiments.The final section concludes this paper.

DS-FMCW Signal Model and Range Transform
2.1.Signal Model and Time Rearrangement Operation.The HF radar system using discontinuous spectrum signal should cooperate with a spectrum monitor which monitors the spectrum environment before the radar starts to transmit a signal.A typical HF spectrum distribution is illustrated in Figure 1.The shading areas below the threshold correspond to the silent frequency bands which can be used by the radar.There are  relative silent segments in the frequency span  0 to  0 +Δ.Let the start frequency and start time of each segment be   and   ,  = 0, 1, . . .,  − 1,  0 = 0.Then, the duration of each segment is   =  +1 −   ,  = 0, 1, . . .,  − 1, and the sweep period is   = ∑ −1 =0   .We use a linear frequency sweep in each segment and fix the sweep slope to  = /  , in which  is the total frequency bandwidth occupied by the radar signal.The ratio of the occupied frequency bandwidth of the radar signal to the whole frequency span is defined as the occupied frequency ratio (OFR) of the signal.Thus, OFR = /Δ.The resulting operation waveform is the so-called discontinuous spectrum-frequency modulated continuous wave (DS-FMCW) signal.The transmitted and received signal frequency time characteristics are shown in Figure 2.
The transmitted signal in a single sweep period is where is a rectangular window function.The received signal from a single target is an attenuated, delayed, and Doppler-shifted version of the transmitted signal.In this section, we only consider a stationary point target and leave the discussion of the Doppler/velocity processing to the next section.Then, the received signal is given by in which  represents the attenuation factor which is mainly related to the target radar cross section (RCS) and  is the round-trip time for the signal to propagate from the radar system and back.A multitarget case can be seen as a superposition of the different echoes at the receiver.Note that the effect of the target delay on the rectangular envelope is omitted since in HF radar applications the maximum potential target delay usually satisfies  max ≪   [19].Similar to the stretch processing of the conventional FMCW signal, the received signal is mixed with the transmitted signal at the receiver.This process is also called deramp processing or dechirp processing.Stretch processing is a technique to reduce the sampling-rate requirements of the following analog-to-digital converter (ADC) [24].The output of the mixer is in which the phase term of the th segment   () is Unlike the mixer output of the conventional FMCW signal, which is a constant-frequency complex sinusoid (beat signal), the corresponding output of the DS-FMCW signal is more analogous to a kind of phase-coded signal with  subpulses of each having a different initial phase related to the target delay.Thus, the range spectrum cannot be obtained by a direct FFT/IFFT operation.However, with a simple transform of (5), the unknown phase hopping quantity converts to a time domain shift related to the waveform parameters which are predefined.Therefore, a rearrangement of the time variable in each segment eliminates the phase hopping.The time vector after the rearrangement operation is For each segment, we have When the receiver samples the deramped signal at a uniform interval   and by letting   , 0 ≤  < ,  = ⌊  /  ⌋ be the sampling instants (⌊⌋ represents the notion of being rounded down), the corresponding sampling instants after the time rearrangement    turn out to be nonuniform (or, more accurately, piecewise uniform).This process of time rearrangement is illustrated in Figure 3. Figure 3(b) shows that the phase hopping is eliminated after the rearrangement.
Accordingly, the sampled deramped signal can be represented as in which  =  2 is a constant phase term and s =  exp () represents a complex amplitude.From ( 9), the sampled time series [] is a constantfrequency complex sinusoid signal sampled at time    .Thus, the range transform problem of the DS-FMCW signal converts to the problem of spectral estimation of nonuniformly sampled data.
We can denote the sampled signal in vector form as where (⋅)  represents the transpose operation,  to target delay , in which  is the considered frequency range, and a(⋅) :  → C ×1 is a known function which solely depends on the spectrum distribution of the DS-FMCW signal; the function output is the matched filter vector (similar to the steering vector in the spatial processing) in the range domain Consider the multitarget situation; the sampled data sequence of the deramped signal is where  ∈ C ×1 is a noise term and  is the unknown number of targets.Equation ( 12) is a common model of spectral analysis of nonuniformly (irregularly) sampled data.
To assist in explaining the spectral analysis algorithms, we extend the model of (12).Let the considered frequency range be  = { :  min ≤  ≤  max }, in which the maximum and the minimum beat frequencies are decided by the considered range window: We then use a uniform grid {  }  =1 to cover , with a specific fine step, and assume that { f } lie on (practically, close to) the grid.This means that there exist  1 , . . .,   such that f =    ,  = 1, . . ., .Also, let Using this notation, we can rewrite ( 12) as 2.2.Periodogram Range Spectrum.The most intuitive method is the classical Fourier transform based periodogram (FP), which is also known as single-frequency least-square (SFLS) method, matched filtering method, or beam-forming method (in array signal processing applications).It uses filter coefficients matched to the received signal structure, which yields the optimum output signal-to-noise ratio (SNR).The periodogram method can also be derived from a solution to the following least-squares (LS) data fitting problem: The solution to the above problem is Then, the spectrum estimation at the frequency grid {  }  =1 is given by The resolution of the periodogram method is proportional to the reciprocal of the observation time interval.
According to (7), in the DS-FMCW situation, the observation interval after a rearrangement of the sampling instants is where Δ is the total frequency span.Thus, the observation interval is stretched from   = / to Δ/ after the rearrangement of the sampling instants.Then, the range resolution is given by Therefore, the range resolution solely depends on the total frequency span (Δ) of the DS-FMCW signal.By adopting the waveform in an appropriate frequency span, we can obtain the desired range resolution.However, the development of the spectral analysis has led to many methods which can resolve targets whose range separation is within Δ.Therefore, we call the resolution of periodogram the inherent resolution.However, the spectrum discontinuity in the frequency span leads to high sidelobes when applying the conventional periodogram method which cannot be suppressed by conventional spectral weighting.The sidelobe levels are related to the waveform occupied frequency ratio (OFR) and the specific spectrum distribution in the frequency span.The larger the OFR and the more uniform the distribution, the lower the sidelobe levels.In this paper, we do not consider extreme cases such as the occupied frequency being at the two sides of the frequency span or being all concentrated into a large continuous block, in which the situation degrades to a continuous spectrum.We assume the occupied frequency is distributed uniformly in the frequency span which is more likely encountered in the practical situation.

IAA Range Spectrum.
The inherent limitations of the periodogram approach motivate a more effective spectral analysis method to be applied here.IAA is a nonparametric and user parameter free weighted least-squares based spectral analysis method first proposed for spatial spectral estimation [22].It is a periodogram related spectral analysis method which eliminates almost completely the leakage problems of the periodogram method in a fully data-adaptive manner.The range spectrum estimation problem is equivalent to a direction of arrival (DOA) estimation problem with single snapshot in the array processing applications.And IAA can properly handle this situation.Consequently, IAA is quite suitable to be applied in the situation of range spectrum estimation of the DS-FMCW signal.
Periodogram uses solution to the LS fitting problem (16) as the estimate in each frequency scan point.It is a dataindependent approach.While estimating a spectrum value of a specific frequency point   , the other components located at frequencies different from   perform as interferences.
Prewhitening of the other frequency components except the current processing one can cancel strong interferences in the data.Therefore, IAA uses a weighted least-squares (WLS) fitting criterion with the weight matrix equal to the inversion of the interferences and noise covariance matrix: in which R +, is the interferences (signals at range cells other than the range cell  of current interest) and noise covariance matrix and square of the weighted ℓ 2 -norm is defined as The iterative solution to the WLS problem ( 21) leads to the IAA algorithm.As the range transform problem of the DS-FMCW signal is equivalent to a general spectral analysis problem of nonuniformly sampled data, the standard IAA algorithm can be directly applied here [22,23].

Extension to Range-Velocity Processing
Conventional FMCW radar uses Fast Fourier Transform (FFT) on the range domain (fast time) and pulse domain (slow time), respectively, to form a range-Doppler spectrum (RD image), thus attaining the target range and Doppler information.In the last section, we have proposed that the range transform of the DS-FMCW signal is equivalent to a spectral estimation problem of a nonuniformly sampled data set.Then, we can replace the FFT adopted by FMCW in the range domain processing with a more effective spectral analysis method.However, the spectrum discontinuity of the waveform has little impact on the Doppler processing.We can still employ FFT in the Doppler processing which is essentially a matched filter bank.But since the carrier frequency of the signal changes at different time segments, a single target's velocity corresponds to different Doppler frequencies.Therefore, we use a matched filter bank in the velocity domain directly which is more intuitive than processing in the Doppler domain and then adjusting different Doppler cells into a single velocity cell according to their carrier frequencies.This is the reason why we use the term range-velocity processing instead of the commonly used term range-Doppler processing.
In this section, we present a scheme for the range-velocity processing of the DS-FMCW signal.We change the sequence of processing steps with applying velocity domain matched filtering first followed by range domain spectral analysis at each velocity channel.This modification has the following benefits: (1) The matched filter coefficients in the velocity domain for each range bin can be easily decided since their carrier frequency is predefined.
(2) After velocity processing, we can compensate the range bias caused by the Doppler shift in a single sweep period by multiplying a correction factor in each velocity channel.This bias is usually omitted in the RD processing of conventional FMCW HF radar.
International Journal of Antennas and Propagation (3) If we choose a spectral analysis method which does not reserve the phase information at the output, we have to do the velocity processing first.
Suppose there are  sweep periods in a coherent processing interval (CPI); the transmitted signal can be expressed as The received signal from a moving target can be expressed as The impact of the target delay on the signal envelope is omitted as before;  represents the attenuation factor.And the delay function is where  0 is the initial range of the target,  0 is the corresponding delay, and V is the radial velocity of the target towards the radar.After mixing the received signal with the transmitted signal, we get the deramped output Let  1 be the time vector in a single sweep period; then the time vector in the th sweep period is  =  1 +   .Thus, the first phase term is which is the same as the stationary model ( 6) in Section 2.
After a time rearrangement, we can get the range spectrum from this term using a spectral analysis method.
We next consider the third term of  , () ≜  (1)  3(,) () +  (2)  3(,) () +  (3)  3(,) () . (29) We have three contributions to the phase: a constant, a linear term in time , and a quadratic term in time,  2 .Firstly, consider the phase variation of the quadratic term in a CPI.Assume that the target's distance of movement in a CPI is less than 3 km; then, V/ < 10 −5 .And the sweep slope of a HF FMCW signal is usually at the order of (10 6 ); thus, the variation of the quadratic term  (3)  3(,) () in a CPI is at the order of 2 × 10 −4 , which can be omitted.Then, consider the linear term  (2)  3(,) () = −2 0 (2V/), which is proportional to the target velocity.As we will see in the following analysis,  2(,) (see (33)) is the so-called Doppler term and has the same form as  (2)  3(,) .Since  0 ≪  0 , the phase term  (3)   3(,) can also be omitted.In conclusion, where Ψ 3 is a constant.For the second term of  (,) (), in which the first phase term is (1,1) 2(,) represents the quadratic disturbance term caused by the target velocity during a sweep period.Its phase variation during a sweep period is at the order of (2 × 10 −4 ).( (1,2)  2(,) ) has the same form as the intersweep Doppler term (the term ( (2,2)  2(,) ()) in (33)), and (( −   )) is comparable to signal carrier frequency; thus, we should consider it in matched filtering in the velocity domain.
In conclusion, the terms influenced by target velocity are Note that the first term and the last term in (34) also exist in conventional FMCW radar, but, in HF applications, they are usually omitted as their impact is minor.However, the neglect of the third term will lead to both a velocity bias and a range bias and the neglect of the first term will lead to a range bias.Thus, we choose to compensate both terms in our system.
As the compensation of the third term is handled along the slow time domain and the compensation of the first term is handled along the fast time domain, we call them intersweep compensation and in-sweep compensation, respectively.Intersweep compensation can be achieved by merging the compensation coefficients into the velocity domain matched filter coefficients.Let the th sampling point in the fast time domain be  1 [] = [ +   ] −   ; then, we can get the velocity domain matched filter coefficients at the point  1 []: and   in a  (V) have to be chosen according to current fast time point.Let the considered velocity range be  = {V : V min ≤ V ≤ V max }, and, similar to the range domain, we use a uniform grid {V  }  =1 to cover .Since the th velocity channel after the velocity domain matched filtering represents the fast time echoes from all the targets with the specific velocity V  , we can compensate the range bias caused by the target velocity according to the first term of  V(,) at each velocity channel.The in-sweep compensation coefficient B in the th fast time sampling point of the th velocity channel is in which   and   are decided from the current fast time sampling point  1 [𝑛].
The block diagram of the DS-FMCW radar system is shown in Figure 4. We employ a stretch processing based receiver and signal processor.The signal processor part follows our proposed range-velocity processing scheme with the details shown in Algorithm 1.Consider a series of  sweeps; after the postmixing sampling, we have the sampled sequence [],  = 0, 1, . . ., −1.With each sweep placed as a row, we get × data matrix Y.The samples in each row of Y are successive samples of the returns in a single sweep, that is, successive range bins.Each column represents a series of measurements from the same range bin over successive sweep periods.We denote the th row of a matrix Y by Y ,: and the th column of a matrix Y by Y :, .We also denote the velocity matched filtering matrix at the range bin  by A  , which includes all the vectors at different velocity grids In Figure 4 and Algorithm 1, Z ∈ C × is the velocity domain matched filtering output.After the velocity compensation in each velocity channel, we get the compensated output Hadamard (element-wise) product and (⋅) * denotes the conjugate operator.Finally, the spectral analysis in each velocity channel gives the final range-velocity spectrum Ξ ∈ C × , in which  is the range bin number that coincides with the range domain model in Section 2.

Numerical Experiments
In the first experiment, we demonstrate the determination of the postmixing sampling rate by using the spectral window method [25].In the second experiment, we illustrate the range spectrum formed by the periodogram and the IAA algorithm.In the third experiment, we evaluate the performance of IAA when applying it in range domain spectral analysis as the OFR of the DS-FMCW signal varies by several Monte Carlo simulations.The second and third experiments only consider the range domain processing while the fourth experiment gives a complete range-velocity processing result in accordance with the processing scheme shown in Figure 4 and Algorithm 1.
In Experiments 1, 2, and 4, we use the same DS-FMCW signal designed according to Canada 1999 (summer) HF band spectrum monitoring result [7]; the available frequency band distribution is shown in Table 1.
The designed DS-FMCW signal utilizes all the clear bands between the start frequency  0 = 5.11 MHz and the end

Experiment 1: Determination of Postmixing Sampling
Rate.In a stretch processing receiver, the sampling process comes after the mixer, by which the sampling rate can be much less than the signal bandwidth.In a conventional FMCW radar, the minimum postmixing sampling rate without any aliasing problem is proportional to the considered range window extent: where  win =  max −  min .This is defined through the Nyquist-Shannon sampling theorem essentially.However, for the DS-FMCW signal, the sampling instants after the time rearrangement operation become nonuniform (Figure 6), and there is no well accepted definition of the nonaliasing sampling rate for the nonuniformly sampled data sequence.
In this paper, we adopt the Nyquist or rollover frequency definition for nonuniform samples described in [25], which is called the spectral window method.Given  nonuniform samples with sampling instants {  }  =1 , the spectral window at any frequency  is defined as We can easily verify that (0) = 1, and the frequency range without any aliasing problem can be determined by examining the largest range for  in which the only peak with height equal (or close) to 1 is at  = 0.
For our designed DS-FMCW signal, at a postmixing sampling rate of   = 4 kHz, the sampling instants after the time rearrangement operation (   ) are shown in Figure 6. Figure 6 shows that the postrearrangement sampling instants    are nonuniformly (more precisely, piecewise uniformly) distributed and total observation interval is stretched from the sweep period   = 128 ms to an equivalent virtual interval of Δ/ = 402 ms.And, on the other hand, comparing with Figure 5, we can see that two figures have an identical distribution.This illustrates the role of the time rearrangement which is converting the discontinuity of the signal spectrum to the discontinuity of the sampling instants.
In fact, we can derive that the periodogram spectrum is just a superposition of spectral windows with different shift amount related to the target range.Thus, the spectral window and the periodogram spectrum have the same spectral shape.We use a range scan step  that is five times smaller than the traditional range resolution according to (20): The spectral window of the designed DS-FMCW signal is shown in Figure 7, and the spectral window of a full occupied FMCW signal in the same frequency span ( 0 ∼  0 + Δ) is plotted in the same figure for comparison.In Figure 7, the spectral window of the FMCW signal is aliased in the range window 0∼200 km while the DS-FMCW signal is not.And the maximum range window for the FMCW signal is 111.3 km which is consistent with (38).Therefore, the spectral window can be used to ensure that the range spectrum of the DS-FMCW signal is unambiguous under the selected sampling rate.
Figure 7 also indicates that, compared with the conventional FMCW signal within the same frequency span, the range spectrum of the DS-FMCW signal has the same main lobe width (/(2Δ)) but has much higher sidelobe levels.Thus, other methods rather than periodogram should be considered to be applied in the range spectrum estimation of the DS-FMCW signal.

Experiment 2:
Periodogram Range Spectrum and IAA Range Spectrum.In this experiment, we illustrate the improvement of the sidelobe levels and range resolution provided by IAA when applied in the range spectrum analysis.We consider two targets located at the range of 120 km and 120.3 km whose separation is slightly larger than the inherent range resolution and with the input SNR of 0 dB and −5 dB, respectively.The range profiles generated by periodogram and IAA are shown in Figure 8.It can be inferred from Figure 8 that the periodogram method suffers from heavy leakage problems which make the peak of weaker target almost masked by the high sidelobes of a strong target nearby.
And the high sidelobe level also causes false alarms and contaminates the detection output.IAA, on the other hand, resolves two closely separated targets clearly and has a much lower sidelobe level.The sidelobe level of the IAA spectrum can hardly be determined from Figure 8 as it is even lower than the noise power.In fact, the spectrum output of IAA is quite close to a line spectrum.There are sharp peaks at the locations of the targets and near-zero values at the other locations.

Experiment 3: Performance of the IAA Range Spectrum.
This subsection designs two Monte Carlo experiments to investigate the performance of IAA when applied in the range processing of the DS-FMCW signal.In these two experiments, we use a similar waveform parameter of the DS-FMCW signal.We fix the total frequency span Δ to 600 kHz.As stated at the end of Section 2, we do not consider the extreme situations and when we consider a specific OFR, we generate randomly a uniformly distributed discontinuous spectrum with total free bandwidth  = Δ ⋅ OFR at each Monte Carlo trial.Firstly, we examine the peak sidelobe level (PSL) of the IAA range spectrum of the DS-FMCW signal as the OFR varies.PSL is defined as the ratio of the maximum sidelobe level to the main lobe level.In this experiment, we fix the SNR to 10 dB. Figure 9 shows the PSL of periodogram and IAA range spectrum via 500 Monte Carlo trials when the OFR varies from 20% to 100%.We can infer from the figure that   as the PSL of IAA spectrum can be lower than −20 dB with an OFR larger than 40%, the PSL of periodogram is larger than −13 dB at any OFR and cannot be suppressed via spectral weighting when the OFR is less than 1.
The second Monte Carlo simulation examines the capability of resolving closely spaced targets by IAA.Two targets are placed at ranges 120 km and 120.2 km, whose separation is slightly less than the inherent range resolution.By definition, two targets are resolved in a given run if both , where R and   denote the estimated and true range for the th target, respectively.Figure 10 shows that, at a SNR higher than 0 dB and an OFR higher than 40%, IAA can properly resolve the targets.

Experiment 4: A Complete Range-Velocity Processing
Example of DS-FMCW Signal.In this subsection, we demonstrate the whole range-velocity process shown in Figure 4 and Algorithm 1.The velocity domain matched filtering is done first followed by correcting in each velocity channel the phase offset caused by the target velocity in a single sweep period (in-sweep compensation), and the range spectrum analysis is then done in each velocity channel to output the final rangevelocity spectrum.
We consider five targets whose information is shown in Table 2.The first to the fourth targets are closely spaced both in range and in velocity domain, and the fifth target    Figure 12 is the range-velocity spectrum formed by IAA.It shows that IAA significantly outperforms periodogram both in main lobe width and in sidelobe levels which result in a higher resolution and a much better detection performance.
Similar to Figure 11, we also lay out the processing result with and without the velocity compensation.It shows that with velocity domain processing coming first, followed by range domain processing in each velocity channel.For the range domain processing, we propose a time rearrangement operation which converts the range transform problem to a general spectral estimation problem of a nonuniformly sampled data sequence.We have presented the periodogram range spectrum and the IAA range spectrum.IAA is shown to outperform the periodogram method both in resolution and in sidelobe levels.And IAA has a quite ideal spectrum output with sharp peaks at the target positions and nearzero values at other positions.The Monte Carlo simulation shows that, with a waveform OFR above 40% and a signal SNR above 0 dB, the range spectrum of the DS-FMCW signal generated by IAA can achieve quite a low sidelobe level and an acceptable resolving ability.For the velocity processing, we propose a strategy of velocity compensation which includes the intersweep and in-sweep compensation.Finally, a complete design example showed the indispensability of the velocity compensation and validated the effectiveness of the proposed range-velocity processing scheme.

Figure 5 :
Figure 5: Time-frequency characteristic of the transmit signal.

Figure 6 :
Figure 6: Sampling instants in a single sweep period after the time rearrangement operation.

Figure 7 :
Figure 7: Spectral window of the designed DS-FMCW signal and a continuous FMCW signal within the same frequency span.

Figure 8 :
Figure 8: Range spectrum generated by periodogram and IAA.

Figure 10 :
Figure 10: IAA resolution probability versus SNR for closely located targets at different OFR.

Figure 11 is
the range-velocity spectrum generated by periodogram, where the circles represent the true positions of the targets in the range-velocity space.The heavy leakage problem in the range domain is clearly illustrated in the figure.We also show the indispensability of the intersweep and in-sweep compensation as mentioned in Section 3. Figure 11(a) is the result with neither compensation applied.The data cursors mark two of the spectrum peaks of the moving targets.The peaks of the four moving targets have several bins shifted from the truth in both range and velocity domain while peak position of the stationary one is accurate.After applying both the intersweep and in-sweep compensation, as shown in Figure 11(b), the peak positions of all the five targets turn out to be accurate.