Skip to content

定义Definition

率失真函数

Rate-distortion function

允许平均失真D时的最小互信息,并以高斯平方误差证明下界、最优测试信道及源信道资源比例。

形式陈述 ​

先固定本页的单字母模型:X1,X2,… 是分布为 PX 的离散无记忆源,源字母表与重构字母表有限;d(x,x^)∈[0,∞) 是单字母失真,预算取 D≥0;若约束下没有测试信道,约定下确界为 +∞。块失真取可分平均

dn(xn,x^n)=1n∑i=1nd(xi,x^i).

率失真函数把选择测试信道写成一个优化问题:遍历从源字母到重构字母的条件分布 PX^∣X,单字母率失真函数定义为

R(D)=infPX^∣X:E[d(X,X^)]≤DI(X;X^),

其中期望与互信息都按 PXPX^∣X 形成的联合分布计算。对固定 PX 和 d,R(D) 随 D 非增且为凸函数;改变源分布、重构字母表或失真准则,就改变了这条曲线。

率失真定理把这个信息量优化与长块编码联系起来。一个块长为 n、码本大小为 Mn 的有损码先把 Xn 编成索引,再由索引重构 X^n;它的每符号码率是 n−1log2⁡Mn。在上述有限字母表模型中,若要求

lim supn→∞E[dn(Xn,X^n)]≤D,

则可达的最小渐近每符号码率是 R(D):高于它的速率可由充分长的块码达到,低于它则不能维持该平均失真约束。这里的操作结论依赖 IID 源、可分单字母失真和渐近块长,不能只看单字母公式便推广到所有源。

本页另行固定的连续模型 ​

以下高斯部分改取实字母表:X∼N(0,σ2)、0<σ2<∞,失真为 (X−X^)2。遍历从实数Borel空间到自身的所有概率核 PX^∣X,仍以联合分布相对边缘乘积的KL定义互信息。重构分布可以离散、连续或奇异,不要求具有密度。这个模型的答案为

R(D)={+∞,D=0,12log2⁡σ2D,0<D<σ2,0,D≥σ2.

对IID高斯源、可分平方误差和期望块失真,连续字母表率失真编码定理赋予该式相同的渐近操作意义。此处源有有限二阶矩,固定重构 0 的平均失真为 σ2,满足相应可积性条件;这不是由前面的有限字母表定理直接推出的。下面先完整求出信息量优化,再说明编码与传输如何接入。

直觉

率失真函数问:若允许平均重构误差不超过 D,编码器至少要为每个源符号保留多少 bit 的信息。每个候选测试信道 PX^∣X 都描述一种随机重构机制;它与源分布共同生成 (X,X^),互信息则衡量重构仍携带多少源信息。约束更宽时,可选机制更多,所以最小互信息不会上升。

测试信道不是要求实际 codec 逐符号随机输出。它刻画最优长块码所诱导的单字母统计关系,证明再通过典型序列或随机码本把这种关系转成块编码。这个视角既保留了“质量—码率折衷”的图像,也解释了为什么一次标量量化通常达不到理论边界。

高斯下界:比较联合分布,不拆未必存在的熵 ​

固定 D>0,任取可行联合分布 J=PXX^,记实际失真 δ=E(X−X^)2≤D。设 ϕv 是 N(0,v) 的密度,构造两个参考分布

M(dx,dx^)=ϕσ2(x)dxPX^(dx^),QD(dx,dx^)=ϕD(x−x^)dxPX^(dx^).

M 是边缘乘积;QD 保留同一个重构边缘,但把给定重构后的源改成均值 x^、方差 D 的高斯。两者相对于 dx⊗PX^ 都有严格正密度,因此彼此等价,即使 PX^ 没有Lebesgue密度也没有问题。

若 I(X;X^)=+∞,所需下界已成立。否则 J≪M≪QD,代入一元正态分布的密度公式得到的对数比

log2⁡dMdQD(x,x^)=12log2⁡Dσ2+(x−x^)22Dln⁡2−x22σ2ln⁡2

在 J 下绝对可积,因为 δ<∞、EX2=σ2。有限KL的对数密度比也绝对可积:其负部积分由 1/(eln⁡2) 控制,正部则由KL有限性控制。因此可以合法相加并取期望,得到

I(X;X^)=D2(J‖QD)+12log2⁡σ2D+D−δ2Dln⁡2.

KL非负性与 δ≤D 给出 I≥12log2⁡(σ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。此时联合分布就是 QD;它与边缘乘积的对数比为上面二次式的相反数,在联合分布下绝对可积,故互信息有限。精确差额中的两项均为零,于是

I(X;X^)=12log2⁡σ2D.

虽然构造从重构向源抽样,但联合分布确定了合法的正向核。由二维高斯分布的条件均值和方差,它等价于

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−H2(D)。D=0 恢复无损压缩率 1 bit/符号,D=1/2 时无需发送信息也可随机猜测达到允许失真。

方差4、失真1:直接加噪声为什么多耗信息 ​

最优正向核是 X^=34X+V,其中 V∼N(0,3/4) 独立于 X∼N(0,4)。逐项复算为

E(X−X^)2=(14)24+34=1,Var(X^)=(34)24+34=3,I(X;X^)=12log2⁡33/4=1 bit/源符号.

作为比较,取 X~=X+W,W∼N(0,1) 独立于源,也有失真 1,但

I(X;X~)=12log2⁡51≈1.160964 bit/源符号.

两个互信息都可用有限高斯密度的积分计算。最优核在加入正向噪声前先收缩信号;只检查噪声方差等于目标失真,无法判断互信息是否最小。

连续源的两个端点 ​

D≥σ2 时恒取 X^=0,失真为 σ2,互信息为零。D=0 时必有 X^=X 几乎处处:联合分布把全部质量放在对角线,而非原子的高斯边缘乘积给对角线质量零。联合分布不绝对连续于边缘乘积,所以互信息为无穷。连续源精确重构的端点与有限字母表Hamming零失真不同。

前面的正向公式用于 0<D<σ2;D=σ2 应直接使用常量重构,不能把退化的条件方差继续代入高斯密度的对数比。

率失真定理是长块、平均失真意义下的极限,不保证每个样本都低于 D,也不提供有限块长、有限时延或低复杂度编码器。实际 codec 还受模型族、算力和缓冲时延约束;若失真函数没有表达人的感知差异,数学上的最优重构也未必主观最好。

有记忆源或一般源通常不再由这个单字母式完整描述。此时需要研究 n 维分布与块失真的多字母优化,再取归一化极限;更一般的非平稳情形还可能使用信息谱表述。连续字母表也需要可测性、矩条件和失真可积性等假设,不能把有限字母表定理删去条件后直接沿用。

推论与应用

率失真理论给出图像、音频和有损压缩的基准极限,并连接量化、感知质量和信息瓶颈。互信息 是单字母优化目标,熵 给出离散无记忆源在零 Hamming 失真处的无损端点;与信道容量结合后,它还能判断给定源与信道是否存在满足目标失真的分离式传输方案。

把每源符号与每信道使用的单位对齐 ​

设 Xn 是方差 σ2 的IID高斯源;编码器把它映成 m 个实信道输入 Um,通过 Yj=Uj+Zj,其中 N>0,噪声 Zj∼N(0,N) 独立同分布并独立于源与编码器随机性。无反馈,接收端由 Ym 重构 X^n。固定 0≤P<∞,功率约束现在按实际源分布取期望:

1m∑j=1mEUj2≤P.

只需考察平均失真有限的方案,否则以下失真下界自动成立。令 δi=E(Xi−X^i)2,δ¯=n−1∑iδi。率失真的凸性、定义、源的独立性、数据处理不等式以及高斯信道容量的块逆界依次给出

nR(δ¯)≤∑i=1nR(δi)≤∑i=1nI(Xi;X^i)≤I(Xn;X^n)≤I(Um;Ym)≤mC(P).

其中容易看反的第三步可由链式法则核对:因为 Xi 与 Xi−1 独立,

I(Xn;X^n)=∑iI(Xi;X^n,Xi−1)≥∑iI(Xi;X^i).

块信道的上界有限,也保证这里相关的互信息有限。代入两个闭式公式,若 δ¯<σ2,得

δ¯≥σ2(1+PN)−m/n.

若 δ¯≥σ2,该不等式自动成立。固定资源比

α=limn→∞mn∈(0,∞)(信道使用次数/源符号),

最优渐近期望失真因而至少为 D∗(α)=σ2(1+P/N)−α。

可达性、严格余量与无界失真 ​

当目标满足 R(D)<αC(P) 且 0<D<σ2 时,由连续性先取略小的 D′<D,仍使 R(D′)<αC(P),为拼接误差留下失真余量。再选压缩速率略大于 R(D′)、每信道使用的码率略小于 C(P),拼接率失真码和信道码。源索引未必均匀,所以信道码使用最大消息错误概率保证;每个信道码字都满足能量约束,故也满足实际源分布下的功率要求。

平方误差无界,还需控制信道译码出错时的失真。连续源编码/分离定理允许本高斯模型选取重构码字能量统一不超过 nK 的码本,其中 K 与块长无关。若信道最大错误概率为 ϵ,记错误事件为 E,由 Pr(E∣Xn=xn)≤ϵ 和平方三角不等式,

1nE[‖Xn−X^n‖21E]≤2nE[‖Xn‖21E]+2KPr(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σ2X,Y=U+Z,X^=Pσ2P+NY.

逐符号便有 EU2=P,而

X−X^=NP+NX−Pσ2P+NZ,E(X−X^)2=N2σ2+Pσ2N(P+N)2=σ2NP+N=D∗(1).

它诱导的正向系数 P/(P+N) 与噪声方差 Pσ2N/(P+N)2,恰好对应前面取 D=σ2N/(P+N) 的最优测试信道。数值例 P=3,N=1,σ2=4 中,发送和接收均乘 3/2,合成后为 X^=34X+32Z,失真正好为1。

这种高斯匹配解释了一个特殊的有限延时最优方案。其功率按源分布取期望;当 P>0 时,高斯源幅度无界,所以它不满足“每个可能源块都对应一个能量至多 mP 的发送块”。把功率量词改强后,不能继续声称该逐符号方案原样可行。

若译码器另有相关观察,Wyner–Ziv 编码在辅助变量 U−X−Y 的限制下最小化条件互信息,明确区分编码器是否也看到边信息。压缩转发进一步把这种描述送过一条中继信道;描述所需速率必须不超过实际可用链路,不能仅算压缩质量而省掉传送费用。

参考资料
  • 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:渐近闭包口径及高斯源的逐符号最优缩放方案;其资源比率为本页 α 的倒数。
关系图谱24 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系