Skip to main content

Drone Detection via Low Complexity Zadoff-Chu Sequence Root Estimation


Summary

Based on 
Yeung, 2025: Yeung, C.K.A., Lo, B.F. and Torborg, S. Drone detection via low complexity zadoff-chu sequence root estimation. In 2020 IEEE 17th Annual Consumer Communications & Networking Conference (CCNC) (pp. 1-4). IEEE, 2020, January.
 

The rise in drone usage—from agriculture and delivery to surveillance and racing—has introduced major privacy and security challenges. Modern drones often use OFDM (Orthogonal Frequency Division Multiplexing) with Zadoff-Chu (ZC) sequences for synchronization. While powerful, detecting these sequences blindly (without knowing their parameters) remains a challenge.

Aim

This article presents a low-complexity solution to blindly detect ZC sequences used by unknown drones. The approach uses a novel double differential method that works without large correlation banks, making it efficient and real-time capable.

ZC Sequence Fundamentals

A ZC sequence of prime length P and root u is defined by:

x_u[n] = exp(-jπu·n(n+1)/P),  0 ≤ n < P
  • P is an odd prime (e.g., 599, 1021).
  • Each root u creates a distinct ZC sequence.
  • ZC sequences have zero autocorrelation and low cross-correlation—ideal for synchronization.

Time Domain Detection with Unknown Root

The paper proposes using single and double differential operators to estimate the ZC root:


y_d[n]   = y[n] * conj(y[n-1])
y_dd[n]  = y_d[n] * conj(y_d[n-1])
  

These steps eliminate frequency offset and isolate the root information. The detection metric is:


γ_td[n] = ∑ y_dd[n] / sqrt(∑ |y_dd[n]|²)
û = round(P × angle(γ_td[n*]) / 2π)
  

SNR Bounds

The effective SNR of this method is:

SNR_td = (P - 2) · |h|⁸ / E[|w_dd[n]|²]

Asymptotic bounds:

High SNR: SNR_td ≥ (P - 2)/(8 + 2√10)
Low SNR: SNR_td ≥ (P - 2)·|h|⁸ / (2σ⁸)

Frequency Domain ZC Without Guard Band

When ZC sequences occupy all subcarriers, the detection is similar. Using inverse DFT properties and root inverses:


γ_fd1[n] = ∑ Y_dd[n] / sqrt(∑ |Y_dd[n]|²)
û⁻¹ = round(P × angle(γ_fd1[n*]) / 2π)
û = modular_inverse(û⁻¹, P)
  

SNR bounds remain the same as time domain.

Frequency Domain ZC with Guard Bands

In practical systems, edge subcarriers are unused. ZC detection must then work with only the middle subcarriers:


Y[k] = x_u[k - (P-1)/2] · exp(-j2πΔk/N) + W[k]
Y_dd[k] = exp(-2jπu/P)
γ_fd2 = ∑ Y_dd[k]
û = round(P × angle(γ_fd2[(P-1)/2]) / 2π)
  

CFO (Carrier Frequency Offset) has a stronger negative effect here than in time domain methods.

Simulation Results

Monte Carlo simulations (30,000 runs) verify theoretical predictions:

  • Time and frequency domain methods reach 1% error rate at 15 dB SNR.
  • CFO has minimal impact in time domain and full-bandwidth frequency domain.
  • Guard-banded frequency domain suffers significantly under CFO.

Conclusion

The proposed double differential blind ZC detection method is a major step toward low-cost, real-time drone detection for consumer applications. It's lightweight, scalable, and resilient to synchronization challenges—ideal for future smart cities and surveillance tools.



Further Reading

People are good at skipping over material they already know!

View Related Topics to







Admin & Author: Salim

s

  Website: www.salimwireless.com
  Interests: Signal Processing, Telecommunication, 5G Technology, Present & Future Wireless Technologies, Digital Signal Processing, Computer Networks, Millimeter Wave Band Channel, Web Development
  Seeking an opportunity in the Teaching or Electronics & Telecommunication domains.
  Possess M.Tech in Electronic Communication Systems.


Contact Us

Name

Email *

Message *

Popular Posts

Constellation Diagrams of M-ary QAM | M-ary Modulation

📘 Overview of QAM 🧮 MATLAB Code for m-ary QAM (4-QAM, 16-QAM, 32-QAM, ...) 🧮 Online Simulator for M-ary QAM Constellations 📚 Further Reading 📂 Other Topics on Constellation Diagrams of QAM configurations ... 🧮 MATLAB Code for 4-QAM 🧮 MATLAB Code for 16-QAM 🧮 MATLAB Code for m-ary QAM (4-QAM, 16-QAM, 32-QAM, ...) 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM 🧮 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK QAM Unlike this, the M-ary PSK signal is modulated with a different phase-shifted version of the carrier signal and varying amplitude levels. Let me give an example for better comprehension. QAM = ASK +...

Theoretical vs. simulated BER vs. SNR for ASK, FSK, and PSK

📘 Overview 🧮 Simulator for calculating BER 🧮 MATLAB Codes for calculating theoretical BER 🧮 MATLAB Codes for calculating simulated BER 📚 Further Reading   BER vs. SNR denotes how many bits in error are received in a communication process for a particular Signal-to-noise (SNR) ratio. In most cases, SNR is measured in decibel (dB). For a typical communication system, a signal is often affected by two types of noises 1. Additive White Gaussian Noise (AWGN) 2. Rayleigh Fading In the case of additive white Gaussian noise (AWGN), random magnitude is added to the transmitted signal. On the other hand, Rayleigh fading (due to multipath) attenuates the different frequency components of a signal differently. A good signal-to-noise ratio tries to mitigate the effect of noise.  Simulator for calculating BER vs SNR for binary ASK, FSK, and PSK Calculate BER for Binary ASK Modulation The theoretical BER for binary ASK (BASK) in an AWGN channel is...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

📘 Overview of BER and SNR 🧮 Online Simulator for BER calculation of m-ary QAM and m-ary PSK 🧮 MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... 📚 Further Reading 📂 View Other Topics on M-ary QAM, M-ary PSK, QPSK ... 🧮 Online Simulator for Constellation Diagram of m-ary QAM 🧮 Online Simulator for Constellation Diagram of m-ary PSK 🧮 MATLAB Code for BER calculation of ASK, FSK, and PSK 🧮 MATLAB Code for BER calculation of Alamouti Scheme 🧮 Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for bit error rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. It is defined as,  In mathematics, BER = (number of bits received in error / total number of transmitted bits)  On the other hand, SNR ...

Hybrid Beamforming | Page 2

Beamforming Techniques Hybrid Beamforming... Page 1 | Page 2 | clear all; close all; clc; Nt = 64; Nr = 16; NtRF = 4; NrRF = 4; At both the transmitter and receiver ends, there are four RF chains only for a hybrid beamforming system. Alternatively, every 16 antenna elements on the transmitter side is connected to a single RF chain, while every 4 antenna elements on the receiver side are connected to a single RF chain. Mixers, amplifiers, and other critical wireless communication components make up the RF chain. Now, in the case of hybrid beamforming, there can be four different data streams between the transmitter and receiver, as both sides have four RF chains, each of which is accountable for a separate data stream. For Analog Beamforming: All 64 Tx antenna elements create a beam or focus the resultant correlated signal spread from adjacent antennas to a particular direction. Similarly, it may be used for beam...

Difference between AWGN and Rayleigh Fading

📘 Introduction, AWGN, and Rayleigh Fading 🧮 Simulator for the effect of AWGN and Rayleigh Fading on a BPSK Signal 🧮 MATLAB Codes 📚 Further Reading Wireless Signal Processing Gaussian and Rayleigh Distribution Difference between AWGN and Rayleigh Fading 1. Introduction Rayleigh fading coefficients and AWGN, or additive white gaussian noise [↗] , are two distinct factors that affect a wireless communication channel. In mathematics, we can express it in that way.  Fig: Rayleigh Fading due to multi-paths Let's explore wireless communication under two common noise scenarios: AWGN (Additive White Gaussian Noise) and Rayleigh fading. y = h*x + n ... (i) Symbol '*' represents convolution. The transmitted signal  x  is multiplied by the channel coefficient or channel impulse response (h)  in the equation above, and the symbol  "n"  stands for the white Gaussian noise that is added to the si...

Comparisons among ASK, PSK, and FSK | And the definitions of each

📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator for calculating Bandwidth of ASK, FSK, and PSK 🧮 MATLAB Code for BER vs. SNR Analysis of ASK, FSK, and PSK 📚 Further Reading 📂 View Other Topics on Comparisons among ASK, PSK, and FSK ... 🧮 Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. 🧮 MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK 🧮 Online Simulator for ASK, FSK, and PSK Generation 🧮 Online Simulator for ASK, FSK, and PSK Constellation 🧮 Some Questions and Answers Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK    Comparisons among ASK, PSK, and FSK   Simulator for Calculating Bandwidth of ASK, FSK, and PSK The baud rate represents the number of symbols transmitted per second. Both baud rate and bit rate a...

Comparing Baseband and Passband Implementations of ASK, FSK, and PSK

📘 Overview 🧮 Baseband and Passband Implementations of ASK, FSK, and PSK 🧮 Difference betwen baseband and passband 📚 Further Reading 📂 Other Topics on Baseband and Passband ... 🧮 Baseband modulation techniques 🧮 Passband modulation techniques   Baseband modulation techniques are methods used to encode information signals onto a baseband signal (a signal with frequencies close to zero), allowing for efficient transmission over a communication channel. These techniques are fundamental in various communication systems, including wired and wireless communication. Here are some common baseband modulation techniques: Amplitude Shift Keying (ASK) [↗] : In ASK, the amplitude of the baseband signal is varied to represent different symbols. Binary ASK (BASK) is a common implementation where two different amplitudes represent binary values (0 and 1). ASK is simple but susceptible to noise...

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSk, BPSK, ...

📘 Overview 🧮 Simulator for m-ary QAM and m-ary PSK 🧮 MATLAB Code for BPSK, M-ary PSK, and M-ary QAM Together 🧮 MATLAB Code for M-ary QAM 🧮 MATLAB Code for M-ary PSK 📚 Further Reading   MATLAB Script for  BER vs. SNR for M-QAM, M-PSK, QPSk, BPSK %Written by Salim Wireless %Visit www.salimwireless.com for study materials on wireless communication %or, if you want to learn how to code in MATLAB clc; clear; close all; % Parameters num_symbols = 1e5; % Number of symbols snr_db = -20:2:20; % Range of SNR values in dB % PSK and QAM orders to be tested psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; % Initialize BER arrays ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); % BER calculation for each PSK order and SNR value for i = 1:length(psk_orders) psk_order = psk_orders(i); for j = 1:length(snr_db) % Generate random symbols ...