Waveforms and the Ambiguity Function

Why Waveform Design Matters for Imaging

The transmitted waveform determines what the radar can resolve. A waveform with duration TT and bandwidth WW cannot simultaneously achieve arbitrarily fine range and Doppler resolution β€” there is a fundamental trade-off encoded in the ambiguity function.

For imaging, the ambiguity function is the point-spread function (PSF) of the matched-filter estimator. Its shape directly determines the sidelobe structure we will encounter in Chapter 13 when we study backpropagation imaging.

Definition:

The Ambiguity Function

For a transmitted waveform s(t)s(t) with unit energy (∫∣s(t)∣2 dt=1\int |s(t)|^2\,dt = 1), the ambiguity function is:

Ο‡A(Ο„,Ξ½)=βˆ«βˆ’βˆžβˆžs(t) sβˆ—(tβˆ’Ο„) ej2πνt dt,\chi_A(\tau, \nu) = \int_{-\infty}^{\infty} s(t)\,s^*(t - \tau)\,e^{j2\pi\nu t}\,dt,

where Ο„\tau is the delay (range) variable and Ξ½\nu is the Doppler (velocity) variable.

The squared magnitude βˆ£Ο‡A(Ο„,Ξ½)∣2|\chi_A(\tau, \nu)|^2 is the ambiguity diagram β€” it measures how well the radar can distinguish a target at delay Ο„\tau and Doppler Ξ½\nu from a target at the origin.

At the origin, βˆ£Ο‡A(0,0)∣2=1|\chi_A(0, 0)|^2 = 1 (perfect self-correlation). A target at (Ο„,Ξ½)(\tau, \nu) produces a response proportional to βˆ£Ο‡A(Ο„,Ξ½)∣2|\chi_A(\tau, \nu)|^2 in the matched-filter output β€” this is the range-Doppler sidelobe level.

Theorem: The Radar Uncertainty Principle

For any unit-energy waveform s(t)s(t), the total volume of the ambiguity function is constant:

βˆ«βˆ’βˆžβˆžβˆ«βˆ’βˆžβˆžβˆ£Ο‡A(Ο„,Ξ½)∣2 dτ dΞ½=1.\int_{-\infty}^{\infty}\int_{-\infty}^{\infty} |\chi_A(\tau, \nu)|^2\,d\tau\,d\nu = 1.

Consequently, it is impossible to design a waveform with βˆ£Ο‡A(Ο„,Ξ½)∣2β‰ˆ0|\chi_A(\tau, \nu)|^2 \approx 0 for all (Ο„,Ξ½)β‰ (0,0)(\tau, \nu) \neq (0,0). Compressing the ambiguity function in one dimension necessarily spreads it in another.

This is the radar analog of the Heisenberg uncertainty principle in quantum mechanics: energy is conserved and must appear somewhere in the delay-Doppler plane. A narrow mainlobe requires sidelobes elsewhere.

Definition:

Linear Frequency Modulation (LFM) Chirp

The LFM chirp (or linear FM, "chirp") waveform is:

s(t)=1T rect ⁣(tT)ejπμt2,s(t) = \frac{1}{\sqrt{T}}\,\text{rect}\!\left(\frac{t}{T}\right) e^{j\pi \mu t^2},

where TT is the pulse duration, ΞΌ=W/T\mu = W/T is the chirp rate, and WW is the swept bandwidth.

The instantaneous frequency is fi(t)=ΞΌtf_i(t) = \mu t, sweeping linearly from βˆ’W/2-W/2 to +W/2+W/2 over [βˆ’T/2,T/2][{-T/2}, {T/2}].

The time-bandwidth product WTW T determines the pulse compression ratio: a chirp with WT=1000W T = 1000 achieves 30 dB processing gain.

Theorem: Ambiguity Function of the LFM Chirp

The ambiguity function of an LFM chirp with bandwidth WW and duration TT is:

βˆ£Ο‡A(Ο„,Ξ½)βˆ£β‰ˆβˆ£sinc ⁣[(WΟ„βˆ’Ξ½T)(1βˆ’βˆ£Ο„βˆ£T)]∣(1βˆ’βˆ£Ο„βˆ£T),|\chi_A(\tau, \nu)| \approx \left|\text{sinc}\!\left[(W\tau - \nu T)\left(1 - \frac{|\tau|}{T}\right)\right]\right|\left(1 - \frac{|\tau|}{T}\right),

for βˆ£Ο„βˆ£<T|\tau| < T. The zero-delay cut gives βˆ£Ο‡A(0,Ξ½)∣=∣sinc(Ξ½T)∣|\chi_A(0, \nu)| = |\text{sinc}(\nu T)| (Doppler resolution 1/T1/T). The zero-Doppler cut gives βˆ£Ο‡A(Ο„,0)βˆ£β‰ˆβˆ£sinc(WΟ„)∣|\chi_A(\tau, 0)| \approx |\text{sinc}(W\tau)| (range resolution 1/W1/W, i.e., Ξ”R=c/(2W)\Delta R = c/(2W)).

The ambiguity diagram has a sheared ridge along the line Ξ½=ΞΌΟ„=(W/T)Ο„\nu = \mu\tau = (W/T)\tau, reflecting range-Doppler coupling.

The LFM trades a clean thumbtack ambiguity for a ridge-shaped one. The ridge means that a Doppler shift appears as a range shift in the matched filter output β€” this coupling must be corrected in radar processing.

Example: Designing an LFM Chirp for Automotive Radar

An automotive radar at f0=77f_0 = 77 GHz requires range resolution Ξ”R=0.15\Delta R = 0.15 m and unambiguous range Rua=200R_{\text{ua}} = 200 m. Design the LFM waveform parameters.

Definition:

Phase-Coded Waveforms

A phase-coded waveform consists of NN sub-pulses of duration TcT_c, each modulated by a phase Ο•n\phi_n:

s(t)=1Nβˆ‘n=0Nβˆ’1ejΟ•n rect ⁣(tβˆ’nTcTc).s(t) = \frac{1}{\sqrt{N}} \sum_{n=0}^{N-1} e^{j\phi_n}\,\text{rect}\!\left(\frac{t - nT_c}{T_c}\right).

The total duration is T=NTcT = NT_c and the bandwidth is Wβ‰ˆ1/TcW \approx 1/T_c, giving time-bandwidth product WT=NW T = N.

Barker codes achieve the minimum possible peak sidelobe level 1/N1/N among binary (Ο•n∈{0,Ο€}\phi_n \in \{0, \pi\}) codes. They exist only for N≀13N \leq 13.

Polyphase codes (Frank, Zadoff-Chu) allow larger NN with near-ideal ambiguity properties.

Ambiguity Function Explorer

Compare the ambiguity functions of LFM, Barker-13, and rectangular pulse waveforms. The LFM shows the characteristic sheared ridge; the Barker code approaches a thumbtack with low sidelobes; the rectangular pulse has a diamond shape with no range resolution.

Parameters
100

Ambiguity Function as the Point-Spread Function

The connection to imaging (Ch 07.2, Ch 13) is direct: the matched-filter image of a point scatterer is precisely βˆ£Ο‡A(Ο„,Ξ½)∣|\chi_A(\tau, \nu)| centered at the scatterer's delay and Doppler. The PSF of the sensing operator AHA\mathbf{A}^{H}\mathbf{A} in the range-Doppler plane is a product of the waveform's ambiguity function and the array's spatial beampattern.

This means that waveform design is not just a radar engineering concern β€” it directly determines the sidelobe structure of the imaging operator, which in turn determines whether sparse recovery algorithms (Ch 14) can separate closely spaced targets.

Comparison of Radar Waveform Types

PropertyRectangular PulseLFM ChirpBarker CodeZadoff-Chu
WTW T product1≫1\gg 1N≀13N \leq 13NN (arbitrary)
Range resolutioncT/2cT/2 (poor)c/(2W)c/(2W) (fine)cTc/2cT_c/2cTc/2cT_c/2
Doppler resolution1/T1/T (good)1/T1/T (good)1/(NTc)1/(NT_c)1/(NTc)1/(NT_c)
Range sidelobesN/Aβˆ’13.3-13.3 dB (unweighted)βˆ’22.3-22.3 dB (N=13N=13)β‰ˆβˆ’13\approx -13 dB
Ambiguity shapeDiamondSheared ridgeNear thumbtackNear thumbtack
Range-Doppler couplingNoneYes (Ξ½=ΞΌΟ„\nu = \mu\tau)MinimalMinimal
Doppler toleranceGoodGood (with correction)Poor for large NNModerate
Hardware complexityLowLow (simple chirp)Moderate (phase switching)Moderate

Common Mistake: Attempting to Design a Thumbtack Ambiguity Function

Mistake:

Trying to design a waveform with simultaneously narrow mainlobe and no sidelobes in both range and Doppler.

Correction:

The volume constraint βˆ¬βˆ£Ο‡A∣2=1\iint |\chi_A|^2 = 1 makes this impossible. A thumbtack-like ambiguity (e.g., Barker code) suppresses close-in sidelobes but necessarily has sidelobes at larger delay/Doppler offsets. Waveform design is about managing where the sidelobe energy goes, not eliminating it.

Historical Note: The Invention of Pulse Compression

1950s

Pulse compression via LFM was invented independently by several groups in the early 1950s. The key insight β€” that a long waveform with internal modulation can achieve the energy of a long pulse and the resolution of a short one β€” was first described in a classified 1951 patent by Sidney Darlington of Bell Labs. The name "chirp" comes from the similarity of the LFM waveform to a bird's call with rising pitch.

Barker codes were introduced by R.H. Barker in 1953 as binary sequences with optimal autocorrelation properties. Despite 70+ years of searching, no Barker code longer than 13 has been found, and it is widely conjectured that none exists.

Ambiguity Function

The cross-correlation of a waveform with a delayed-and-Doppler-shifted version of itself: Ο‡A(Ο„,Ξ½)=∫s(t)sβˆ—(tβˆ’Ο„)ej2πνt dt\chi_A(\tau, \nu) = \int s(t) s^*(t-\tau) e^{j2\pi\nu t}\,dt. Characterizes the joint range-Doppler resolution of the waveform.

Key Takeaway

The ambiguity function βˆ£Ο‡A(Ο„,Ξ½)∣2|\chi_A(\tau, \nu)|^2 is the fundamental tool for waveform design. Its volume is constant (the radar uncertainty principle), so the designer chooses where to place sidelobe energy. For imaging, the ambiguity function is the range-Doppler component of the PSF of AHA\mathbf{A}^{H}\mathbf{A}.

LFM Chirp Pulse Compression

An LFM (linear frequency modulated) chirp waveform is transmitted and passes through the matched filter. The animation shows the instantaneous frequency sweeping across the pulse bandwidth WW, followed by the matched filter output contracting into a narrow sinc-like peak with range resolution Ξ”R=c/(2W)\Delta R = c/(2W). Sidelobe structure and the effect of windowing are highlighted.