Waveform Design for ISAC

Waveform Design for Joint Radar and Communication

The choice of waveform is the most fundamental design decision in an ISAC system. The waveform must simultaneously provide good radar sensing performance (sharp ambiguity function, low sidelobes) and efficient communication (high data rate, low bit error rate). This section analyses two dominant waveform families β€” OFDM and FMCW β€” from a dual-function perspective and discusses joint optimisation criteria.

Orthogonal frequency-division multiplexing (OFDM) is the dominant modulation for modern communication systems (4G LTE, 5G NR, WiFi). Its multicarrier structure also makes it a natural radar waveform, with a particularly elegant connection between the subcarrier and symbol dimensions and the range and Doppler domains.

Consider an OFDM frame with KK subcarriers, subcarrier spacing Ξ”f\Delta f, and MM OFDM symbols. The total bandwidth is B=KΞ”fB = K \Delta f and the frame duration (excluding cyclic prefix overhead) is Tframe=MTsymT_{\text{frame}} = M T_{\text{sym}}, where Tsym=1/Ξ”f+TCPT_{\text{sym}} = 1/\Delta f + T_{\text{CP}} is the total OFDM symbol duration including the cyclic prefix TCPT_{\text{CP}}.

The transmitted baseband signal for the mm-th OFDM symbol is:

sm(t)=βˆ‘k=0Kβˆ’1Xm,k ej2Ο€kΞ”f(tβˆ’mTsym)s_m(t) = \sum_{k=0}^{K-1} X_{m,k} \, e^{j2\pi k \Delta f (t - m T_{\text{sym}})}

where Xm,kX_{m,k} is the data symbol on the kk-th subcarrier of the mm-th OFDM symbol. After reflection from a target at range RR with velocity vv, the received signal on subcarrier kk of symbol mm (after CP removal and FFT demodulation) is:

Ym,k=α Xm,k eβˆ’j2Ο€kΞ”fτ ej2Ο€mTsymfd+Nm,kY_{m,k} = \alpha \, X_{m,k} \, e^{-j2\pi k \Delta f \tau} \, e^{j2\pi m T_{\text{sym}} f_d} + N_{m,k}

where Ο„=2R/c\tau = 2R/c and fd=2v/Ξ»f_d = 2v/\lambda.

Definition:

OFDM Radar Processing

OFDM radar processing extracts target range and velocity from the phase progression of the received signal across subcarriers and OFDM symbols. After standard OFDM demodulation (CP removal and KK-point FFT), the received frequency-domain symbols are:

Ym,k=α Xm,k eβˆ’j2Ο€kΞ”fτ ej2Ο€mTsymfd+Nm,kY_{m,k} = \alpha \, X_{m,k} \, e^{-j2\pi k \Delta f \tau} \, e^{j2\pi m T_{\text{sym}} f_d} + N_{m,k}

Element-wise division by the known transmitted symbols Xm,kX_{m,k} removes the data modulation:

Y~m,k=Ym,kXm,k=α eβˆ’j2Ο€kΞ”fτ ej2Ο€mTsymfd+N~m,k\tilde{Y}_{m,k} = \frac{Y_{m,k}}{X_{m,k}} = \alpha \, e^{-j2\pi k \Delta f \tau} \, e^{j2\pi m T_{\text{sym}} f_d} + \tilde{N}_{m,k}

The compensated matrix Y~∈CMΓ—K\tilde{\mathbf{Y}} \in \mathbb{C}^{M \times K} contains pure range-Doppler phase information:

  • Along subcarriers (columns): The phase eβˆ’j2Ο€kΞ”fΟ„e^{-j2\pi k \Delta f \tau} is a complex sinusoid in kk with frequency proportional to the target delay Ο„\tau. An IFFT across kk yields the range profile.

  • Along symbols (rows): The phase ej2Ο€mTsymfde^{j2\pi m T_{\text{sym}} f_d} is a complex sinusoid in mm with frequency proportional to the Doppler shift fdf_d. An FFT across mm yields the Doppler profile.

The range resolution is Ξ”r=c/(2KΞ”f)=c/(2B)\Delta r = c/(2K\Delta f) = c/(2B) and the velocity resolution is Ξ”v=Ξ»/(2MTsym)\Delta v = \lambda/(2 M T_{\text{sym}}).

A critical advantage of OFDM radar is that the communication data symbols Xm,kX_{m,k} are known at the ISAC transmitter and can be divided out in the radar processing chain. This makes OFDM a natural DFRC waveform: the same signal carries data to the communication receiver and provides radar illumination, with no loss in sensing performance due to the data modulation.

2D-FFT Range-Doppler Processing for OFDM Radar

Input: Received frequency-domain matrix Y∈CMΓ—K\mathbf{Y} \in \mathbb{C}^{M \times K};
transmitted data matrix X∈CMΓ—K\mathbf{X} \in \mathbb{C}^{M \times K};
subcarrier spacing Ξ”f\Delta f; symbol duration TsymT_{\text{sym}};
carrier frequency f0f_0; detection threshold Ξ³\gamma
Output: Range-Doppler map Z\mathbf{Z}; detected target list {(Ri,vi)}\{(R_i, v_i)\}
1. Data removal (element-wise division):
For m=0,…,Mβˆ’1m = 0, \ldots, M-1 and k=0,…,Kβˆ’1k = 0, \ldots, K-1:
Y~m,k=Ym,k/Xm,k\tilde{Y}_{m,k} = Y_{m,k} / X_{m,k}
2. Range processing (IFFT along subcarrier dimension):
For each symbol m=0,…,Mβˆ’1m = 0, \ldots, M-1:
rm=IFFTK[Y~m,0,Y~m,1,…,Y~m,Kβˆ’1]\mathbf{r}_m = \mathrm{IFFT}_K[\tilde{Y}_{m,0}, \tilde{Y}_{m,1}, \ldots, \tilde{Y}_{m,K-1}]
Form range-compressed matrix R∈CMΓ—K\mathbf{R} \in \mathbb{C}^{M \times K}
3. Doppler processing (FFT along symbol dimension):
For each range bin p=0,…,Kβˆ’1p = 0, \ldots, K-1:
dp=FFTM[R0,p,R1,p,…,RMβˆ’1,p]\mathbf{d}_p = \mathrm{FFT}_M[R_{0,p}, R_{1,p}, \ldots, R_{M-1,p}]
Form range-Doppler map Z∈CKΓ—M\mathbf{Z} \in \mathbb{C}^{K \times M}
(equivalently, Z=FFTM[IFFTK[Y~T]]T\mathbf{Z} = \mathrm{FFT}_M[\mathrm{IFFT}_K[\tilde{\mathbf{Y}}^T]]^T)
4. Target detection (CFAR or fixed threshold):
For each bin (p,q)(p, q): if ∣Zp,q∣2>γ|Z_{p,q}|^2 > \gamma, report target at
range R=pβ‹…c/(2KΞ”f)R = p \cdot c / (2K\Delta f) and velocity
v=qβ‹…Ξ»/(2MTsym)v = q \cdot \lambda / (2 M T_{\text{sym}})
Complexity: O(KMlog⁑K+KMlog⁑M)=O(KMlog⁑(KM))O(KM\log K + KM\log M) = O(KM\log(KM)).

FMCW Chirp Radar

The frequency-modulated continuous wave (FMCW) waveform is a linear chirp that sweeps the instantaneous frequency linearly from f0βˆ’B/2f_0 - B/2 to f0+B/2f_0 + B/2 over the chirp duration TcT_c:

s(t)=ejπμt2,0≀t≀Tcs(t) = e^{j\pi \mu t^2}, \quad 0 \leq t \leq T_c

where ΞΌ=B/Tc\mu = B/T_c is the chirp rate (Hz/s). Upon reflection from a target at range RR with velocity vv, the received signal is mixed with the transmitted signal (de-chirping), producing a beat signal:

yb(t)=α ej2Ο€(fbt+fdt)y_b(t) = \alpha \, e^{j2\pi (f_b t + f_d t)}

where the beat frequency is:

fb=ΞΌΟ„=BTcβ‹…2Rcf_b = \mu \tau = \frac{B}{T_c} \cdot \frac{2R}{c}

and fd=2vf0/cf_d = 2v f_0 / c is the Doppler shift. The range is extracted from the beat frequency via FFT:

R=fbβ‹…cβ‹…Tc2B=fbβ‹…c2ΞΌR = \frac{f_b \cdot c \cdot T_c}{2B} = \frac{f_b \cdot c}{2\mu}

FMCW is the dominant waveform in automotive radar (77 GHz) due to its simple hardware implementation (a single mixer performs de-chirping) and its ability to achieve wide bandwidth with low-rate ADCs (the beat frequency is much lower than the signal bandwidth).

Criterion OFDM FMCW
Communication Native multicarrier modulation Requires additional modulation overlay
Range processing IFFT of subcarrier phase FFT of beat frequency
Doppler processing FFT across symbols FFT across chirps
Peak-to-average power ratio High (requires linear PA) Low (constant envelope)
Hardware complexity Standard comm. transceiver Simple mixer-based front-end
Waveform flexibility Adaptive subcarrier allocation Limited to chirp parameters
Sidelobe control Window functions, null subcarriers Window functions, chirp shaping

For ISAC systems, OFDM is generally preferred when the primary platform is a communication system (5G NR, WiFi) because the radar processing can be added as a software module with no hardware changes. FMCW is preferred when the primary platform is a radar system (automotive) and communication is a secondary function.

Joint Waveform Optimisation Metrics

When designing a waveform to serve both radar and communication, two fundamental metrics characterise the sensing performance:

Mutual information (MI) for sensing: The radar channel between the transmitter and the target can be modelled as a Gaussian channel with a frequency-dependent response determined by the target's scattering characteristics. The mutual information between the transmitted waveform and the target impulse response is:

Isense=∫0Blog⁑2 ⁣(1+∣S(f)∣2Οƒh2(f)Οƒn2)dfI_{\text{sense}} = \int_0^B \log_2\!\left(1 + \frac{|S(f)|^2 \sigma_h^2(f)}{\sigma_n^2}\right) df

where ∣S(f)∣2|S(f)|^2 is the power spectral density of the waveform, Οƒh2(f)\sigma_h^2(f) is the target's spectral response, and Οƒn2\sigma_n^2 is the noise power spectral density. Maximising IsenseI_{\text{sense}} leads to water-filling over the target spectrum.

Cram'{e}r-Rao bound (CRB) for parameter estimation: The CRB provides a lower bound on the variance of any unbiased estimator of the target parameters Ξ·=[Ο„,fd,ΞΈ]T\boldsymbol{\eta} = [\tau, f_d, \theta]^T. For delay estimation:

CRB(Ο„)=18Ο€2β‹…SNRβ‹…Ξ²2\mathrm{CRB}(\tau) = \frac{1}{8\pi^2 \cdot \text{SNR} \cdot \beta^2}

where β2=∫f2∣S(f)∣2df/∫∣S(f)∣2df\beta^2 = \int f^2 |S(f)|^2 df / \int |S(f)|^2 df is the effective (RMS) bandwidth of the waveform. Minimising the CRB favours waveforms that spread energy to the spectral edges.

Communication rate: The standard Shannon capacity applies:

Rcomm=∫0Blog⁑2 ⁣(1+∣S(f)∣2∣Hc(f)∣2Οƒn2)dfR_{\text{comm}} = \int_0^B \log_2\!\left(1 + \frac{|S(f)|^2 |H_c(f)|^2}{\sigma_n^2}\right) df

The joint optimisation problem is:

max⁑∣S(f)∣2β€…β€ŠRcomms.t.CRB(Ο„)≀ϡ,β€…β€Šβ€…β€Šβˆ«βˆ£S(f)∣2df≀P\max_{|S(f)|^2} \; R_{\text{comm}} \quad \text{s.t.} \quad \mathrm{CRB}(\tau) \leq \epsilon, \;\; \int |S(f)|^2 df \leq P

This formulation reveals a fundamental trade-off: the communication channel favours water-filling over ∣Hc(f)∣2|H_c(f)|^2 (allocating power to strong subcarriers), while the CRB favours maximising the RMS bandwidth (allocating power to edge subcarriers). The Pareto frontier of this trade-off characterises the achievable rate-accuracy region of the ISAC system.

Why OFDM Is a Natural ISAC Waveform

The elegance of OFDM for ISAC stems from three structural properties:

  1. Data transparency. The element-wise division Y~m,k=Ym,k/Xm,k\tilde{Y}_{m,k} = Y_{m,k}/X_{m,k} removes the random data modulation without affecting the sensing performance. This means that the same OFDM frame that carries user data also serves as a perfect radar waveform β€” there is no need to sacrifice subcarriers for radar-specific pilots.

  2. Resolution decoupling. Range resolution depends on the total bandwidth B=KΞ”fB = K\Delta f (subcarrier dimension), while velocity resolution depends on the frame duration T=MTsymT = MT_{\text{sym}} (symbol dimension). These can be tuned independently through numerology selection (Ξ”f\Delta f, KK, MM).

  3. Compatibility with existing infrastructure. Any 5G NR or WiFi device already has the OFDM modulator, FFT engine, and channel estimation pipeline. Adding radar functionality requires only software modifications for the 2D-FFT range-Doppler processing and target detection β€” no hardware changes are needed.

The main limitation of OFDM for radar is the high peak-to-average power ratio (PAPR), which reduces the efficiency of the power amplifier and limits the effective sensing range. PAPR reduction techniques (clipping, selected mapping, tone reservation) from the communication literature can be applied, but they may affect the ambiguity function sidelobe structure.

⚠️Engineering Note

PAPR Limitations of OFDM Radar

The high peak-to-average power ratio (PAPR) of OFDM is a critical practical constraint for ISAC systems. The power amplifier (PA) must operate with sufficient backoff to avoid nonlinear distortion, which reduces the effective isotropic radiated power (EIRP) and hence the maximum sensing range.

Typical PAPR values:

  • OFDM with random QAM data: PAPR β‰ˆ10\approx 10-12 dB for 99.9% complementary CDF
  • FMCW (constant envelope): PAPR = 0 dB

Impact on sensing range: The radar equation gives detection range Rmax⁑∝Pavg1/4R_{\max} \propto P_{\text{avg}}^{1/4}. A 6 dB PA backoff (typical for OFDM) reduces PavgP_{\text{avg}} by 6 dB, cutting the sensing range by a factor of 106/(4Γ—10)=1.4110^{6/(4 \times 10)} = 1.41 β€” a 29% range reduction compared to a constant-envelope waveform.

5G NR considerations: The 3GPP ISAC study item (Rel-19, TR 22.837) identifies PAPR as a key challenge. Proposed mitigations include:

  • DFT-spread OFDM (SC-FDMA) for uplink sensing: reduces PAPR to 3-4 dB
  • Dedicated sensing signals with low-PAPR sequences
  • Hybrid waveforms with OFDM for communication and FMCW for sensing
Practical Constraints
  • β€’

    OFDM PAPR: 10-12 dB (random data), 3-4 dB (DFT-spread)

  • β€’

    PA backoff reduces effective sensing range by ~30% vs constant-envelope

  • β€’

    Nonlinear distortion raises range-Doppler sidelobes

πŸ“‹ Ref: 3GPP TR 22.837 (ISAC study, Rel-19)
πŸ”§Engineering Note

5G NR Numerology for ISAC

The OFDM numerology of 5G NR directly determines the sensing parameters. The following table maps NR numerology to radar resolution for FR1 and FR2:

Parameter FR1 (sub-6 GHz) FR2 (mmWave)
Subcarrier spacing Ξ”f\Delta f 30 kHz 120 kHz
Max bandwidth 100 MHz 400 MHz
Range resolution c/(2B)c/(2B) 1.5 m 0.375 m
Max unamb. range c/(2Ξ”f)c/(2\Delta f) 5 km 1.25 km
Slot duration (14 symbols) 0.5 ms 0.125 ms
Velocity res. (1 slot) Ξ»/(2Γ—0.5ms)\lambda/(2 \times 0.5\text{ms}) Ξ»/(2Γ—0.125ms)\lambda/(2 \times 0.125\text{ms})
At f0=3.5f_0 = 3.5 GHz / 28 GHz 85.7 mm / 0.5 ms β†’\to 85.7 m/s 10.7 mm / 0.125 ms β†’\to 42.9 m/s

Key insight: FR2 provides excellent range resolution (37.5 cm with 400 MHz) but limited unambiguous range (1.25 km) due to the wide subcarrier spacing. FR1 has better unambiguous range (5 km) but coarser range resolution (1.5 m). Multi-slot processing extends the Doppler resolution at the cost of longer CPI.

For automotive ISAC at 77 GHz, the 4 GHz bandwidth achievable with dedicated radar waveforms provides 3.75 cm resolution β€” far beyond what 5G NR can achieve with its 400 MHz maximum.

Practical Constraints
  • β€’

    FR1 max BW: 100 MHz (range res. 1.5 m)

  • β€’

    FR2 max BW: 400 MHz (range res. 37.5 cm)

  • β€’

    Max unambiguous range limited by subcarrier spacing

πŸ“‹ Ref: 3GPP TS 38.104 (NR base station radio transmission and reception)

OFDM Radar Resolution Parameters

Explore how OFDM numerology parameters affect radar sensing performance. The subcarrier spacing Ξ”f\Delta f and number of subcarriers K=B/Ξ”fK = B/\Delta f determine the range resolution Ξ”r=c/(2B)\Delta r = c/(2B) and the maximum unambiguous range Rmax⁑=c/(2Ξ”f)R_{\max} = c/(2\Delta f). The number of OFDM symbols MM and the symbol duration Tsymβ‰ˆ1/Ξ”fT_{\text{sym}} \approx 1/\Delta f determine the velocity resolution Ξ”v=Ξ»/(2MTsym)\Delta v = \lambda/(2 M T_{\text{sym}}) and the maximum unambiguous velocity vmax⁑=λΔf/2v_{\max} = \lambda \Delta f / 2. Adjust the parameters to see the resulting range-velocity resolution cell and the unambiguous range-velocity window.

Parameters
120
100
14

Common Mistake: Division by Zero in OFDM Radar Data Removal

Mistake:

Performing element-wise division Y~m,k=Ym,k/Xm,k\tilde{Y}_{m,k} = Y_{m,k}/X_{m,k} without checking for null subcarriers where Xm,k=0X_{m,k} = 0.

Correction:

OFDM systems reserve some subcarriers as null (guard subcarriers at band edges, DC subcarrier). On these subcarriers, Xm,k=0X_{m,k} = 0 and division is undefined. In OFDM radar processing:

  • Exclude null subcarriers from the range-Doppler FFT, or
  • Set Y~m,k=0\tilde{Y}_{m,k} = 0 for null subcarriers (zero-padding)
  • Interpolate across null subcarriers from adjacent data

Null subcarrier gaps create spectral leakage and raise range sidelobes. In 5G NR, the DC subcarrier and ~10% guard subcarriers are null, causing range sidelobes approximately 13 dB below the peak without windowing.

OFDM Radar

A radar technique that uses an OFDM communication waveform for simultaneous sensing. Range information is extracted from the phase progression across subcarriers (IFFT), and Doppler information from the phase progression across OFDM symbols (FFT). The communication data is removed by element-wise division.

Related: ISAC, DFRC

Quick Check

In OFDM radar processing, after element-wise division by the known data symbols Xm,kX_{m,k}, the compensated received signal Y~m,k=Ξ±eβˆ’j2Ο€kΞ”fΟ„ej2Ο€mTsymfd+N~m,k\tilde{Y}_{m,k} = \alpha e^{-j2\pi k \Delta f \tau} e^{j2\pi m T_{\text{sym}} f_d} + \tilde{N}_{m,k}. Which FFT operation extracts the target range?

FFT across the symbol index mm (slow-time dimension)

IFFT across the subcarrier index kk (frequency dimension), because the range information is encoded in the phase progression across subcarriers

2D-FFT applied simultaneously across both dimensions

Correlation of the received signal with a reference chirp