Model Mismatches and Robustness

Model Mismatches — When AtrueAassumed\mathbf{A}_{\text{true}} \neq \mathbf{A}_{\text{assumed}}

Every reconstruction algorithm assumes a specific forward model y=Ac+w\mathbf{y} = \mathbf{A}\mathbf{c} + \mathbf{w}. In practice, the true model always deviates:

y=(A+ΔA)c+weff,\mathbf{y} = (\mathbf{A} + \Delta\mathbf{A})\mathbf{c} + \mathbf{w}_{\text{eff}},

where ΔA\Delta\mathbf{A} captures model mismatch and weff\mathbf{w}_{\text{eff}} may include clutter and other non-Gaussian components. This section catalogs the dominant sources of mismatch and strategies for robustness.

Definition:

Born Approximation Breakdown

The Born approximation assumes single scattering: each point in the scene interacts only with the incident field, not with fields scattered by other points. This fails when:

  1. High contrast: χ(p)=εr11|\chi(\mathbf{p})| = |\varepsilon_r - 1| \gg 1 (e.g., metal objects, human body at RF).
  2. Electrically large objects: κa1\kappa a \gg 1, where aa is the object dimension.
  3. Dense scenes: Multiple scatterers interact (multiple scattering).

The Born validity criterion (approximate rule of thumb):

χκa<C,|\chi| \cdot \kappa a < C,

where C0.3C \approx 0.3-11 depending on geometry. When violated, the true received signal includes multiple-scattering terms that the linear model y=Ac+w\mathbf{y} = \mathbf{A}\mathbf{c} + \mathbf{w} cannot capture.

Alternative models for strong scatterers:

  • Rytov approximation: Replaces Born's additive error with a multiplicative one; better for smooth media.
  • Distorted Born iterative method (DBIM): Iteratively updates the background field to include multiple scattering.
  • Contrast source inversion (CSI): Jointly solves for the contrast and the total field.
,

Theorem: Reconstruction Error Under Model Mismatch

Let c\mathbf{c}^* be the true reflectivity and c^\hat{\mathbf{c}} be the Tikhonov-regularized estimate using the mismatched sensing matrix A\mathbf{A} (instead of the true A+ΔA\mathbf{A} + \Delta\mathbf{A}):

c^=(AHA+λI)1AHy.\hat{\mathbf{c}} = (\mathbf{A}^{H}\mathbf{A} + \lambda\mathbf{I})^{-1}\mathbf{A}^{H}\mathbf{y}.

The reconstruction error satisfies:

c^c2ΔA2λc2mismatch error+σ2λA21noise amplification+λc2regularization bias.\|\hat{\mathbf{c}} - \mathbf{c}^*\|_2 \leq \underbrace{\frac{\|\Delta\mathbf{A}\|_2}{\lambda} \cdot \|\mathbf{c}^*\|_2}_{\text{mismatch error}} + \underbrace{\frac{\sigma^2}{\lambda} \cdot \|\mathbf{A}\|_2^{-1}}_{\text{noise amplification}} + \underbrace{\lambda\,\|\mathbf{c}^*\|_2}_{\text{regularization bias}}.

The regularization parameter λ\lambda trades off noise amplification against mismatch and bias: larger λ\lambda reduces mismatch sensitivity but increases bias.

Over-regularization makes the reconstruction less sensitive to model errors (because the regularizer dominates over the data fidelity term), at the cost of a blurred, biased estimate. This is the fundamental tradeoff in model-mismatched imaging.

Definition:

Gridding Error (Basis Mismatch)

When targets do not lie on the assumed grid, the sensing matrix A\mathbf{A} does not perfectly model the measurements. The gridding error has two components:

  1. Spectral leakage: Energy from an off-grid target leaks to neighboring grid points, creating a sinc-like artifact pattern. The leakage energy is:

    Eleak1sinc2(δ/Δx),E_{\text{leak}} \approx 1 - \text{sinc}^2(\delta/\Delta x),

    where δ\delta is the offset from the nearest grid point and Δx\Delta x is the grid spacing.

  2. Phase mismatch: The off-grid target has a different round-trip phase than any grid point, creating a structured model error:

    ΔA(i,j,k),q=[A](i,j,k),q(ejκi,j,kTδp1).\Delta\mathbf{A}_{(i,j,k),q} = [\mathbf{A}]_{(i,j,k),q}\, (e^{-j\boldsymbol{\kappa}_{i,j,k}^T \delta\mathbf{p}} - 1).

For compressed sensing, the gridding error destroys the sparsity assumption. A truly sparse scene (few point targets) becomes non-sparse in the grid basis due to leakage.

Definition:

Calibration Errors

Calibration errors arise from imperfect knowledge of the system parameters used to construct A\mathbf{A}:

Error source Model effect Typical magnitude
Antenna position Phase error in steering vector <λ/20< \lambda/20 with GPS/INS
Mutual coupling Effective pattern distortion 15-15 to 25-25 dB (element-dependent)
Phase noise Random phase on each measurement 80-80 to 100-100 dBc/Hz at 1 MHz
I/Q imbalance Image at conjugate frequency 1-5% gain, 1-5 degrees phase
Frequency offset Range shift PPM-level with modern oscillators

These errors are modeled as a perturbation to the sensing matrix:

Aactual=Aideal+ΔAcal,\mathbf{A}_{\text{actual}} = \mathbf{A}_{\text{ideal}} + \Delta\mathbf{A}_{\text{cal}},

and reconstruction quality depends on ΔAcal/A\|\Delta\mathbf{A}_{\text{cal}}\| / \|\mathbf{A}\|.

Example: Born Approximation Validity for a Dielectric Cylinder

A dielectric cylinder of radius a=λ/4a = \lambda/4 and contrast χ=εr1\chi = \varepsilon_r - 1 is illuminated at frequency f0f_0.

(a) Compute the Born validity parameter χκa|\chi| \cdot \kappa a for εr{1.5,3,6,10}\varepsilon_r \in \{1.5, 3, 6, 10\}.

(b) Using the rule of thumb χκa<0.5|\chi|\kappa a < 0.5, determine which materials are within the Born regime.

(c) For concrete (εr=6\varepsilon_r = 6), at what maximum object size aa does the Born approximation remain valid at f0=5f_0 = 5 GHz?

,

Born Approximation Validity

Explore the Born approximation error as a function of contrast strength χ|\chi| and electrical size κa\kappa a. The color map shows the relative error between the Born approximation and a reference multi-scattering solution.

Parameters
5
3

Reconstruction Robustness vs. Model Mismatch

Compare the reconstruction error of back-projection, Tikhonov, and 1\ell_1 minimization as the model mismatch ΔAF/AF\|\Delta\mathbf{A}\|_F / \|\mathbf{A}\|_F increases. Over-regularized Tikhonov trades resolution for robustness.

Parameters
0.1
1

Strategies for Robustness to Model Mismatch

  1. Over-regularization: Increase λ\lambda beyond the noise-optimal value. Reduces sensitivity to ΔA\Delta\mathbf{A} at the cost of resolution (bias).

  2. Robust optimization: Minimax formulations that optimize for the worst-case ΔA\Delta\mathbf{A} within a bounded set: mincmaxΔAϵy(A+ΔA)c2\min_{\mathbf{c}} \max_{\|\Delta\mathbf{A}\| \leq \epsilon} \|\mathbf{y} - (\mathbf{A} + \Delta\mathbf{A})\mathbf{c}\|^2.

  3. Joint estimation: Simultaneously estimate c\mathbf{c} and ΔA\Delta\mathbf{A} (auto-calibration, autofocus in SAR).

  4. Learned robustness: Train neural networks on data with realistic model mismatches (Part IV). The network implicitly learns to be robust.

Understanding the magnitude and structure of model mismatch is essential for choosing reconstruction algorithms — a theme that pervades Part III.

🎓CommIT Contribution(2026)

Caire's Unified Illumination and Sensing Model

G. CaireCommIT Research Note

Caire's unified framework derives the forward model y=Ac+w\mathbf{y} = \mathbf{A}\mathbf{c} + \mathbf{w} from first principles (diffraction tomography), connecting the illumination pattern to the sensing operator through wavenumber-domain analysis. The derivation in this chapter (Sections s01-s02) follows Caire's framework, which provides the Kronecker-structured operator, the far-field Taylor approximation, and the link budget normalization that unifies diverse sensing geometries under a single formalism.

rf-imagingforward-modelsensing-operator
⚠️Engineering Note

Calibration — Closing the Model-Reality Gap

Calibration measures and corrects hardware impairments so that AactualAideal\mathbf{A}_{\text{actual}} \approx \mathbf{A}_{\text{ideal}}:

Calibration type What it corrects Method
Antenna calibration Pattern errors, coupling Anechoic chamber measurement
Channel calibration Gain/phase across channels Known target (corner reflector)
Motion calibration Platform trajectory errors INS + GPS + autofocus
System calibration End-to-end transfer function Loop-back or reference target

For MIMO radar arrays, online calibration techniques estimate the coupling matrix and channel gain/phase errors from the radar data itself, analogous to self-calibration in radio astronomy.

Practical Constraints
  • Anechoic chamber measurements are expensive and one-time

  • Online calibration requires known reference targets in the scene

  • Autofocus adds computational cost to reconstruction

Common Mistake: The Inverse Crime

Mistake:

Testing a reconstruction algorithm using synthetic data generated from the same forward model used for reconstruction (same grid, same A\mathbf{A}, same discretization).

Correction:

Always generate test data from a finer grid or a different forward model than the one used for reconstruction. Otherwise, the results are overly optimistic because the model mismatch that dominates real-world performance is absent. This is called the "inverse crime" in the inverse problems literature.

Quick Check

Increasing the regularization parameter λ\lambda in Tikhonov regularization makes the reconstruction:

More sensitive to model mismatch but sharper

Less sensitive to model mismatch but more blurred

Both sharper and more robust

Key Takeaway

Model mismatch is the dominant practical challenge in RF imaging. Born approximation breakdown, gridding errors, and calibration imperfections all introduce structured perturbations ΔA\Delta\mathbf{A}. Over-regularization trades resolution for robustness. The inverse crime — testing with the same model used for reconstruction — hides the true impact of mismatch.