28. The Rate for a Source Relative to a Fidelity Evaluation
Defining the Rate-Distortion Function
Let \(R(D)\) = minimum rate (bits per source symbol) required to represent a source with average distortion \(\leq D\).
Formally:
Minimization over all conditional distributions \(p(y|x)\) (test channels) that achieve distortion at most \(D\).
Properties of \(R(D)\)
- Monotonically decreasing: Lower distortion requires higher rate
- Convex: \(R(\lambda D_1 + (1-\lambda) D_2) \leq \lambda R(D_1) + (1-\lambda) R(D_2)\)
- \(R(0) = H(x)\): Lossless rate equals source entropy
- \(R(D_{\max}) = 0\): At maximum tolerable distortion, send nothing
graph LR
subgraph "R(D) Curve"
direction TB
A["R(0) = H"] -->
B["R(D) decreases"] -->
C["R(D_max) = 0"]
endThe Gaussian Source (Squared Error)
For a Gaussian source with variance \(\sigma^2\) and MSE distortion:
Shannon's reverse water-filling: For a Gaussian, the rate-distortion function has a closed form. At distortion \(D\):
- Allocate \(R = \frac{1}{2}\log(\sigma^2/D)\) bits
- Reconstruction error variance = \(D\)
- The test channel is additive Gaussian noise: \(y = x + z\) where \(z \sim \mathcal{N}(0, D)\)
The Binary Source (Hamming Error)
For a Bernoulli(\(p\)) source with Hamming distortion:
for \(0 \leq D \leq \min(p, 1-p)\).
At \(D = 0\): \(R(0) = H_2(p)\) (lossless rate).
At \(D = p\): \(R(p) = 0\) (just guess the more likely symbol).
Operational Meaning
\(R(D)\) is both:
1. Information-theoretic: minimum mutual information over test channels
2. Operational: there exist codes achieving any rate \(\u003e R(D)\) with distortion \(\leq D\)
This is the lossy analog of the source coding theorem.
Application to Modern Compression
| Standard | Source | Distortion | Typical Operating Point |
|---|---|---|---|
| JPEG | Images | DCT-weighted MSE | \(R \approx 0.5\)–\(2\) bpp |
| MP3 | Audio | Perceptual (masking) | \(R \approx 128\) kbps |
| H.264/AVC | Video | Motion-compensated MSE | \(R \approx 2\)–\(10\) Mbps |
| JPEG 2000 | Images | Wavelet MSE | \(R \approx 0.25\)–\(4\) bpp |
None achieve the theoretical \(R(D)\) exactly, but modern codecs get within 1–3 dB (in rate) of optimal.