Signal Processing Foundations for Sensing

Radar Sensing Fundamentals

The fundamental task of radar sensing is to extract information about targets β€” their range (distance), radial velocity, and angle β€” from the echoes of a transmitted waveform. The key tool is the matched filter, which maximises the output signal-to-noise ratio (SNR) for a known waveform in additive white Gaussian noise.

Consider a transmitted baseband signal s(t)s(t) with bandwidth BB and duration TT. A point target at range RR with radial velocity vv produces an echo:

y(t)=α s(tβˆ’Ο„) ej2Ο€fdt+n(t)y(t) = \alpha \, s(t - \tau) \, e^{j2\pi f_d t} + n(t)

where Ξ±\alpha is the complex target reflectivity (incorporating path loss and radar cross section), Ο„=2R/c\tau = 2R/c is the round-trip delay, fd=2v/Ξ»f_d = 2v/\lambda is the Doppler frequency shift, and n(t)∼CN(0,Οƒ2)n(t) \sim \mathcal{CN}(0, \sigma^2) is additive white Gaussian noise. The factor of 2 in both Ο„\tau and fdf_d arises from the two-way (transmit and receive) propagation in monostatic radar.

The matched filter for delay estimation correlates the received signal with a time-shifted copy of the transmitted waveform:

z(Ο„^)=βˆ«βˆ’βˆžβˆžy(t) sβˆ—(tβˆ’Ο„^) dtz(\hat{\tau}) = \int_{-\infty}^{\infty} y(t) \, s^*(t - \hat{\tau}) \, dt

This correlation peaks at Ο„^=Ο„\hat{\tau} = \tau, yielding the target range R=cΟ„/2R = c\tau/2. The width of the correlation peak determines the range resolution.

Definition:

Ambiguity Function

The ambiguity function of a waveform s(t)s(t) is defined as:

Ο‡(Ο„,fd)=βˆ«βˆ’βˆžβˆžs(t) sβˆ—(tβˆ’Ο„) ej2Ο€fdt dt\chi(\tau, f_d) = \int_{-\infty}^{\infty} s(t) \, s^*(t - \tau) \, e^{j2\pi f_d t} \, dt

The squared magnitude βˆ£Ο‡(Ο„,fd)∣2|\chi(\tau, f_d)|^2 is called the ambiguity surface and characterises the joint range-Doppler resolution of the waveform. It satisfies the following properties:

  1. Peak normalisation: βˆ£Ο‡(0,0)∣2=Es2|\chi(0, 0)|^2 = E_s^{2}, where Es=∫∣s(t)∣2 dtE_s = \int |s(t)|^2 \, dt is the waveform energy.

  2. Symmetry: βˆ£Ο‡(Ο„,fd)∣=βˆ£Ο‡(βˆ’Ο„,βˆ’fd)∣|\chi(\tau, f_d)| = |\chi(-\tau, -f_d)|.

  3. Volume conservation (Moyal's identity): βˆ¬βˆ£Ο‡(Ο„,fd)∣2 dτ dfd=Es2\iint |\chi(\tau, f_d)|^2 \, d\tau \, df_d = E_s^{2} This fundamental constraint means that the ambiguity volume is fixed: narrowing the main lobe in one dimension necessarily raises sidelobes elsewhere. The ambiguity surface cannot be made arbitrarily narrow in both delay and Doppler simultaneously.

  4. Marginals: The zero-Doppler cut Ο‡(Ο„,0)\chi(\tau, 0) is the autocorrelation of s(t)s(t); the zero-delay cut Ο‡(0,fd)\chi(0, f_d) is related to the spectral autocorrelation.

The ambiguity function is the central tool for waveform design in radar signal processing. An ideal "thumbtack" ambiguity function β€” a sharp peak at the origin with uniformly low sidelobes β€” cannot exist due to Moyal's identity. In practice, waveform design aims to shape the ambiguity function to concentrate energy near the origin with acceptable sidelobe levels.

Definition:

Range-Doppler Map

A range-Doppler map is a two-dimensional representation of the radar scene obtained by evaluating the matched filter output across a grid of delay-Doppler hypotheses. For a sequence of MM received pulses (or OFDM symbols), the range-Doppler map is:

Z(p,q)=βˆ‘m=0Mβˆ’1βˆ‘k=0Kβˆ’1Ym,k eβˆ’j2Ο€kp/K ej2Ο€mq/M\mathbf{Z}(p, q) = \sum_{m=0}^{M-1} \sum_{k=0}^{K-1} Y_{m,k} \, e^{-j2\pi k p / K} \, e^{j2\pi m q / M}

where Ym,kY_{m,k} is the received signal on the kk-th subcarrier (or range bin) of the mm-th symbol (or pulse), pp indexes the range bins, and qq indexes the Doppler bins. The map displays peaks at the range-Doppler coordinates of each target. In matrix form, this is a 2D-FFT applied first across the subcarrier (fast-time) dimension for range processing, then across the symbol (slow-time) dimension for Doppler processing.

The range-Doppler map is the radar analogue of the time-frequency representation used in communication signal processing. Each bin in the map corresponds to a hypothesis about a specific delay-Doppler pair, and targets appear as peaks above the noise floor. Target detection is performed by applying a threshold (e.g., CFAR β€” constant false alarm rate) to the map.

Theorem: Radar Range and Velocity Resolution

For a waveform with bandwidth BB and coherent observation time TT:

  1. Range resolution: Ξ”r=c2B\Delta r = \frac{c}{2B} Two targets separated by less than Ξ”r\Delta r in range cannot be resolved. This follows from the width of the autocorrelation (zero-Doppler cut of the ambiguity function), which is approximately 1/B1/B in delay.

  2. Velocity resolution: Ξ”v=Ξ»2T\Delta v = \frac{\lambda}{2T} equivalently, the Doppler resolution is Ξ”fd=1/T\Delta f_d = 1/T. Two targets separated by less than Ξ”v\Delta v in radial velocity cannot be resolved.

  3. Maximum unambiguous range (for pulsed radar with PRI TpT_p): Rmax⁑=cTp2R_{\max} = \frac{c T_p}{2}

  4. Maximum unambiguous velocity (for pulsed radar with PRI TpT_p): vmax⁑=λ4Tpv_{\max} = \frac{\lambda}{4 T_p}

The time-bandwidth product BTBT characterises the waveform's degrees of freedom: a large BTBT enables simultaneously fine range and velocity resolution. The resolution cell volume in the range-velocity plane is Ξ”rβ‹…Ξ”v=cΞ»/(4BT)\Delta r \cdot \Delta v = c\lambda/(4BT), which shrinks as BTBT increases.

Example: Resolution of a 77 GHz Automotive Radar

A 77 GHz automotive radar has bandwidth B=150B = 150 MHz and a coherent processing interval of T=10T = 10 ms (corresponding to M=128M = 128 chirps at a PRI of Tp=78.125β€…β€ŠΞΌsT_p = 78.125\;\mu\text{s}). Compute the range resolution, velocity resolution, maximum unambiguous range and velocity, and the time-bandwidth product.

Ambiguity Function Geometry

Ambiguity Function Geometry
The ambiguity function βˆ£Ο‡(Ο„,fd)∣2|\chi(\tau, f_d)|^2 for three representative waveforms: (a) a linear frequency modulated (LFM) chirp showing the characteristic ridge tilted at 45Β°45Β° in the delay-Doppler plane, (b) an unweighted rectangular pulse with a "bed of nails" ambiguity along the Doppler axis, and (c) a random phase-coded waveform with a thumbtack-like shape (sharp main lobe, uniformly distributed sidelobes). Moyal's identity guarantees that the total volume under each surface is identical; only the distribution of energy differs.

Ambiguity Function: Bandwidth vs Resolution

Watch the ambiguity function main lobe narrow as the signal bandwidth sweeps from 10 MHz to 500 MHz. The narrower the main lobe, the finer the range resolution (Ξ”r=c/(2B)\Delta r = c/(2B)). This animation illustrates the fundamental bandwidth-resolution relationship.
As bandwidth increases, the zero-Doppler cut of the ambiguity function narrows, improving range resolution.

Ambiguity Function Surface

Explore the ambiguity function βˆ£Ο‡(Ο„,fd)∣2|\chi(\tau, f_d)|^2 for different waveform types. Select among three waveforms: (1) a linear frequency modulated (LFM) chirp, which produces a ridge-shaped ambiguity tilted in the delay-Doppler plane; (2) an OFDM waveform with random data, yielding a near-thumbtack shape with data-dependent sidelobes; and (3) a random phase-coded waveform with low autocorrelation sidelobes. Adjust the bandwidth and duration to observe how these parameters control the width of the main lobe in the delay and Doppler dimensions, respectively. Note that Moyal's identity ensures the total volume is conserved: narrowing the main lobe in one dimension spreads energy to sidelobes.

Parameters
1
100
10

Range-Doppler Processing via 2D-FFT

For a coherent processing interval (CPI) consisting of MM pulses (or OFDM symbols), each sampled at KK range bins, the received data can be arranged into a KΓ—MK \times M matrix Y\mathbf{Y} where the (k,m)(k, m) entry is the signal at the kk-th range sample of the mm-th pulse. A single point target at range bin p0p_0 with Doppler bin q0q_0 contributes:

Yk,m=α eβˆ’j2Ο€kp0/K ej2Ο€mq0/MY_{k,m} = \alpha \, e^{-j2\pi k p_0 / K} \, e^{j2\pi m q_0 / M}

The range-Doppler map is obtained by applying:

  1. An IFFT across rows (subcarrier/fast-time dimension) for range compression;
  2. An FFT across columns (symbol/slow-time dimension) for Doppler processing.

This 2D-FFT has computational complexity O(KMlog⁑(KM))O(KM \log(KM)) and produces a KΓ—MK \times M range-Doppler map where each target appears as a peak at coordinates (p0,q0)(p_0, q_0). The peak height is proportional to ∣α∣KM|\alpha| KM (coherent processing gain), while the noise floor is proportional to ΟƒKM\sigma\sqrt{KM}, yielding an output SNR gain of KMKM over a single sample.

Range-Doppler Map with Targets

Visualise the range-Doppler map produced by 2D-FFT processing of a simulated OFDM radar signal. Place multiple targets at random range-velocity positions and observe how they appear as peaks in the map. Adjust the number of targets to see how closely spaced targets interact. Vary the SNR to observe the transition from clear peaks to noise-limited detection. Increase the bandwidth to sharpen the range resolution (narrower peaks along the range axis).

Parameters
3
20
100

Ambiguity Function

The function Ο‡(Ο„,fd)=∫s(t)sβˆ—(tβˆ’Ο„)ej2Ο€fdtdt\chi(\tau, f_d) = \int s(t) s^*(t-\tau) e^{j2\pi f_d t} dt that characterises the joint range-Doppler resolution of a radar waveform. Its squared magnitude βˆ£Ο‡(Ο„,fd)∣2|\chi(\tau,f_d)|^2 is the ambiguity surface; Moyal's identity constrains its total volume to be constant.

Related: ISAC, Range-Doppler Map

Range-Doppler Map

A two-dimensional representation of the radar scene produced by applying a 2D-FFT to the received signal matrix. Each pixel corresponds to a range-velocity hypothesis; targets appear as peaks above the noise floor.

Related: Ambiguity Function

Common Mistake: Confusing Resolution with Maximum Unambiguous Range/Velocity

Mistake:

Assuming that finer range resolution (larger bandwidth) also increases the maximum unambiguous range.

Correction:

Range resolution Ξ”r=c/(2B)\Delta r = c/(2B) and maximum unambiguous range Rmax⁑=cTp/2R_{\max} = cT_p/2 depend on different parameters. Increasing bandwidth improves resolution but does not affect Rmax⁑R_{\max}, which depends on the pulse repetition interval TpT_p. Similarly, velocity resolution Ξ”v=Ξ»/(2T)\Delta v = \lambda/(2T) depends on the coherent processing interval, while vmax⁑=Ξ»/(4Tp)v_{\max} = \lambda/(4T_p) depends on the PRI.

In OFDM radar: Ξ”r\Delta r depends on total bandwidth B=KΞ”fB = K\Delta f, while Rmax⁑=c/(2Ξ”f)R_{\max} = c/(2\Delta f) depends on subcarrier spacing. Increasing KK (more subcarriers) improves resolution without affecting the unambiguous range.

Quick Check

A radar system operates at a carrier frequency of 28 GHz with a signal bandwidth of B=200B = 200 MHz and a coherent processing interval of T=5T = 5 ms. What are the range resolution Ξ”r\Delta r and the velocity resolution Ξ”v\Delta v?

Ξ”r=0.75\Delta r = 0.75 m, Ξ”v=1.07\Delta v = 1.07 m/s

Ξ”r=1.5\Delta r = 1.5 m, Ξ”v=0.54\Delta v = 0.54 m/s

Ξ”r=0.75\Delta r = 0.75 m, Ξ”v=2.14\Delta v = 2.14 m/s

Ξ”r=1.5\Delta r = 1.5 m, Ξ”v=1.07\Delta v = 1.07 m/s