Signal-Space Detection and the Minimum-Distance Decoder

Why Signal Space?

Continuous-time passband signals {sm(t)}m=0Mβˆ’1\{s_m(t)\}_{m=0}^{M-1} live in L2(R)L^2(\mathbb{R}) β€” an infinite-dimensional Hilbert space. Detection in such a space looks hopelessly abstract. The key structural result of this section is that only an N≀MN \leq M-dimensional subspace matters. By projecting the received waveform onto an orthonormal basis of this subspace, we reduce detection of waveforms in function space to detection of vectors in RN\mathbb{R}^N (or CN\mathbb{C}^N) β€” the latter is exactly the setting of Section 3.1. This reduction is the theoretical backbone of every digital demodulator.

Definition:

M-ary Signal Detection in AWGN

Let {sm(t)}m=0Mβˆ’1\{s_m(t)\}_{m=0}^{M-1} be a finite set of finite-energy real (or complex baseband) waveforms on [0,T][0,T] with Em=βˆ₯smβˆ₯2=∫0T∣sm(t)∣2dtE_m = \|s_m\|^2 = \int_0^T |s_m(t)|^2 dt. Under Hm\mathcal{H}_m the received signal is Y(t)β€…β€Š=β€…β€Šsm(t)+Z(t),t∈[0,T],Y(t) \;=\; s_m(t) + Z(t), \qquad t\in[0,T], where Z(t)Z(t) is white Gaussian noise with two-sided PSD N0/2N_0/2. The average symbol energy is Es=1Mβˆ‘m=0Mβˆ’1EmE_s = \frac{1}{M}\sum_{m=0}^{M-1} E_m and the SNR is SNR=Es/N0\text{SNR} = E_s/N_0.

The signals sms_m need not be orthogonal and need not have equal energy; we will handle the general case below.

Theorem: Existence of an Orthonormal Signal-Space Basis

For any set of MM finite-energy waveforms {sm(t)}\{s_m(t)\} there exist orthonormal functions {Ο•n(t)}n=1N\{\phi_n(t)\}_{n=1}^{N} with N≀MN \leq M such that each sms_m can be written as sm(t)β€…β€Š=β€…β€Šβˆ‘n=1Nxm,n ϕn(t),xm,n=∫0Tsm(t)Ο•n(t)‾ dt,s_m(t) \;=\; \sum_{n=1}^{N} x_{m,n}\, \phi_n(t), \qquad x_{m,n} = \int_0^T s_m(t)\overline{\phi_n(t)}\,dt, and the map sm↦xm=(xm,1,…,xm,N)⊀s_m \mapsto \mathbf{x}_m = (x_{m,1},\ldots,x_{m,N})^\top is an isometry: βˆ₯smβˆ’sjβˆ₯2=βˆ₯xmβˆ’xjβˆ₯2\|s_m - s_j\|^2 = \|\mathbf{x}_m - \mathbf{x}_j\|^2. The dimension NN equals the rank of the Gram matrix Gmj=⟨sm,sj⟩G_{mj} = \langle s_m, s_j\rangle.

This is exactly the Gram-Schmidt procedure from finite-dimensional linear algebra, applied in the Hilbert space L2[0,T]L^2[0,T]. The isometric embedding means that inner products, norms, and distances all survive the reduction from waveforms to vectors.

Theorem: The Projection Vector is a Sufficient Statistic

Define the projection coefficients Yn=⟨Y,Ο•n⟩Y_n = \langle Y, \phi_n\rangle for n=1,…,Nn = 1,\ldots,N and the "tail" process Y~(t)=Y(t)βˆ’βˆ‘n=1NYnΟ•n(t)\tilde Y(t) = Y(t) - \sum_{n=1}^N Y_n \phi_n(t). Then Y=(Y1,…,YN)⊀=xm+Z,Z∼N(0,N02IN)β€…β€Šβ€…β€Š(realΒ case)\mathbf{Y} = (Y_1,\ldots,Y_N)^\top = \mathbf{x}_m + \mathbf{Z},\qquad \mathbf{Z}\sim \mathcal{N}(\mathbf{0}, \tfrac{N_0}{2}\mathbf{I}_N)\;\;(\text{real case}) and Y~\tilde Y is independent of Y\mathbf{Y} and has the same distribution under every Hm\mathcal{H}_m. Consequently Y\mathbf{Y} is a sufficient statistic for detection: no decision rule using the full waveform Y(β‹…)Y(\cdot) can outperform the optimal rule based on Y\mathbf{Y}.

White noise has independent components in any orthonormal basis. The NN "in-subspace" components carry all the signal information; the infinitely many "out-of-subspace" components carry only noise whose distribution is identical under every hypothesis β€” hence irrelevant.

Definition:

Minimum-Distance (ML) Decoder

Under the sufficient statistic Y=xm+Z\mathbf{Y} = \mathbf{x}_m + \mathbf{Z} with Z∼N(0,N02IN)\mathbf{Z} \sim \mathcal{N}(\mathbf{0}, \tfrac{N_0}{2}\mathbf{I}_N), the ML decision rule is gML(y)β€…β€Š=β€…β€Šarg⁑min⁑m∈{0,…,Mβˆ’1}βˆ₯yβˆ’xmβˆ₯2g_{\text{ML}}(\mathbf{y}) \;=\; \arg\min_{m\in\{0,\ldots,M-1\}} \|\mathbf{y} - \mathbf{x}_m\|^2 β€” i.e., decide in favor of the constellation point closest in Euclidean distance to the received vector. With equiprobable signals of equal energy Em=EE_m = E, this is also the MAP rule.

Expanding the squared distance gives an equivalent correlator form: arg⁑max⁑m(Re⟨y,xmβŸ©βˆ’12βˆ₯xmβˆ₯2)\arg\max_m \bigl( \text{Re}\langle \mathbf{y},\mathbf{x}_m\rangle - \tfrac{1}{2}\|\mathbf{x}_m\|^2 \bigr). When energies are unequal the βˆ’12βˆ₯xmβˆ₯2-\tfrac{1}{2}\|\mathbf{x}_m\|^2 term is the energy correction.

Example: QPSK: Decision Regions and Equivalent Scalar Statistic

Apply the minimum-distance decoder to QPSK with xm=Es(cos⁑θm,sin⁑θm)\mathbf{x}_m = \sqrt{E_s}(\cos\theta_m, \sin\theta_m) and ΞΈm=Ο€/4+mΟ€/2\theta_m = \pi/4 + m\pi/2 for m=0,1,2,3m = 0,1,2,3. What are the decision regions, and how does the decoder simplify?

Voronoi Decision Regions for 16-QAM

Voronoi Decision Regions for 16-QAM
The Voronoi cells for a standard 16-QAM constellation: decision regions are axis-aligned squares (boundaries at half-integer multiples of the grid spacing). Corner points have two neighbors, edge points have three, interior points have four.
⚠️Engineering Note

QAM Demapping Complexity and the I/Q Shortcut

A naive ML demapper for MM-QAM requires MM squared-distance computations per received symbol. For 1024-QAM (used in 5G NR with high-rank MIMO) this becomes prohibitive β€” 10610^6 multiplications per symbol at hundreds of millions of symbols per second. The key optimization: a square QAM constellation is the Cartesian product of two independent PAM constellations, so I and Q can be demapped separately in O(M)O(\sqrt M) operations each, bringing the total to O(M)O(\sqrt M) rather than O(M)O(M). Modern soft-output demappers additionally use the max-log approximation to avoid exponentials.

Practical Constraints
  • β€’

    3GPP NR supports constellations up to 1024-QAM (Rel-17) β€” demapping must fit in microseconds per transport block

  • β€’

    Soft LLR output requires the second-nearest point per bit, not only the nearest

πŸ“‹ Ref: 3GPP TS 38.211, Section 5.1

Received Samples and Minimum-Distance Decoder

Transmit random constellation symbols through an AWGN channel, color received samples by their ML decision, and compare against the true (transmitted) label to count errors. Adjust the SNR to see how the noise cloud shrinks and errors become rare.

Parameters
15
500

Common Mistake: Forgetting the Energy-Correction Term

Mistake:

A student writes the ML decoder as m^=arg⁑max⁑m⟨y,xm⟩\hat m = \arg\max_m \langle \mathbf{y}, \mathbf{x}_m\rangle β€” correlation β€” and applies it to a non-uniform constellation like a PSK subset with amplitude-modulated symbols.

Correction:

Pure correlation is optimal only when βˆ₯xmβˆ₯2\|\mathbf{x}_m\|^2 is constant across mm. For unequal-energy constellations (e.g. APSK, amplitude-and-phase-shift keying) you must include the energy correction: m^=arg⁑max⁑m(Re⟨y,xmβŸ©βˆ’12βˆ₯xmβˆ₯2)\hat m = \arg\max_m \bigl(\text{Re}\langle \mathbf{y}, \mathbf{x}_m\rangle - \tfrac{1}{2}\|\mathbf{x}_m\|^2\bigr). The geometric intuition: correlation finds the point with the largest projection, which biases decisions toward high-energy symbols.

Voronoi Cell

For a finite point set {xm}\{\mathbf{x}_m\} in Euclidean space, the Voronoi cell of xm\mathbf{x}_m is Vm={y:βˆ₯yβˆ’xmβˆ₯≀βˆ₯yβˆ’xjβˆ₯β€‰βˆ€jβ‰ m}\mathcal{V}_m = \{\mathbf{y} : \|\mathbf{y} - \mathbf{x}_m\| \leq \|\mathbf{y} - \mathbf{x}_j\|\,\forall j \neq m\} β€” the set of points closer to xm\mathbf{x}_m than to any other point in the set. ML decision regions for AWGN with equiprobable, equal-energy symbols are exactly Voronoi cells.

Related: Decision Region, Minimum Distance Decoder

Quick Check

Three orthogonal waveforms s0,s1,s2s_0, s_1, s_2 of unit energy are used for 3-ary signaling in AWGN. What is the dimension NN of the signal space, and what is the Euclidean distance between every pair?

N=3N = 3, d=2d = \sqrt{2}

N=3N = 3, d=2d = 2

N=2N = 2, d=2d = \sqrt{2}

N=2N = 2, d=1d = 1

Why This Matters: Signal-Space View of OFDM and CDMA

OFDM and CDMA are both instances of the signal-space framework with different basis choices. OFDM picks Ο•n(t)=ej2Ο€nt/T1[0≀t≀T]\phi_n(t) = e^{j2\pi n t/T} \mathbb{1}[0\leq t \leq T] β€” complex exponentials separated in frequency. CDMA picks Ο•n(t)\phi_n(t) as time-translates of a PN sequence β€” separated in a code domain. The minimum-distance decoder of this section is then matched-filter demodulation in each case, and the analysis below applies directly.

See full treatment in Chapter 14