Estimation in ISAC Systems

From Either-Or to Both

Radar and communications historically compete for the same spectrum, waveforms, and hardware. Integrated Sensing and Communications (ISAC) refuses to choose: one waveform, transmitted from a multi-antenna base station, should simultaneously convey data to a user and illuminate targets for the base station to estimate their range, angle, and velocity. This dual purpose introduces a fundamental tradeoff --- waveforms optimised for data (typically near-Gaussian, high-entropy) are good communicators but poor radars; waveforms optimised for sensing (deterministic, low-entropy, tightly designed) are good radars but poor communicators.

The right way to quantify this tradeoff is via the Cramer-Rao bound for the sensing parameters --- delay, angle, Doppler --- combined with the achievable communication rate. For each choice of transmit sample covariance Rs\mathbf{R}_s, one gets a rate R(Rs)R(\mathbf{R}_s) and a CRB matrix CRB(Rs)\mathrm{CRB}(\mathbf{R}_s). The boundary of the achievable (rate, sensing-accuracy) region is traced by sweeping a Lagrange multiplier that trades one against the other. CommIT and adjacent research groups have spent the last five years mapping this region for progressively richer system models; the section that follows presents the core tools.

Definition:

ISAC Signal Model (Monostatic)

A base station with NtN_t transmit antennas and NrN_r co-located receive antennas transmits TT samples of a matrix signal SCNt×T\mathbf{S} \in \mathbb{C}^{N_t\times T}. A single point target with complex reflectivity α\alpha, azimuth θ\theta, delay τ\tau, and Doppler ν\nu reflects the signal, producing received samples Y  =  αar(θ)at(θ)TDνSτ+W,\mathbf{Y} \;=\; \alpha\,\mathbf{a}_r(\theta)\mathbf{a}_t(\theta)^T\, \mathbf{D}_\nu \,\mathbf{S}_\tau + \mathbf{W}, where at,ar\mathbf{a}_t,\mathbf{a}_r are transmit/receive steering vectors, Dν\mathbf{D}_\nu is a diagonal Doppler phase ramp, Sτ\mathbf{S}_\tau is the delayed signal, and W\mathbf{W} is AWGN with per-entry variance σ2\sigma^2. Meanwhile, a communication user with a single antenna receives yt(c)  =  hHst+wt(c),y^{(c)}_t \;=\; \mathbf{h}^H \mathbf{s}_t + w^{(c)}_t, where h\mathbf{h} is the downlink channel to the user and st\mathbf{s}_t is the tt-th column of S\mathbf{S}.

The transmitter's only lever is the joint distribution of s1,,sT\mathbf{s}_1,\ldots,\mathbf{s}_T, parametrised by the sample covariance Rs=1TtE[ststH]\mathbf{R}_s = \tfrac{1}{T}\sum_t \mathbb{E}[\mathbf{s}_t \mathbf{s}_t^H], which controls both the beampattern and the communication rate.

Theorem: CRB for Joint Delay-Angle-Doppler

Under the ISAC signal model with a single target and AWGN, assuming α\alpha is known and we estimate η=(θ,τ,ν)\boldsymbol\eta = (\theta,\tau,\nu), the Fisher information matrix is J(η)  =  2α2σ2Re ⁣{ηGHηG},\mathbf{J}(\boldsymbol\eta) \;=\; \frac{2|\alpha|^2}{\sigma^2}\, \mathrm{Re}\!\left\{\partial_{\boldsymbol\eta}\mathbf{G}^H \partial_{\boldsymbol\eta}\mathbf{G}\right\}, where G(η)=ar(θ)at(θ)TDνSτ\mathbf{G}(\boldsymbol\eta) = \mathbf{a}_r(\theta)\mathbf{a}_t(\theta)^T \mathbf{D}_\nu \mathbf{S}_\tau is the noise-free received template, and η\partial_{\boldsymbol\eta} stacks the partial derivatives with respect to angle, delay, and Doppler. The CRB for each parameter is [J1]ii[\mathbf{J}^{-1}]_{ii}, and the overall sensing error is governed by tr(J1)\mathrm{tr}(\mathbf{J}^{-1}).

The FIM is quadratic in the transmit covariance Rs\mathbf{R}_s through the waveform S\mathbf{S}. Angular precision is driven by Rs\mathbf{R}_s aligning energy with at(θ)\mathbf{a}_t(\theta); delay precision by the effective bandwidth of the waveform; Doppler precision by the effective duration. The three parameters couple through the off-diagonal entries of the FIM.

Definition:

Rate-CRB Achievability Region

The rate-CRB region is the set of achievable (rate, sensing accuracy) pairs RISAC  =  Rs0,tr(Rs)Pt{(R,ϵ2):Rlog ⁣(1+hHRsh/σ2),\mathcal{R}_{\text{ISAC}} \;=\; \bigcup_{\mathbf{R}_s \succeq \mathbf{0},\, \mathrm{tr}(\mathbf{R}_s) \leq P_t} \Big\{(R, \epsilon^2) : R \leq \log\!\left(1 + \mathbf{h}^H\mathbf{R}_s\mathbf{h}/\sigma^2\right), ϵ2tr(J1(Rs))}.\quad \epsilon^2 \geq \mathrm{tr}(\mathbf{J}^{-1}(\mathbf{R}_s))\Big\}. The Pareto boundary of this region is traced by the beampattern design problem: maximise RR subject to a CRB constraint, or minimise tr(J1)\mathrm{tr}(\mathbf{J}^{-1}) subject to a rate constraint.

Both problems are non-convex in Rs\mathbf{R}_s because the FIM is quadratic in S\mathbf{S}, and hence the CRB involves the inverse of a matrix affine in Rs\mathbf{R}_s. Semidefinite-relaxation and majorisation-minimisation methods are the standard workhorses.

Theorem: Convex Relaxation of the Rate-CRB Problem

Consider the ISAC problem minRs0  tr(J1(Rs))\min_{\mathbf{R}_s \succeq \mathbf{0}} \;\mathrm{tr}(\mathbf{J}^{-1}(\mathbf{R}_s)) subject to hHRshγσ2\mathbf{h}^H \mathbf{R}_s \mathbf{h} \geq \gamma\, \sigma^2 (rate constraint) and tr(Rs)Pt\mathrm{tr}(\mathbf{R}_s) \leq P_t (power constraint). Introducing a slack matrix UJ1(Rs)\mathbf{U}\succeq \mathbf{J}^{-1}(\mathbf{R}_s) and applying the Schur complement, the problem becomes the SDP min  tr(U)s.t.[UIIJ(Rs)]0,  hHRshγσ2,  tr(Rs)Pt.\min \;\mathrm{tr}(\mathbf{U}) \quad \text{s.t.}\quad \begin{bmatrix} \mathbf{U} & \mathbf{I}\\ \mathbf{I} & \mathbf{J}(\mathbf{R}_s) \end{bmatrix} \succeq \mathbf{0},\; \mathbf{h}^H \mathbf{R}_s \mathbf{h} \geq \gamma\,\sigma^2,\; \mathrm{tr}(\mathbf{R}_s) \leq P_t. This is a convex program (the FIM is linear in Rs\mathbf{R}_s, the Schur condition is LMI, the rate and power constraints are linear).

The non-convex CRB minimisation is transformed into a convex SDP by the Schur complement. The price is an auxiliary variable U\mathbf{U} whose trace upper-bounds tr(J1)\mathrm{tr}(\mathbf{J}^{-1}). At the optimum, the Schur inequality holds with equality and tr(U)\mathrm{tr}(\mathbf{U}) equals the sensing cost.

,

SDP-Based ISAC Beampattern Design

Complexity: Each SDP solves in O((ntx^2)^3) time by interior-point; K SDPs total
Input: steering vectors a_t(theta), a_r(theta), delay/Doppler operators,
user channel h, noise variance sigma2, power budget P,
grid of rate constraints {R_1,...,R_K}
Output: achievable (rate, CRB) Pareto frontier
Compute partial derivatives of G w.r.t. theta, tau, nu at nominal point
Build the linear maps A_{ij}(R_s) = partial_i G^H partial_j G
for each R_k in rate grid:
Solve SDP:
minimize tr(U)
subject to [[U, I], [I, FIM(R_s)]] >= 0
h^H R_s h >= (exp(R_k) - 1) * sigma2
tr(R_s) <= P
R_s >= 0
record crb_k = tr(U) from the optimiser
(optional) extract transmit covariance, draw N samples s_t ~ CN(0, R_s)
return {(R_k, crb_k)}_{k=1..K}

In practice, warm-starting each SDP with the previous Rs\mathbf{R}_s dramatically accelerates the sweep. For very large arrays, ADMM or block-coordinate methods replace CVX/Mosek.

Rate-CRB Tradeoff Region

Sweep the Lagrange multiplier λ\lambda between a pure communication objective (λ=0\lambda=0, maximise rate, ignore sensing) and a pure sensing objective (λ=\lambda=\infty, minimise CRB, ignore rate). The resulting Pareto frontier shows the achievable (rate, angular-CRB) pairs for a 16-element ULA with a known user channel and target direction.

Parameters
16

Number of transmit antennas

10
5
45

Transmit Beampattern on the Pareto Frontier

At each Lagrange multiplier, the optimal Rs\mathbf{R}_s produces a specific transmit beampattern P(θ)=at(θ)HRsat(θ)P(\theta) = \mathbf{a}_t(\theta)^H \mathbf{R}_s \mathbf{a}_t(\theta). Move the slider to see how the pattern morphs from a pencil beam toward the user (pure communication) to a pattern with energy at both the user and the target (ISAC) to a pencil beam toward the target (pure sensing).

Parameters
16
0.5

Sensing/comm weight (0 = pure comm, 1 = pure sensing)

-20
30

Morphing from Comm-Only to Sensing-Only Beampattern

Animated sweep of the Lagrange weight λ\lambda from 00 (pure communication pencil beam toward the user) through 0.50.5 (ISAC pattern illuminating both user and target) to 11 (pure sensing pencil beam toward the target).
⚠️Engineering Note

ISAC in 6G Standardisation

ISAC is one of the 2024 ITU IMT-2030 usage scenarios for 6G. The standardisation problem is not waveform optimisation per se --- the air interface is largely inherited from 5G --- but resource allocation: which subcarriers, time slots, and spatial streams to dedicate to sensing, and under what QoS guarantees for the communication users. The rate-CRB region is the theoretical benchmark against which practical schedulers are measured, and the SDP formulation above is a building block of every published ISAC scheduler.

Practical Constraints
  • Per-user rate constraints turn the SDP into a multi-constraint problem

  • Clutter and multi-target scenarios require extending the FIM to block-diagonal form

  • Robust ISAC with imperfect CSI replaces fixed h\mathbf{h} with a worst-case set

📋 Ref: ITU-R Recommendation M.2160 (IMT-2030)
🎓CommIT Contribution(2023)

CRB-Rate Region for Multi-User ISAC

G. Caire, Y. Xiong, F. LiuIEEE Transactions on Information Theory

This CommIT contribution extends the single-user rate-CRB region to the multi-user broadcast channel. Each communication user demands a minimum rate, and the radar task is to estimate the parameters of multiple targets, each with its own reflectivity. The fundamental result is that the rate-CRB region has the same convex structure as the single-user case, but the Pareto boundary decomposes into an angular-allocation sub-problem (which users and which targets to serve simultaneously) and an SDP beampattern problem per angular group. The paper gives closed-form expressions for two-user, two-target cases and numerical sweeps for larger instances.

isacmulti-usercrb-rate-tradeoff

Common Mistake: Random vs. Deterministic Signal in the FIM

Mistake:

Treating S\mathbf{S} as deterministic in the FIM derivation while simultaneously treating it as random for the communication rate.

Correction:

The rate depends on Rs=E[ststH]\mathbf{R}_s = \mathbb{E}[\mathbf{s}_t \mathbf{s}_t^H], a second-moment quantity; the FIM, as written here, assumes deterministic S\mathbf{S} with empirical covariance Rs\mathbf{R}_s. A coherent treatment averages the FIM over the random S\mathbf{S} (expected FIM) --- the distinction matters when TT is small. For TT large enough that empirical and ensemble covariances agree, the two formulations coincide.

Historical Note: A Research Explosion

2018-present

ISAC as a distinct research agenda crystallised around 2018-2020 through a sequence of papers by Fan Liu, Christos Masouros, Bruno Clerckx and collaborators, who reframed long-standing radar-comm cohabitation questions as joint waveform design problems. The Caire-Kobayashi 2018 paper on information-theoretic ISAC fundamentals gave the field its canonical rate-distortion formulation. By 2024 the topic was one of the top-three submission tracks at IEEE Transactions on Signal Processing.

Quick Check

On the Pareto frontier of the rate-CRB region, what happens to the transmit beampattern as the rate constraint is loosened (i.e., as we allow lower rate in exchange for better sensing)?

It remains a pencil beam toward the user; only the power scales.

It shifts monotonically from the user toward the target.

It becomes isotropic.

It collapses to the target direction with zero side-lobe.

ISAC (Integrated Sensing and Communications)

A system paradigm in which a single waveform is transmitted from a multi-antenna transmitter to simultaneously convey data and probe the environment for sensing parameters (angle, delay, Doppler).

Related: Rate-CRB Achievability Region, Beampattern Design

Rate-CRB Region

The closure of achievable pairs (rate, sensing CRB) over all feasible transmit covariances in an ISAC system. Its Pareto boundary is traced by SDP-based beampattern design.

Related: ISAC Signal Model (Monostatic), Cramer (1946), Rao (1945), and a Near-Simultaneous Discovery