Error Probability: Union Bounds and Exact Formulas

Why Bounds Matter

Exact symbol-error probability expressions exist only for a handful of highly symmetric constellations (orthogonal, simplex, M-PSK, square QAM). For the vast majority of practical signal sets β€” and certainly for any lattice- or code-based constellation β€” closed-form PeP_e is out of reach. Bounds that are provably tight for high SNR (union bound, nearest-neighbor bound) and simple enough to evaluate by hand fill this gap. They are the daily working tools of every communication engineer.

Definition:

Pairwise Error Probability

Given constellation points xm,xj\mathbf{x}_m, \mathbf{x}_j in AWGN, the pairwise error probability (PEP) is the probability that a minimum-distance decoder β€” restricted to just these two signals β€” decides jj given that mm was transmitted: P(mβ†’j)β€…β€Š=β€…β€ŠPr⁑(βˆ₯Yβˆ’xjβˆ₯2<βˆ₯Yβˆ’xmβˆ₯2β€‰βˆ£β€‰Hm).P(m\to j) \;=\; \Pr\bigl(\|\mathbf{Y} - \mathbf{x}_j\|^2 < \|\mathbf{Y}-\mathbf{x}_m\|^2 \,\bigm|\, \mathcal{H}_m\bigr). For Z∼N(0,N02IN)\mathbf{Z} \sim \mathcal{N}(\mathbf{0}, \tfrac{N_0}{2}\mathbf{I}_N) this simplifies to P(mβ†’j)β€…β€Š=β€…β€ŠQ ⁣(βˆ₯xmβˆ’xjβˆ₯2N0)β€…β€Š=β€…β€ŠQ ⁣(dm,j22N0).P(m\to j) \;=\; Q\!\left(\frac{\|\mathbf{x}_m - \mathbf{x}_j\|}{\sqrt{2 N_0}}\right) \;=\; Q\!\left(\sqrt{\frac{d_{m,j}^2}{2 N_0}}\right).

The PEP depends only on the distance dm,j=βˆ₯xmβˆ’xjβˆ₯d_{m,j} = \|\mathbf{x}_m - \mathbf{x}_j\| β€” a consequence of the isotropy of the AWGN noise vector.

Theorem: Union Bound on Symbol Error Probability

For a constellation with equiprobable symbols in AWGN, Peβ€…β€Šβ‰€β€…β€Š1Mβˆ‘m=0Mβˆ’1βˆ‘jβ‰ mP(mβ†’j)β€…β€Š=β€…β€Š1Mβˆ‘m=0Mβˆ’1βˆ‘jβ‰ mQ ⁣(dm,j2N0).P_e \;\leq\; \frac{1}{M}\sum_{m=0}^{M-1}\sum_{j\neq m} P(m\to j) \;=\; \frac{1}{M}\sum_{m=0}^{M-1}\sum_{j\neq m} Q\!\Bigl(\tfrac{d_{m,j}}{\sqrt{2 N_0}}\Bigr). A looser β€” but useful β€” universal bound is obtained by replacing every pairwise distance with the minimum distance: Peβ€…β€Šβ‰€β€…β€Š(Mβˆ’1) Q ⁣(dmin⁑2N0),dmin⁑=min⁑mβ‰ jdm,j.P_e \;\leq\; (M-1)\, Q\!\Bigl(\tfrac{d_{\min}}{\sqrt{2 N_0}}\Bigr), \qquad d_{\min} = \min_{m\neq j} d_{m,j}.

The event "the ML detector makes an error" is the union ⋃jβ‰ m{βˆ₯yβˆ’xjβˆ₯<βˆ₯yβˆ’xmβˆ₯}\bigcup_{j\neq m} \{\|\mathbf{y}-\mathbf{x}_j\|<\|\mathbf{y}-\mathbf{x}_m\|\}. Replacing the union probability by the sum of individual PEPs (the union bound) yields a bound that is tight when the component events are approximately disjoint β€” which is the case at high SNR.

Theorem: Nearest-Neighbor (High-SNR) Bound

Let Nnn(m)N_{\text{nn}}(m) be the number of constellation points at minimum distance from xm\mathbf{x}_m, and NΛ‰nn=1Mβˆ‘mNnn(m)\bar N_{\text{nn}} = \tfrac{1}{M}\sum_m N_{\text{nn}}(m). Then at high SNR Peβ€…β€Šβ‰ˆβ€…β€ŠNΛ‰nnβ‹…Q ⁣(dmin⁑2N0).P_e \;\approx\; \bar N_{\text{nn}} \cdot Q\!\left(\frac{d_{\min}}{\sqrt{2 N_0}}\right). More precisely, the ratio of PeP_e to this expression tends to 1 as SNRβ†’βˆž\text{SNR} \to \infty.

At high SNR, errors to points beyond the nearest neighbors are exponentially less likely than errors to the nearest neighbors (the QQ decays as eβˆ’d2/(4N0)e^{-d^2/(4 N_0)}). The dominant contribution to PeP_e comes from the NΛ‰nn\bar N_{\text{nn}} minimum-distance neighbors, so the union bound is asymptotically tight with the correct constant.

Example: Exact PeP_e for BPSK and QPSK

Derive the exact symbol-error probability for BPSK (x0,1=Β±Es\mathbf{x}_{0,1} = \pm\sqrt{E_s}) and QPSK in AWGN.

Theorem: Exact SER for M-PSK

For equi-energy MM-PSK with xm=Es(cos⁑(2Ο€m/M),sin⁑(2Ο€m/M))\mathbf{x}_m = \sqrt{E_s}(\cos(2\pi m/M),\sin(2\pi m/M)) in AWGN, the exact symbol-error probability is PeMPSKβ€…β€Š=β€…β€Š1βˆ’βˆ«βˆ’Ο€/MΟ€/MpΘ(ΞΈ) dΞΈ,P_e^{\text{MPSK}} \;=\; 1 - \int_{-\pi/M}^{\pi/M} p_\Theta(\theta)\,d\theta, where pΘp_\Theta is the density of the angle of Y\mathbf{Y} given x0\mathbf{x}_0 was sent. For large SNR\text{SNR} this integrates to Peβ‰ˆ2 Q ⁣(2SNRsin⁑(Ο€/M))P_e \approx 2\,Q\!\bigl(\sqrt{2\text{SNR}}\sin(\pi/M)\bigr).

By rotational symmetry, assume x0\mathbf{x}_0 was sent; the received point is correctly decoded iff its phase lies in the wedge ∣θ∣<Ο€/M|\theta| < \pi/M. The exact angular density has a closed form in terms of erfc and is given in the proof.

Example: Exact SER for Square M-QAM

Compute the exact symbol-error probability for square MM-QAM with M=L2M = L^2 (so L=ML = \sqrt M is an integer) and grid spacing Ξ”\Delta.

Union Bound vs. Exact Monte-Carlo SER

Compare three curves for a selected constellation: (i) the exact analytical PeP_e where available, (ii) the loose union bound (Mβˆ’1)Q(dmin⁑/2N0)(M-1)Q(d_{\min}/\sqrt{2N_0}), and (iii) the nearest-neighbor tight bound NΛ‰nnQ(dmin⁑/2N0)\bar N_{\text{nn}}Q(d_{\min}/\sqrt{2N_0}). At high SNR, the nearest-neighbor curve closes the gap to exact.

Parameters
0
25

SER Curves for M-PSK

Compute the symbol-error probability of MM-PSK for M∈{2,4,8,16}M \in \{2,4,8,16\} over a range of symbol SNRs, using the exact angular-integral formula. The curves illustrate the 3-dB penalty incurred on each doubling of MM beyond QPSK.

Parameters
0
25

SER Curves for Square M-QAM

Exact symbol-error curves for MM-QAM with M∈{4,16,64,256}M \in \{4,16,64,256\}. The 6-dB penalty per quadrupling of MM (at fixed EsE_s) reflects the shrinking of dmin⁑2=6Es/(Mβˆ’1)d_{\min}^2 = 6E_s/(M-1).

Parameters
0
30
⚠️Engineering Note

BER Is Not SER: The Gray-Coding Gift

Link-budget engineers often quote "BER = 10βˆ’610^{-6}" but the curves in standards documents show SER. Under Gray coding β€” where nearest constellation neighbors differ in exactly one bit β€” a single symbol error causes on average one bit error out of log⁑2M\log_2 M, so BERβ‰ˆSER/log⁑2M\text{BER} \approx \text{SER}/\log_2 M at high SNR. For 64-QAM (log⁑2M=6\log_2 M = 6), that is a factor of 6 reduction. Without Gray coding (e.g. natural binary labeling) a single symbol error can flip Θ(log⁑2M)\Theta(\log_2 M) bits, and the high-SNR SER β‰ˆ\approx BER slope is lost. Every modern standard (LTE, NR, DVB, Wi-Fi) uses Gray-labeled QAM.

Practical Constraints
  • β€’

    Gray labeling exists only for constellations with 'neighbor-of-neighbor' graphs that are bipartite (PSK, PAM, square QAM all work)

  • β€’

    Non-square QAM (e.g. 32-QAM cross) has no perfect Gray labeling β€” standards use near-Gray approximations

πŸ“‹ Ref: 3GPP TS 38.211, Section 5.1.3; ITU-T G.9960

High-SNR SER Approximations for Common Constellations

Constellationdmin⁑2/Esd_{\min}^2/E_sHigh-SNR SERSpectral Eff. (b/s/Hz)
BPSK44Q(2SNR)Q(\sqrt{2\text{SNR}})1
QPSK222 Q(SNR)2\,Q(\sqrt{\text{SNR}}) (approx.)2
8-PSK2(1βˆ’cos⁑π/4)2(1-\cos\pi/4)2 Q(2SNRsin⁑π/8)2\,Q(\sqrt{2\text{SNR}}\sin\pi/8)3
16-QAM0.40.43 Q(SNR/5)3\,Q(\sqrt{\text{SNR}/5})4
64-QAMβ‰ˆ0.095\approx 0.0953.5 Q(SNR/21)3.5\,Q(\sqrt{\text{SNR}/21})6
256-QAMβ‰ˆ0.0235\approx 0.02353.75 Q(SNR/85)3.75\,Q(\sqrt{\text{SNR}/85})8

Common Mistake: The Loose Union Bound Exceeds 1 at Low SNR

Mistake:

Plotting (Mβˆ’1)Q(dmin⁑/2N0)(M-1)Q(d_{\min}/\sqrt{2N_0}) against SNR on log-log axes, a student reports "error probability = 5" at low SNR and concludes the simulation is broken.

Correction:

A probability bound can exceed 1 β€” it is simply uninformative in that regime. Always cap bounds at 1: Pe≀min⁑(1,(Mβˆ’1)Q(β‹…))P_e \leq \min\bigl(1, (M-1)Q(\cdot)\bigr). At very low SNR, PeP_e approaches 1βˆ’1/M1 - 1/M (random guessing) while the raw union bound grows linearly with MM, so the cap is essential.

Quick Check

For 16-QAM with symbol energy EsE_s, what is dmin⁑2d_{\min}^2?

Es/5E_s/5

2Es/52E_s/5

Es/4E_s/4

6Es6E_s

Pairwise Error Probability (PEP)

The probability P(m→j)P(m\to j) that a minimum-distance decoder restricted to two candidate symbols {xm,xj}\{\mathbf{x}_m, \mathbf{x}_j\} decides jj given that mm was sent. In AWGN with isotropic noise of variance N0/2N_0/2, P(m→j)=Q(dm,j/2N0)P(m\to j) = Q(d_{m,j}/\sqrt{2N_0}).

Related: Union Bound on Symbol Error Probability, Qfn