# Radar-Assisted UAV Detection and Identification Based on 5G in the Internet of Things

Authors: Jingcheng Zhao; Xinru Fu; Zongkai Yang; Fengtong Xu
Journal: Wireless Communications and Mobile Computing
 (2019)

---

## Abstract

Unmanned aerial vehicles (UAVs) have broad application potential for the Internet of Things (IoT) due to their small size, low cost, and flexible control. At present, the main positioning method for UAVs is the use of GPS. However, GPS positioning may be affected by stronger electromagnetic signals from spoofing attacks. In this study, a radar-assisted positioning method based on 5G millimeter waves is proposed. In 5G end-to-end network slices, the rotors of UAVs can be detected and identified by deploying 5G millimeter wave radar. High-resolution range profile (HRRP) is used to obtain the UAV location in the detection zone. Micro-Doppler characteristics are used to identify the UAVs and the cepstrum method is used to extract the number and speed information of the UAV rotor. The sinusoidal frequency modulation (SFM) parameter optimization method is used to separate multiple UAVs. The proposed method provides information on the number of UAVs, the position of the UAV, the number of rotors, and the rotation speed of each rotor. The simulation results show that the proposed radar detection method is well suited for UAV detection and identification and provides a valid GPS-independent method for UAV tracking.


---

## Body

## 1. Introduction

With the maturity of unmanned aerial vehicle (UAV) technology and the improvement of relevant laws and regulations, UAVs are increasingly being used for the development of the Internet of Things (IoT). For example, UAVs have been broadly used in the military IoT, smart agriculture, and smart cities to obtain and transmit geospatial information, sensor data information, and controlling information. In smart agriculture, UAV systems are used to gather near real-time remote sensing data for precision farming. In such applications, UAV systems need to be positioned precisely to obtain data from sensors. In smart cities, almost all aspects of the city are combined with the IoT, a task that requires large amounts of data transmission. This data volume requires a large number of base stations for IoT data transmission. Researchers have equipped UAVs with various sensors, such as high-definition cameras, as well as temperature, humidity, and air pollution sensors. In addition, todays UAVs are equipped with high quality wireless communication functions, including 5G, Wifi, Bluetooth, radio-frequency identification (RFID), and other communication means [1]. This type of equipment means that UAVs have a very important role in the IoT and can act as a mobile air base station, providing reliable downlink and uplink communications to terrestrial users, thereby increasing the capacity of wireless networks [2]. In addition, UAVs can be used as airborne wireless relays, as well as signal boosters and mobile nodes for the IoT and sensor networks [3]. If they are used as reliable air base stations, UAVs need to be positioned precisely. Another application of UAVs is for product delivery, in which case the UAVs location has to be monitored constantly. In all the aforementioned applications, the positioning of UAVs is very important. However, due to the mobility of UAVs, in certain cases, accurate positioning remains an unsolved problem. In particular, high-precision positioning is very important for remote sensing and automatic waypoint navigation of UAVs [4]. If the UAV cannot be accurately positioned, the UAV may collide with objects or other aircraft or may be deceived, captured, or interfered with.At present, the positioning technology of UAVs is mainly achieved by using the GPS system but this method is sensitive to certain environments. Theoretically, the GPS system only requires three satellites to obtain positioning data but in practice, the UAV needs to communicate with at least 10 satellites to obtain a stable GPS system. When flying indoors, through tunnels, or dense construction areas, the GPS signal is not reliable and UAV positioning based on the GPS system will cause difficulties [5]. In addition, if UAVs are used for transporting goods, if criminals want to capture UAVs, it is difficult to prevent these attacks using GPS systems. Since the civilian GPS signal is not encrypted, the signal structure and data of UAV positioning system used in the civil IoT is predictable. This characteristic makes the civilian Global Navigation Satellite System (GNSS) signals a simple target for spoofing attacks. Moreover, after the UAV is captured, the criminal can take control of the captured UAV and fly away from its original flight path, while the original control source still receives data that appear correct and no alert is sent. This results in great risk in the use of UAVs. Many different technologies that assist in the positioning of UAVs have been developed, such as visual aids. However, visual aids are inevitably prone to drifting during long-range surveys without a priori feature maps and are only suitable for good weather and lighting conditions [6].In this study, the precise positioning of UAVs is achieved through radar-assisted detection and identification of UAV technology. The method prevents many problems associated with traditional GPS positioning technology, such as GPS signal failure indoors, satellite occlusion areas, and interference with GPS signals. The radar-assisted detection and identification method has the advantage of working under all weather conditions and is capable of extracting additional features of the UAV due to the micro-Doppler features [7]. According to the current technology development plan, 5G will adopt the previous generations of wireless solutions (3G and 4G) and introduce new promising mmWave technology so that a 5G base station not only can serve as a communication site for 5G networks but can also provide many powerful radar sensors in a functional device network, allowing for the detection of drones in urban environments [8–10]. The background of the UAV flight area is usually complex with other aircraft and birds; therefore, advanced clutter suppression and target detection algorithms are required. Following target detection, it is determined whether the detected target is a UAV and its category is identified. For the purposes of classification and identification, several features unique to UAVs are extracted [11]. In Doppler radar, the Doppler effect is used to determine the velocity of distant objects and accurate measurements can be obtained of the radial component of a targets velocity relative to the radar. Stepped frequency waveforms (SFW) are used to produce a high-resolution range profile (HRRP) of the target. With a wide-band Doppler radar, the HRRP and Doppler information of the target can be obtained to detect and identify UAVs.The rest of this paper is organized as follows. The second part introduces the system model of the radar-assisted UAV positioning method in the IoT. In the third part, the scattering model of the UAV and rotors is established, the HRRP of a rotor-type UAV is introduced, and the micro-Doppler effect of the UAV is analyzed. The cepstrum method is used to extract the speed of the UAV and the number of rotors and the method for identifying multiple UAVs is described. In the fourth part, the simulation results and one-dimensional imaging of a multitarget model of two UAVs are provided. The identification of the UAV based on the micro-Doppler information is provided and the method for extracting the number and speed of the rotors are described. The conclusions are provided in the final section.


## 2. System Model

5G mmWave technology and distributed networks will facilitate radar detection and identification of small UAVs. The narrow steerable beam enables the base station to function as a millimeter wave radar system so that the base station is able to detect and identify drones in urban environments [12, 13]. In addition, through end-to-end fragmentation technology, the 5G network provides a logically isolated virtual private segmentation network for different services to implement a manual flight-specific traffic channel network [14].The model for detecting and identifying the UAV in the 5G network is shown in Figure1. The detection network channel of the UAV is established in the 5G network and the millimeter wave base station is established. The electromagnetic signal is transmitted by the millimeter wave base station and beamforming technology is used to reduce the interference of the clutter signal [15, 16]. The reflected echo of the target is obtained by multipath reflection, and the location of UAV is transmitted through the data processing center [17, 18]. Accurate identification and positioning are achieved in a timely manner.Figure 1
UAV detection based on 5G network.


### 2.1. Radar-Assisted Detection and Identification Method

The flowchart of the process of radar-assisted detection and identification of UAVs is shown in Figure2. The 5G millimeter wave radar transmits a wideband signal to the target area and receives the echo reflected from the target. The HRRP provides the number of targets in the detection area. For a single target, the Doppler method is used to determine whether the target is a UAV and the cepstrum is used to extract the rotor parameters of the UAV. For multiple targets, after the Doppler and cepstrum analyses, the sinusoidal frequency modulation (SFM) parameter optimization method is used to obtain the position of each rotor so that multiple targets can be separated and identified.Figure 2
Flowchart of radar-based identification.


### 2.2. Stepped Frequency Waveforms

High bandwidth, low latency, and good beamforming of 5G networks provide a reliable basis for resolution because future 5G mmWave systems cannot only support very high carrier frequencies but also wider bandwidths (effective operation up to 2 GHz). Wider bandwidth facilitates the detection of UAVs flying in close formation. The SFW is a type of 5G mmWave signal and produces an HRRP of the target. A typical SFW burst is shown in Figure3; it consists of N narrow-band pulses. The frequency from pulse to pulse is stepped by a fixed frequency and each group of pulses is called a pulse train [19].Figure 3
SFW burst.Within a pulse widthτ′, the transmitted signal of the nth pulse is defined in (1)sn,t=ρncos⁡2πfnt+θn,nT≤t≤nT+τ′whereθn is the initial phase of the pulse and ρn is the amplitude constant.The echo signal of the target is defined in(2)srn,t=ρn′cos⁡2πfnt-τt+θnwhereρn′ is the amplitude of the echo signal and τt is the echo delay of the signal; τt is expressed as(3)τt=2Rtcwherec is the speed of light and Rt is the distance between the target and the radar at time t.The received signal is split in the quadrature phase detector and enters the mixer of the two synchronous detectors. In the first mixer, the reference signal and the received signal are mixed and the 90° phase-shifted echo signal is mixed with the reference signal in the second mixer. After combining the outputs of the two mixers and low pass filtering, the complex signal is(4)srn,t=ρnexp⁡-j2πf0+nΔf2RtcIf we assume that the targets’ moving distance is less than one range resolution cell for each group of pulse signals, thensrn,t is expressed as(5)srn=ρnexp⁡-j2πf0+nΔf2RcThe inverse discrete Fourier transform (IDFT) ofsrn is(6)sk=1N∑n=0N-1srnexp⁡j2πknNBy normalizing [19], (6) yields(7)sk=∑n=0N-1exp⁡j2πnNk-2RΔfcwhich can be simplified as(8)sk=sin⁡πχsin⁡πχ/Nexp⁡jN-12πχNwhereχ=-2NRΔf/c+k. The peak of sk appears at k=2NΔfR0/c±mN,m∈Z. The distance resolution is c/(2BW), where BW is the signal bandwidth.


### 2.3. Micro-Doppler Effect in Radar

In traditional radar applications, the antenna illuminates the target with a microwave signal and receives the echo reflected by the target. The echo signal includes the target characteristics of interest. For example, if the transmit signal hits a moving target, the carrier frequency of the echo signal shifts, which is called the Doppler effect. The Doppler shift reflects the moving speed of the target. Mechanical vibration or rotation of a structure on the target can also cause additional frequency shifts in the returned radar signal; this is known as the micro-Doppler effect. The micro-Doppler effect allows us to determine additional properties of the target [20].In many cases, some components of the target may have rotations or vibrations in addition to the main target translation, such as a rotor on a UAV or the flapping wing of a bird. The motion dynamics of the rotating rotor or wing will frequency-modulate the backscattered signals and cause additional Doppler shifts near the center of the translational Doppler shift, which is known as the micro-Doppler shift. Therefore, the micro-Doppler shift can be considered a unique signature of targets with vibrations, rotations, or other nonuniform motions. The micro-Doppler shift depends on the signal wavelength, the vibration or rotation speed of the target, and the incident angle of the radar line of sight (LOS) relative to the rotation or vibration plane. The reflected radar signal has the largest micro-Doppler shift when the radar LOS is parallel to the rotation or vibration plane. When the radar LOS is perpendicular to the rotation or vibration plane, the reflected signal has the smallest micro-Doppler shift. The micro-Doppler shift has been widely used to classify targets with rotations or vibrations.Almost all UAVs have at least one or more rotating rotors. For example, a four-rotor UAV has two rotors that rotate clockwise and two that rotate counterclockwise about a vertical axis. The UAV can take off vertically, hover, and fly forward, backward, and sideways. Other multirotor UAVs, such as UAVs with six or eight rotors, have matching rotor sets that rotate in opposite directions. Therefore, it is possible to detect and identify UAVs by using the micro-Doppler effect of the rotors. In this study, we extract the characteristics of UAVs based on the micro-Doppler effect [21].


## 2.1. Radar-Assisted Detection and Identification Method

The flowchart of the process of radar-assisted detection and identification of UAVs is shown in Figure2. The 5G millimeter wave radar transmits a wideband signal to the target area and receives the echo reflected from the target. The HRRP provides the number of targets in the detection area. For a single target, the Doppler method is used to determine whether the target is a UAV and the cepstrum is used to extract the rotor parameters of the UAV. For multiple targets, after the Doppler and cepstrum analyses, the sinusoidal frequency modulation (SFM) parameter optimization method is used to obtain the position of each rotor so that multiple targets can be separated and identified.Figure 2
Flowchart of radar-based identification.


## 2.2. Stepped Frequency Waveforms

High bandwidth, low latency, and good beamforming of 5G networks provide a reliable basis for resolution because future 5G mmWave systems cannot only support very high carrier frequencies but also wider bandwidths (effective operation up to 2 GHz). Wider bandwidth facilitates the detection of UAVs flying in close formation. The SFW is a type of 5G mmWave signal and produces an HRRP of the target. A typical SFW burst is shown in Figure3; it consists of N narrow-band pulses. The frequency from pulse to pulse is stepped by a fixed frequency and each group of pulses is called a pulse train [19].Figure 3
SFW burst.Within a pulse widthτ′, the transmitted signal of the nth pulse is defined in (1)sn,t=ρncos⁡2πfnt+θn,nT≤t≤nT+τ′whereθn is the initial phase of the pulse and ρn is the amplitude constant.The echo signal of the target is defined in(2)srn,t=ρn′cos⁡2πfnt-τt+θnwhereρn′ is the amplitude of the echo signal and τt is the echo delay of the signal; τt is expressed as(3)τt=2Rtcwherec is the speed of light and Rt is the distance between the target and the radar at time t.The received signal is split in the quadrature phase detector and enters the mixer of the two synchronous detectors. In the first mixer, the reference signal and the received signal are mixed and the 90° phase-shifted echo signal is mixed with the reference signal in the second mixer. After combining the outputs of the two mixers and low pass filtering, the complex signal is(4)srn,t=ρnexp⁡-j2πf0+nΔf2RtcIf we assume that the targets’ moving distance is less than one range resolution cell for each group of pulse signals, thensrn,t is expressed as(5)srn=ρnexp⁡-j2πf0+nΔf2RcThe inverse discrete Fourier transform (IDFT) ofsrn is(6)sk=1N∑n=0N-1srnexp⁡j2πknNBy normalizing [19], (6) yields(7)sk=∑n=0N-1exp⁡j2πnNk-2RΔfcwhich can be simplified as(8)sk=sin⁡πχsin⁡πχ/Nexp⁡jN-12πχNwhereχ=-2NRΔf/c+k. The peak of sk appears at k=2NΔfR0/c±mN,m∈Z. The distance resolution is c/(2BW), where BW is the signal bandwidth.


## 2.3. Micro-Doppler Effect in Radar

In traditional radar applications, the antenna illuminates the target with a microwave signal and receives the echo reflected by the target. The echo signal includes the target characteristics of interest. For example, if the transmit signal hits a moving target, the carrier frequency of the echo signal shifts, which is called the Doppler effect. The Doppler shift reflects the moving speed of the target. Mechanical vibration or rotation of a structure on the target can also cause additional frequency shifts in the returned radar signal; this is known as the micro-Doppler effect. The micro-Doppler effect allows us to determine additional properties of the target [20].In many cases, some components of the target may have rotations or vibrations in addition to the main target translation, such as a rotor on a UAV or the flapping wing of a bird. The motion dynamics of the rotating rotor or wing will frequency-modulate the backscattered signals and cause additional Doppler shifts near the center of the translational Doppler shift, which is known as the micro-Doppler shift. Therefore, the micro-Doppler shift can be considered a unique signature of targets with vibrations, rotations, or other nonuniform motions. The micro-Doppler shift depends on the signal wavelength, the vibration or rotation speed of the target, and the incident angle of the radar line of sight (LOS) relative to the rotation or vibration plane. The reflected radar signal has the largest micro-Doppler shift when the radar LOS is parallel to the rotation or vibration plane. When the radar LOS is perpendicular to the rotation or vibration plane, the reflected signal has the smallest micro-Doppler shift. The micro-Doppler shift has been widely used to classify targets with rotations or vibrations.Almost all UAVs have at least one or more rotating rotors. For example, a four-rotor UAV has two rotors that rotate clockwise and two that rotate counterclockwise about a vertical axis. The UAV can take off vertically, hover, and fly forward, backward, and sideways. Other multirotor UAVs, such as UAVs with six or eight rotors, have matching rotor sets that rotate in opposite directions. Therefore, it is possible to detect and identify UAVs by using the micro-Doppler effect of the rotors. In this study, we extract the characteristics of UAVs based on the micro-Doppler effect [21].


## 3. Detection and Identification Method of UAV

### 3.1. Scattering Model of UAV with Rotors

A UAV is a slow small target that can be simplified to a point target in traditional radar; however, since we plan to identify the UAV by using information on the rotors, we have to develop a model of a UAV in which the rotor is separate from the fuselage. Figure4 shows the scattering model of a rotor-type UAV. The distance between the center of the rotor and the radar is R0 and the elevation angle of the radar observation is β.Figure 4
Scattering model of a rotor-type UAV.The fuselage scattering of the UAV can be obtained by (4). For a rotating point with a center of rotation of l and a rotational speed of Ω, the position relative to the radar at time t can be expressed as(9)Rt=R0+lcos⁡βcos⁡Ωt+φ0Then the echo signal can be defined as(10)srn,t=ρexp⁡-j2πf0+nΔf×2R0+lcos⁡βcos⁡Ωt+φ0cWe assume that the blade length is L as shown in Figure3. Each blade can be thought of as a uniform, rigid rectangular plate that consists of arrays of rectangular facets with the scattering center of the surface of each facet assigned to its tip. The baseband signal originating from the rotors of a UAV can be expressed as(11)srn,t=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+nΔfc×R0i,k+li,m,kcosβi,kcos⁡Ωit+φ0i,mwhereNT is the total number of rotors, NM is the number of blades per rotor, and Nk represents the number of facets.


### 3.2. HRRP of Rotor-Type UAV

The movement of the UAV can be divided into two parts, namely, the movement of the main body of the UAV and the rotation of the rotor. Since the time of a one-step sweep measurement is about hundreds of milliseconds and the UAV is a small slow target, its displacement during the one-step sweep measurement is negligible. Therefore, the UAV can be regarded as a stationary target during the one-step sweep measuring time. In this study, we do not consider the movement of the fuselage.(12)srt=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+PRFΔftc×R0i,k+li,m,kcosβi,kcos⁡Ωit+φ0i,kIt is evident from (12) that the echo signal form of the rotor is an SFM signal and the equation can be expanded according to the following Bessel function:(13)ejAsin⁡θ=∑n=-∞∞JnAejnθwhere Jn· is the first type of an n-order Bessel function.The IDFT of the expanded rotor echo signal based on Eq. (13) can be obtained as [22]:(14)Srf=∑i=1NT∑m=1NM∑k=1Nk∑n=-N1N1bi,m,k,n×δf-fR,i-NMnfΩ,iwhere fR,i is the frequency offset caused by the position of the rotor center, fΩ,i is the micro-Doppler frequency caused by the rotor speed, and b is the Bessel coefficient. The signal obtained by (14) is a superposition of the pulse signals of different amplitudes with a period of NmfΩ, which is a comb-like function whose amplitude is modulated.According to the characteristics of the Bessel function:(15)∑m=-∞∞Jm2A=1The energy of the frequency-domain signal produces a divergence effect. The magnitude of the harmonic amplitude is modulated by the Bessel function. According to the Carson criterion [23], the number of harmonics is(16)N1≈2A+3The width of spectral broadening caused by the rotation of the rotor is [22](17)B≈8πfΩLcos⁡βNMλ0The rotational motion appears in the HRRP as a symmetrical range distribution centered on the position of the fuselage. When the rotor scattering is smaller than the fuselage scattering or slightly larger than the fuselage scattering, the amplitude of the rotor in the HRRP is much smaller than the amplitude of the fuselage because the rotor rotational motion is dispersed after the Bessel series expansion. For the purpose of detection, the influence of the rotor rotation on the HRRP can be neglected and the position information of the UAV can be obtained by directly reading the position of the main scattering from the fuselage. In most cases, the fuselage scattering is greater than the rotor scattering and the effects of the rotor rotation do not have to be considered for obtaining the UAV position.When the rotor scattering is much larger than the fuselage scattering, the main body scattering is submerged in the distance harmonic caused by the rotation of the rotor. In this case, the harmonic boundary caused by the rotation can be estimated and the center position of the two boundary positions can be used to estimate the target position. Alternatively, the SFM parameter optimization method described below can be used to compensate for the effect of rotor rotation on the HRRP.


### 3.3. UAV Identification

#### 3.3.1. Micro-Doppler Effect of UAV with Rotors

The Doppler shift of the received signal can be expressed as(18)fD=12πdφtdt=2cf0+PRF×Δft-ΩLcosβsin⁡Ωt+φ0+ΔfLcosβcos⁡Ωt+φ0+R0×PRF=2cLcosβΩf0+ΩPRF×Δft2+Δf×PRF2·cosΩt+θ+ΔfR0×PRFThe micro-Doppler frequency shift changes over time and the joint time-domain analysis method is the most intuitive means to observe the target’s micro-Doppler effect. Short-time Fourier transform (STFT) [24] is a commonly used joint time domain analysis methods. The basic process of the STFT is to decompose a long time signal into shorter segments of the same length, calculate the Fourier transform on each short segment, and visualize the variation of the echo frequency over time.


#### 3.3.2. Micro-Doppler Parameter Extraction Algorithm

The extraction of the parameters of the rotor, such as the number of rotors and the speed, is used to identify the type of UAV. For single-rotor UAVs, the STFT can be used for the identification. For multirotor UAVs, the micro-Doppler frequencies of different rotors are aliased together, which makes the target difficult to resolve and requires a new extraction method.It is evident in (14) that the echo signals are superimposed on the pulse signals of different amplitudes in the frequency domain with NmfΩ as the period. The period of the signal is related to the rotor speed and the rotors of different speeds produce comb-like functions of different periods. The number of rotors can be determined by extracting the number of different periodic signals and the rotational speed of each rotor of the UAV is estimated by extracting the period of the signal.Cepstrum is an effective method used in speech signal processing [25, 26] for the periodic detection of the pitch and is now also commonly used for mechanical state detection and fault diagnosis. In this study, the power cepstrum is used for parameter extraction. By multiplying the signal by its own conjugate signal, the value of the micro-Doppler can be doubled [11], which contributes to providing better periodic signal extraction. We set gt to(19)gt=st∗s∗-tA Fourier transform is performed ongt and the result is expressed as Gf. Gf can be regarded as a comb function Cf that is multiplied by some form of amplitude modulation Af. The results of Gf are shown in Figure 5.Figure 5
Function diagram ofGf.By calculating the log spectrum, the multiplicative relationship between the comb function and amplitude modulation is converted into a summation. Then, the inverse Fourier transform is performed and the square of the amplitude is used to obtain the cepstrum of the function.The specific equation is as follows:(20)Cq=IFFTlog⁡FFTgt2The cepstrum domain is shown in Figure6; the modulated function Af is a nonperiodic function whose peak appears at the zero value of the cepstrum and Gf appears as a peak at its periodic frequency. By using a high-pass filter, the influence of Gf can be removed and the rotation period of the UAV is obtained.Figure 6
Cepstrum domain.


### 3.4. Estimation of Multi-UAV Positions

The rotor rotation widens the HRRP distance. If the distance between two targets is large enough so that they can be separated in the HRRP, each target can be separately processed to determine whether it is a UAV. If the distance between the targets is very small due to distance expansion by rotation, a multi-SFM signal separation method based on parameter optimization is used to compensate for the rotational motion.


#### 3.4.1. Multi-SFM Signal Separation Based on Parameter Optimization

We discretize (12):(21)srn=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+nΔfc×R0i,k+li,m,kcosβi,kcos⁡ΩinPRF+φ0i,kTable1 shows the simplified variables in (21):Table 1
Variable simplification.
Initial parameter Simplified parameters △ f 2 R 0 c f c   f 0 + n △ f 2 l cos ⁡ β cos ⁡ Ω n / P R F + φ 0 c m f i   f 0 2 R 0 c θ   Ω P R F f m iEquation (21) can be expressed in the form of multiple SFMs:(22)xn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θiThe unknown parameters in (22) include the center frequency fci, frequency modulation factor mfi, modulation period fmi, initial phase of the SFM signal ϕi, initial phase of the FM signal θi, and the amplitude Ai. fci is related to the target position and thus we mainly focus on the estimation of fci from different signals. The parameter optimization method estimates the parameter estimation order of fci.The modulation period of the signal can be obtained by the cepstrum method and will not be described here.(a) Estimation of the modulation periodfmi.(b) Estimation of the initial phaseφi and the frequency modulation factor mfi.We assumeφiϵ0,2π, mfiϵm1,m2. The initial estimated value φi=0, mfi=m1. The reference signal can be expressed as(23)x~mi,φin=exp⁡j2πm^ficos⁡2πf^min+φ^iEquation (23) is multiplied by (22):(24)xnx~m,φn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^iFor the ith SFM component,(25)xinx~mi,φin=Aiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^i=Aiexp⁡-j2πfcin+θiexp⁡-j2π-2mfisin⁡2πfmin+φi+φ^i2sinφi-φ^i2-m^fi-mficos⁡2πf^min+φ^i(c) Estimation off^ci: when φi=φ^i, the spectrum of the signal is concentrated in a band of width 2m^fi-mfif^mi. When m^fi=mfi, the rotational motion is compensated. The original signal becomes(26)Aiexp⁡-j2πfcin+θiFast Fourier transform (FFT) is performed on (26); the position corresponding to the maximum value after the transformation is f^ci, which can be expressed as(27)f^ci=argmaxf⁡FFTxinx~mi,φin


#### 3.4.2. Rotor Position Estimation

The steps for identifying close-range multi-UAVs using the SFM parameter optimization method are as follows:(a) The cepstrum method is used to calculate the rotation period of each rotating target.(b) Since the fuselage scattering will affect the SFM parameter optimization method, the influence of the subject scattering must be removed prior to using the SFM parameter optimization method. First, two subsequent frames of signals are measured and subtracted from each other for the cancellation of the fuselage scattering. Second, HRRP is used to determine the maximum value of the rotor scattering. Finally, the maximum value is set as a threshold. If the single frame in the HRRP is larger than the threshold, the single signal is set as the mean value of the subtracted HRRP.(c) The SFM parameter optimization method is used for the rotation periods of the rotors obtained in step (a) to determine their position.


## 3.1. Scattering Model of UAV with Rotors

A UAV is a slow small target that can be simplified to a point target in traditional radar; however, since we plan to identify the UAV by using information on the rotors, we have to develop a model of a UAV in which the rotor is separate from the fuselage. Figure4 shows the scattering model of a rotor-type UAV. The distance between the center of the rotor and the radar is R0 and the elevation angle of the radar observation is β.Figure 4
Scattering model of a rotor-type UAV.The fuselage scattering of the UAV can be obtained by (4). For a rotating point with a center of rotation of l and a rotational speed of Ω, the position relative to the radar at time t can be expressed as(9)Rt=R0+lcos⁡βcos⁡Ωt+φ0Then the echo signal can be defined as(10)srn,t=ρexp⁡-j2πf0+nΔf×2R0+lcos⁡βcos⁡Ωt+φ0cWe assume that the blade length is L as shown in Figure3. Each blade can be thought of as a uniform, rigid rectangular plate that consists of arrays of rectangular facets with the scattering center of the surface of each facet assigned to its tip. The baseband signal originating from the rotors of a UAV can be expressed as(11)srn,t=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+nΔfc×R0i,k+li,m,kcosβi,kcos⁡Ωit+φ0i,mwhereNT is the total number of rotors, NM is the number of blades per rotor, and Nk represents the number of facets.


## 3.2. HRRP of Rotor-Type UAV

The movement of the UAV can be divided into two parts, namely, the movement of the main body of the UAV and the rotation of the rotor. Since the time of a one-step sweep measurement is about hundreds of milliseconds and the UAV is a small slow target, its displacement during the one-step sweep measurement is negligible. Therefore, the UAV can be regarded as a stationary target during the one-step sweep measuring time. In this study, we do not consider the movement of the fuselage.(12)srt=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+PRFΔftc×R0i,k+li,m,kcosβi,kcos⁡Ωit+φ0i,kIt is evident from (12) that the echo signal form of the rotor is an SFM signal and the equation can be expanded according to the following Bessel function:(13)ejAsin⁡θ=∑n=-∞∞JnAejnθwhere Jn· is the first type of an n-order Bessel function.The IDFT of the expanded rotor echo signal based on Eq. (13) can be obtained as [22]:(14)Srf=∑i=1NT∑m=1NM∑k=1Nk∑n=-N1N1bi,m,k,n×δf-fR,i-NMnfΩ,iwhere fR,i is the frequency offset caused by the position of the rotor center, fΩ,i is the micro-Doppler frequency caused by the rotor speed, and b is the Bessel coefficient. The signal obtained by (14) is a superposition of the pulse signals of different amplitudes with a period of NmfΩ, which is a comb-like function whose amplitude is modulated.According to the characteristics of the Bessel function:(15)∑m=-∞∞Jm2A=1The energy of the frequency-domain signal produces a divergence effect. The magnitude of the harmonic amplitude is modulated by the Bessel function. According to the Carson criterion [23], the number of harmonics is(16)N1≈2A+3The width of spectral broadening caused by the rotation of the rotor is [22](17)B≈8πfΩLcos⁡βNMλ0The rotational motion appears in the HRRP as a symmetrical range distribution centered on the position of the fuselage. When the rotor scattering is smaller than the fuselage scattering or slightly larger than the fuselage scattering, the amplitude of the rotor in the HRRP is much smaller than the amplitude of the fuselage because the rotor rotational motion is dispersed after the Bessel series expansion. For the purpose of detection, the influence of the rotor rotation on the HRRP can be neglected and the position information of the UAV can be obtained by directly reading the position of the main scattering from the fuselage. In most cases, the fuselage scattering is greater than the rotor scattering and the effects of the rotor rotation do not have to be considered for obtaining the UAV position.When the rotor scattering is much larger than the fuselage scattering, the main body scattering is submerged in the distance harmonic caused by the rotation of the rotor. In this case, the harmonic boundary caused by the rotation can be estimated and the center position of the two boundary positions can be used to estimate the target position. Alternatively, the SFM parameter optimization method described below can be used to compensate for the effect of rotor rotation on the HRRP.


## 3.3. UAV Identification

### 3.3.1. Micro-Doppler Effect of UAV with Rotors

The Doppler shift of the received signal can be expressed as(18)fD=12πdφtdt=2cf0+PRF×Δft-ΩLcosβsin⁡Ωt+φ0+ΔfLcosβcos⁡Ωt+φ0+R0×PRF=2cLcosβΩf0+ΩPRF×Δft2+Δf×PRF2·cosΩt+θ+ΔfR0×PRFThe micro-Doppler frequency shift changes over time and the joint time-domain analysis method is the most intuitive means to observe the target’s micro-Doppler effect. Short-time Fourier transform (STFT) [24] is a commonly used joint time domain analysis methods. The basic process of the STFT is to decompose a long time signal into shorter segments of the same length, calculate the Fourier transform on each short segment, and visualize the variation of the echo frequency over time.


### 3.3.2. Micro-Doppler Parameter Extraction Algorithm

The extraction of the parameters of the rotor, such as the number of rotors and the speed, is used to identify the type of UAV. For single-rotor UAVs, the STFT can be used for the identification. For multirotor UAVs, the micro-Doppler frequencies of different rotors are aliased together, which makes the target difficult to resolve and requires a new extraction method.It is evident in (14) that the echo signals are superimposed on the pulse signals of different amplitudes in the frequency domain with NmfΩ as the period. The period of the signal is related to the rotor speed and the rotors of different speeds produce comb-like functions of different periods. The number of rotors can be determined by extracting the number of different periodic signals and the rotational speed of each rotor of the UAV is estimated by extracting the period of the signal.Cepstrum is an effective method used in speech signal processing [25, 26] for the periodic detection of the pitch and is now also commonly used for mechanical state detection and fault diagnosis. In this study, the power cepstrum is used for parameter extraction. By multiplying the signal by its own conjugate signal, the value of the micro-Doppler can be doubled [11], which contributes to providing better periodic signal extraction. We set gt to(19)gt=st∗s∗-tA Fourier transform is performed ongt and the result is expressed as Gf. Gf can be regarded as a comb function Cf that is multiplied by some form of amplitude modulation Af. The results of Gf are shown in Figure 5.Figure 5
Function diagram ofGf.By calculating the log spectrum, the multiplicative relationship between the comb function and amplitude modulation is converted into a summation. Then, the inverse Fourier transform is performed and the square of the amplitude is used to obtain the cepstrum of the function.The specific equation is as follows:(20)Cq=IFFTlog⁡FFTgt2The cepstrum domain is shown in Figure6; the modulated function Af is a nonperiodic function whose peak appears at the zero value of the cepstrum and Gf appears as a peak at its periodic frequency. By using a high-pass filter, the influence of Gf can be removed and the rotation period of the UAV is obtained.Figure 6
Cepstrum domain.


## 3.3.1. Micro-Doppler Effect of UAV with Rotors

The Doppler shift of the received signal can be expressed as(18)fD=12πdφtdt=2cf0+PRF×Δft-ΩLcosβsin⁡Ωt+φ0+ΔfLcosβcos⁡Ωt+φ0+R0×PRF=2cLcosβΩf0+ΩPRF×Δft2+Δf×PRF2·cosΩt+θ+ΔfR0×PRFThe micro-Doppler frequency shift changes over time and the joint time-domain analysis method is the most intuitive means to observe the target’s micro-Doppler effect. Short-time Fourier transform (STFT) [24] is a commonly used joint time domain analysis methods. The basic process of the STFT is to decompose a long time signal into shorter segments of the same length, calculate the Fourier transform on each short segment, and visualize the variation of the echo frequency over time.


## 3.3.2. Micro-Doppler Parameter Extraction Algorithm

The extraction of the parameters of the rotor, such as the number of rotors and the speed, is used to identify the type of UAV. For single-rotor UAVs, the STFT can be used for the identification. For multirotor UAVs, the micro-Doppler frequencies of different rotors are aliased together, which makes the target difficult to resolve and requires a new extraction method.It is evident in (14) that the echo signals are superimposed on the pulse signals of different amplitudes in the frequency domain with NmfΩ as the period. The period of the signal is related to the rotor speed and the rotors of different speeds produce comb-like functions of different periods. The number of rotors can be determined by extracting the number of different periodic signals and the rotational speed of each rotor of the UAV is estimated by extracting the period of the signal.Cepstrum is an effective method used in speech signal processing [25, 26] for the periodic detection of the pitch and is now also commonly used for mechanical state detection and fault diagnosis. In this study, the power cepstrum is used for parameter extraction. By multiplying the signal by its own conjugate signal, the value of the micro-Doppler can be doubled [11], which contributes to providing better periodic signal extraction. We set gt to(19)gt=st∗s∗-tA Fourier transform is performed ongt and the result is expressed as Gf. Gf can be regarded as a comb function Cf that is multiplied by some form of amplitude modulation Af. The results of Gf are shown in Figure 5.Figure 5
Function diagram ofGf.By calculating the log spectrum, the multiplicative relationship between the comb function and amplitude modulation is converted into a summation. Then, the inverse Fourier transform is performed and the square of the amplitude is used to obtain the cepstrum of the function.The specific equation is as follows:(20)Cq=IFFTlog⁡FFTgt2The cepstrum domain is shown in Figure6; the modulated function Af is a nonperiodic function whose peak appears at the zero value of the cepstrum and Gf appears as a peak at its periodic frequency. By using a high-pass filter, the influence of Gf can be removed and the rotation period of the UAV is obtained.Figure 6
Cepstrum domain.


## 3.4. Estimation of Multi-UAV Positions

The rotor rotation widens the HRRP distance. If the distance between two targets is large enough so that they can be separated in the HRRP, each target can be separately processed to determine whether it is a UAV. If the distance between the targets is very small due to distance expansion by rotation, a multi-SFM signal separation method based on parameter optimization is used to compensate for the rotational motion.


### 3.4.1. Multi-SFM Signal Separation Based on Parameter Optimization

We discretize (12):(21)srn=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+nΔfc×R0i,k+li,m,kcosβi,kcos⁡ΩinPRF+φ0i,kTable1 shows the simplified variables in (21):Table 1
Variable simplification.
Initial parameter Simplified parameters △ f 2 R 0 c f c   f 0 + n △ f 2 l cos ⁡ β cos ⁡ Ω n / P R F + φ 0 c m f i   f 0 2 R 0 c θ   Ω P R F f m iEquation (21) can be expressed in the form of multiple SFMs:(22)xn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θiThe unknown parameters in (22) include the center frequency fci, frequency modulation factor mfi, modulation period fmi, initial phase of the SFM signal ϕi, initial phase of the FM signal θi, and the amplitude Ai. fci is related to the target position and thus we mainly focus on the estimation of fci from different signals. The parameter optimization method estimates the parameter estimation order of fci.The modulation period of the signal can be obtained by the cepstrum method and will not be described here.(a) Estimation of the modulation periodfmi.(b) Estimation of the initial phaseφi and the frequency modulation factor mfi.We assumeφiϵ0,2π, mfiϵm1,m2. The initial estimated value φi=0, mfi=m1. The reference signal can be expressed as(23)x~mi,φin=exp⁡j2πm^ficos⁡2πf^min+φ^iEquation (23) is multiplied by (22):(24)xnx~m,φn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^iFor the ith SFM component,(25)xinx~mi,φin=Aiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^i=Aiexp⁡-j2πfcin+θiexp⁡-j2π-2mfisin⁡2πfmin+φi+φ^i2sinφi-φ^i2-m^fi-mficos⁡2πf^min+φ^i(c) Estimation off^ci: when φi=φ^i, the spectrum of the signal is concentrated in a band of width 2m^fi-mfif^mi. When m^fi=mfi, the rotational motion is compensated. The original signal becomes(26)Aiexp⁡-j2πfcin+θiFast Fourier transform (FFT) is performed on (26); the position corresponding to the maximum value after the transformation is f^ci, which can be expressed as(27)f^ci=argmaxf⁡FFTxinx~mi,φin


### 3.4.2. Rotor Position Estimation

The steps for identifying close-range multi-UAVs using the SFM parameter optimization method are as follows:(a) The cepstrum method is used to calculate the rotation period of each rotating target.(b) Since the fuselage scattering will affect the SFM parameter optimization method, the influence of the subject scattering must be removed prior to using the SFM parameter optimization method. First, two subsequent frames of signals are measured and subtracted from each other for the cancellation of the fuselage scattering. Second, HRRP is used to determine the maximum value of the rotor scattering. Finally, the maximum value is set as a threshold. If the single frame in the HRRP is larger than the threshold, the single signal is set as the mean value of the subtracted HRRP.(c) The SFM parameter optimization method is used for the rotation periods of the rotors obtained in step (a) to determine their position.


## 3.4.1. Multi-SFM Signal Separation Based on Parameter Optimization

We discretize (12):(21)srn=∑i=1NT∑m=1NM∑k=1Nkρi,m,kexp⁡-j4πf0+nΔfc×R0i,k+li,m,kcosβi,kcos⁡ΩinPRF+φ0i,kTable1 shows the simplified variables in (21):Table 1
Variable simplification.
Initial parameter Simplified parameters △ f 2 R 0 c f c   f 0 + n △ f 2 l cos ⁡ β cos ⁡ Ω n / P R F + φ 0 c m f i   f 0 2 R 0 c θ   Ω P R F f m iEquation (21) can be expressed in the form of multiple SFMs:(22)xn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θiThe unknown parameters in (22) include the center frequency fci, frequency modulation factor mfi, modulation period fmi, initial phase of the SFM signal ϕi, initial phase of the FM signal θi, and the amplitude Ai. fci is related to the target position and thus we mainly focus on the estimation of fci from different signals. The parameter optimization method estimates the parameter estimation order of fci.The modulation period of the signal can be obtained by the cepstrum method and will not be described here.(a) Estimation of the modulation periodfmi.(b) Estimation of the initial phaseφi and the frequency modulation factor mfi.We assumeφiϵ0,2π, mfiϵm1,m2. The initial estimated value φi=0, mfi=m1. The reference signal can be expressed as(23)x~mi,φin=exp⁡j2πm^ficos⁡2πf^min+φ^iEquation (23) is multiplied by (22):(24)xnx~m,φn=∑i=1MAiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^iFor the ith SFM component,(25)xinx~mi,φin=Aiexp⁡-j2πfcin+mficos⁡2πfmin+φi+θi+j2πm^ficos⁡2πf^min+φ^i=Aiexp⁡-j2πfcin+θiexp⁡-j2π-2mfisin⁡2πfmin+φi+φ^i2sinφi-φ^i2-m^fi-mficos⁡2πf^min+φ^i(c) Estimation off^ci: when φi=φ^i, the spectrum of the signal is concentrated in a band of width 2m^fi-mfif^mi. When m^fi=mfi, the rotational motion is compensated. The original signal becomes(26)Aiexp⁡-j2πfcin+θiFast Fourier transform (FFT) is performed on (26); the position corresponding to the maximum value after the transformation is f^ci, which can be expressed as(27)f^ci=argmaxf⁡FFTxinx~mi,φin


## 3.4.2. Rotor Position Estimation

The steps for identifying close-range multi-UAVs using the SFM parameter optimization method are as follows:(a) The cepstrum method is used to calculate the rotation period of each rotating target.(b) Since the fuselage scattering will affect the SFM parameter optimization method, the influence of the subject scattering must be removed prior to using the SFM parameter optimization method. First, two subsequent frames of signals are measured and subtracted from each other for the cancellation of the fuselage scattering. Second, HRRP is used to determine the maximum value of the rotor scattering. Finally, the maximum value is set as a threshold. If the single frame in the HRRP is larger than the threshold, the single signal is set as the mean value of the subtracted HRRP.(c) The SFM parameter optimization method is used for the rotation periods of the rotors obtained in step (a) to determine their position.


## 4. Simulation

### 4.1. Simulation Model

Figure7 shows a schematic of the simulation model. There are three targets in the observation area, i.e., a quad-rotor UAV at a distance of 300 m from the radar, a non-UAV target 800 m from the radar, and a single-rotor UAV 850 m from the radar. The rotor length of the quad-rotor UAV is 0.12 m and the rotational speeds are 158 r/s, 150 r/s, 142 r/s, and 139 r/s, respectively. The rotor length of the single-rotor UAV is 0.15 m and the rotational speed is 88 r/s. The parameters of the radar are as follows: starting frequency F0=25 GHz, which is in the 5G band; signal bandwidth BW = 500 MHz, the resolution cell is 0.3m; sampling points N = 8001; PRF=150 kHz.Figure 7
Schematic of the simulation model.


### 4.2. Simulation Results

The HRRP of the echo signal is shown in Figure8. Three targets are detected at distances of 300 m, 799.2 m, and 850.2 m from the radar. Target 2 and target 3 cannot be separated in the HRRP because they are in close proximity and have to be analyzed as one target. The SFM optimization method can be used to determine the UAV locations. Target 1 is far from target 2 and target 3. Therefore, the time-domain echo signal of target 1 can be separated and identified by using a distance-domain window and FFT.Figure 8
HRRP of the target.


#### 4.2.1. Identification of Target 1

Figure9 shows the result of the micro-Doppler characteristics of target 1. It can be inferred that target 1 is a rotor-type UAV. Due to the presence of the stepping frequency, the echo signal appears as a fixed frequency value that does not change over time in the time-frequency image. The micro-Doppler frequency caused by the rotor changes sinusoidally over time. The difference in the amplitude of the micro-Doppler at different moments is due to the different contributions of rotor scattering to echoes at different radar observation angles.Figure 9
Micro-Doppler of target 1.There are multiple rotors in this target; therefore, the micro-Doppler signal is aliased and the type parameters of the UAV rotor cannot be directly obtained. The cepstrum method is used to determine the rotor type.Figure10 shows the cepstrum extraction results. It can be seen that target 1 is a UAV with 4 rotors. The rotational speed of the rotors is 157.9 r/s, 150.0 r/s, 142.0 r/s, and 138.9 r/s. The speed obtained by the cepstrum method is nearly the same as the predetermined speed of the rotors.Figure 10
Cepstrum extraction of target 1.In order to verify the robustness of the cepstrum method, noise with different signal-to-noise ratios (SNRs) is added to the echo signal of target 1. Figure11 shows the simulation results with SNR = -5, SNR = 0, and SNR = 5. It is evident from Figure 11 that after adding the noise, the peak energy of the position corresponding to the rotation period of the rotor is slightly lower. However, there is still a significant difference between the rotor values and the values of other positions and the addition of the noise has no influence on the identification of the number of rotors and the rotational speed.Figure 11
Cepstrum results with the addition of different SNRs.
(a)
SNR= −5dB (b)
SNR= 0dB (c)
SNR= 5dB


#### 4.2.2. Identifications of Target 2 and Target 3

Target 2 and target 3 are analyzed as a single target. Figure12 shows the micro-Doppler characteristics of the two targets.Figure 12
Micro-Doppler of target 2 and target 3.As shown in Figure12, there is one rotor for target 2 and target 3. Figure 13 shows that the rotational speed of this rotor obtained by the cepstrum method is 88 r/s.Figure 13
Cepstrum result of target 2 and target 3.The position of the single rotor is estimated by the SFM parameter optimization method after the two target body scatterings are removed. The estimated results are shown in Figure14. The parameters of the three targets in the detection area are shown in Table 2.Table 2
Target parameters.
Target type Position (m) Rotor speed (r/s) Target 1 Quad-rotor UAV 300 157.9, 150.0, 142.0, 138.9  Target 2 Other targets 799.2 -  Target 3 One-rotor UAV 850.2 88Figure 14
Rotor position estimation.Table2 shows the good accuracy of the detection of the UAV and its position. Figure 15 shows a comparison of the UAVs predetermined route based on the GNSS signal and the actual route based on radar tracking. In the case of a GNSS spoofing attack, the captured UAV has flown away from its predetermined route but the original control source was not alerted because it received data that appeared to be correct from the attackers. In this case, the radar-assisted detection and identification method can be used for real-time monitoring of the UAV independent from the GNSS signal; the deviation between the actual route and predetermined route can be determined and an alert can be sent. It is evident that there is a considerable deviation between the GPS-based positions and the radar positions. The latter method provides accurate positioning of the UAV and allows for effective monitoring of UAVs in the IoT.Figure 15
UAV flight route.


## 4.1. Simulation Model

Figure7 shows a schematic of the simulation model. There are three targets in the observation area, i.e., a quad-rotor UAV at a distance of 300 m from the radar, a non-UAV target 800 m from the radar, and a single-rotor UAV 850 m from the radar. The rotor length of the quad-rotor UAV is 0.12 m and the rotational speeds are 158 r/s, 150 r/s, 142 r/s, and 139 r/s, respectively. The rotor length of the single-rotor UAV is 0.15 m and the rotational speed is 88 r/s. The parameters of the radar are as follows: starting frequency F0=25 GHz, which is in the 5G band; signal bandwidth BW = 500 MHz, the resolution cell is 0.3m; sampling points N = 8001; PRF=150 kHz.Figure 7
Schematic of the simulation model.


## 4.2. Simulation Results

The HRRP of the echo signal is shown in Figure8. Three targets are detected at distances of 300 m, 799.2 m, and 850.2 m from the radar. Target 2 and target 3 cannot be separated in the HRRP because they are in close proximity and have to be analyzed as one target. The SFM optimization method can be used to determine the UAV locations. Target 1 is far from target 2 and target 3. Therefore, the time-domain echo signal of target 1 can be separated and identified by using a distance-domain window and FFT.Figure 8
HRRP of the target.


### 4.2.1. Identification of Target 1

Figure9 shows the result of the micro-Doppler characteristics of target 1. It can be inferred that target 1 is a rotor-type UAV. Due to the presence of the stepping frequency, the echo signal appears as a fixed frequency value that does not change over time in the time-frequency image. The micro-Doppler frequency caused by the rotor changes sinusoidally over time. The difference in the amplitude of the micro-Doppler at different moments is due to the different contributions of rotor scattering to echoes at different radar observation angles.Figure 9
Micro-Doppler of target 1.There are multiple rotors in this target; therefore, the micro-Doppler signal is aliased and the type parameters of the UAV rotor cannot be directly obtained. The cepstrum method is used to determine the rotor type.Figure10 shows the cepstrum extraction results. It can be seen that target 1 is a UAV with 4 rotors. The rotational speed of the rotors is 157.9 r/s, 150.0 r/s, 142.0 r/s, and 138.9 r/s. The speed obtained by the cepstrum method is nearly the same as the predetermined speed of the rotors.Figure 10
Cepstrum extraction of target 1.In order to verify the robustness of the cepstrum method, noise with different signal-to-noise ratios (SNRs) is added to the echo signal of target 1. Figure11 shows the simulation results with SNR = -5, SNR = 0, and SNR = 5. It is evident from Figure 11 that after adding the noise, the peak energy of the position corresponding to the rotation period of the rotor is slightly lower. However, there is still a significant difference between the rotor values and the values of other positions and the addition of the noise has no influence on the identification of the number of rotors and the rotational speed.Figure 11
Cepstrum results with the addition of different SNRs.
(a)
SNR= −5dB (b)
SNR= 0dB (c)
SNR= 5dB


### 4.2.2. Identifications of Target 2 and Target 3

Target 2 and target 3 are analyzed as a single target. Figure12 shows the micro-Doppler characteristics of the two targets.Figure 12
Micro-Doppler of target 2 and target 3.As shown in Figure12, there is one rotor for target 2 and target 3. Figure 13 shows that the rotational speed of this rotor obtained by the cepstrum method is 88 r/s.Figure 13
Cepstrum result of target 2 and target 3.The position of the single rotor is estimated by the SFM parameter optimization method after the two target body scatterings are removed. The estimated results are shown in Figure14. The parameters of the three targets in the detection area are shown in Table 2.Table 2
Target parameters.
Target type Position (m) Rotor speed (r/s) Target 1 Quad-rotor UAV 300 157.9, 150.0, 142.0, 138.9  Target 2 Other targets 799.2 -  Target 3 One-rotor UAV 850.2 88Figure 14
Rotor position estimation.Table2 shows the good accuracy of the detection of the UAV and its position. Figure 15 shows a comparison of the UAVs predetermined route based on the GNSS signal and the actual route based on radar tracking. In the case of a GNSS spoofing attack, the captured UAV has flown away from its predetermined route but the original control source was not alerted because it received data that appeared to be correct from the attackers. In this case, the radar-assisted detection and identification method can be used for real-time monitoring of the UAV independent from the GNSS signal; the deviation between the actual route and predetermined route can be determined and an alert can be sent. It is evident that there is a considerable deviation between the GPS-based positions and the radar positions. The latter method provides accurate positioning of the UAV and allows for effective monitoring of UAVs in the IoT.Figure 15
UAV flight route.


## 4.2.1. Identification of Target 1

Figure9 shows the result of the micro-Doppler characteristics of target 1. It can be inferred that target 1 is a rotor-type UAV. Due to the presence of the stepping frequency, the echo signal appears as a fixed frequency value that does not change over time in the time-frequency image. The micro-Doppler frequency caused by the rotor changes sinusoidally over time. The difference in the amplitude of the micro-Doppler at different moments is due to the different contributions of rotor scattering to echoes at different radar observation angles.Figure 9
Micro-Doppler of target 1.There are multiple rotors in this target; therefore, the micro-Doppler signal is aliased and the type parameters of the UAV rotor cannot be directly obtained. The cepstrum method is used to determine the rotor type.Figure10 shows the cepstrum extraction results. It can be seen that target 1 is a UAV with 4 rotors. The rotational speed of the rotors is 157.9 r/s, 150.0 r/s, 142.0 r/s, and 138.9 r/s. The speed obtained by the cepstrum method is nearly the same as the predetermined speed of the rotors.Figure 10
Cepstrum extraction of target 1.In order to verify the robustness of the cepstrum method, noise with different signal-to-noise ratios (SNRs) is added to the echo signal of target 1. Figure11 shows the simulation results with SNR = -5, SNR = 0, and SNR = 5. It is evident from Figure 11 that after adding the noise, the peak energy of the position corresponding to the rotation period of the rotor is slightly lower. However, there is still a significant difference between the rotor values and the values of other positions and the addition of the noise has no influence on the identification of the number of rotors and the rotational speed.Figure 11
Cepstrum results with the addition of different SNRs.
(a)
SNR= −5dB (b)
SNR= 0dB (c)
SNR= 5dB


## 4.2.2. Identifications of Target 2 and Target 3

Target 2 and target 3 are analyzed as a single target. Figure12 shows the micro-Doppler characteristics of the two targets.Figure 12
Micro-Doppler of target 2 and target 3.As shown in Figure12, there is one rotor for target 2 and target 3. Figure 13 shows that the rotational speed of this rotor obtained by the cepstrum method is 88 r/s.Figure 13
Cepstrum result of target 2 and target 3.The position of the single rotor is estimated by the SFM parameter optimization method after the two target body scatterings are removed. The estimated results are shown in Figure14. The parameters of the three targets in the detection area are shown in Table 2.Table 2
Target parameters.
Target type Position (m) Rotor speed (r/s) Target 1 Quad-rotor UAV 300 157.9, 150.0, 142.0, 138.9  Target 2 Other targets 799.2 -  Target 3 One-rotor UAV 850.2 88Figure 14
Rotor position estimation.Table2 shows the good accuracy of the detection of the UAV and its position. Figure 15 shows a comparison of the UAVs predetermined route based on the GNSS signal and the actual route based on radar tracking. In the case of a GNSS spoofing attack, the captured UAV has flown away from its predetermined route but the original control source was not alerted because it received data that appeared to be correct from the attackers. In this case, the radar-assisted detection and identification method can be used for real-time monitoring of the UAV independent from the GNSS signal; the deviation between the actual route and predetermined route can be determined and an alert can be sent. It is evident that there is a considerable deviation between the GPS-based positions and the radar positions. The latter method provides accurate positioning of the UAV and allows for effective monitoring of UAVs in the IoT.Figure 15
UAV flight route.


## 5. Conclusion

In this paper, a radar-assisted UAV detection and identification method that is independent of a GPS system was proposed. The HRRP technology of wideband radar was used for UAV detection and positioning and the micro-Doppler signal, which is capable of detecting rotating targets, was used for UAV identification. Cepstrum analysis was used to extract the number and rotation speed of the UAV rotors. The simulation results demonstrated the good robustness of this method. An SFM parameter optimization algorithm was used to compensate for the rotational motion and estimate the rotor positions; two UAVs with aliasing signals in the HRRP were separated effectively using the proposed method. The simulation results showed that the UAVs can be identified and the number of UAVs, the number of rotors on each UAV, and the rotation speeds of the rotors can be determined. In addition, the proposed radar-assisted UAV detection and identification method can provide alerts if UAVs deviate from their routes during GNSS spoofing attacks.



---

*Source: 2850263-2019-07-14.xml*