Reading
Stories Mode

Applications — Radar, Sonar & Communications

~14 min read Lesson 4 of Module 10

Bringing It All Together

Throughout Module 10 we have built a complete theoretical toolkit: autocorrelation and cross-correlation for measuring signal similarity, the Neyman–Pearson detection framework for making binary decisions in noise, and the matched filter for extracting the maximum signal-to-noise ratio from a noisy observation. Now we deploy these tools in three of the most important application domains in electrical engineering: radar, sonar, and digital communications.

Each domain poses the same fundamental problem — detecting and characterizing a signal of known shape embedded in noise — but with different physical channels, different time and frequency scales, and different figures of merit. Seeing how the same mathematical machinery operates across domains deepens intuition and shows why the theory developed in this module is not academic abstraction but engineering bedrock.

Unifying Theme

All three systems transmit or receive a known waveform, apply a matched filter (or correlator), compare the output to a threshold, and read off time delay, frequency shift, or symbol identity. The underlying algorithm is identical; only the physical layer differs.

Radar: Range, Velocity, and Resolution

Radar (Radio Detection And Ranging) transmits a radio-frequency pulse and listens for the echo reflected from a target. The round-trip propagation delay τ determines the range R = cτ/2, where c is the speed of light. Detecting the echo is exactly a binary hypothesis test: H0 (no target, noise only) versus H1 (target present, signal + noise). The NP-optimal test statistic is the cross-correlation of the received waveform with the transmitted waveform — i.e., the matched filter output.

The matched filter does two things simultaneously in radar. First, it maximizes the peak SNR of the echo, enabling detection at the greatest possible range. Second, when the transmitted waveform is a chirp (LFM), the matched filter compresses the echo into a narrow spike whose width is approximately 1/B, where B is the swept bandwidth. This pulse compression decouples range resolution (controlled by B) from transmitted energy (controlled by pulse duration T), allowing modern radar to achieve both long detection range and fine range resolution.

Range Equation
R = \frac{c\,\tau}{2}
Round-trip delay τ gives range R. Factor of 2 accounts for two-way propagation.
Range Resolution
\Delta R = \frac{c}{2B}
After pulse compression, the range resolution is set by the signal bandwidth B, not the pulse length T.

A moving target introduces a Doppler frequency shift fd = 2v/λ, where v is the radial velocity and λ is the wavelength. This shifts the echo spectrum away from the transmitted spectrum, reducing the matched filter output. To measure velocity, modern radars use a bank of matched filters each tuned to a different Doppler shift, or equivalently, a 2-D correlation over both delay and Doppler — the ambiguity function, which is the cross-correlation of the transmitted signal with a Doppler-shifted copy of itself.

Ambiguity Function
\chi(\tau, f_d) = \int s(t)\, s^*(t-\tau)\, e^{j2\pi f_d t}\, dt
The ambiguity function χ(τ, fd) is the matched filter output as a function of time delay and Doppler shift. Its shape determines the radar’s joint range-velocity resolution.
c/2B
Range Resolution
2v/λ
Doppler Shift
10 log TB
Pulse Compression Gain (dB)

Sonar: Underwater Acoustics

Sonar (Sound Navigation And Ranging) is the acoustic analogue of radar. Instead of electromagnetic waves it uses sound waves propagating through water at approximately 1500 m/s — about 200,000 times slower than light, but with far lower attenuation at the frequencies used. The mathematics is identical: transmit a known waveform, apply a matched filter to the received echo, and compare the output to a threshold.

Active sonar transmits a ping and listens for the return. Passive sonar only listens, using cross-correlation between sensors to determine the direction and range of a sound source. In passive sonar, the “known signal” is replaced by the cross-correlation between two spatially separated hydrophones: a peak in the cross-correlation at lag τ indicates that the same wavefront arrived at the two sensors with a time difference of τ, giving the angle of arrival.

The key challenge in sonar is multipath propagation. Sound bounces off the ocean surface, the seafloor, and thermal layers, creating multiple delayed and attenuated copies of the transmitted signal at the receiver. The matched filter output shows multiple peaks — one for each propagation path — and the receiver must identify the direct-path return to avoid range errors. This is closely related to the channel equalization problem in digital communications.

Passive Sonar Principle

Two hydrophones separated by distance d receive the same wavefront at slightly different times. The cross-correlation peak at lag τ gives the time difference of arrival (TDOA): τ = d sinθ/c, where θ is the bearing angle and c is the speed of sound in water.

τTDOA = d sinθ / c  —  bearing from time delay

Digital Communications: Demodulation and Channel Estimation

In a digital communications system, the transmitter sends a stream of symbols, each represented by a pulse shape g[n]. The channel adds noise, and possibly multipath reflections. The receiver must decide which symbol was sent at each time instant, minimizing the probability of a wrong decision. This is again a hypothesis test: for binary signaling, H0 (bit 0, pulse −g[n]) versus H1 (bit 1, pulse +g[n]).

The matched filter for this problem has impulse response h[n] = g*[N−1−n], the time-reversed conjugate of the pulse shape. Sampling the output at the end of each symbol period gives the sufficient statistic for the optimal decision. For binary antipodal signaling in AWGN, the resulting bit error probability is:

Bit Error Probability (BER)
P_b = Q\!\left(\sqrt{\frac{2E_b}{N_0}}\right)
Q is the Q-function (tail probability of the standard normal). Eb/N0 is the energy-per-bit to noise-density ratio. The matched filter achieves this minimum BER.

When the channel introduces multipath (inter-symbol interference, ISI), the matched filter alone is insufficient. The received signal is the convolution of the transmitted sequence with the channel impulse response, and overlapping pulses cause each symbol to interfere with its neighbors. Modern receivers combine matched filtering with an equalizer (e.g., a minimum mean square error (MMSE) equalizer or a Viterbi decoder on a trellis) to remove ISI before making decisions.

Channel estimation uses cross-correlation directly. The transmitter sends a known pilot sequence p[n] whose autocorrelation is impulsive (like a PN sequence). The receiver correlates the received signal with the known pilot to estimate the channel impulse response h[n]: the cross-correlation peak at delay k gives the channel coefficient at tap k. This is exactly the matched-filter approach applied to the identification problem rather than the detection problem.

Pilot-Based Channel Estimation
\hat{h}[k] = \frac{1}{E_p}\sum_{n} r[n]\, p^*[n-k] \approx h[k] + \text{noise}
Correlating the received signal r[n] with the known pilot p[n] recovers the channel h[k]. Works well when the pilot autocorrelation Rpp[k] ≈ δ[k].

GPS Acquisition: Correlation in Two Dimensions

GPS acquisition is a striking case where correlation must be performed over two unknowns simultaneously: the code phase (which chip offset to use) and the carrier Doppler offset (the satellite is moving relative to the receiver). The receiver maintains a 2-D search grid and computes the correlation of the incoming signal with every combination of code phase and Doppler hypothesis.

Modern GPS receivers use the FFT to accelerate this search. For a fixed Doppler hypothesis, the correlation over all code phases is a circular cross-correlation and can be computed in O(N log N) operations using the FFT, instead of O(N²) for a direct search. Once the correct cell (code phase, Doppler) is identified by a peak exceeding the detection threshold, the receiver enters tracking mode, using phase-locked loops (PLLs) and delay-locked loops (DLLs) to maintain synchronization.

FFT-Based Acquisition

GPS acquisition correlates a 1023-chip PRN code over all possible code phases. Direct search: O(1023²) multiplications per Doppler bin. FFT-based: O(1023 log 1023) per Doppler bin — roughly 100× speedup. Modern SoCs run hundreds of Doppler bins in parallel, acquiring a lock in under a second.

Medical Ultrasound: Correlation Imaging

Medical ultrasound is a close cousin of sonar. A transducer emits a short acoustic pulse (typically 2–15 MHz) and records the echoes returning from tissue boundaries. The time-of-flight gives the depth of each reflector, and the amplitude of the matched-filter output gives the reflection coefficient. Scanning the beam electronically produces a 2-D cross-sectional image (B-mode scan) in real time.

Doppler ultrasound measures blood flow velocity by tracking the Doppler shift of echoes from moving red blood cells. This is exactly the radar Doppler problem scaled down to acoustic frequencies and centimeter ranges. The same ambiguity function trade-off applies: a long pulse resolves velocity well but provides poor axial (depth) resolution; a short pulse gives sharp depth resolution but poor velocity resolution.

The Correlation Engine: A Common Architecture

Stepping back, every application in this lesson shares an identical processing chain:

1. Transmit (or observe) a known waveform s[n].
2. Receive r[n] = s[n − τ] * h[n] + w[n] (delayed, convolved with channel, plus noise).
3. Correlate: compute Rrs[k] = Σ r[n] s*[n−k].
4. Threshold: declare detection if max|Rrs[k]| > γ.
5. Estimate: read off delay, Doppler, or symbol from the peak location.

This architecture is the matched filter / correlator receiver, and it is optimal in the Neyman–Pearson and minimum-BER senses. The only degrees of freedom are the choice of waveform s[n] — which shapes the ambiguity function and the autocorrelation sidelobes — and the threshold γ, which trades off Pfa against Pd.

Key Takeaways
  • Radar, sonar, GPS, communications, and medical ultrasound all solve the same detection/estimation problem using matched filtering and cross-correlation.
  • Radar range is determined by round-trip delay R = cτ/2; range resolution after pulse compression is ΔR = c/(2B).
  • Moving targets introduce a Doppler shift fd = 2v/λ. The ambiguity function describes joint delay-Doppler resolution for any waveform.
  • Passive sonar uses cross-correlation between two sensors to find the time difference of arrival (TDOA), giving the bearing angle of a sound source.
  • In digital communications, the matched filter achieves the minimum BER = Q(√(2Eb/N0)) for binary antipodal signaling in AWGN.
  • Channel estimation uses pilot sequences with impulsive autocorrelation. Correlating the received signal with the pilot directly recovers the channel impulse response.
  • GPS acquisition is a 2-D correlation search over code phase and Doppler. FFT-based methods reduce complexity from O(N²) to O(N log N) per Doppler bin.
  • All these systems share the same processing chain: transmit known waveform → correlate received signal → threshold → estimate delay/Doppler/symbol.
Previous Matched Filtering Module Overview Next Module Decimation (Downsampling)