Design and Processing Method for Doppler-Tolerant Stepped-Frequency Waveform Using Staggered PRF

Stepped-frequency waveform may be used to synthesize a wideband signal with several narrow-band pulses and achieve a high-resolution range profile without increasing the instantaneous bandwidth. Nevertheless, the conventional stepped-frequency waveform is Doppler sensitive, which greatly limits its application to moving targets. For this reason, this paper proposes a waveform design method using a staggered pulse repetition frequency to improve the Doppler tolerance effectively. First, a generalized echo model of the stepped-frequency waveform is constructed in order to analyze the Doppler sensitivity. Then, waveform design is carried out in the stepped-frequency waveform by using a staggered pulse repetition frequency so as to eliminate the high-order phase component that is caused by the target’s velocity. Further, the waveform design method is extended to the sparse stepped-frequency waveform, and we also propose corresponding methods for high-resolution range profile synthesis and motion compensation. Finally, experiments with electromagnetic data verify the high Doppler tolerance of the proposed waveform.


Introduction
The wideband imaging ability is one of the most important functions of modern radar that may provide high-resolution information about targets. Wideband signals may be roughly divided into two categories, i.e., instantaneous wideband signals and the synthetic ones [1]. The instantaneous wideband signal achieves a wide bandwidth within a single pulse, which has the advantages of a short observation duration and a simple signal processing method. However, it has high requirements for hardware, which raises the cost of the radar system or even makes it unrealizable. The synthetic wideband signal, such as the stepped-frequency waveform (SFW), achieves a wide bandwidth by a sequence of pulses that is called a burst. Each burst has different carrier frequencies, and a wide bandwidth may be achieved by the synthetic processing method. In general, a synthetic wideband signal can be realized more easily for engineering practice due to its instantaneous narrow bandwidth.
The SFW has been used in radar systems [2,3], and different intra-and inter-sub-pulse modulation methods are applied. For instance, intra-sub-pulse modulation includes a single-frequency signal, a linear frequency modulation signal, a phase, or a frequencycoded signal [4,5], and inter sub-pulse modulation includes a linear frequency-stepping signal and a random or coded frequency-hopping signal [6,7]. In recent years, many studies on the SFW have been presented to improve radar imaging quality and anti-jamming ability, such as the compensation methods of system error [8,9], grating lobe suppression [10,11], and velocity compensation [12,13]. However, the duration of the SFW burst is too long, which leads to a reduction in detection efficiency and Doppler sensitivity, limiting its application to various scenes [14]. Therefore, to reduce the total duration, the sparse SFW (SSFW) has been designed for sparse target scenes [15]. Normally, the SSFW may reduce the Moreover, the echo of the SPRF-SFW is analyzed. Suppose the initial distance between target and radar is R 0 . The radial velocity of the target is defined as v, which is positive when the target is away from the radar. The velocity of the electromagnetic wave is c. Therefore, the two-way path delay under the initial distance is τ 0 = 2R 0 /c and under the ith sub-pulse is τ i = τ 0 + (2v/c)T i . Therefore, the base-band echo signal of the SPRF-SFW is expressed as: The model above is actually the generalized echo model of the SFW, which is composed of intra-pulse and inter-pulse terms. After constructing the model, the processing methods of the synthetic HRRP and the reasons of the Doppler sensitivity for CPRF waveforms are analyzed based on it.

Analysis of Doppler Sensitivity for CPRF Waveforms
The common signal processing method of the synthetic range profile is called twice pulse compression: one is intra-pulse compression, and the other is inter-pulse compression. The intra-pulse compression is to perform matched filtering on each sub-pulse, and the result of it is related to the waveform structure of a sub-pulse that can be a single-frequency, chirp and phase, or frequency-coded signal. Without loss of generality, we deduce the signal processing method with chirp, and the waveform structure of sub-pulse is expressed as: where rect t T p = 1 0 , , −T p /2 ≤ t ≤ T p /2 others is the rectangular window function and γ = B p /T p is the slope of sub-pulse chirp. Thus, according to Equation (3), the base-band echo is expressed as: Further, the obtained signal after intra-pulse compression is expressed as the structure of amplitude and phase, namely: where sinc(x) = (sin x)/x and the inter-pulse phase term is ϕ = −2π f i τ i . We find that the velocity of the target makes the peak of the intra-pulse compression shift from τ i to τ i + (2v f i )/(cγ) and makes the peak value from 1 to 1 − (2v f i )/ cB p . The results of intra-pulse compression are different in different waveform structures. Nevertheless, if the sub-pulse is a single-frequency signal, the bandwidth will be the reciprocal of the time width, which is not favorable for increasing the bandwidth and time width at the same time. Similarly, if the sub-pulse is a phase-coded signal, the result of pulse compression will be so sensitive to velocity that motion compensation would be difficult. Thus, the chirp signal is usually used in sub-pulse.
When the spectrum is complete, the phase term of the CPRF-SFW is: The phase term, which includes the stationary phase caused by the radial distance of static target, is ϕ r0 = −2π( f 0 + i∆ f )τ 0 , including only the primary phase, but the motion phase caused by the velocity of the moving target is ϕ v = −2π( f 0 + i∆ f )(2v/c)iT r , including the primary phase and the quadratic phase. As a result, when inverse fast Fourier transform (IFFT) is used during inter-pulse compression, imaging defocusing occurs due to the existence of the quadratic phase term. In addition, the position of the target is r p = cτ 0 /2 + v f 0 T r /∆ f and the broadening is r w = 2vMT r which is obtained by the use of the stationary phase method.
Further, when the spectrum is incomplete, the phase term of the CPRF-SSFW is: The stationary phase is Thus, when the CS algorithm is adopted during inter-pulse compression, the two phase terms above cannot be expressed by a unified basis, which leads to the defocusing of the synthetic HRRP.
In conclusion, the reason for Doppler sensitivity for the conventional CPRF waveform is the high-order phase component and filter mismatch caused by a moving target.

Waveform Design Method
After analyzing the Doppler sensitivity of CPRF waveforms, if the carrier frequencies of sub-pulses are determined, the Doppler tolerance may be improved by staggering the PRF of the sub-pulse.
By analyzing Equation (8), the design of the SPRF-SFW with high Doppler tolerance should meet the following two conditions. Condition 1: The PRF of the sub-pulse is modulated so that the high-order phase is eliminated in ϕ v . Condition 2: The pulse repetition interval (PRI), which is the reciprocal of the PRF, should have physical meaning that should be greater than zero. Thus, the condition 1 is expressed as: where k 0 is a constant, which represents the weighting coefficient. Thus, Equation (10) is solved as: At this time, the modulation mode of T i has met condition 1. Further, we give it physical meaning and connect it with the constant T r of the CPRF-SFW. Thus, we make k 0 = G 0 T r , and the transmitting time of each sub-pulse T i is rewritten as: where G 0 is a waveform parameter, which determines the value range of the PRI. As for condition 2, it is expressed as: where T rmin is the minimum value of the PRI. We substitute the transmission time of each sub-pulse T i into Equation (13), and condition 2 is rewritten as: The SPRF-SFW under a negative frequency modulation slope is also designed by this method. Then, they are simplified to the following expression: In addition, further analysis shows that when the waveform parameter G 0 < 0, the PRI of each sub-pulse is T ri = which is a monotonic increasing function. When the waveform parameter G 0 > 0, the PRI of each sub-pulse is which is a monotonic decreasing function, as shown in Figure 2.   Table 1.

Parameters Symbols Value
Initial carrier frequency   Table 1.
Further, when the spectrum is incomplete, the design of the SPRF-SSFW with high Doppler tolerance should also meet two conditions. Condition 2 remains unchanged, but Condition 1 changes to the following: the PRF of the sub-pulse is modulated so that the motion phase ϕ v and the stationary phase ϕ r0 are expressed by a unified basis. Then, the transmitting time of each sub-pulse T i is expressed as: where c(i) is the coding law of the transmission frequency. When the waveform parameter G 0 > 0, c(i) > c(i − 1), which is equal to the spectrum selection of positive frequency modulation coding. When the waveform parameter G 0 < 0, c(i) < c(i − 1), which is equal to the spectrum selection of negative frequency modulation coding. The timefrequency profiles of waveforms designed by the proposed method using the SPRF are shown in Figure 3.  In conclusion, the waveform design method of SPRF waveforms with high Doppler tolerance is proposed. At the same time, as one of the staggered PRF waveforms, the SPRF-SSFW also has the characteristics of improving blind speed: the minimum blind velocity is increased from the constant PRF to the least-common-multiple PRF of each sub-pulse.
In addition, Doppler tolerance is defined as the corresponding value when the pulse compression gain drops to the minimum allowable level, which is expressed as: where peak A is the peak amplitude after pulse compression and DT G is the minimum allowable gain. Thus, the imaging results for different velocities and the curve of pulse compression gains versus velocities are shown in Figure 4.   Table 1.
Furthermore, by analyzing the echo model of the CPRF-SFW and SPRF-SFW, we get the ambiguity function of the synthetic wideband signal by twice pulse compression as: The contour maps of the ambiguity function for the two kinds of SFW, which have the same waveform parameters, after wideband synthesis are shown in Figure 5. Figure  5a is the ambiguity function of the CPRF-SFW, and Figure 5b is the ambiguity function of the SPRF-SFW. Both of them have range velocity coupling. However, when the velocity deviates, the broadening of the ambiguity function of the CPRF-SFW is obvious, which is the reason for the defocusing synthetic HRRP. In contrast, the ambiguity function of the In conclusion, the waveform design method of SPRF waveforms with high Doppler tolerance is proposed. At the same time, as one of the staggered PRF waveforms, the SPRF-SSFW also has the characteristics of improving blind speed: the minimum blind velocity is increased from the constant PRF to the least-common-multiple PRF of each sub-pulse.
In addition, Doppler tolerance is defined as the corresponding value when the pulse compression gain drops to the minimum allowable level, which is expressed as: where A peak is the peak amplitude after pulse compression and G DT is the minimum allowable gain. Thus, the imaging results for different velocities and the curve of pulse compression gains versus velocities are shown in Figure 4.  In conclusion, the waveform design method of SPRF waveforms with high Doppler tolerance is proposed. At the same time, as one of the staggered PRF waveforms, the SPRF-SSFW also has the characteristics of improving blind speed: the minimum blind velocity is increased from the constant PRF to the least-common-multiple PRF of each sub-pulse.
In addition, Doppler tolerance is defined as the corresponding value when the pulse compression gain drops to the minimum allowable level, which is expressed as: where peak A is the peak amplitude after pulse compression and DT G is the minimum allowable gain. Thus, the imaging results for different velocities and the curve of pulse compression gains versus velocities are shown in Figure 4.   Table 1.
Furthermore, by analyzing the echo model of the CPRF-SFW and SPRF-SFW, we get the ambiguity function of the synthetic wideband signal by twice pulse compression as: The contour maps of the ambiguity function for the two kinds of SFW, which have the same waveform parameters, after wideband synthesis are shown in Figure 5. Figure  5a is the ambiguity function of the CPRF-SFW, and Figure 5b is the ambiguity function of the SPRF-SFW. Both of them have range velocity coupling. However, when the velocity deviates, the broadening of the ambiguity function of the CPRF-SFW is obvious, which is the reason for the defocusing synthetic HRRP. In contrast, the ambiguity function of the  Table 1.
Furthermore, by analyzing the echo model of the CPRF-SFW and SPRF-SFW, we get the ambiguity function of the synthetic wideband signal by twice pulse compression as: The contour maps of the ambiguity function for the two kinds of SFW, which have the same waveform parameters, after wideband synthesis are shown in Figure 5. Figure 5a is the ambiguity function of the CPRF-SFW, and Figure 5b is the ambiguity function of the SPRF-SFW. Both of them have range velocity coupling. However, when the velocity deviates, the broadening of the ambiguity function of the CPRF-SFW is obvious, which is the reason for the defocusing synthetic HRRP. In contrast, the ambiguity function of the proposed SPRF-SFW has no broadening, which means high Doppler tolerance.

Synthetic HRRP and Motion Compensation Process
The synthetic HRRP for the proposed SPRF-SSFW with high Doppler tole also processed by twice pulse compression, as shown in Figure 6. Because the SPR one inter-pulse modulation way, the intra-pulse compression is still the same f waveforms.

Synthetic HRRP and Motion Compensation Process
The synthetic HRRP for the proposed SPRF-SSFW with high Doppler tolerance is also processed by twice pulse compression, as shown in Figure 6. Because the SPRF is only one inter-pulse modulation way, the intra-pulse compression is still the same for CPRF waveforms.

Synthetic HRRP and Motion Compensation Process
The synthetic HRRP for the proposed SPRF-SSFW with high Doppler tolerance is also processed by twice pulse compression, as shown in Figure 6. Because the SPRF is only one inter-pulse modulation way, the intra-pulse compression is still the same for CPRF waveforms. First, when the spectrum is complete, the phase term of the proposed SPRF-SFW is expressed as: First, when the spectrum is complete, the phase term of the proposed SPRF-SFW is expressed as: The phase term includes the stationary phase ϕ r0 = −2π( f 0 + i∆ f )τ 0 and the motion phase ϕ v = −2π(i∆ f )(2v/c)G 0 T r . Thus, none of them contains a quadratic or high-order phase term. The inter-pulse compression by IFFT is expressed as: (20) In this case, the range profiles without defocusing are obtained, and the target position is Then, when the spectrum is incomplete, the phase term of the proposed SPRF-SSFW is expressed as: The phase term includes the stationary phase Thus, the two phase-terms above are expressed by a unified basis.
The sparse characteristics of the target are used in processing of the SPRF-SSFW, and the CS algorithm is proposed to transform the target imaging into a parameter estimation optimization problem. The measurement process is expressed as: where the matrix is expressed in bold, Y = [y 1 , y 2 , · · · , y M c ] T is the measurement value vector, and y i is the echo of the i-th sub-pulse. X = x 1 , x 2 , · · · , x M d T represents the target to be reconstructed, and x j is the amplitude of the i-th resolution unit. A is the sensing matrix with M c × M d , Φ is the measurement matrix with M c × M c , and Ψ denotes the basis matrix with M c × M d . The target is sparse, which means that most of the elements in matrix X are zero or minimum. Thus, the solution of the measurement process is transformed into the l 0 norm optimization problem as follows: The l 0 norm optimization problem is a nondeterministic polynomial problem. Thus, in the process of optimization, it is usually relaxed to an l 1 norm optimization problem. At the same time, considering the influence of weak scattering points and noise, the final optimization problem is expressed as follows: where ε represents one minimum value, the sensing matrix A = ΦΨ, and the measurement matrix Φ is usually set as a Gaussian random matrix, which expresses the influence of thermal noise in a single channel of the radar receiver. According to Equation (21), the elements of the basis matrix designed for the SPRF-SSFW, Ψ, are as follows: where r x = c/(2Ω) is the synthetic range resolution and Ω is the synthetic bandwidth.
As for the optimization method, greedy algorithms, such as orthogonal matching pursuit, are used to solve this problem. CS methods have made great progress, and to improve the reconstruction efficiency, we use the joint-block sparse method for reconstruction [37,38]. The flowchart of the joint-block sparse CS algorithm is shown in Figure 7.
= Ω is the synthetic range resolution and Ω is the synthetic bandwidth.
As for the optimization method, greedy algorithms, such as orthogonal matching pursuit, are used to solve this problem. CS methods have made great progress, and to improve the reconstruction efficiency, we use the joint-block sparse method for reconstruction [37,38]. The flowchart of the joint-block sparse CS algorithm is shown in Figure  7. According to the analysis, the proposed SPRF waveform has high Doppler tolerance and has a focusing range profile of the moving target. However, the migration of the range profile still appears, which is not conducive to the accumulation of multi-frame imaging. Therefore, it is necessary to carry on the correlation velocity measurement and motion compensation.
The echo motion compensation of the SPRF waveform is similar to that of the CPRF [21,22], which usually consists of two steps: intra-pulse range alignment and inter-pulse phase compensation, as shown in Figure 6. It is necessary to go on the intra-pulse range alignment when the echo envelope movement cannot be ignored in a single frame. It According to the analysis, the proposed SPRF waveform has high Doppler tolerance and has a focusing range profile of the moving target. However, the migration of the range profile still appears, which is not conducive to the accumulation of multi-frame imaging. Therefore, it is necessary to carry on the correlation velocity measurement and motion compensation.
The echo motion compensation of the SPRF waveform is similar to that of the CPRF [21,22], which usually consists of two steps: intra-pulse range alignment and interpulse phase compensation, as shown in Figure 6. It is necessary to go on the intra-pulse range alignment when the echo envelope movement cannot be ignored in a single frame. It effectively accumulates the intra-pulse compression results in the burst, reduces the loss of amplitude, and obtains the synthetic HRRP effectively. The envelope alignment step should ensure that the envelope movement caused by the residual velocity error is less than half of the sub-pulse resolution, which is expressed as |v e T M | < c/4B p , namely: The requirement of intra-pulse range alignment is not high, which is called coarse motion compensation. Generally, intra-pulse range alignment is attained by system measurement or intra-pulse correlation velocity measurement [21]. After intra-pulse envelope alignment, the target is located in the same range resolution cell, and inter-pulse compression is performed.
In inter-pulse phase compensation, we obtain the motion parameters by the difference in the range profile. In the radar system, the difference in the synthetic range profile is obtained by the adjacent correlation method, and the velocity of the target is further obtained as v = ∆R p /T M . The velocity measurement error should meet the requirement that the range difference is not greater than the resolution of the synthetic range profile c/2Ω, and then the velocity measurement error is expressed as: The requirement of inter-pulse phase compensation is high, which is called fine-motion compensation [22]. Furthermore, the difference in the proposed SPRF waveform is more accurate than the CPRF waveform due to the synthetic range profile without defocusing. After processing, the estimated motion parameters of the target are obtained.
Furthermore, it is necessary to consider the influence of the target's acceleration on motion compensation [23]. If the acceleration of the target is defined as a and the two-way path delay is τ i = τ 0 + (2v/c)T i + (a/c)T i 2 , the base-band signal with uniform acceleration is presented as: Further, the motion compensation term is expressed as: where v p represents the estimated velocity, a p represents the estimated acceleration, and τ p denotes the peak position after intra-pulse compression, used as the reference time delay. Thus, the SPRF-SSFW base-band echo after accurate motion compensation is expressed as: where v e = v p − v 0 represents the residual velocity, a e = a p − a 0 represents the residual acceleration after motion compensation, and τ e = τ p − τ i is the two-way delay of the difference between the reference point and each scattering point. The result shows that to ensure that the defocusing caused by acceleration is ignored, it is necessary to meet the following requirements: In addition, τ e is small and the phase term related to τ e is approximately constant, so the results of the SPRF-SSFW base-band echo after motion compensation are consistent with the results for a stationary target.

Experiment and Discussion
In this section, three simulation experiments combined with the electromagnetic data of the UAV are presented to verify the effectiveness of the proposed waveform design method. First, the proposed SPRF-SFW is compared with the conventional CPRF-SFW in order to verify the high-Doppler tolerance of the staggered PRF waveform. Second, the proposed SPRF-SSFW is designed and processed. On the one hand, it is compared with the CPRF-SSFW to verify the high Doppler tolerance. On the other hand, it verifies the effectiveness of the processing algorithm based on CS. Third, the velocities of target and waveform parameters are changed to verify the effectiveness of the waveform design method for high Doppler tolerance under different conditions.
In addition, the ideal HRRP and ISAR images obtained from the electromagnetic data of the UAV at a certain angle are shown in Figure 8 and are used as benchmarks to compare with the imaging results of the designed waveform.
In addition, e τ is small and the phase term related to e τ is approximately constant, so the results of the SPRF-SSFW base-band echo after motion compensation are consistent with the results for a stationary target.

Experiment and Discussion
In this section, three simulation experiments combined with the electromagnetic data of the UAV are presented to verify the effectiveness of the proposed waveform design method. First, the proposed SPRF-SFW is compared with the conventional CPRF-SFW in order to verify the high-Doppler tolerance of the staggered PRF waveform. Second, the proposed SPRF-SSFW is designed and processed. On the one hand, it is compared with the CPRF-SSFW to verify the high Doppler tolerance. On the other hand, it verifies the effectiveness of the processing algorithm based on CS. Third, the velocities of target and waveform parameters are changed to verify the effectiveness of the waveform design method for high Doppler tolerance under different conditions.
In addition, the ideal HRRP and ISAR images obtained from the electromagnetic data of the UAV at a certain angle are shown in Figure 8 and are used as benchmarks to compare with the imaging results of the designed waveform. Further, the evaluation factors are set as the signal-to-noise ratio (SNR) and the main lobe width W ML for the HRRP image and the image entropy H IE and the image contrast H IC for the ISAR image. The SNR is defined as the ratio of the peak power of the pulsecompression-result-to-noise variance, which is expressed as: where A pot is the peak amplitude, A st is the amplitude of the starting point, A f h is the amplitude of the ending point, N p is the number of the protecting points around the target, and var represents the variance function. Thus, the higher the SNR is, the higher the imaging quality is. The main lobe width was defined as the corresponding broadening at half of the peak, which is expressed as: where A(x i ) is the amplitude of the range profile. Thus, the narrower the main lobe width is, the higher the imaging quality is.
The image entropy is expressed as: where P(x i ) is the amplitude of the imaging point. Thus, the lower the image entropy is, the higher the imaging quality is. In addition, according to [45], the image contrast is expressed as: Thus, the higher the image contrast is, the higher the imaging quality is.

Experiment I: SPRF-SFW
This subsection shows the HRRP and ISAR images to verify that the proposed SPRF-SFW has high Doppler tolerance. The design of the SPRF-SFW is carried out under the condition that the carrier frequencies obey the linear law with f i = f 0 + i∆ f , and the same parameters of the SPRF-SFW and CPRF-SFW are shown in the Table 1. In addition, the target is the UAV electromagnetic data.
We further obtain the range resolution as 0.075 m, and the azimuth resolution is 0.2 m. To compare the SPRF-SFW with the CPRF-SFW, their burst durations are equal, that is, T M = MT r . Then, G 0 ∆ f f 0 +M∆ f MT r = MT r is satisfied, and the waveform parameter G 0 is represented as: Therefore, the PRI of the CPRF-SFW is 200 µs, and the PRI of the SPRF-SFW is obtained according to Equation (15), among which the maximum repetition interval is 245.92 µs, the minimum repetition interval is 165.29 µs, and other PRIs are shown in Figure 2. The time-frequency curves of the CPRF-SFW and the proposed SPRF-SFW under the designed parameters are presented in Figure 9. The base-band signal of the SRPF-SFW with acceleration is discussed in the above waveform parameters. According to Equation (31), to ensure that the defocusing caused by acceleration is ignored, the residual acceleration after motion compensation should meet . It is easy to be satisfied for aerodynamic targets. The residual velocity after coarse motion compensation is set as 20 m/s, which is far away from the radar direction, in order to better reflect the advantage of the proposed waveform. The received signal is contaminated by additive white complex Gaussian noise. Then, the SNR after intra-pulse compression is set as 16 dB, and the theoretical SNR is 26 dB after inter-pulse compression with 100 pulses. Therefore, the synthetic HRRPs of the conventional CPRF-SFW and the proposed SPRF-SFW are shown in Figure 10, and the ISAR images are shown in Figure 11.  The base-band signal of the SRPF-SFW with acceleration is discussed in the above waveform parameters. According to Equation (31), to ensure that the defocusing caused by acceleration is ignored, the residual acceleration after motion compensation should meet |a e | ≤ 37.5 m/s 2 . It is easy to be satisfied for aerodynamic targets. The residual velocity after coarse motion compensation is set as 20 m/s, which is far away from the radar direction, in order to better reflect the advantage of the proposed waveform. The received signal is contaminated by additive white complex Gaussian noise. Then, the SNR after intra-pulse compression is set as 16 dB, and the theoretical SNR is 26 dB after interpulse compression with 100 pulses. Therefore, the synthetic HRRPs of the conventional CPRF-SFW and the proposed SPRF-SFW are shown in Figure 10, and the ISAR images are shown in Figure 11. ity after coarse motion compensation is set as 20 m/s, which is far away from the radar direction, in order to better reflect the advantage of the proposed waveform. The received signal is contaminated by additive white complex Gaussian noise. Then, the SNR after intra-pulse compression is set as 16 dB, and the theoretical SNR is 26 dB after inter-pulse compression with 100 pulses. Therefore, the synthetic HRRPs of the conventional CPRF-SFW and the proposed SPRF-SFW are shown in Figure 10, and the ISAR images are shown in Figure 11.  As shown in Figure 10, the blue solid line is the synthesis HRRP of the CP and the red solid line is the synthesis HRRP of the SPRF-SFW. The imaging off SPRF-SFW caused by the target motion is 2 m, which is consistent with the th analysis. Additionally, compared with Figure 8a, it is observed that the CPRF-SF ing has both offset and defocusing caused by the target motion, with the SNR as and the main lobe width as 0.71 m, while the SPRF-SFW imaging only has offset broadening, with the SNR as 25.93 dB and the main lobe width as 0.08 m. For com the SNR for the HRRP of the SPRF-SFW is higher and the main lobe width is n than that of the CPRF-SFW, which means that the imaging quality of the SPR higher than that of the CPRF-SFW for the same moving target.
The ISAR images of the CPRF-SFW and the SPRF-SFW are shown in Figu respectively. Compared with Figure 8b, the target motion causes defocusing o SFW imaging, with the image entropy as 2484.08 bits and the image contrast as 2. the proposed SPRF-SFW imaging is focused, with the image entropy as 1811.15 the image contrast as 3.13. For comparison, the image entropy for ISAR imag SPRF-SFW is lower and the image contrast is higher than that of the CPRF-SFW means that the imaging quality of the SPRF-SFW is higher than that of the CPRF the same moving target.
Consequently, the above simulation results verify the high Doppler toleran proposed SPRF-SFW, which makes up for the Doppler sensitivity of the CPRF-S has more practical value in complex scenes, especially in an imaging scene with velocity targets such as target separation due to the difficulty in accurate motion sation.

Experiment II: SPRF-SSFW
In this subsection, the following simulation experiments are carried out w spectrum is incomplete to verify that the SPRF-SSFW has high Doppler toleranc tains the wideband synthetic image by the CS algorithm. Therefore, the degr transmission spectrum integrity is set to 0.6. As a result, when compared with si experiment I, the number of sub-pulses is reduced from 100 to 60, the burst du As shown in Figure 10, the blue solid line is the synthesis HRRP of the CPRF-SFW and the red solid line is the synthesis HRRP of the SPRF-SFW. The imaging offset of the SPRF-SFW caused by the target motion is 2 m, which is consistent with the theoretical analysis. Additionally, compared with Figure 8a, it is observed that the CPRF-SFW imaging has both offset and defocusing caused by the target motion, with the SNR as 16.58 dB and the main lobe width as 0.71 m, while the SPRF-SFW imaging only has offset without broadening, with the SNR as 25.93 dB and the main lobe width as 0.08 m. For comparison, the SNR for the HRRP of the SPRF-SFW is higher and the main lobe width is narrower than that of the CPRF-SFW, which means that the imaging quality of the SPRF-SFW is higher than that of the CPRF-SFW for the same moving target.
The ISAR images of the CPRF-SFW and the SPRF-SFW are shown in Figure 11a,b, respectively. Compared with Figure 8b, the target motion causes defocusing of CPRF-SFW imaging, with the image entropy as 2484.08 bits and the image contrast as 2.32, while the proposed SPRF-SFW imaging is focused, with the image entropy as 1811.15 bits and the image contrast as 3.13. For comparison, the image entropy for ISAR images of the SPRF-SFW is lower and the image contrast is higher than that of the CPRF-SFW, which means that the imaging quality of the SPRF-SFW is higher than that of the CPRF-SFW for the same moving target.
Consequently, the above simulation results verify the high Doppler tolerance of the proposed SPRF-SFW, which makes up for the Doppler sensitivity of the CPRF-SFW, and has more practical value in complex scenes, especially in an imaging scene with multiple velocity targets such as target separation due to the difficulty in accurate motion compensation.

Experiment II: SPRF-SSFW
In this subsection, the following simulation experiments are carried out when the spectrum is incomplete to verify that the SPRF-SSFW has high Doppler tolerance and obtains the wideband synthetic image by the CS algorithm. Therefore, the degree of the transmission spectrum integrity is set to 0.6. As a result, when compared with simulation experiment I, the number of sub-pulses is reduced from 100 to 60, the burst duration is reduced from 20 ms to 12 ms, and the carrier frequency of each sub-pulse is selected by random coding. Further, the PRI of the sub-pulse is calculated according to Equation (16), and other parameters are the same as those given in Table 1. The time-frequency curves of the CPRF-SSFW [15] and the proposed SPRF-SSFW under the designed parameters are presented in Figure 12. Due to the incomplete spectrum, the burst duration of the SRPF-SSFW is lower than that of the SRPF-SFW in experiment I. Thus, the acceleration compensation tolerance of the SPRF-SSFW is higher. According to Equation (31), the residual acceleration after motion compensation should meet   Due to the incomplete spectrum, the burst duration of the SRPF-SSFW is lower than that of the SRPF-SFW in experiment I. Thus, the acceleration compensation tolerance of the SPRF-SSFW is higher. According to Equation (31), the residual acceleration after motion compensation should meet |a e | ≤ 104.17 m/s 2 . The residual velocity after coarse motion compensation is also set as 20 m/s. Meanwhile, the synthetic HRRP and ISAR images of the CPRF-SSFW and the SPRF-SSFW are shown in Figures 13 and 14, respectively. Due to the incomplete spectrum, the burst duration of the SRPF-SSFW is lower than that of the SRPF-SFW in experiment I. Thus, the acceleration compensation tolerance of the SPRF-SSFW is higher. According to Equation (31), the residual acceleration after motion compensation should meet  As shown in Figure 13, the blue solid line is the HRRP of the CPRF-SSFW and the red solid line is the HRRP of the SPRF-SSFW. The imaging offset caused by the target motion is 1.2 m, which is consistent with the theoretical analysis. Moreover, comparing these results with Figure 8a, we observe that the CPRF-SSFW imaging has both offset and defocusing caused by the target motion, with the SNR as 15.98 dB and the main lobe width as 0.39 m, while the SPRF-SSFW imaging only has offset without broadening, with the SNR as 23.95 dB and the main lobe width as 0.08 m. For comparison, the SNR for the HRRP of the SPRF-SSFW is higher and the main lobe width is narrower than that of the CPRF-SSFW, which means that the imaging quality of the SPRF-SSFW is higher than that of the CPRF-SSFW for the same moving target.  The ISAR images of the CPRF-SSFW and the SPRF-SSFW are shown in Figure 14a,b, respectively. Compared with Figure 8b, the target motion causes serious defocusing of CPRF-SSFW images, which is caused by the mismatch of the sensing basis, with the image entropy as 2712.49 bits and the image contrast as 2.05, while the SPRF-SSFW imaging has better focusing, and the sensing basis matches well, with the image entropy as 1947.65 bits and the image contrast as 2.96. For comparison, the image entropy for ISAR images of the SPRF-SSFW is lower and the image contrast is higher than that of the CPRF-SSFW, which also means that the imaging quality of the SPRF-SSFW is higher than that of the CPRF-SSFW for the same moving target.
Consequently, the above simulation results verify that the proposed SPRF-SSFW ha high Doppler tolerance and is processed by the CS algorithm well. Clearly, the SPRF-SSFW has a shorter burst duration and it is more suitable for a complex electromagnetic environment.

Experiment III: Different Conditions
In this subsection, the velocities of target and waveform parameters are changed, and the evaluation factors for ISAR images of conventional CPRF waveforms [15] (blue star line), optimized CPRF waveforms [11,17] (green cross line), and proposed SPRF waveforms (red diamond line) are further compared under different conditions in order to verify the effectiveness of the waveform design method for high Doppler tolerance.
First, when the spectrum is complete, the waveform parameters are fixed, as shown in experiment I, and the residual velocity of target changes. The imaging qualities of the conventional CPRF-SFW, optimized CPRF-SFW [11], and proposed SPRF-SFW are compared in Figure 15. The image entropy of the conventional CPRF-SFW increases and the image contrast decreases with velocity, which corroborates the decrease in imaging quality. The optimized CPRF-SFW increases the imaging quality in a certain velocity tolerance but decreases significantly when the velocity is too high. In addition, all evaluation factors of the proposed SPRF-SFW change little, which verifies the high Doppler tolerance for different velocities.
Then, when the degree of spectrum integrity is set to 0.6, the waveform parameters are fixed, as shown in experiment II, and the residual velocity of target also changes. The imaging qualities of the conventional CPRF-SSFW [15], optimized CPRF-SSFW [17], and proposed SPRF-SFW are compared in Figure 16. The image entropy of the conventional CPRF-SSFW increases and the image contrast decreases with velocity, which corroborates the decrease in imaging quality. The change is more rapid than that in Figure 15, because the imaging quality decreases more obviously when the echo of the moving target does not match the CS basis. The optimized CPRF-SSFW is similar to the conventional one, and the imaging focuses only when the velocity is low. In addition, all evaluation factors of the proposed SPRF-SSFW change little, and the imaging quality is stable, which also verifies the high Doppler tolerance for different velocities. pared in Figure 15. The image entropy of the conventional CPRF-SFW increases image contrast decreases with velocity, which corroborates the decrease in imagin ity. The optimized CPRF-SFW increases the imaging quality in a certain velocity to but decreases significantly when the velocity is too high. In addition, all evaluation of the proposed SPRF-SFW change little, which verifies the high Doppler tolera different velocities. Then, when the degree of spectrum integrity is set to 0.6, the waveform para are fixed, as shown in experiment II, and the residual velocity of target also chang imaging qualities of the conventional CPRF-SSFW [15], optimized CPRF-SSFW [1 proposed SPRF-SFW are compared in Figure 16. The image entropy of the conve CPRF-SSFW increases and the image contrast decreases with velocity, which corro the decrease in imaging quality. The change is more rapid than that in Figure 15, b the imaging quality decreases more obviously when the echo of the moving targ not match the CS basis. The optimized CPRF-SSFW is similar to the conventional o the imaging focuses only when the velocity is low. In addition, all evaluation fa the proposed SPRF-SSFW change little, and the imaging quality is stable, which a ifies the high Doppler tolerance for different velocities. Next, when the spectrum is complete, the residual velocity of target is fixed as and the burst duration changes. The imaging qualities of the conventional CPR optimized CPRF-SFW [11], and proposed SPRF-SFW are compared in Figure 17. age entropy of the conventional CPRF-SFW increases and the image contrast de with the burst duration, which verifies that the longer the burst duration, the lo Doppler tolerance of the conventional CPRF-SFW. The optimized CPRF-SFW in the imaging quality in a certain burst duration but decreases significantly when th duration is too long. In addition, all evaluation factors of the proposed SPRF-SFW little, which verifies the high Doppler tolerance for different burst durations.  Next, when the spectrum is complete, the residual velocity of target is fixed as 20 m/s and the burst duration changes. The imaging qualities of the conventional CPRF-SFW, optimized CPRF-SFW [11], and proposed SPRF-SFW are compared in Figure 17. The image entropy of the conventional CPRF-SFW increases and the image contrast decreases with the burst duration, which verifies that the longer the burst duration, the lower the Doppler tolerance of the conventional CPRF-SFW. The optimized CPRF-SFW increases the imaging quality in a certain burst duration but decreases significantly when the burst duration is too long. In addition, all evaluation factors of the proposed SPRF-SFW change little, which verifies the high Doppler tolerance for different burst durations.
Finally, when the spectrum is incomplete, the residual velocity of the target is fixed as 20 m/s and the degree of spectrum integrity changes. The imaging qualities of the conventional CPRF-SSFW [15], optimized CPRF-SSFW [17], and proposed SPRF-SFW are compared in Figure 18. The image entropy of the conventional CPRF-SSFW is always higher and the image contrast is always lower, which indicates that the images of the moving targets defocus under various degrees of spectrum integrity. The optimized CPRF-SSFW has higher imaging quality than the conventional CPRF-SSFW but also has imaging defocusing. As a comparison, the image entropy of the proposed SPRF-SSFW increases slowly and the image contrast decreases slowly when the degree of spectrum integrity is bigger than 0.5. The main reason for the decrease in imaging quality is the decrease in the SNR, which is caused by the decrease in the accumulated pulse number. However, when the degree of spectrum integrity is smaller than 0.5, image entropy increases and image contrast decreases rapidly, which indicates that the degree of spectrum integrity should not be too low, otherwise the imaging quality will be affected. The experiment verifies the high Doppler tolerance of the proposed SPRF-SSFW for a high degree of spectrum integrity. age entropy of the conventional CPRF-SFW increases and the image contrast de with the burst duration, which verifies that the longer the burst duration, the lo Doppler tolerance of the conventional CPRF-SFW. The optimized CPRF-SFW in the imaging quality in a certain burst duration but decreases significantly when th duration is too long. In addition, all evaluation factors of the proposed SPRF-SFW little, which verifies the high Doppler tolerance for different burst durations.  Finally, when the spectrum is incomplete, the residual velocity of the target as 20 m/s and the degree of spectrum integrity changes. The imaging qualities of t ventional CPRF-SSFW [15], optimized CPRF-SSFW [17], and proposed SPRF-S compared in Figure 18. The image entropy of the conventional CPRF-SSFW is higher and the image contrast is always lower, which indicates that the image moving targets defocus under various degrees of spectrum integrity. The op CPRF-SSFW has higher imaging quality than the conventional CPRF-SSFW but a imaging defocusing. As a comparison, the image entropy of the proposed SPR increases slowly and the image contrast decreases slowly when the degree of sp integrity is bigger than 0.5. The main reason for the decrease in imaging qualit decrease in the SNR, which is caused by the decrease in the accumulated pulse n However, when the degree of spectrum integrity is smaller than 0.5, image entr creases and image contrast decreases rapidly, which indicates that the degree of sp integrity should not be too low, otherwise the imaging quality will be affected. The iment verifies the high Doppler tolerance of the proposed SPRF-SSFW for a high of spectrum integrity.

Conclusions
In this paper, we proposed a waveform design and a corresponding pro method using the SPRF to improve the Doppler tolerance of the SFW and the SSFW oretical analysis and experiments verified that both the SPRF-SFW and the SPR designed by the proposed method have high Doppler tolerance and a focusing sy range profile. In addition, the SPRF-SSFW processed by CS has not only high D tolerance but also a short burst duration and low probability of intercept, which im the anti-jamming ability. Further work can expand the sub-pulse optimization de improve the adaptability to maneuvering targets in a complex electromagnetic e ment.

Conclusions
In this paper, we proposed a waveform design and a corresponding processing method using the SPRF to improve the Doppler tolerance of the SFW and the SSFW. Theoretical analysis and experiments verified that both the SPRF-SFW and the SPRF-SSFW designed by the proposed method have high Doppler tolerance and a focusing synthetic range profile. In addition, the SPRF-SSFW processed by CS has not only high Doppler tolerance but also a short burst duration and low probability of intercept, which improves the antijamming ability. Further work can expand the sub-pulse optimization design to improve the adaptability to maneuvering targets in a complex electromagnetic environment.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.