DennyQi's Log

01 Information Entropy 2

“信息(information)”是什么呢?提到信息,我们会联想到填表时的个人信息表格,联想到计算机中数字的二进制表示。在日常语言中,我们会这样使用“信息”这一词汇:这一句话是信息量很大的;这是一句废话(信息量几乎为0);……在自然语言中,“信息”一词和“文字”、“符号”、“句子”、“语言”、“特征”等词语相关联,但并没有精确回答“什么是信息”。

二十世纪人类最伟大的成就之一就是在数学上定义了什么是信息。这一定义极大地改变了人类的思维方式,创造出了许多前所未有的高效工具,使得许多概念得到精确地重新定义。

所有关于“信息”的现象都和人与人(主体与主体)之间的通信(communication)有关。信息产生于人与人的交流——一方将自己知道的告知(inform)另一方。一方要将“已知信息”告知另一方时,必须将这种“已知信息”通过某种物理过程作用在接收方身上。这种物理作用可以是呼喊(声波,也就是介质的振动),电磁波等等。更奇妙的是,物理过程本身并不重要,而是那些能够从物理过程中抽象出来的模式(pattern)构成了信息。只要发送方和接收方达成了某一统一意见(称为“协议(protocol)”),发送方就可以把想要发送的内容“编码(encode)”为某一物理模式,接收方可以通过“解码(decode)”把信息从物理模式中提取出来。这样一来,关于“信息”的讨论就可以完全独立于物理过程,而在一个抽象的数学语言中得到描述。

“计算机”恰好就是处理信息的最好工具,它本质上是对一系列存储在内部的0和1的物理信号的处理。这就是为什么,当我们提到“信息科学”经常指的是“计算机科学”,事实上这二者还是有区别。对信息科学的研究很多时候可以独立于物理学的研究,这是不同层级学科之间的独立性的一个最显著案例。物理学的发展会促进计算机的存储空间、计算效率、体系结构(传统图灵机还是量子计算机)等。同时,任何学科都可以采用信息科学地思维方式。

可能出人意料的是,“信息”这一概念的精确定义和“概率”紧密联系在一起。信息的“获得”就是把我们认知中不确定性的东西转化为确定的东西。我最初听到这一说法时并不理解。令我印象深刻的是,有一天我坐在地铁上,前面站着一个背对我的人,我想象那个人的相貌,产生了各种可能性,而就在这一瞬间那个人转过了头来,一瞬间我领悟了“什么是信息”——原本我对于那个人的相貌存在各种想象,而当那个人转过头来时我终于确定了那个人的相貌究竟是什么,也即“我获得了关于那个人的容貌的信息”。我们生活中也会有这样的体验,当预感到即将“获知信息”时我们会感到紧张:当我们等待高考成绩发送到手机短信时,我们会紧张,因为此时我们不知道我们会收到令人满意的高考成绩还是令人沮丧的高考成绩,当手机接收到短信的一瞬间,未知转化为了已知。这种“把不确定转化为确定”“把未知转化为已知”的过程是信息的本质信息的本质是这一转化的过程,而不是完成这个转化的符号。比如说,我们收到了高考总分是650分。我们可以说,我们收到的信息就是数字“650”吗?事实上重要的不仅是“数字是650”。“数字不是651”、“数字不是700”、“数字不是600”都很重要。在我们收到成绩前,我们脑中可能会有一个关于分数预计的概率分布,但是在收到成绩以后我们得到了样本空间中一个具体的值,不再是一个概率分布了(坍缩了!)。这就是“信息”和“概率”的紧密联系:信息是概率空间的缩小。

信息与随机事件

有了以上认知以后,我们可以用数学语言构建“信息”这一概念。

考虑一个最简单的“通信”的场景:A和B分别位于不同的房间内。A会在它的房间进行一次抛均匀硬币的实验。由于B在另一个房间,除非A把实验结果“告知”B,B始终处在“硬币朝上、硬币朝下”这二者的不确定性之间。在实验前,A和B达成协议,用一张只能写0或1的字符串的纸条传递实验结果。如果正面朝上,纸条上写0;如果反面朝上,纸条上写1。通过传递纸条这一过程,A和B就硬币的抛掷结果完成了通信。在此通信期间,纸条上的数字0或1传递了实验结果的信息。在这个实验中,无论如何,我们只需要一个二进制位(bit)就能传递“抛均匀硬币的实验结果”这一信息。

如果场景稍微复杂一些,光靠一个二进制位就不足以完成信息传递了。假设A在房间内抛一个八面体的均匀骰子,那么就有八种可能结果,因此要把信息传递给B,至少需要33个二进制位(23=82^3=8)。如果少于3个二进制位,B永远都不可能推断出A的实验结果究竟是什么。

抛硬币只需要11个二进制位就足够传递信息,而抛骰子却需要33个二进制位。所以我们说,“抛骰子的结果”的信息量要比“抛硬币的结果”的信息量更大。为什么后者必须要用更多的位数呢?因为抛骰子这一随机事件比抛硬币这一随机事件的“不确定性(uncertainty)”更大。抛硬币时,B事先知道结果只有两种,A在通知它结果以后只帮B排除了一种情况;抛骰子时,B事先知道结果可能有8种,A通知它以后立即帮B排除了其余的7种情况。抛骰子这一随机事件的结果的不确定性比抛硬币更大。得知一个不确定性更大的事件的结果意味着得知更多的信息。由此可见,信息量的大小来源于所消除的不确定性的大小。

按照这样的思路,我们很容易把信息量的定义推广到一般情形。对于任何一个随机事件,其结果的可能情况越多,得知其结果时所获得的信息量就越大。假如一个随机事件有nn种可能的结果,那么描述其结果至少需要log2n\log_2 n个二进制位。

熵(Entropy)

然而我们认识到,以上的对信息的描述方式并不精确。按照以上说法,信息似乎只与随机变量的“取值个数”有关,而并不与具体的“分布”有关。比如,如果我们抛的是一个带有倾向性的硬币——正面概率为60%,反面概率为40%。按照以上说法,在得知抛这枚带有倾向性的硬币的结果时,我们所获得的信息量是和抛均匀硬币时相同的。当我们把这种“倾向性”推广到极端,就会立刻意识到不对劲:假如某一硬币正面向上的概率为99.99%,反面向上的概率为0.01%,那么我们在得知其结果的时候会如何反应呢?如果我们得知结果为“正面向上”,我们会觉得这几乎是一句废话,信息量几乎为0;但是如果我们得知结果为“反面向上”,我们会非常惊讶,好像发生了一件非常神奇的事情,换言之这一结果包含很大的信息量。就好像,对于一个常年干旱的地区,如果你告诉那里的居民“明天不会下雨”,居民会觉得你这句话没任何价值;但是如果你说“明天会下雨”,居民就会认为这句话有很大的信息量。

我们应当定义一个关于分布的函数,作为对信息量的数学描述。我们对这一描述会有一些基本的要求:

  • 首先,对于均匀分布,信息量应当等于可能情况的对数。
  • 其次,从上面的例子不难想到,这应当是一个连续的函数——这正是我们用“极端法”论证的基础,既然100%正面向上的硬币的结果毫无信息可言,那么99%正面向上的硬币的结果信息量也不会太大。
  • 最后,信息量还应当满足链式法则:设分布{p1,,pn}\{p_1,\cdots,p_n\}对应的信息量为H(p1,,pn)H(p_1,\cdots,p_n),考虑一个不均匀分布的随机变量XXX=aX=a的概率为1/21/2X=bX=b的概率为1/31/3X=cX=c的概率为1/61/6。那么,XX的分布为{1/2,1/3,1/6}\{1/2,1/3,1/6\},这个分布对应的信息量为H(1/2,1/3,1/6)H(1/2,1/3,1/6),我们可以这样理解信息量:当我们想要获取XX的取值时,首先询问是否成立X=aX=a,这个问题的答案的信息量为H(1/2,1/2)H(1/2,1/2);如果答案为否(这一事件发生的概率为1/21/2),那么我们再次询问是否成立X=bX=b,这个问题的答案的信息量为H(2/3,1/3)H(2/3,1/3)。所以,HH应当满足以下分解性质:H(1/2,1/3,1/6)=H(1/2,1/2)+1/2H(2/3,1/3)H(1/2,1/3,1/6)=H(1/2,1/2)+1/2H(2/3,1/3)

综上所述,我们想找的函数H(p1,,pn)H(p_1,\cdots,p_n)要满足以下三条基本性质:

  1. H(1/n,,1/n)=log2nH(1/n,\cdots,1/n)=\log_2 n
  2. HH是连续函数;
  3. H(p1,,pn1,αpn,(1α)pn)=H(p1,,pn1,pn)+pnH(α,1α)H(p_1,\cdots,p_{n-1},\alpha p_n,(1-\alpha)p_n)=H(p_1,\cdots,p_{n-1},p_n)+p_nH(\alpha,1-\alpha)

我们来推导这三条性质的必要条件:对于任何一个离散的分布{p1,,pN}\{p_1,\cdots,p_N\},如果pip_i都是有理数,我们总可以设存在整数n1,,nNn_1,\cdots,n_N使得pi=nin1++nNp_i=\dfrac{n_i}{n_1+\cdots+n_N}。设S=n1++nNS=n_1+\cdots+n_N,所以分布{p1,,pN}\{p_1,\cdots,p_N\}可以分解为SS个样本的均匀分布,并且该分解满足性质3。这意味着,log2(S)=H(p1,,pn)+pilog2(ni)\log_2(S)=H(p_1,\cdots,p_n)+\sum p_i\log_2(n_i)。这就导出了H(p1,,pn)=pilog2(ni)+log2(S)H(p_1,\cdots,p_n)=-\sum p_i\log_2(n_i)+\log_2(S) =pilog2(ni)+pilog2(S)=-\sum p_i\log_2(n_i)+\sum p_i\log_2(S) =pilog2(niS)=-\sum p_i\log_2(\dfrac{n_i}{S}) =pilog2pi=-\sum p_i\log_2 p_i

于是,我们找到了这个著名的表达式

H(p1,,pn)=pilog2piH(p_1,\cdots,p_n)=-\sum p_i\log_2 p_i

经检验,它满足我们提出的三个要求。这一结果最早出现在香农(Shannon)在1948年发表的名为A Mathematical Theory of Communication的论文中。文中把这一函数和统计力学中的Boltzmann定理联系起来,把这个函数命名为entropy(熵)。由于这是用于刻画信息量的熵,所以称为信息熵。

因为熵是关于离散概率分布的函数,所以对于任何离散随机变量XX,我们可以定义随机变量的信息熵。假设XX概率质量函数为p(x)p(x),其中p(x)=Pr[X=x]p(x)=\Pr[X=x],则定义XX的熵为H(X)=xXp(x)logp(x)H(X)=-\sum\limits_{x\in\mathcal{X}}p(x)\log p(x)。其中X\mathcal{X}是所有p(x)>0p(x)>0xx构成的集合。

现在我们来看我们导出的这一表达式的含义。我们可以把熵写成这样的形式:H(p1,,pn)=pilog1piH(p_1,\cdots,p_n)=\sum p_i \log \dfrac{1}{p_i}。这很容易写作期望的形式:H(p1,,pn)=E[log1pi]H(p_1,\cdots,p_n)=\mathbb{E}[\log \dfrac{1}{p_i}]log1pi\log \dfrac{1}{p_i}这一项的含义是什么呢?根据之前举的例子,这就可以看作是随机事件的每种结果带给我们的“惊讶程度”:当我们抛一枚99%正面朝上的硬币时,正面朝上带给我们的惊讶程度为log10099\log \dfrac{100}{99},这是一个很接近0的小量,我们对这个结果丝毫不感到惊讶;反面朝上带给我们的惊讶程度为log10016.64\log \dfrac{100}{1}\approx 6.64。最终,整个事件的信息量是各个情况的惊讶程度的期望,也就是我们的平均惊讶程度。尽管6.646.64的惊讶值很高,但只有11%的概率会发生。最终,我们的平均惊讶程度为0.01×6.64+0.99×0.010.07630.01\times 6.64+0.99\times 0.01\approx 0.0763

对于一个随机变量XX,当其均匀分布时,它就具有最大的熵。只要分布稍不均匀,熵就会降低。我们可以证明这一点:均匀分布时,H(X)=logXH(X)=\log |\mathcal{X}|。而对于H(X)=xXp(x)log1p(x)H(X)=\sum\limits_{x\in \mathcal{X}} p(x)\log \dfrac{1}{p(x)},我们把它看作在上凸函数logx\log x上分别以p(x)p(x)的权重选取xi=1p(x)x_i=\dfrac{1}{p(x)},由Jensen不等式可得H(X)log(xXp(x)1p(x))=logXH(X) \leq \log\left(\sum\limits_{x\in \mathcal{X}}p(x) \cdot \dfrac{1}{p(x)}\right)=\log |\mathcal{X}|。当且仅当p(x)=1Xp(x)=\dfrac{1}{| \mathcal{X}|}时取到等号。由此可见,离散的随机变量的熵始终满足0H(X)logX0\leq H(X)\leq \log |\mathcal{X}|,均匀分布时取到最大值。

联合熵(Joint Entropy)

两个随机变量的联合分布可以导出“联合熵”。这是很自然的,因为熵是一个仅仅关于分布的函数,只需要一系列离散的概率密度就可以定义。设X,YX,Y有联合分布的密度函数p(x,y)p(x,y),那么定义X,YX,Y的联合熵为H(X,Y)=xXyYp(x,y)logp(x,y)H(X,Y)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in \mathcal{Y}}p(x,y)\log p(x,y)。从期望的角度,H(X,Y)=E[logp(X,Y)]H(X,Y)=-\mathbb{E}[\log p(X,Y)]。事实上,我们可以把(X,Y)(X,Y)看作一个整体(一个随机向量),那么X,YX,Y联合分布的概率密度实际就是这单个随机向量的概率分布,它衡量这个随机向量(另一个新的随机变量)的不确定性。从对称性容易看出,H(X,Y)=H(Y,X)H(X,Y)=H(Y,X)

容易验证,如果X=YX=Y,那么p(x,y)>0p(x,y)>0当且仅当x=yx=yp(x,x)=p(x)p(x,x)=p(x),代入定义式可得H(X,X)=xXyXp(x,y)logp(x,y)=xXp(x,x)logp(x,x)H(X,X)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in \mathcal{X}}p(x,y)\log p(x,y)=-\sum\limits_{x \in \mathcal{X}}p(x,x)\log p(x,x) xXp(x)logp(x)=H(X)-\sum\limits_{x\in\mathcal{X}}p(x)\log p(x)=H(X)。所以,两个相同的随机变量的联合熵就等于单个随机变量的熵。从信息量的角度,增加一个相同的随机变量并没有增加信息量。

如果XXYY的函数,也即YY确定时XX会被唯一确定,那么yy确定时使得p(x,y)>0p(x,y)>0的只有唯一的xx,因此p(x,y)=p(y)p(x,y)=p(y)。于是代入定义可得H(X,Y)=xXyYp(x,y)logp(x,y)H(X,Y)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in \mathcal{Y}}p(x,y)\log p(x,y) =yYp(y)logp(y)=H(Y)=-\sum\limits_{y \in \mathcal{Y}}p(y)\log p(y)=H(Y)XX的信息完全被包含在YY以内,因此增加XX并不能带来更多的信息。

如果X,YX,Y独立,那么p(x,y)=p(x)p(y)p(x,y)=p(x)p(y)。那么H(X,Y)=xXyYp(x)p(y)[logp(x)+logp(y)]H(X,Y)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in \mathcal{Y}}p(x)p(y)[\log p(x)+\log p(y)] =xXp(x)logp(x)yYp(y)xXp(x)yYp(y)logp(y)=-\sum\limits_{x \in \mathcal{X}}p(x)\log p(x)\sum\limits_{y \in \mathcal{Y}}p(y)-\sum\limits_{x \in \mathcal{X}}p(x)\sum\limits_{y \in \mathcal{Y}}p(y)\log p(y) =H(X)+H(Y)=H(X)+H(Y)。两个独立的随机变量的熵恰好是它们熵的和。X,YX,Y中并没有互相重叠的信息。

联合熵可以继续推广到多元:定义H(X1,,Xn)=p(x1,,xn)logp(x1,,xn)H(X_1,\cdots,X_n)=-\sum p(x_1,\cdots,x_n)\log p(x_1,\cdots,x_n) =E[logp(X1,,Xn)]=-\mathbb{E}[\log p(X_1,\cdots,X_n)]

条件熵(Conditional Entropy)

由随机变量的条件分布可以导出条件熵。对于两个离散随机变量X,YX,Yp(YX=x)p(Y\mid X=x)依然是一个概率分布,由此定义H(YX=x)=yYp(yX=x)logp(yX=x)H(Y\mid X=x)=-\sum\limits_{y\in \mathcal{Y}}p(y\mid X=x)\log p(y\mid X=x)。从期望的角度,可以写作E[logp(yX=x)]-\mathbb{E}[\log p(y\mid X=x)]。基于H(YX=x)H(Y\mid X=x),定义X,YX,Y的条件熵H(YX)=xXp(x)H(YX=x)H(Y \mid X)=\sum\limits_{x \in \mathcal{X}}p(x)H(Y\mid X=x)。它表示已知XXYY的不确定性,而“已知XX”是期望意义下的已知。展开H(YX=x)H(Y\mid X=x)这一项,得到H(YX)=xXyYp(x)p(yx)logp(yx)H(Y\mid X)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in\mathcal{Y}}p(x)p(y\mid x)\log p(y\mid x)。而p(x)p(yx)=p(x,y)p(x)p(y\mid x)=p(x,y),因此得到条件熵的一般表达式H(YX)=xXyYp(x,y)logp(yx)H(Y\mid X)=-\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in\mathcal{Y}}p(x,y)\log p(y\mid x) =E[logp(YX)]=-\mathbb{E}[\log p(Y\mid X)]

注意,H(XY)H(X\mid Y)一般不等于H(YX)H(Y \mid X)。但可以证明:H(XY)+H(Y)=H(YX)+H(X)=H(X,Y)H(X\mid Y)+H(Y)=H(Y\mid X)+H(X)=H(X,Y)。这称为熵的计算的链式法则。这可以从概率的链式法则p(x,y)=p(xy)p(y)p(x,y)=p(x\mid y)p(y)直接导出:从期望的角度,H(X,Y)=E[logp(X,Y)]=E[logp(XY)+logp(Y)]H(X,Y)=-\mathbb{E}[\log p(X,Y)]=-\mathbb{E}[\log p(X\mid Y)+\log p(Y)] =H(XY)+H(Y)=H(X\mid Y)+H(Y)。另一个是对称的。推广到nn元情形:H(X1,,Xn)=i=1nH(XiXi1,,X1)H(X_1,\cdots,X_n)=\sum\limits_{i=1}^{n}H(X_i\mid X_{i-1},\cdots,X_1)

同样的,根据条件概率的定义容易验证p(x,yz)=p(xz)p(yx,z)p(x,y\mid z)=p(x\mid z)\cdot p(y\mid x,z)。用同样的方法可以证明H(X,YZ)=H(XZ)+H(YX,Z)H(X,Y\mid Z)=H(X\mid Z)+H(Y\mid X,Z)

XXYY的函数时,H(X,Y)=H(Y)H(X,Y)=H(Y)。而H(X,Y)=H(Y)+H(XY)H(X,Y)=H(Y)+H(X\mid Y),可见此时H(XY)=0H(X\mid Y)=0YY已知时XX没有任何不确定性。而反过来,如果H(XY)=0H(X\mid Y)=0,那么xXyYp(x,y)logp(xy)=0\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in\mathcal{Y}}p(x,y)\log p(x\mid y)=0,这当且仅当p(xy)p(x\mid y)恒等于1,也即yy确定xx确定,XXYY的函数。综上我们得到,X=f(y)    H(XY)=0X=f(y)\iff H(X\mid Y)=0

Mutual Information(互信息)

比较H(YX)H(Y\mid X)H(Y)H(Y)的大小,从直观上,“X已知”本身提供了信息,这一信息势必会使得YY的不确定性降低,或至少不会让YY变得更不确定。因此应当成立不等式H(YX)H(Y)H(Y \mid X)\leq H(Y)。什么时候成立等号呢?代入H(YX)=H(X,Y)H(X)H(Y\mid X)=H(X,Y)-H(X),等号成立时H(X,Y)=H(X)+H(Y)H(X,Y)=H(X)+H(Y)。我们先前验证了,如果X,YX,Y是独立的,那么这个等式就成立。直观上,这个不等式(也即差值H(Y)H(YX)H(Y)-H(Y\mid X))在衡量随机变量X,YX,Y之间距离独立还有多远。我们定义这个差值为X,YX,Y的互信息I(X;Y)=H(Y)H(YX)I(X;Y)=H(Y)-H(Y\mid X)(或对称的I(X;Y)=H(X)H(XY)I(X;Y)=H(X)-H(X\mid Y))。代入化简可得I(X;Y)=xXyYp(x,y)logp(x,y)p(x)p(y)I(X;Y)=\sum\limits_{x \in \mathcal{X}}\sum\limits_{y \in \mathcal{Y}}p(x,y)\log \dfrac{p(x,y)}{p(x)p(y)}。互信息具有对称性:I(X;Y)=I(Y;X)I(X;Y)=I(Y;X)。从信息的角度,它描述X,YX,Y之间有多少共同的信息。如果没有共同的信息(独立),那么互信息为0。

事实上,表达式xXp(x)logp(x)q(x)\sum\limits_{x \in \mathcal{X}}p(x)\log \dfrac{p(x)}{q(x)}是一种用来衡量分布之间“距离”的一般方式,它称为Kullback-Leibler距离,记为D(p(x)q(x))D(p(x)||q(x)),又称为分布分别为p,qp,q的两个随机变量的相对熵(Relative Entropy)。互信息可以用KL距离写作I(X;Y)=D(p(x,y)p(x)p(y))I(X;Y)=D(p(x,y)||p(x)p(y))。(注意,Kullback-Leibler距离是不具有对称性的)

下面我们证明,始终成立D(p(x)q(x))0D(p(x)||q(x))\geq 0。这是信息论中最重要的不等式之一,称为信息不等式(Information Inequality)。根据定义,D(p(x)q(x))=xXp(x)logp(x)q(x)=xXp(x)logq(x)p(x)D(p(x)||q(x))=\sum\limits_{x \in \mathcal{X}}p(x)\log \dfrac{p(x)}{q(x)}=-\sum\limits_{x \in \mathcal{X}}p(x)\log \dfrac{q(x)}{p(x)}。由于log\log是上凸函数,根据Jensen不等式有xXp(x)logq(x)p(x)log(xXp(x)q(x)p(x))\sum\limits_{x \in \mathcal{X}}p(x)\log \dfrac{q(x)}{p(x)}\leq\log \left(\sum\limits_{x \in \mathcal{X}} p(x)\cdot \dfrac{q(x)}{p(x)}\right) =log1=0=\log 1 = 0。因此D(pq)0D(p||q)\geq 0。由于Jensen不等式只在所有点都重合时取等,因此当且仅当p,qp,q为同一分布时D(pq)=0D(p||q)=0。由I(X;Y)=D(p(x,y)p(x)p(y))I(X;Y)=D(p(x,y)||p(x)p(y)),可得I(X;Y)0I(X;Y)\geq 0。信息不等式表明,互信息始终是非负的!

作为例子,我们取qq为均匀分布,也即q(x)1Xq(x)\equiv \dfrac{1}{|\mathcal{X}|},那么D(pq)=xp(x)logp(x)+xp(x)logXD(p||q)=\sum\limits_{x}p(x)\log p(x)+\sum\limits_xp(x)\log |\mathcal{X}| =logXH(X)=\log|\mathcal{X}|-H(X)。由于D(pq)0D(p||q)\geq 0,这再次表明H(X)H(X)只能在均匀分布时取到最大值logX\log |\mathcal{X}|。后续在微分熵中,这是更普适的证明方法。

信息图(The Information Diagram)

X,YX,Y的熵、联合熵、互信息始终满足H(X,Y)=H(X)+H(Y)I(X;Y)H(X,Y)=H(X)+H(Y)-I(X;Y)。这意味着,我们可以用韦恩图来理解熵与互信息的关系:H(X),H(Y)H(X),H(Y)是单个圆的面积,H(X,Y)H(X,Y)是并集的面积,而I(X;Y)I(X;Y)是交集的面积。H(XY)H(X\mid Y)YY去掉XX部分的面积,H(YX)H(Y\mid X)XX去掉YY部分的面积。

image-20240304165448878

我们可以把信息图推广到以下的三元情形:

image-20240304165358837

在熵与联合熵中,只会涉及逗号与竖线,其中逗号表示对两块面积取并,竖线表示去除对应部分的面积,逗号的优先级高于竖线。在互信息中,会出现分号,其中分号表示对两块面积取交,分号的优先级高于竖线,低于逗号。综合起来,优先级从高到低为,>;>, > ;>|。我们可以验证,根据信息图做恒等变形始终是成立的。其中,为I(X;YZ)=H(XZ)H(XY,Z)I(X;Y\mid Z)=H(X\mid Z)-H(X\mid Y,Z)称为条件互信息(Conditional Mutual Information)。条件互信息也具有非负性。I(X1;X2;X3)I(X_1;X_2;X_3)仅仅是一个形式上的记号,它并不是互信息,不具有非负性(它也是信息图中唯一可能取负值的一片区域。我们可以证明,当X=Y=ZX=Y=ZI(X;Y;Z)>0I(X;Y;Z)>0,而Z=X+YZ=X+YI(X;Y;Z)<0I(X;Y;Z)<0。)

信息图为我们完整描述了三元以内的所有熵与互信息之间的等式关系。而对于三元以上的信息图,我们不能找到一个把它在平面上画出来的简单方式。但是可以验证,以下基于信息图的理解经过验证在三元以上的情形也是正确的:

①熵的链式法则:H(X1,,Xn)=i=1nH(XiX1,,Xi1)H(X_1,\cdots,X_n)=\sum\limits_{i=1}^{n}H(X_i\mid X_1,\cdots,X_{i-1})(一系列面积取并,等价于每次累加一个新面积去除已经计算过的所有面积);

②互信息的链式法则:I(X1,,Xn;Y)=i=1nI(Xi;YXi1,,X1)I(X_1,\cdots,X_n;Y)=\sum\limits_{i=1}^{n}I(X_i;Y\mid X_{i-1},\cdots,X_1)(一系列面积与另一个面积取交,等价于每次累加一个面积与它的交去除所有已经计算过的部分);

③互信息与熵的转化:I(X1,,Xn;Y)=H(X1,,Xn)H(X1,,XnY)I(X_1,\cdots,X_n;Y)=H(X_1,\cdots,X_n)-H(X_1,\cdots,X_n\mid Y)(将(X1,,Xn)(X_1,\cdots,X_n)看作一个随机向量);

H(XY)H(X)H(X\mid Y)\leq H(X)中的X,YX,Y看作随机向量推广到多元,可以验证不等式依然成立。那么基于熵的链式法则H(X1,,Xn)=i=1nH(XiX1,,Xi1)H(X_1,\cdots,X_n)=\sum\limits_{i=1}^{n}H(X_i\mid X_1,\cdots,X_{i-1}),可以得到以下不等式,称为The Independence Bound: H(X1,,Xn)i=1nH(Xi)H(X_1,\cdots,X_n)\leq\sum\limits_{i=1}^{n}H(X_i)。这直观上表明nn个随机变量联合熵总是不超过各自熵的和。这种系统间的相互影响(重叠的信息)而造成的。如果nn个变量全都互相独立,那么恰好取到等号。这个不等式可以看作信息不等式的一个推论。信息图中真正本质的不等关系只有信息不等式一个(而它的本质是Jensen不等式)。

同样的,基于D(pq)=Ep[logp(x)q(x)]D(p||q)=\mathbb{E}_p\left[\log \dfrac{p(x)}{q(x)}\right],可以定义条件相对熵(Conditional Relative Entropy) D(p(yx)q(yx))=Ep(x,y)[logp(YX)q(YX)]D(p(y\mid x)||q(y\mid x))=\mathbb{E}_{p(x,y)}\left[\log \dfrac{p(Y\mid X)}{q(Y\mid X)}\right] =xyp(x,y)logp(yx)q(yx)=\sum\limits_{x}\sum\limits_{y}p(x,y)\log \dfrac{p(y\mid x)}{q(y\mid x)}

对于D(p(x,y)q(x,y))D(p(x,y)||q(x,y)),会出现logp(x,y)q(x,y)\log \dfrac{p(x,y)}{q(x,y)}一项,根据条件概率可以展开为logp(x)p(yx)q(x)q(yx)=logp(x)q(x)+logp(yx)q(yx)\log \dfrac{p(x)p(y\mid x)}{q(x)q(y\mid x)}=\log \dfrac{p(x)}{q(x)}+\log \dfrac{p(y\mid x)}{q(y\mid x)}。因此D(p(x,y)q(x,y))=D(p(x)q(x))+D(p(yx)q(yx))D(p(x,y)||q(x,y))=D(p(x)||q(x))+D(p(y\mid x)||q(y\mid x))。这是相对熵的链式法则。

马尔科夫链(Markov Chain)

一般来说,根据链式法则,三个随机变量的分布满足p(x,y,z)=p(x)p(yx)p(zx,y)p(x,y,z)=p(x)p(y\mid x)p(z\mid x,y)。假如我们发现分布可以进一步满足p(x,y,z)=p(x)p(yx)p(zy)p(x,y,z)=p(x)p(y\mid x)p(z\mid y),也即在X,Y,ZX,Y,Z的联合分布中ZZ总是只依赖于YY而不依赖于XX,就称这三个随机变量形成了马尔可夫链XYZX\to Y \to Z。一个很常见的情形是,Z=f(Y)Z=f(Y),此时自然有XYf(Y)X \to Y \to f(Y)

对于XYZX\to Y\to Z,根据马尔可夫链定义,由p(x,zy)=p(x,y,z)p(y)p(x,z\mid y)=\dfrac{p(x,y,z)}{p(y)} =p(x)p(yx)p(zy)p(y)=\dfrac{p(x)p(y\mid x)p(z\mid y)}{p(y)} =p(x,y)p(zy)p(y)=p(xy)p(zy)=\dfrac{p(x,y)p(z\mid y)}{p(y)}=p(x\mid y)p(z\mid y)。这说明,马尔可夫链等价于在中间随机变量的条件概率意义下,前后的两个事件是独立的。而这样的定义是对称的,因此XYZX\to Y\to Z一定同时意味着ZYXZ\to Y \to X.

对于马尔可夫链,变量之间的互信息满足以下重要的不等式,称为数据处理不等式(Data-processing inequality):如果XYZX\to Y \to Z,那么I(X;Y)I(X;Z)I(X;Y)\geq I(X;Z)。 它表明,在马尔可夫链中相距更近的两个变量之间的关联一定比更远的变量更紧密。仅仅通过处理YY的数据来得到的变量ZZ不可能帮助我们获得更多信息。证明如下:由于XYZX\to Y\to Z,因此在YY的条件下X,ZX,Z独立,那么有I(X;ZY)=0I(X;Z\mid Y)=0。根据链式法则,I(X;Y,Z)=I(X;Z)+I(X;YZ)I(X;Y,Z)=I(X;Z)+I(X; Y\mid Z),对称的也有I(X;Y,Z)=I(X;Y)+I(X;ZY)=I(X;Y)I(X;Y,Z)=I(X;Y)+I(X;Z\mid Y)=I(X;Y)。因为I(X;YZ)0I(X;Y\mid Z)\geq 0,因此I(X;Z)I(X;Y)I(X;Z)\leq I(X;Y)

根据I(X;Y)=H(X)H(XY)I(X;Y)=H(X)-H(X\mid Y)I(X;Z)=H(X)H(XZ)I(X;Z)=H(X)-H(X\mid Z),数据处理不等式也可以等价地写为H(XY)H(XZ)H(X\mid Y)\leq H(X\mid Z)。这说明给定一个马尔可夫链上间隔越远的已知条件,对不确定性的约束效果更弱。

在马尔可夫链中,由于I(X;ZY)I(X;Z\mid Y)恒为0,因此我们不需要韦恩图中(XZ)Y(X\cap Z)\setminus Y那一片区域,因此可以把韦恩图画成三个山峰的形式。多元的马尔可夫链的韦恩图也是类似的(要保证X1X_1XnX_n相交)。

image-20240304165319733
image-20240304165336326