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







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

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 Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude Frequency ...

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 ...

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 ...

Constellation Diagrams of ASK, PSK, and FSK

📘 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...

MATLAB Code for ASK, FSK, and PSK

📘 Overview & Theory 🧮 MATLAB Code for ASK 🧮 MATLAB Code for FSK 🧮 MATLAB Code for PSK 🧮 Simulator for binary ASK, FSK, and PSK Modulations 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for ASK Modulation and Demodulation % The code is written by SalimWireless.Com % Clear previous data and plots clc; clear all; close all; % Parameters Tb = 1; % Bit duration (s) fc = 10; % Carrier frequency (Hz) N_bits = 10; % Number of bits Fs = 100 * fc; % Sampling frequency (ensure at least 2*fc, more for better representation) Ts = 1/Fs; % Sampling interval samples_per_bit = Fs * Tb; % Number of samples per bit duration % Generate random binary data rng(10); % Set random seed for reproducibility binary_data = randi([0, 1], 1, N_bits); % Generate random binary data (0 or 1) % Initialize arrays for continuous signals t_overall = 0:Ts:(N_bits...

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

🧮 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 clc; clear; close all; num_symbols = 1e5; snr_db = -20:2:20; psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); for i = 1:length(psk_orders) psk_order = psk_orders(i); for j = 1:length(snr_db) data_symbols = randi([0, psk_order-1], 1, num_symbols); modulated_signal = pskmod(data_symbols, psk_order, pi/psk_order); received_signal = awgn(modulated_signal, snr_db(j), 'measured'); demodulated_symbols = pskdemod(received_signal, psk_order, pi/psk_order); ber_psk_results(i, j) = sum(data_symbols ~= demodulated_symbols) / num_symbols; end end for ...

MATLAB code for GMSK

📘 Overview & Theory 🧮 MATLAB Codes for GMSK 🧮 Online Simulator for GMSK 🧮 Simulation Results for GMSK 📚 Further Reading   Copy the MATLAB code from here  % The code is developed by SalimWireless.com clc; clear; close all; % Parameters samples_per_bit = 36; bit_duration = 1; num_bits = 20; sample_interval = bit_duration / samples_per_bit; time_vector = 0:sample_interval:(num_bits * bit_duration); time_vector(end) = []; % Generate and modulate binary data binary_data = randi([0, 1], 1, num_bits); modulated_bits = 2 * binary_data - 1; upsampled_signal = kron(modulated_bits, ones(1, samples_per_bit)); figure; plot(time_vector, upsampled_signal); title('Message Signal'); % Apply Gaussian filter filtered_signal = conv(GMSK_gaussian_filter1(bit_duration, samples_per_bit), upsampled_signal); filtered_signal = [filtered_signal, filtered_signal(end)]; figure; plot(filtered_signal); title('Filtered Signal'); % Integration ...

Simulation of ASK, FSK, and PSK using MATLAB Simulink

📘 Overview 🧮 How to use MATLAB Simulink 🧮 Simulation of ASK using MATLAB Simulink 🧮 Simulation of FSK using MATLAB Simulink 🧮 Simulation of PSK using MATLAB Simulink 🧮 Simulator for ASK, FSK, and PSK 🧮 Digital Signal Processing Simulator 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Simulation Simulation of Amplitude Shift Keying (ASK) using MATLAB Simulink      In Simulink, we pick different components/elements from MATLAB Simulink Library. Then we connect the components and perform a particular operation.  Result A sine wave source, a pulse generator, a product block, a mux, and a scope are shown in the diagram above. The pulse generator generates the '1' and '0' bit sequences. Sine wave sources produce a specific amplitude and frequency. The scope displays the modulated signal as well as the original bit sequence created by the pulse generator. Mux is a tool for displaying b...