Rotor UAV’s Micro-Doppler Signal Detection and Parameter Estimation Based on FRFT-FSST

The micro-Doppler signal generated by the rotors of an Unmanned Aerial Vehicle (UAV) contains the structural features and motion information of the target, which can be used for detection and classification of the target, however, the standard STFT has the problems such as the lower time-frequency resolution and larger error in rotor parameter estimation, an FRFT (Fractional Fourier Transform)-FSST (STFT based synchrosqueezing)-based method for micro-Doppler signal detection and parameter estimation is proposed in this paper. Firstly, the FRFT is used in the proposed method to eliminate the influence of the velocity and acceleration of the target on the time-frequency features of the echo signal from the rotors. Secondly, the higher time-frequency resolution of FSST is used to extract the time-frequency features of micro-Doppler signals. Moreover, the specific solution methodologies for the selection of window length in STFT and the estimation of rotor parameters are given in the proposed method. Finally, the effectiveness and accuracy of the proposed method for target detection and rotor parameter estimation are verified through simulation and measured data.


Introduction
In recent years, the incidents such as UAV crashing, injuring people, or damaging properties, etc. that civil UAVs endanger public safety have frequently occurred, for which the implementation of safety supervision and prevention has been a "new normal" in [1,2].
Since most of the civil UAVs have rotors, the micro-Doppler signals thus produced contain information such as the structural features and motion status of the target, which can be used to improve the accuracy of the classification and recognition of the target in [3,4]. For the analysis and recognition of the micro-Doppler in [5][6][7], the micro-Doppler features caused by the UAV rotors were studied, where Victor C. Chen pointed out that the micro-Doppler from the radial velocity was a very important feature for the target with multiple rotors while the micro-Doppler spectra caused by multiple rotors might affect each other, for which the high-precision time-frequency analysis techniques should be further studied. In [8], the micro-Doppler effect caused by small UAVs was studied, and the effectiveness of the micro-Doppler signal for target detection was verified by using simulation and measurement data. In [9], the STFT was used to analyze the parameters of the micro-Doppler signal and the window function was designed according to the instantaneous period of the data set, improving the time-frequency resolution of the signal. In [10], the FMCW radar was used to detect Nano-UAV (less than 5 cm), and the micro-Doppler features were extracted from experimental data. In [11], an IRNN-based algorithm for UAV detection and classification recognition was studied based on the micro-Doppler features with the recognition rate reaching more than 90%. In [12], a time-frequency analysis method for micro-Doppler signal based on STFRFT (Short-Time Fractional Fourier Transform) was studied, significantly improving the time-frequency resolution compared with STFT. In [13][14][15][16], the adaptive detection method of targets was studied in a clutter environment, which could be used to improve the detection performance of small UAV targets. In [17,18], a target detection method based on acoustic array signals was studied, and the experimental results indicated that it could be used as an effective means for low-altitude UAV detection. In [19][20][21], the automatic classification and detection methods of UAVs based on the micro-Doppler spectrum were studied, and the research indicated that the introduction of the neural network could improve the target recognition rate. In [22][23][24], a variety of methods for UAV detection and classification were studied, and the research indicated that the target recognition rate could be improved by introducing multimodal and multi-dimensional deep learning methods. In [25][26][27], the time-frequency focus performance of the methods such as SST (Synchrosqueezing Transform), SST2 (Secondorder Synchrosqueezing Transform), HSST (High-order Synchrosqueezing Transform), etc. to the time-varying signals was studied, which could be used to improve the time-frequency resolution of the micro-Doppler signal. In [28], it was pointed out that the micro-Doppler effect was an important feature for identifying small UAVs, and it was further proved that SST had a higher time-frequency resolution with measured data.
The above analysis indicates that the micro-Doppler signal plays an important role in the detectionand identification of UAVs, but the existing methods often fail to consider the influence of velocity and acceleration on the micro-Doppler signal. Meanwhile, when the number of rotors increases, the time-frequency resolution by the STFT method is too low to accurately judge the rotor parameters. To this end, starting from the influence of velocity, acceleration, and multi-rotor parameters on the time-frequency features of the micro-Doppler signal, the FRFT-based method for velocity and acceleration estimation has been studied in this paper, together with the study on the method for extraction of the multi-component micro-Doppler signal features based on FSST, combined with the feature of higher time-frequency resolution of SST, and the high-precision estimation of rotor parameters has been realized.
An echo model of the rotor UAV is given in the Section 2 of this paper; a proposed algorithm is given in the Section 3; a method for the parameter estimation of the proposed algorithm is indicated in the Section 4; the influence of the rotor parameters on the focus performance of micro-Doppler signal is analyzed in the Section 5; the effectiveness of the algorithm is verified by using simulation and measured data in the Section 6.

Signal Model
The geometric configuration relationship between the UAV and the radar is shown in Figure 1, where the radar is located at the origin of the coordinate system (X,Y), the reference coordinate system (x,y) is the translation of the radar coordinate system, whose origin O is the geometric center of the UAV. In order to simplify the analysis, here, it is assumed that the azimuth angle and the elevation angle of the UAV relative to the radar are both 0. Figure 1. Geometric configuration between radar and UAV rotors.
Suppose 0 P is a scattering point on the rotor, which rotates around the center p o,with the rotation frequency being r f . The distance from 0 P to point o is l , and the  Suppose P 0 is a scattering point on the rotor, which rotates around the center point o, with the rotation frequency being f r . The distance from P 0 to point o is l, and the distance from radar to point o is R 0 , where R 0 is a function of time t. Here, let t = 0, the initial rotation angle of P 0 is θ 0 , then point P 0 rotates to point P t when going through time t. Since it is assumed that the radar and the rotor are on the same plane (x, y, z = 0), the distance from point P t to the radar can be expressed as [4,29]: where, V R = v 0 + at is the radial velocity of the UAV, v 0 is the initial velocity, and a is the acceleration, usually (l/R 0 ) 2 → 0 in the far field. Therefore, the echo signal of the scattering point P t received by the radar can be simplified as: where, δ m is the reflection coefficient of the rotor and the phase function is Φ P (t) = 4πR P (t)/λ. Performing down-conversion processing on (2) and ignoring the influence of the constant term, the baseband echo signal of scattering point P t can be obtained as: Suppose the length of the blade of the rotor is L, when it is integrated according to (3), the baseband echo signal of the entire blade is: According to (4), the total received signal of Z blades can be expressed as: where, θ k = θ 0 + 2kπ/Z(k = 0, 1, 2 · · · Z − 1) the initial rotation angle of Z different blades, φ k (t) represents its phase function, expressed as: When β = 0, Formula (5) can be turned as: where, φ k (t) = 4π λ V R t + 4π λ L 2 cos β sin(2π f r t + θ k ). When considering the echo signal from the body, the echo signal of the UAV can be obtained as: where, S m (t) = s Σ (t) is the echo signal of the rotor, and S D (t) is the echo signal of the body, which can be expressed as:

Doppler Features of UAV
Performing derivative operation for the phase function of S D (t) in (8), and the Doppler frequency caused by the body can be obtained as: Performing derivative operation for the phase function of S m (t) in (8), and the Doppler frequency caused by the kth rotor can be obtained as: Formula (11) shows that the micro-Doppler frequency of the signal is modulated by the rotation rate f r through the functions sin 2π f r t and cos 2π f r t.

Velocity and Acceleration Compensation Based on FRFT
It can be seen from (10) and (11) that velocity and acceleration will cause the translation and tilt of the time-frequency curve time-frequency, even folding in severe cases, making it impossible to estimate the relevant parameters. The velocity and acceleration would cause the phase of the target echo to turn into the approximate form of a linear FM signal, and the FRFT has advantages in LFM parameter detection. Therefore, the FRFT has been used to estimate the velocity and acceleration, and further construct the penalty function so as to eliminate its influence on the time-frequency features of the echo.
Example 1: Suppose that the radar carrier frequency is 5 GHz, the baseband sampling rate 10 KHz, the sampling time 1 s, the number of target blades 2, the rotation frequency 2 Hz/s, the length 6 m, the elevation angle of the blades relative to the radar β = 0 0 , not considering the body echo in the analysis, the signal-to-noise ratio of the blade echo is −3 dB. Figure 2a,b indicate the influence of velocity and acceleration on the time-frequency curve; Figure 2c shows the time-frequency curve when processed by FRFT compensation (at velocity of 19.6 m/s, acceleration of 60.4 m/s 2 ). It can be seen from Figure 2 that when the velocity and acceleration increase, the time-frequency curve will be folded, which will lead to an increase of errors in the subsequent parameter estimation. When the compensation process is carried out by using (17), the influence of the velocity and acceleration on the accuracy of parameter estimation can be eliminated, as shown in Figure 2c.

Time-Frequency Analysis Based on FSST
STFT had been widely used in the analysis of micro-Doppler signals [34,35], but the selection of window length was very important for the analysis of echo signals at different

Time-Frequency Analysis Based on FSST
STFT had been widely used in the analysis of micro-Doppler signals [34,35], but the selection of window length was very important for the analysis of echo signals at different rotation frequencies; therefore, a selection method is presented in this paper by the experiential knowledge to realize the adaptive selection of window length; at the same time, in order to further improve the aggregation performance of the time-frequency curve, the SST method is introduced to improve the time-frequency resolution of the signal.
STFT is performed to (17) to obtain: where, The time-frequency resolution performance in (18) is directly related to the selection of window length. Here, according to the requirement that the micro-Doppler bandwidth of the signal to be analyzed and the length of the analyzed signal shall be less than a quarter period, the selected window length g W is as: where, F is the signal sampling frequency, and f mBW is the micro-Doppler spectrum bandwidth caused by the tip of the blade, it can be expressed from the (11) as follows: In order to further improve the time-frequency concentration performance of STFT, SST is used in this section for subsequent processing.
Performing SST on Formula (18), the time-frequency diagram after rearrangement can be expressed as [25,36]: where,ŵ f (τ, w) is the instantaneous frequency, which can be expressed aŝ It can be seen by comparing Equations (21) and (18) that in Equation (21), the timefrequency focus processing has been performed, which is conducive to the time-frequency resolution and parameter estimation of multi-component micro-Doppler signals.

Parameter Estimation
In order to simplify the analysis, assuming that the velocity and acceleration have been compensated, the form of the numerical dispersion of (17) is rewritten as: where, n = 0, 1, 2 · · · N is the number of signal sampling points, and T s = 1/F is the interval of signal sampling time.

Estimation of the Number of Rotors
The multiple time points τ Q on the time-frequency image T Y [τ, w] in (21) are randomly extracted to estimateẑ, the number of UAV rotors, namely: where, q i represents the sum of the number of points with the amplitude greater than the threshold in the extracted ith point, which can be expressed as: where η ∈ [0.6, 0.9] is the threshold coefficient. The higher coefficient is generally suitable for high SNR signals. The lower coefficient is conducive to the detection of low SNR signals, but it will also lead to the increase of false targets.

Estimation of Rotation Frequency f r
When performing FFT processing on (17), the frequency interval ∆ f of adjacent spectral lines is obtained, and the rotation frequency f r is obtained according to ∆ f : where, n a , n a+1 is the number of frequency points corresponding to the adjacent spectral line of the maximum peak point in (17). In addition, by performing FFT on T Y [τ, w] by the time dimension in (21), the rotation frequency f r can also be obtained, namely: where, N is the point number of signal analysis, n a and n a+1 are respectively the frequency points corresponding to the main lobe and the first side lobe in FFT{T Y [τ, w]} τ spectrum.

Estimation of Rotor LengthL
According to (20), the expression ofL can be obtained as [37]: where, the numerical expression off mBW iŝ where, a 0+ , a 0− represents the number of frequency points corresponding to the maximum value on the left and right sides of the origin in FFT{S L (nT s )} spectrum. Figure 3 shows the parameter estimation process of the proposed algorithm.The steps of the proposed algorithm are as follows: (1) FRFT is performed on the echo signal to estimate the velocity and acceleration of the target, and then a penalty function is constructed to eliminate the influence of the velocity and acceleration on the time-frequency curve to obtain the function S L (t); (2) FFT is performed on S L (t) to obtain ∆ f ,the interval of adjacent spectral lines,f mBW , the micro-Doppler bandwidth, and then g W ,the estimated window length is given; (3) g W is used to perform STFT and SST processing on the echo signal to obtain a highresolution time-spectrogram;

Influence of Rotor Number and Rotation Frequency on Time-Frequency Features
Suppose that the simulation parameters are the same as those in Example 1. Figure 4 shows the results of the proposed algorithm when the number of rotors is 2 and 4 respectively, of which Table 1 shows the configuration of each parameter in Figure 4. Figure 4a,b indicate that the time-frequency resolution will not decrease as the number of rotors increases when the rotation frequency is 2 and the number of rotors is changed. Figure 4a,c show that the time-frequency resolution will be significantly decreased when the number of rotors is 2, the rotation frequency is changed, and the window length remains the same. Figure 4c,d show that the time-frequency resolution will be significantly improved when the rotation frequency is doubled compared to Figure 4a and the window length is reduced. The above results indicate that the time-frequency resolution of the proposed algorithm is mainly related to the window length of the window function determined by the rotation frequency. The larger the rotation frequency, the shorter the window length, and vice versa.

Influence of Rotor Number and Rotation Frequency on Time-Frequency Features
Suppose that the simulation parameters are the same as those in Example 1. Figure 4 shows the results of the proposed algorithm when the number of rotors is 2 and 4 respectively, of which Table 1 shows the configuration of each parameter in Figure 4. Figure 4a,b indicate that the time-frequency resolution will not decrease as the number of rotors increases when the rotation frequency is 2 and the number of rotors is changed. Figure 4a,c show that the time-frequency resolution will be significantly decreased when the number of rotors is 2, the rotation frequency is changed, and the window length remains the same. Figure 4c,d show that the time-frequency resolution will be significantly improved when the rotation frequency is doubled compared to Figure 4a and the window length is reduced. The above results indicate that the time-frequency resolution of the proposed algorithm is mainly related to the window length of the window function determined by the rotation frequency. The larger the rotation frequency, the shorter the window length, and vice versa. Table 1. Configuration of each parameter in Figure 4.

Figure 4ẑ
f rfmBW g W (Theoretical) g W (Practical)    Figure 5 shows the processing results of the proposed algorithm when r f Δ is at 1 Hz/s and 2 Hz/s, respectively, of which Table 2 shows the parameter configuration used in each sub-graph in Figure 5. Comparing Figure 5a,b, it can be seen that as r f Δ increases, the time-frequency focus performance of the proposed algorithm becomes worse.

Influence of Time-Varying Rotation Frequency on Time-Frequency Features
In order to solve the above problems, Equation (29) can be utilized to estimate BW m fˆ for the signal segmentation. Figure 6 gives out the FFT processing results of different segmented signals when r f Δ is at 1 Hz/s and 2 Hz/s respectively; meanwhile, Table 3 gives out BW m fˆ corresponding to each sub-graph in Figure 6, from which it can be seen that the is, and the larger the variation range of W g is. Therefore, for the echo signals with variable rotation frequencies, no fixed window length can be used for time-frequency matching. Figure 7 shows the processing results of the proposed algorithm when r f Δ is at 2 Hz/s, of which Items (c) and (d) in Table 3 are the parameter configurations used in Figure 7.
Comparing Figures 5b and 7, it can be seen that the time-frequency focus performance of the proposed method is obviously improved when the window length W g is selected by segmentation; the above results show that for the echo signals with variable   Table 1 in Section 5.1 show that whenẑ and f r are fixed, it can obtain better time-frequency resolution to use g W as the window length of the FSST; but in the actual UAV echo signal processing, f r may vary as ∆ f r does, for which the adaptability of the proposed algorithm has been analyzed under different ∆ f r . Figure 5 shows the processing results of the proposed algorithm when ∆ f r is at 1 Hz/s and 2 Hz/s, respectively, of which Table 2 shows the parameter configuration used in each sub-graph in Figure 5. Comparing Figure 5a,b, it can be seen that as ∆ f r increases, the time-frequency focus performance of the proposed algorithm becomes worse. rotation frequencies, it can improve the time-frequency resolution of the signal when the method for signal segmentation is used to select the window length. In addition, further in-depth research shall be conducted on how to reasonably segment the signal. Table 2. Configuration of each parameter in Figure 5.    Table 3. Calculation of each parameter in Figure 6.  Table 2. Configuration of each parameter in Figure 5. In order to solve the above problems, Equation (29) can be utilized to estimatef m BW for the signal segmentation. Figure 6 gives out the FFT processing results of different segmented signals when ∆ f r is at 1 Hz/s and 2 Hz/s respectively; meanwhile, Table 3 gives outf m BW corresponding to each sub-graph in Figure 6, from which it can be seen that the larger ∆ f r is, the wider the broadening range off m BW is, and the larger the variation range of g W is. Therefore, for the echo signals with variable rotation frequencies, no fixed window length can be used for time-frequency matching. Figure 7 shows the processing results of the proposed algorithm when ∆ f r is at 2 Hz/s, of which Items (c) and (d) in Table 3 are the parameter configurations used in Figure 7.    Table 3. Calculation of each parameter in Figure 6.

Analysis on Algorithm Performance
In this section, the simulation data and measured data are used for verification in order to further validate the effectiveness and adaptability of the proposed algorithm, of which the simulation data are the small UAV echo signals collected by a W-band radar, which are used in the experiment to verify, on the one hand, the higher time-frequency resolution of the proposed algorithm than that of STFT, and, on the other hand, the correctness of the proposed algorithm to the parameter estimation such as the parameter compensation, window length selection W g , micro-Doppler bandwidth BW m fˆ, and rotation frequency; and the measured data are the echo data from a model helicopter collected by a P-band radar, which are used in the experiment to verify, on the one hand, the micro-Doppler features of the UAV under low-frequency radar detection, and on the other hand, it can be seen from Equation (20) that the micro-Doppler frequency bandwidth is proportional to the length of the rotor and the radar frequency, which means increasing the length of the rotor can obtain the micro-Doppler signal equivalent to that generated with the high-frequency radar detecting UAV, based on which to further verify the performance of the proposed algorithm in terms of target detection and parameter estimation. Comparing Figures 5b and 7, it can be seen that the time-frequency focus performance of the proposed method is obviously improved when the window length g W is selected by segmentation; the above results show that for the echo signals with variable rotation frequencies, it can improve the time-frequency resolution of the signal when the method for signal segmentation is used to select the window length. In addition, further in-depth research shall be conducted on how to reasonably segment the signal.

Analysis on Algorithm Performance
In this section, the simulation data and measured data are used for verification in order to further validate the effectiveness and adaptability of the proposed algorithm, of which the simulation data are the small UAV echo signals collected by a W-band radar, which are used in the experiment to verify, on the one hand, the higher time-frequency resolution of the proposed algorithm than that of STFT, and, on the other hand, the correctness of the proposed algorithm to the parameter estimation such as the parameter compensation, window length selection g W , micro-Doppler bandwidthf mBW , and rotation frequency; and the measured data are the echo data from a model helicopter collected by a P-band radar, which are used in the experiment to verify, on the one hand, the micro-Doppler features of the UAV under low-frequency radar detection, and on the other hand, it can be seen from Equation (20) that the micro-Doppler frequency bandwidth is proportional to the length of the rotor and the radar frequency, which means increasing the length of the rotor can obtain the micro-Doppler signal equivalent to that generated with the high-frequency radar detecting UAV, based on which to further verify the performance of the proposed algorithm in terms of target detection and parameter estimation.

Simulation Experiment
Assuming that the radar is provided with a CW Regime, the radar parameters and target parameters are shown in Table 4; the signal-to-noise ratio(SNR) of the UAV fuselage is 5 dB higher than that of the rotor echo, and the SNR ratio between the rotor echo and Gaussian noise is 5 dB. The results of the proposed algorithm processed according to the flow path in Figure 3 are shown in Figure 8, of which Figure 8a is the results from FFT, it can be obtained by Equation (29) thatf mBW is at about 20 KHz, quite different from 4 KHz, the theoretical value calculated by Equation (20), and the experimental results indicate that the existence of the velocity and acceleration of the target has a greater influence on the estimation of micro-Doppler bandwidth, especially the spectrum broadening caused by acceleration has an even greater influence on the estimation off mBW ; Figure 8b is the results from performing FRFT compensation processing (with the velocity at 4.97 m/s, and the acceleration at 59.96 m/s 2 ), it can be obtained by Equation (29) thatf mBW is approximately at 3.88 KHz, basically consistent with the theoretical value of 4 KHz, verifying the effectiveness of FRFT compensation processing; Figure 8c is the amplified results of Figure 8b, from which it can be obtained that the adjacent frequency spectrum interval is 40 Hz. Figure 8d is the STFT processing results of Figure 8b, showing that the timefrequency focus performance is poor; Figure 8e is the FSST processing results of Figure 8b, whose time-frequency focus performance has been greatly improved compared with that of Figure 8d; Figure 8f is the estimated result of the number of rotors in Figure 8b when processed by Equation (20), showing that the number of rotors is 4; based on the above results, combined with Figure 8c (with adjacent spectrum interval at 40 Hz), Equations (27) and (28), it can be obtained that the rotation frequency is 10 Hz/s and the rotor length is 9.7 cm, basically consistent with the theoretical value, and verifying the correctness of the parameter estimation of the established algorithm; Figure 8g is the FSST results without velocity and acceleration compensation processing, whose time-frequency curve has a serious deviation compared with Figure 8e, making it difficult to complete the parameter estimation accurately.
is the results from performing FRFT compensation processing (with the velocity at 4.97 m/s, and the acceleration at 59.96 m/s 2 ), it can be obtained by Equation (29) that BW m fˆ is approximately at 3.88 KHz, basically consistent with the theoretical value of 4KHz,verifying the effectiveness of FRFT compensation processing; Figure 8c is the amplified results of Figure 8b, from which it can be obtained that the adjacent frequency spectrum interval is 40 Hz. Figure 8d is the STFT processing results of Figure 8b, showing that the timefrequency focus performance is poor; Figure 8e is the FSST processing results of Figure  8b, whose time-frequency focus performance has been greatly improved compared with that of Figure 8d; Figure 8f is the estimated result of the number of rotors in Figure 8b when processed by Equation (20), showing that the number of rotors is 4; based on the above results, combined with Figure 8c (with adjacent spectrum interval at 40 Hz), Equations (27) and (28), it can be obtained that the rotation frequency is 10 Hz/s and the rotor length is 9.7 cm, basically consistent with the theoretical value, and verifying the correctness of the parameter estimation of the established algorithm; Figure 8g is the FSST results without velocity and acceleration compensation processing, whose time-frequency curve has a serious deviation compared with Figure 8e, making it difficult to complete the parameter estimation accurately.

Measured Data
In order to verify the performance of the algorithm for the detection and parameter estimation of the actual target, a P-band radar is used in the experiment to conduct data acquisition and algorithm verification on the model helicopter. The system configuration is shown in Figure 9, and the system parameters are shown in Table 5.

Measured Data
In order to verify the performance of the algorithm for the detection and parameter estimation of the actual target, a P-band radar is used in the experiment to conduct data acquisition and algorithm verification on the model helicopter. The system configuration is shown in Figure 9, and the system parameters are shown in Table 5.   For the echo signals in Figure 9, the processing results of the proposed algorithm are displayed in Figure 10. Figure 10a is the FFT of the echo signal, from which it can be seen that the SNR of the target fuselage is 5dB higher than that of the rotor. Figure 10b is the FFT results of Figure 10a when the fuselage echo is removed by using the CLEAN [38,39], that's because, usually, the echo signal noise of the fuselage is relatively large, and the fuselage echo can be removed by the CLEAN method so as to facilitate the subsequent analysis of the micro-Doppler signal. In Figure 10b (26) and (29); Figure 10c is the STFT results (wherein the window length W g obtained according to Equation (19) is 40), indicating that the timefrequency curves produced by the two rotors have poorer focus performance; Figure 10d is the FSST-processed results of the proposed algorithm, showing that the time-frequency focus is better than that in Figure 10c, and the number of rotors is 2. Substituting the above parameters into Equations (26) and (28), it can be obtained that the rotation frequency is about 37.5 Hz, and the blade length is about 680 mm, basically consistent with the actual values, and verifying the correctness of the parameter estimation by the proposed algorithm.  For the echo signals in Figure 9, the processing results of the proposed algorithm are displayed in Figure 10. Figure 10a is the FFT of the echo signal, from which it can be seen that the SNR of the target fuselage is 5 dB higher than that of the rotor. Figure 10b is the FFT results of Figure 10a when the fuselage echo is removed by using the CLEAN [38,39], that's because, usually, the echo signal noise of the fuselage is relatively large, and the fuselage echo can be removed by the CLEAN method so as to facilitate the subsequent analysis of the micro-Doppler signal. In Figure 10b,f mBW ≈ 950 Hz and ∆ f ≈ 75 Hz can be obtained by using Equations (26) and (29); Figure 10c is the STFT results (wherein the window length g W obtained according to Equation (19) is 40), indicating that the time-frequency curves produced by the two rotors have poorer focus performance; Figure 10d is the FSSTprocessed results of the proposed algorithm, showing that the time-frequency focus is better than that in Figure 10c, and the number of rotors is 2. Substituting the above parameters into Equations (26) and (28), it can be obtained that the rotation frequency is about 37.5 Hz, and the blade length is about 680 mm, basically consistent with the actual values, and verifying the correctness of the parameter estimation by the proposed algorithm.
is the FSST-processed results of the proposed algorithm, showing that the time-frequ focus is better than that in Figure 10c, and the number of rotors is 2. Substituting the parameters into Equations (26) and (28), it can be obtained that the rotation freque about 37.5 Hz, and the blade length is about 680 mm, basically consistent with the a values, and verifying the correctness of the parameter estimation by the proposed rithm.

Analysis of Algorithm Complexity
Assuming that the signal length is N , the window width is W g , the compen number of the fractional-order is ' k , the complex multiplication operational quan standard STFT in Equation (18)  , taking Section 6.1 as ample, the complexity of the proposed algorithm is about one order of magnitude than that of standard STFT while improving performance.

Conclusions
In this paper, the echo model of the rotor UAV was analyzed and FRFT was u perform the compensation processing on the velocity and acceleration of the target, improved the estimation accuracy of the micro-Doppler signal parameters of the Simultaneously, FSST was used to carry out the time-frequency analysis on the roto signal, the method to estimate the parameters such as the rotor rotation frequency number, and rotor length was given out, and the empirical Formula was given f selection of the window length of Gaussian window function, which improved the frequency resolution and operation efficiency of STFT. Simulation and experiment proved that the proposed method could accurately estimate the motion parameters UAV's rotors. Relevant research results can be applied to the detection, classificatio recognition of the low-altitude UAVs; in addition, for the signals with further re SNR, the algorithm performance can be improved combined with deep learning me

Analysis of Algorithm Complexity
Assuming that the signal length is N, the window width is g W , the compensation number of the fractional-order is k , the complex multiplication operational quantity of standard STFT in Equation (18) is O[ N 2 log 2 (g w )], the complex multiplication operational quantity of FRFT in Equation (12) is O[ N 2 k log 2 (N)], the operational quantity of SST in Equation (21) is O[g W N], and the total complex multiplication operational quantity of the proposed algorithm is approximately O[ N 2 k log 2 (N) + g W N], taking Section 6.1 as an example, the complexity of the proposed algorithm is about one order of magnitude higher than that of standard STFT while improving performance.

Conclusions
In this paper, the echo model of the rotor UAV was analyzed and FRFT was used to perform the compensation processing on the velocity and acceleration of the target, which improved the estimation accuracy of the micro-Doppler signal parameters of the rotors. Simultaneously, FSST was used to carry out the time-frequency analysis on the rotor echo signal, the method to estimate the parameters such as the rotor rotation frequency, rotor number, and rotor length was given out, and the empirical Formula was given for the selection of the window length of Gaussian window function, which improved the time-frequency resolution and operation efficiency of STFT. Simulation and experimental data proved that the proposed method could accurately estimate the motion parameters of the UAV's rotors. Relevant research results can be applied to the detection, classification, and recognition of the low-altitude UAVs; in addition, for the signals with further reduced SNR, the algorithm performance can be improved combined with deep learning methods.