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])
This is equivalent to a conjugate delay multiply. It’s a way to estimate and remove frequency offset by exploiting the phase difference between adjacent samples.
 
 y_dd[n]  = y_d[n] * conj(y_d[n-1])
 Another conjugate delay product — isolates residual modulation, often connected with root index  in ZC sequences.
 

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Ï€)
  

Conjugate Delay Product

The conjugate delay multiplication is defined as:

\( y_d[n] = y[n] \cdot \text{conj}(y[n-1]) \)

Expanding assuming \( y[n] = A e^{j\omega n} \):

\( y_d[n] = A e^{j\omega n} \cdot A e^{-j\omega(n-1)} = A^2 e^{j\omega} \)

Key Idea:

  • The result has a constant phase \( \omega \): this is the frequency.
  • This operation isolates the frequency difference between adjacent samples.

Second Conjugate Delay (Double Differencing)

Applying the same idea again:

\( y_{dd}[n] = y_d[n] \cdot \text{conj}(y_d[n-1]) \)

This Second Differencing:

  • Removes residual modulation — small leftover variations in phase.
  • The remaining structure is often linked to properties like the Zadoff-Chu sequence’s root index \( u \).

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







Contact Us

Name

Email *

Message *

Popular Posts

Constellation Diagrams of ASK, PSK, and FSK with MATLAB Code + Simulator

📘 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 📚 Further Reading 📂 Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...

Fading : Slow & Fast and Large & Small Scale Fading (with MATLAB Code + Simulator)

📘 Overview 📘 LARGE SCALE FADING 📘 SMALL SCALE FADING 📘 SLOW FADING 📘 FAST FADING 🧮 MATLAB Codes 📚 Further Reading LARGE SCALE FADING The term 'Large scale fading' is used to describe variations in received signal power over a long distance, usually just considering shadowing.  Assume that a transmitter (say, a cell tower) and a receiver  (say, your smartphone) are in constant communication. Take into account the fact that you are in a moving vehicle. An obstacle, such as a tall building, comes between your cell tower and your vehicle's line of sight (LOS) path. Then you'll notice a decline in the power of your received signal on the spectrogram. Large-scale fading is the term for this type of phenomenon. SMALL SCALE FADING  Small scale fading is a term that describes rapid fluctuations in the received signal power on a small time scale. This includes multipath propagation effects as well as movement-induced Doppler fr...

Online Simulator for ASK, FSK, and PSK

Try our new Digital Signal Processing Simulator!   Start Simulator for binary ASK Modulation Message Bits (e.g. 1,0,1,0) Carrier Frequency (Hz) Sampling Frequency (Hz) Run Simulation Simulator for binary FSK Modulation Input Bits (e.g. 1,0,1,0) Freq for '1' (Hz) Freq for '0' (Hz) Sampling Rate (Hz) Visualize FSK Signal Simulator for BPSK Modulation ...

Theoretical BER vs SNR for binary ASK, FSK, and PSK with MATLAB Code + Simulator

📘 Overview & Theory 🧮 MATLAB Codes 📚 Further Reading Theoretical BER vs SNR for Amplitude Shift Keying (ASK) The theoretical Bit Error Rate (BER) for binary ASK depends on how binary bits are mapped to signal amplitudes. For typical cases: If bits are mapped to 1 and -1, the BER is: BER = Q(√(2 × SNR)) If bits are mapped to 0 and 1, the BER becomes: BER = Q(√(SNR / 2)) Where: Q(x) is the Q-function: Q(x) = 0.5 × erfc(x / √2) SNR : Signal-to-Noise Ratio N₀ : Noise Power Spectral Density Understanding the Q-Function and BER for ASK Bit '0' transmits noise only Bit '1' transmits signal (1 + noise) Receiver decision threshold is 0.5 BER is given by: P b = Q(0.5 / σ) , where σ = √(N₀ / 2) Using SNR = (0.5)² / N₀, we get: BER = Q(√(SNR / 2)) Theoretical BER vs ...

DFTs-OFDM vs OFDM: Why DFT-Spread OFDM Reduces PAPR Effectively (with MATLAB Code)

DFT-spread OFDM (DFTs-OFDM) has lower Peak-to-Average Power Ratio (PAPR) because it "spreads" the data in the frequency domain before applying IFFT, making the time-domain signal behave more like a single-carrier signal rather than a multi-carrier one like OFDM. Deeper Explanation: Aspect OFDM DFTs-OFDM Signal Type Multi-carrier Single-carrier-like Process IFFT of QAM directly QAM → DFT → IFFT PAPR Level High (due to many carriers adding up constructively) Low (less fluctuation in amplitude) Why PAPR is High Subcarriers can add in phase, causing spikes DFT "pre-spreads" data, smoothing it Used in Wi-Fi, LTE downlink LTE uplink (as SC-FDMA) In OFDM, all subcarriers can...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

📘 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. BER = (number of bits received in error) / (total number of tran...

BER performance of QPSK with BPSK, 4-QAM, 16-QAM, 64-QAM, 256-QAM, etc (MATLAB + Simulator)

📘 Overview 📚 QPSK vs BPSK and QAM: A Comparison of Modulation Schemes in Wireless Communication 📚 Real-World Example 🧮 MATLAB Code 📚 Further Reading   QPSK provides twice the data rate compared to BPSK. However, the bit error rate (BER) is approximately the same as BPSK at low SNR values when gray coding is used. On the other hand, QPSK exhibits similar spectral efficiency to 4-QAM and 16-QAM under low SNR conditions. In very noisy channels, QPSK can sometimes achieve better spectral efficiency than 4-QAM or 16-QAM. In practical wireless communication scenarios, QPSK is commonly used along with QAM techniques, especially where adaptive modulation is applied. Modulation Bits/Symbol Points in Constellation Usage Notes BPSK 1 2 Very robust, used in weak signals QPSK 2 4 Balanced speed & reliability 4-QAM ...

Theoretical vs. simulated BER vs. SNR for ASK, FSK, and PSK (MATLAB Code + Simulator)

📘 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 for a given signal-to-noise ratio, typically measured in dB. Common noise types in wireless systems: 1. Additive White Gaussian Noise (AWGN) 2. Rayleigh Fading AWGN adds random noise; Rayleigh fading attenuates the signal variably. A good SNR helps reduce these effects. Simulator for calculating BER vs SNR for binary ASK, FSK, and PSK Calculate BER for Binary ASK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary FSK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary PSK Modulation Enter SNR (dB): Calculate BER BER vs. SNR Curves MATLAB Code for Theoretical BER % The code is written by SalimWireless.Com clc; clear; close all; % SNR va...