形式陈述
先固定本页的单字母模型:X 1 , X 2 , … 是分布为 P X 的离散无记忆源,源字母表与重构字母表有限;d ( x , x ^ ) ∈ [ 0 , ∞ ) 是单字母失真,预算取 D ≥ 0 ;若约束下没有测试信道,约定下确界为 + ∞ 。块失真取可分平均
d n ( x n , x ^ n ) = 1 n ∑ i = 1 n d ( x i , x ^ i ) . 率失真函数把选择测试信道写成一个优化问题 公理库 优化问题 Optimization problem 在可行解集合上最小化或最大化目标函数的计算问题。 :遍历从源字母到重构字母的条件分布 公理库 条件分布 Conditional distribution · Regular conditional distribution 给定观测值后随机量的概率律,以及它与原联合分布相容的核表示。 P X ^ ∣ X ,单字母率失真函数定义为
R ( D ) = inf P X ^ ∣ X : E [ d ( X , X ^ ) ] ≤ D I ( X ; X ^ ) , 其中期望 公理库 期望 Expectation · Expected value 实值或复值随机变量关于概率测度的 Lebesgue 积分,概括加权平均与总体质量平衡。 与互信息都按 P X P X ^ ∣ X 形成的联合分布计算。对固定 P X 和 d ,R ( D ) 随 D 非增且为凸函数;改变源分布、重构字母表或失真准则,就改变了这条曲线。
率失真定理把这个信息量优化与长块编码联系起来。一个块长为 n 、码本大小为 M n 的有损码先把 X n 编成索引,再由索引重构 X ^ n ;它的每符号码率是 n − 1 log 2 M n 。在上述有限字母表模型中,若要求
lim sup n → ∞ E [ d n ( X n , X ^ n ) ] ≤ D , 则可达的最小渐近每符号码率是 R ( D ) :高于它的速率可由充分长的块码达到,低于它则不能维持该平均失真约束。这里的操作结论依赖 IID 源、可分单字母失真和渐近块长,不能只看单字母公式便推广到所有源。
本页另行固定的连续模型
以下高斯部分改取实字母表:X ∼ N ( 0 , σ 2 ) 、0 < σ 2 < ∞ ,失真为 ( X − X ^ ) 2 。遍历从实数Borel空间到自身的所有概率核 公理库 概率核(Markov 核) Probability kernel · Markov kernel · 转移核 从每个输入状态可测地指定一个输出概率分布的映射。 P X ^ ∣ X ,仍以联合分布相对边缘乘积的KL定义互信息。重构分布可以离散、连续或奇异,不要求具有密度。这个模型的答案为
R ( D ) = { + ∞ , D = 0 , 1 2 log 2 σ 2 D , 0 < D < σ 2 , 0 , D ≥ σ 2 . 对IID高斯源、可分平方误差和期望块失真,连续字母表率失真编码定理赋予该式相同的渐近操作意义。此处源有有限二阶矩,固定重构 0 的平均失真为 σ 2 ,满足相应可积性条件;这不是由前面的有限字母表定理直接推出的。下面先完整求出信息量优化,再说明编码与传输如何接入。
直觉
率失真函数问:若允许平均重构误差不超过 D ,编码器至少要为每个源符号保留多少 bit 的信息。每个候选测试信道 P X ^ ∣ X 都描述一种随机重构机制;它与源分布共同生成 ( X , X ^ ) ,互信息则衡量重构仍携带多少源信息。约束更宽时,可选机制更多,所以最小互信息不会上升。
测试信道不是要求实际 codec 逐符号随机输出。它刻画最优长块码所诱导的单字母统计关系,证明再通过典型序列或随机码本把这种关系转成块编码。这个视角既保留了“质量—码率折衷”的图像,也解释了为什么一次标量量化通常达不到理论边界。
高斯下界:比较联合分布,不拆未必存在的熵
固定 D > 0 ,任取可行联合分布 J = P X X ^ ,记实际失真 δ = E ( X − X ^ ) 2 ≤ D 。设 ϕ v 是 N ( 0 , v ) 的密度,构造两个参考分布
M ( d x , d x ^ ) = ϕ σ 2 ( x ) d x P X ^ ( d x ^ ) , Q D ( d x , d x ^ ) = ϕ D ( x − x ^ ) d x P X ^ ( d x ^ ) . M 是边缘乘积;Q D 保留同一个重构边缘,但把给定重构后的源改成均值 x ^ 、方差 D 的高斯。两者相对于 d x ⊗ P X ^ 都有严格正密度,因此彼此等价,即使 P X ^ 没有Lebesgue密度也没有问题。
若 I ( X ; X ^ ) = + ∞ ,所需下界已成立。否则 J ≪ M ≪ Q D ,代入一元正态分布的密度公式 公理库 正态分布 Normal distribution · Gaussian distribution · 高斯分布 具有指数平方密度、在仿射变换与独立求和下封闭的概率分布族。 得到的对数比
log 2 d M d Q D ( x , x ^ ) = 1 2 log 2 D σ 2 + ( x − x ^ ) 2 2 D ln 2 − x 2 2 σ 2 ln 2 在 J 下绝对可积,因为 δ < ∞ 、E X 2 = σ 2 。有限KL的对数密度比也绝对可积:其负部积分由 1 / ( e ln 2 ) 控制,正部则由KL有限性控制。因此可以合法相加并取期望,得到
I ( X ; X ^ ) = D 2 ( J ‖ Q D ) + 1 2 log 2 σ 2 D + D − δ 2 D ln 2 . KL非负性 公理库 KL 散度 Kullback–Leibler divergence · Relative entropy 同一可测空间上分布 P 相对于 Q 的对数 Radon–Nikodym 导数在 P 下的积分。 与 δ ≤ D 给出 I ≥ 1 2 log 2 ( σ 2 / D ) 。这个证明同时处理了连续源与离散重构,也避免使用可能没有定义的 h ( X ∣ X ^ ) 。
后向测试信道达到下界
当 0 < D < σ 2 时,先独立抽取
X ^ ∼ N ( 0 , σ 2 − D ) , Z ∼ N ( 0 , D ) , X = X ^ + Z . 高斯和保证 X ∼ N ( 0 , σ 2 ) ,而平方误差恰为 D 。此时联合分布就是 Q D ;它与边缘乘积的对数比为上面二次式的相反数,在联合分布下绝对可积,故互信息有限。精确差额中的两项均为零,于是
I ( X ; X ^ ) = 1 2 log 2 σ 2 D . 虽然构造从重构向源抽样,但联合分布确定了合法的正向核。由二维高斯分布 公理库 多元正态分布 Multivariate normal distribution · Multivariate Gaussian distribution · Jointly Gaussian vector 以所有线性组合都正态刻画联合高斯向量,并由特征函数连接线性构造、退化支撑、全维密度与条件分布。 的条件均值和方差,它等价于
X ^ = ( 1 − D σ 2 ) X + V , V ⊥ X , V ∼ N ( 0 , D ( 1 − D σ 2 ) ) . 例如 Cov ( X , X ^ ) = σ 2 − D ,所以回归系数是 ( σ 2 − D ) / σ 2 ;条件方差为 ( σ 2 − D ) − ( σ 2 − D ) 2 / σ 2 。后向误差 Z = X − X ^ 独立于重构 ,一般不独立于源;正向独立噪声是另一个变量 V 。
例子与边界
离散无记忆源配 Hamming 失真,且每个正概率源符号都能在重构字母表中原样输出时,D = 0 要求逐符号准确重构,因此 R ( 0 ) = H ( X ) 。当 D 大到某个固定重构符号已经满足约束时,编码器无需观察输入也能达到目标,码率于是降为 0 。率失真函数作为互信息的下确界始终非负,这一点不依赖连续模型中的差分熵是否可能为负。
对 Bernoulli( 1 / 2 ) 源和 Hamming 失真,在 0 ≤ D ≤ 1 / 2 时有 R ( D ) = 1 − H 2 ( D ) 。D = 0 恢复无损压缩率 1 bit/符号,D = 1 / 2 时无需发送信息也可随机猜测达到允许失真。
方差4、失真1:直接加噪声为什么多耗信息
最优正向核是 X ^ = 3 4 X + V ,其中 V ∼ N ( 0 , 3 / 4 ) 独立于 X ∼ N ( 0 , 4 ) 。逐项复算为
源 符 号 E ( X − X ^ ) 2 = ( 1 4 ) 2 4 + 3 4 = 1 , Var ( X ^ ) = ( 3 4 ) 2 4 + 3 4 = 3 , I ( X ; X ^ ) = 1 2 log 2 3 3 / 4 = 1 bit/源符号 . 作为比较,取 X ~ = X + W ,W ∼ N ( 0 , 1 ) 独立于源,也有失真 1 ,但
源 符 号 I ( X ; X ~ ) = 1 2 log 2 5 1 ≈ 1.160964 bit/源符号 . 两个互信息都可用有限高斯密度的积分计算。最优核在加入正向噪声前先收缩信号;只检查噪声方差等于目标失真,无法判断互信息是否最小。
连续源的两个端点
D ≥ σ 2 时恒取 X ^ = 0 ,失真为 σ 2 ,互信息为零。D = 0 时必有 X ^ = X 几乎处处:联合分布把全部质量放在对角线,而非原子的高斯边缘乘积给对角线质量零。联合分布不绝对连续于边缘乘积,所以互信息为无穷。连续源精确重构的端点与有限字母表Hamming零失真不同。
前面的正向公式用于 0 < D < σ 2 ;D = σ 2 应直接使用常量重构,不能把退化的条件方差继续代入高斯密度的对数比。
率失真定理是长块、平均失真意义下的极限,不保证每个样本都低于 D ,也不提供有限块长、有限时延或低复杂度编码器。实际 codec 还受模型族、算力和缓冲时延约束;若失真函数没有表达人的感知差异,数学上的最优重构也未必主观最好。
有记忆源或一般源通常不再由这个单字母式完整描述。此时需要研究 n 维分布与块失真的多字母优化,再取归一化极限;更一般的非平稳情形还可能使用信息谱表述。连续字母表也需要可测性、矩条件和失真可积性等假设,不能把有限字母表定理删去条件后直接沿用。
推论与应用
率失真理论给出图像、音频和有损压缩的基准极限,并连接量化、感知质量和信息瓶颈。互信息 公理库 互信息 Mutual information 用联合分布相对独立边缘乘积的 KL 散度量化统计依赖。 是单字母优化目标,熵 公理库 Shannon 熵 Shannon entropy · Information entropy 随机变量不确定性的平均信息量,以最优编码所需位数为基本解释。 给出离散无记忆源在零 Hamming 失真处的无损端点;与信道容量 公理库 信道容量 Channel capacity 对输入分布最大化输入与输出互信息所得的每次使用信息率。 结合后,它还能判断给定源与信道是否存在满足目标失真的分离式传输方案。
把每源符号与每信道使用的单位对齐
设 X n 是方差 σ 2 的IID高斯源;编码器把它映成 m 个实信道输入 U m ,通过 Y j = U j + Z j ,其中 N > 0 ,噪声 Z j ∼ N ( 0 , N ) 独立同分布并独立于源与编码器随机性。无反馈,接收端由 Y m 重构 X ^ n 。固定 0 ≤ P < ∞ ,功率约束现在按实际源分布 取期望:
1 m ∑ j = 1 m E U j 2 ≤ P . 只需考察平均失真有限的方案,否则以下失真下界自动成立。令 δ i = E ( X i − X ^ i ) 2 ,δ ¯ = n − 1 ∑ i δ i 。率失真的凸性、定义、源的独立性、数据处理不等式 公理库 数据处理不等式 Data processing inequality 对 Markov 链 X→Y→Z,有 I(X;Z)≤I(X;Y)。 以及高斯信道容量 公理库 高斯信道容量 Gaussian channel capacity · AWGN capacity · 实AWGN信道容量 在平均平方功率约束下,以KL恒等式证明实AWGN互信息极值,并连接块编码逆界与高斯源的失真预算。 的块逆界依次给出
n R ( δ ¯ ) ≤ ∑ i = 1 n R ( δ i ) ≤ ∑ i = 1 n I ( X i ; X ^ i ) ≤ I ( X n ; X ^ n ) ≤ I ( U m ; Y m ) ≤ m C ( P ) . 其中容易看反的第三步可由链式法则核对:因为 X i 与 X i − 1 独立,
I ( X n ; X ^ n ) = ∑ i I ( X i ; X ^ n , X i − 1 ) ≥ ∑ i I ( X i ; X ^ i ) . 块信道的上界有限,也保证这里相关的互信息有限。代入两个闭式公式,若 δ ¯ < σ 2 ,得
δ ¯ ≥ σ 2 ( 1 + P N ) − m / n . 若 δ ¯ ≥ σ 2 ,该不等式自动成立。固定资源比
( 信 道 使 用 次 数 源 符 号 ) α = lim n → ∞ m n ∈ ( 0 , ∞ ) (信道使用次数/源符号) , 最优渐近期望失真因而至少为 D ∗ ( α ) = σ 2 ( 1 + P / N ) − α 。
可达性、严格余量与无界失真
当目标满足 R ( D ) < α C ( P ) 且 0 < D < σ 2 时,由连续性先取略小的 D ′ < D ,仍使 R ( D ′ ) < α C ( P ) ,为拼接误差留下失真余量。再选压缩速率略大于 R ( D ′ ) 、每信道使用的码率略小于 C ( P ) ,拼接率失真码和信道码。源索引未必均匀,所以信道码使用最大消息错误概率保证;每个信道码字都满足能量约束,故也满足实际源分布下的功率要求。
平方误差无界,还需控制信道译码出错时的失真。连续源编码/分离定理允许本高斯模型选取重构码字能量统一不超过 n K 的码本,其中 K 与块长无关。若信道最大错误概率为 ϵ ,记错误事件为 E ,由 Pr ( E ∣ X n = x n ) ≤ ϵ 和平方三角不等式,
1 n E [ ‖ X n − X ^ n ‖ 2 1 E ] ≤ 2 n E [ ‖ X n ‖ 2 1 E ] + 2 K Pr ( E ) ≤ 2 ϵ ( σ 2 + K ) ⟶ 0. 正确译码部分不超过源码自身的平均失真。有限二阶矩和这个坏事件控制,使两种编码能够真正拼接,而不只是把两个公式相除。所用的连续率失真及有损分离结论见MIT讲义Theorems24.2、25.4–25.5及Lemma25.1。
让严格余量趋零,边界 D ∗ ( α ) 在渐近闭包意义下可达:给任意正容差,可取足够长的码使资源比趋近 α 、平均失真至多 D ∗ ( α ) 加该容差。因此它就是最优渐近失真;一般不表示存在有限块码恰好达到边界。P = 0 时常量重构直接给出 D ∗ = σ 2 。
例如 σ 2 = 4 、P / N = 3 ,C = 1 bit/实使用:
每源符号的信道次数 α
每源符号的可靠信息预算 α C
最优渐近失真 D ∗ ( α )
1 / 2
1 / 2 bit
2
1
1 bit
1
2
2 bit
1 / 4
一次对一次的匹配:无需长块便达到边界
当 α = 1 时,取缩放发送与线性重构
U = P σ 2 X , Y = U + Z , X ^ = P σ 2 P + N Y . 逐符号便有 E U 2 = P ,而
X − X ^ = N P + N X − P σ 2 P + N Z , E ( X − X ^ ) 2 = N 2 σ 2 + P σ 2 N ( P + N ) 2 = σ 2 N P + N = D ∗ ( 1 ) . 它诱导的正向系数 P / ( P + N ) 与噪声方差 P σ 2 N / ( P + N ) 2 ,恰好对应前面取 D = σ 2 N / ( P + N ) 的最优测试信道。数值例 P = 3 , N = 1 , σ 2 = 4 中,发送和接收均乘 3 / 2 ,合成后为 X ^ = 3 4 X + 3 2 Z ,失真正好为1。
这种高斯匹配解释了一个特殊的有限延时最优方案。其功率按源分布取期望;当 P > 0 时,高斯源幅度无界,所以它不满足“每个可能源块都对应一个能量至多 m P 的发送块”。把功率量词改强后,不能继续声称该逐符号方案原样可行。
若译码器另有相关观察,Wyner–Ziv 编码 公理库 Wyner–Ziv 带边信息有损编码 Wyner-Ziv coding 在编码器看不到边信息时,用辅助描述与分箱得到译码端边信息下的最小有损编码率。 在辅助变量 U − X − Y 的限制下最小化条件互信息,明确区分编码器是否也看到边信息。压缩转发 公理库 压缩转发中继 Compress-forward relaying 让中继压缩自己的观察供终点联合解释,用描述质量与描述链路容量共同约束可达率。 进一步把这种描述送过一条中继信道;描述所需速率必须不超过实际可用链路,不能仅算压缩质量而省掉传送费用。
参考资料
Thomas M. Cover and Joy A. Thomas, Elements of Information Theory , 2nd ed., Wiley, 2006,Chapter 10 “Rate Distortion Theory”。
Claude E. Shannon, “Coding Theorems for a Discrete Source With a Fidelity Criterion ,” IRE National Convention Record , Part 4, 1959, pp. 142–163。
Yury Polyanskiy and Yihong Wu, Lecture Notes on Information Theory , MIT6.441, Spring2016课程讲义:§24.2、Theorem24.2,印刷页247,连续无记忆源的率失真定理;§25.1.2,pp.256–257,高斯源及KL下界;§25.3.1–25.3.2、Theorems25.3–25.5与Lemma25.1,pp.262–265,分离逆界、最大错误概率与无界失真控制。正文展开联合分布比较与每个不等式的方向。
Tsachy Weissman, EE276 Lecture18: Joint source-channel coding2 , 2024-03-12,pp.1–3:渐近闭包口径及高斯源的逐符号最优缩放方案;其资源比率为本页 α 的倒数。