From Finite to Infinite Dimensions

Why Infinite-Dimensional Spaces for RF Imaging?

The fundamental task in RF imaging is reconstruction: given noisy measurements

y=Ac+w,\mathbf{y} = \mathbf{A}\,\mathbf{c} + \mathbf{w},

we seek to recover the scene c(r)c(\mathbf{r}) — a reflectivity map, a permittivity profile, or a scattering distribution defined at every point rΩ\mathbf{r} \in \Omega.

This is fundamentally different from the finite-dimensional linear algebra of Telecom Ch~1, where all signals live in Cn\mathbb{C}^n. A scene function c(r)c(\mathbf{r}) assigns a complex value to every point in a continuous domain — it lives in an infinite-dimensional space. To judge the quality of an estimate c^\hat{c}, we need c^c\|\hat{c} - c\|, and the choice of norm is far from cosmetic:

  • Minimising L2\|\cdot\|_{L^2} yields smooth, energy-optimal reconstructions.
  • Minimising L1\|\cdot\|_{L^1} promotes sparse, edge-preserving solutions.
  • Sobolev norms Hs\|\cdot\|_{H^s} encode prior smoothness.

The mathematical framework of this section — normed spaces, completeness, Banach spaces — is the bedrock on which every reconstruction algorithm rests.

Historical Note: Stefan Banach and the Birth of Functional Analysis

1920–1932

Stefan Banach (1892–1945), working in Lwów (then Poland), created the theory of complete normed spaces in his 1920 doctoral thesis and the landmark 1932 monograph Théorie des opérations linéaires. The monograph, written in French to reach the widest European audience, collected the Hahn–Banach theorem, the open mapping theorem, and the uniform boundedness principle — results developed in parallel by Banach, Hahn, and Steinhaus — into a unified theory.

Banach's spaces provided the language in which John von Neumann simultaneously developed quantum mechanics (Hilbert spaces) and in which Norbert Wiener grounded signal processing. The LpL^p spaces had been defined by Riesz in 1910; Banach's insight was that completeness — not just a norm — was the essential structural property.

Definition:

Normed Space

A normed space (V,)(\mathcal{V}, \|\cdot\|) is a vector space V\mathcal{V} over a field F{R,C}\mathbb{F} \in \{\mathbb{R}, \mathbb{C}\} together with a function  ⁣:V[0,)\|\cdot\| \colon \mathcal{V} \to [0,\infty) — called a norm — satisfying for all f,gVf, g \in \mathcal{V} and αF\alpha \in \mathbb{F}:

# Axiom Statement
N1 Positive definiteness f=0    f=0\|f\| = 0 \iff f = 0
N2 Absolute homogeneity αf=αf\|\alpha f\| = |\alpha|\,\|f\|
N3 Triangle inequality f+gf+g\|f + g\| \leq \|f\| + \|g\|

Every norm induces a metric d(f,g)=fgd(f, g) = \|f - g\|, so a normed space is automatically a metric space and topological notions (convergence, continuity, open/closed sets) are inherited.

The nontrivial content of N1 is the converse direction: f=0\|f\| = 0 implies f=0f = 0. This rules out "pseudo-norms" that vanish on nonzero elements.

Definition:

LpL^p Spaces

Let ΩRd\Omega \subseteq \mathbb{R}^d be a measurable set. For 1p<1 \leq p < \infty, the Lebesgue space Lp(Ω)L^p(\Omega) consists of all equivalence classes of measurable functions f ⁣:ΩCf \colon \Omega \to \mathbb{C} for which the following norm is finite:

fLp=(Ωf(r)pdr)1/p.\|f\|_{L^p} = \Bigl(\int_\Omega |f(\mathbf{r})|^p \, d\mathbf{r}\Bigr)^{1/p}.

For p=p = \infty:

fL=ess suprΩf(r),\|f\|_{L^\infty} = \operatorname{ess\,sup}_{\mathbf{r} \in \Omega}\, |f(\mathbf{r})|,

and L(Ω)L^\infty(\Omega) consists of all measurable functions with fL<\|f\|_{L^\infty} < \infty (essentially bounded functions).

Two functions are identified if they differ only on a set of measure zero: f=gf = g a.e.

L2(Ω)L^2(\Omega) occupies a privileged position: it is the only LpL^p space (for p2p \neq 2) that is also a Hilbert space (see DHilbert Space). It is the natural home for finite-energy signals and the setting where least-squares reconstruction is most naturally formulated.

Example: LpL^p Norms of a Rectangular Pulse

Let f ⁣:R2Rf \colon \mathbb{R}^2 \to \mathbb{R} be the indicator function of the unit square:

f(r)=1[0,1]2(r)={1r[0,1]2,0otherwise.f(\mathbf{r}) = \mathbf{1}_{[0,1]^2}(\mathbf{r}) = \begin{cases} 1 & \mathbf{r} \in [0,1]^2, \\ 0 & \text{otherwise}. \end{cases}

Compute fL1\|f\|_{L^1}, fL2\|f\|_{L^2}, and fL\|f\|_{L^\infty}.

Definition:

Sequence Spaces p\ell^p

The sequence space p\ell^p is the discrete analogue of LpL^p:

p={x=(x1,x2,):xp=(k=1xkp)1/p<},1p<.\ell^p = \Bigl\{\mathbf{x} = (x_1, x_2, \ldots) : \|\mathbf{x}\|_{\ell^p} = \Bigl(\sum_{k=1}^{\infty} |x_k|^p\Bigr)^{1/p} < \infty\Bigr\}, \quad 1 \leq p < \infty.

For p=p = \infty: ={x:supkxk<}\ell^\infty = \{\mathbf{x} : \sup_k |x_k| < \infty\}.

For finite sequences of length NN, all p\ell^p norms are well-defined and the space is CN\mathbb{C}^N with the p\|\cdot\|_p norm.

When we discretise a scene on a grid of NN voxels, the infinite-dimensional problem in Lp(Ω)L^p(\Omega) becomes a finite-dimensional problem in CN\mathbb{C}^N with the p\ell^p norm. This bridge between continuous theory and discrete computation is central to every practical imaging algorithm.

Definition:

Cauchy Sequence

A sequence {fn}n=1\{f_n\}_{n=1}^\infty in a normed space (V,)(\mathcal{V}, \|\cdot\|) is a Cauchy sequence if for every ε>0\varepsilon > 0 there exists NNN \in \mathbb{N} such that

fmfn<εfor all m,nN.\|f_m - f_n\| < \varepsilon \qquad \text{for all } m, n \geq N.

Every convergent sequence is Cauchy. The converse need not hold in an arbitrary normed space (the rationals Q\mathbb{Q} provide a counterexample).

Definition:

Banach Space

A normed space (V,)(\mathcal{V}, \|\cdot\|) is a Banach space if every Cauchy sequence in V\mathcal{V} converges to an element of V\mathcal{V}:

{fn} Cauchy in V    fV:fnf0.\{f_n\} \text{ Cauchy in } \mathcal{V} \implies \exists\, f \in \mathcal{V} : \|f_n - f\| \to 0.

Completeness is not a luxury — it is essential for the convergence guarantees of iterative reconstruction algorithms. Banach's fixed-point theorem (the cornerstone of ISTA, gradient descent, and proximal algorithms) requires completeness as a hypothesis. Without it, an iterative sequence could satisfy fmfn0\|f_{m} - f_{n}\| \to 0 yet have no limit within the space.

,

Theorem: LpL^p Spaces are Banach Spaces

For 1p1 \leq p \leq \infty, the space Lp(Ω)L^p(\Omega) equipped with the Lp\|\cdot\|_{L^p} norm is a Banach space (the Riesz–Fischer theorem).

Extract a "rapidly convergent" subsequence from a Cauchy sequence, build a candidate limit via dominated or monotone convergence, then show the entire original sequence converges to that limit. The argument exploits the completeness of R\mathbb{R} pointwise (a.e.) and then promotes this to LpL^p convergence.

LpL^p Unit Balls in R2\mathbb{R}^2

The unit ball {xR2:xp1}\{\mathbf{x} \in \mathbb{R}^2 : \|\mathbf{x}\|_p \leq 1\} changes shape dramatically with pp. At p=1p = 1: a diamond whose corners on the axes explain why 1\ell^1 minimisation promotes sparsity. At p=2p = 2: the familiar Euclidean circle. As pp \to \infty: a square [1,1]2[-1,1]^2.

For p<1p < 1 the "ball" is no longer convex — p\|\cdot\|_p fails the triangle inequality and is only a quasi-norm. Observe how the pinching makes the coordinate axes even stronger attractors, linking sub-1\ell^1 norms to compressed sensing.

Parameters
2

Norm order

Definition:

Equivalence of Norms

Two norms a\|\cdot\|_a and b\|\cdot\|_b on a vector space V\mathcal{V} are equivalent if there exist constants c,C>0c, C > 0 such that

cfafbCfafV.c\,\|f\|_a \leq \|f\|_b \leq C\,\|f\|_a \qquad \forall\, f \in \mathcal{V}.

Equivalent norms define the same topology (same convergent sequences, same continuous maps, same open/closed sets).

Theorem: All Norms on Finite-Dimensional Spaces are Equivalent

Let V\mathcal{V} be a finite-dimensional vector space over R\mathbb{R} or C\mathbb{C}. Then any two norms on V\mathcal{V} are equivalent.

The unit sphere of any norm on Cn\mathbb{C}^n is compact (Heine–Borel). A continuous function on a compact set attains its minimum (which is positive, since the norm is positive definite) and maximum, supplying the equivalence constants.

Common Mistake: Norms Are NOT Equivalent in Infinite Dimensions

Mistake:

Assuming that convergence in one LpL^p norm implies convergence in another, as holds in Cn\mathbb{C}^n.

Correction:

In infinite-dimensional spaces different LpL^p norms give genuinely different notions of convergence.

Example. On Ω=[0,1]\Omega = [0,1], let fn(x)=n1[0,1/n](x)f_n(x) = n \cdot \mathbf{1}_{[0,\,1/n]}(x) (a tall thin spike).

  • fnL1=nn1=1\|f_n\|_{L^1} = n \cdot n^{-1} = 1 for all nn — does NOT converge to 00.
  • fnL2=n\|f_n\|_{L^2} = \sqrt{n} \to \infty — diverges.
  • fn(x)0f_n(x) \to 0 pointwise for x>0x > 0, yet neither LpL^p norm converges to zero.

The key point: convergence rates and even convergence itself depend on which norm you use. An algorithm converging in L2L^2 (minimising MSE) may behave very differently from one converging in L1L^1 (promoting sparsity and sharp edges).

Why This Matters: Norms and Image Reconstruction Quality

The choice of norm in

c^=argmincyAc2+λR(c)\hat{c} = \arg\min_{c}\, \bigl\|\mathbf{y} - \mathbf{A}\,\mathbf{c}\bigr\|^2 + \lambda\,\mathcal{R}(c)

directly determines the character of the solution via the regulariser R\mathcal{R}:

  • L2L^2 (Tikhonov): R=L22\mathcal{R} = \|\cdot\|_{L^2}^2 — penalises large values, produces smooth reconstructions, leads to linear systems.
  • L1L^1 (sparsity): R=L1\mathcal{R} = \|\cdot\|_{L^1} — promotes exactly zero pixels, the basis of compressed sensing for radar.
  • Total variation: R=L1\mathcal{R} = \|\nabla\cdot\|_{L^1} — promotes piecewise-constant images with sharp edges, ideal for SAR urban scenes.

The completeness of the underlying space guarantees that minimisers exist and iterative algorithms converge — without it, none of these guarantees hold. We return to this in Chapter 2 (Regularisation Theory).

See full treatment in Chapter 2

Quick Check

Which LpL^p space is also a Hilbert space?

L1(Ω)L^1(\Omega)

L2(Ω)L^2(\Omega)

L(Ω)L^\infty(\Omega)

All LpL^p spaces are Hilbert spaces

Quick Check

In R2\mathbb{R}^2, the 1\ell^1 unit ball {x:x11}\{\mathbf{x} : \|\mathbf{x}\|_1 \leq 1\} is a ___.

Circle

Square with sides parallel to the axes

Diamond (rhombus) with vertices on the axes

Ellipse

Quick Check

True or False: All norms on L2([0,1])L^2([0,1]) are equivalent.

True

False

Normed space

A vector space V\mathcal{V} equipped with a norm  ⁣:V[0,)\|\cdot\| \colon \mathcal{V} \to [0,\infty) satisfying positive definiteness, absolute homogeneity, and the triangle inequality. See DNormed Space.

Related: Normed Space, Banach Space, Equivalence of Norms

Banach space

A normed space that is complete: every Cauchy sequence converges to an element of the space. Named after Stefan Banach. See DBanach Space.

Related: Banach Space, Cauchy Sequence, LpL^p Spaces are Banach Spaces

LpL^p space

The Lebesgue space of measurable functions whose pp-th power is integrable (1p<1 \leq p < \infty), or essentially bounded (p=p = \infty). Lp(Ω)L^p(\Omega) is a Banach space for all p[1,]p \in [1,\infty]; only L2L^2 is also a Hilbert space. See LpL^p Spaces" data-ref-type="definition">DLpL^p Spaces.

Related: LpL^p Spaces, LpL^p Spaces are Banach Spaces, Sequence Spaces p\ell^p

Completeness

A metric or normed space is complete if every Cauchy sequence in the space converges to a limit that belongs to the space. Completeness ensures "limits stay inside" and is essential for well-posedness of iterative algorithms. See DBanach Space.

Related: Banach Space, Cauchy Sequence

Key Takeaway

Four core messages of this section:

(1) Normed spaces generalise Cn\mathbb{C}^n to infinite dimensions. Imaging scenes are functions, not finite vectors, requiring Lp(Ω)L^p(\Omega) spaces with norms that measure "size" in problem-specific ways.

(2) Banach spaces are the right setting for iterative algorithms. Completeness guarantees that convergent iterative sequences have a limit within the space — without it algorithms could "converge to nothing."

(3) The LpL^p family is foundational: L2L^2 for energy-based and MSE reconstruction, L1L^1 for sparsity-promoting and edge-preserving solutions, Sobolev HsH^s for smooth prior models.

(4) Unlike finite dimensions, the choice of norm is consequential. In Cn\mathbb{C}^n all norms are equivalent. In infinite-dimensional function spaces they are not — the norm determines the topology, the notion of convergence, and ultimately the character of reconstruction.

⚠️Engineering Note

From Infinite Dimensions to Finite Grids

In practice, every imaging algorithm operates on a discretised grid. Replacing L2(Ω)L^2(\Omega) by CN\mathbb{C}^N (a grid of NN voxels) introduces discretisation error: the continuous forward operator A\mathcal{A} is approximated by a finite matrix ACM×N\mathbf{A} \in \mathbb{C}^{M \times N}.

The theorems of this section (completeness, Parseval's identity, the projection theorem) have exact finite-dimensional counterparts, so the transition is clean. However, the grid resolution NN must be chosen to satisfy the sampling theorem for the maximum spatial frequency supported by the aperture — under-sampling introduces aliasing that no regularisation can remove.

Practical Constraints
  • Grid spacing Δrλ/(2NA)\Delta r \leq \lambda / (2 \mathrm{NA}) where NA\mathrm{NA} is the numerical aperture — this is the Nyquist criterion for imaging.

  • For 3-D scenes at 10 GHz with 10 cm aperture, a λ/4=7.5\lambda/4 = 7.5 mm grid gives N106N \sim 10^6 voxels — GPU memory becomes the binding constraint.

  • Coarser grids speed up computation but blur the PSF and reduce effective resolution; fine grids require O(N2)O(N^2) storage for the full sensing matrix A\mathbf{A}.