信息论基础

主要内容:信息论基本概念、信源模型、自信息、互信息、信息熵、马尔可夫信源等

信息

信息:所获得的新知识,信息是用来减少随机不确定性的东西。- -信息论创始人香农

信息论中最基本、最重要的概念(消息的同义词)

哲学家:“信息是物质成分和意识成分按完全特殊的方式融合起来的产物”。

数学家:“信息是使概率分布发生改变的东西”。

信息论:信息是事物运动状态或存在方式的不确定性的描述。

  • 运动状态:事物在空间上所展现的形状和态势
  • 存在方式:事物在时间上所呈现出的过程和规律

Ex1:简述消息、信号、信息的区别和联系

信号:是信息的物理表达层,是三个层次中最具体的层次。它是一个物理量,是一个载荷信息的实体,可测量、可描述、可显示(电信号、光信号、声音信号等),是消息的运载工具。

消息:是信息的载体,是用文字、符号、数据、语言、图片、图像等形式,把客观事物运动和主观思维活动的状态表达出来。

信息:它是更高层次哲学上的抽象,是信号与消息的更高表达层次,可以定量的描述。信息、物质和能量是构成一切系统的三大要素。

三个层次中,信号最具体,信息最抽象。它们三者之间的关系是哲学上的内涵与外延的关系。

信号 \Rightarrow 消息 \Rightarrow 信息:对通信系统来说,传输的是信号,信号承载着消息,消息中的不确定
成份是信息。

信息的分类

信息的分类:语法、语义、语用构成语言的三个基本方面,对于中间包含的信息按照性质分类

  • 语法信息:最基本也是最抽象的类型是语法信息,也是迄今为止在理论上研究得最多的类型。
    • 语法学研究符号与符号之间的关系。
  • 语义信息
    • 语义学研究符号与所指事物之间的关系。
  • 语用信息
    • 语用学研究符号与使用者之间的关系。

Ex2:文献调研

有一个情报部门,其主要任务是对经济情报进行收集、整理与分析以提供给决策机构。

该部门设三个组:

  1. 信息收集组将收集到的资料按中文、英文或其他文字、明文、密文进行分类;
  2. 信息处理组根据资料的性质进行翻译或破译得到资料的含义
  3. 信息分析组根据资料从中挑选出有价值的情报提交给决策机构。

Ex3:爱因斯坦方程 E=mc2E=mc^2

公式的语法信息:英文字母和数学符号之间的特定排列

公式的语义信息:EE 代表能量,mm 代表质量,cc 代表光速

公式的语用信息:公式的含义,试验的验证证实逻辑合理性——通过改变原子核的质量状态来获得巨大的原子核能。

信息的量化

信息的度量(信息量):和不确定性消除的程度有关,消除多少不确定性,就获得多少信息量。

不确定性:随机性,用概率论和随机过程来测度不确定性的大小。

  • 概率 \downarrow 不确定性 \uparrow 信息量 \uparrow
  • 概率 \uparrow 不确定性 \downarrow 信息量 \downarrow

概率信息/窄义信息:从事物的客观性出发讨论问题,而与事物的主观性无缘。

哈特莱早在 20 世纪 20 年代就提出用对数作为信息量的测度,信息量的计算式可表示为:

I(xk)=log21pk=log2pkI\left(x_k\right)=\log _2 \frac{1}{p_k}=-\log _2 p_k

  • pkp_k:事件 xkx_k 发生的概率

Ex4:有八只灯泡,只知其中有一只灯丝已断,用一节电池来测,问最少需测几次可找到坏灯泡(解除其不定度)?

解:

88 个灯泡断丝的概率相同,等概率事件,且:p(x)=1/8p(x)=1 / 8

log2p(x)=log28=3 bit -\log _2 p(x)=\log _2 8=3 \text { bit }

二元判断:每一次判断后可得到一个是/否信息。

二元信息。用来解除不定度的最小单位,也即信息量的最小单位——比特(bit)

log22=1 bit \log _2 2=1 \text { bit }

每次测量可以解除 1 bit1 \ \text{bit} 不确定度,至少需要测量 33 次。

信息科学与信息论

在人类文明史中,信息传输和传播手段经历了五次重大变革:

①语言的产生。

②文字的产生。

③印刷术的发明。

④电报、电话的发明。

⑤近现代,计算机技术与通信技术相结合,带来了信息技术革命性发展。

信息科学是以信息作为主要研究对象,以信息过程的运动规律作为主要研究内容,以信息科学方法论作为主要研究方法,以扩展人的信息功能(全部信息功能形成的有机整体就是智力功能)作为主要研究目标的一门科学。其主要内容包括:

  • 阐明信息的概念和本质(哲学信息论);
  • 探讨信息的度量和变换(基本信息论);
  • 研究信息的提取方法(识别信息论);
  • 澄清信息的传递规律(通信理论);
  • 探明信息的处理机制(智能理论);
  • 探究信息的再生理论(决策理论);
  • 阐明信息的调节原则(控制理论);
  • 完善信息的组织理论(系统理论)
  • \cdots \cdots

历史:

  • 1924 年,奈奎斯特(Harry Nyquist)(1889-1976)解释了信号带宽和信息速率之间的关系

  • 1928 年,哈特莱(Hartley)首先提出了用对数度量信息的概念。

  • 20 世纪30 年代,新的调制方式,如调频、调相、单边带调制、脉冲编码调制和增量调制的出现,使人们对信息能量、带宽和干扰的关系有了进一步的认识。

  • 1936 年,阿姆斯特朗(Edwin·Armstrong)指出增大带宽可以使抗干扰能力加强,并根据这一思想提出了宽频移的频率调制方法

  • 1939 年,达得利(Homer Dudley)发明了带通声码器,指出通信所需带宽至少同待传送消息的带宽应该一样声码器是最早的语音数据压缩系统。这一时期还诞生了无线电广播和电视广播

  • 1941~1944 年,香农(Shannon)在对通信和密码进行深入研究,用概率论和数理统计的方法系统地讨论了通信的基本问题,得出了几个重要而带有普遍意义的结论:

    • 阐明通信系统传递的对象;
    • 提出了信息熵的概念;
    • 指出通信系统的中心问题;
    • 指明了解决问题的方法。

    以上这些成果 1948 年以“通信的数学理论”(A mathematical theory of communication)为题公开发表,标志着信息论的正式诞生。

  • 维纳(Wiener)在研究火控系统和人体神经系统时,提出了在干扰作用下的信息最佳滤波理论,成为信息论的一个重要分支。

  • 50 年代,信息论在学术界引起了巨大反响。1951年,美国无线电工程师协会(IRE)成立了信息论组,并于 1955 年正式出版了信息论汇刊

  • 1959 年,香农发表了“保真度准则下的离散信源编码定理”(Coding theorems for a discrete source with a fidelity criterion)系统地提出了信息率失真理论(rate-distortion theory)。为信源压缩编码的研究奠定了理论基础。

  • 60 年代,信道编码技术有了较大发展,尤其,以 Viterbi 译码为代表的译码方法被美国卫星通信系统采用后,使它成为信息论的又一重要分支。

  • 1961 年,香农的重要论文“双路通信信道”开拓了网络信息论的研究。

  • 1970 年以来,随着卫星通信、计算机通信网的迅速发展,网络信息理论的研究成为当前信息论的中心研究课题之一。

信息论

信息论研究对象:信息论是一门应用概率论、随机过程、数理统计和近世代数的方法,来研究信息的传输、提取和处理系统中一般规律的工程学科。

 

信息论研究目的:提高信息系统的可靠性、有效性和安全性以便达到系统最优化。

信息论研究内容:主要是香农信息论或狭义信息论。

香农信息论或狭义信息论:应用近代概率统计方法研究信息的基本性质及度量方法,研究信息传输、处理等一般规律的学科。

  • 通信过程作为一个系统进行考察,提出了“通信系统的随机模型”。把有许多复杂的通信机构和过程简化为由信源、编码、信道、噪声、译码及信宿组成的一个信息的发送、传递、加工、接收系统。

  • 统计和概率观点引入通信理论,重新定义了信息和信息量,找到信息传输过程中的共同规律,帮助传递信息,而不是理解信息。

一般信息论:主要是研究信息传输和处理问题。除了香农基本理论之外,还包括噪声理论、信号滤波和预测、统计检测与估计理论、调制理论(维纳)。

广义信息论:概括说来,凡是能够用广义通信系统模型描述的过程或系统,都能用信息基本理论来研究。

信息论研究和学习意义:

  • 一种理论方法:帮助工程师从全局的观点观察和设计信息传输系统
  • 一种指导原则:提供一系列支持信息传输实践的原则和理论界,指导信息传输系统最优化设计

Ex5:13 个外观完全一样的小球,其中有 1 个小球重量与其余 12 个不同,有一架
天枰,问最少称几次,可以找到这个小球,并判断其是轻还是重。

解:

1313 个小球异常的概率相同,等概率事件,且:p(x)=1/13p(x)=1 / 13

判断轻重,需要增加一个二元判断

log2p(x)+1=log213+log22=log213+log22=log226 bit -\log _2 p(x)+1=\log _2 13+\log _2 2=\log _2 13+\log _2 2=\log _2 26 \text { bit }

天枰有三种状态(平衡,左倾,右倾),每称一次可以解除不确定度 log23  bit \log _2 3 \ \text { bit }

至少需要测量 log226 bit log23  bit =log326\frac{\log _2 26 \text { bit }}{\log _2 3 \ \text { bit }}=\log _3 26 次。

log39<log326<log327\log _3 9<\log _3 26<\log _3 27

只需称三次一定可以完全解除不确定度,少于三次不可能求解(可行性评估)

三门问题(蒙提霍尔(Monty Hall)问题):在你眼前有 33 扇巨大的关闭的门,编号分别是 AABBCC。站在旁边的主持人告诉你,其中一扇门的后面摆着极为诱人的大奖,而另外两扇门的后面各站着一头羊,你需要在这 33 扇门中选择 11 扇门,并获得那扇门后面的奖品。你经过深思熟虑,选择了编号为 AA 的门,在你紧张兮兮正准备打开时,主持人说慢着,然后他打开了编号为 CC 的门,后面正好是一头山羊,然后他问你:现在再给你一次选择的机会,你是坚持选择现在的门 AA,还是更换成门 BB?给出你的理由。

Shannon

Claude El-wood Shannon(1916-2001)美国科学家, 信息论创始人。

教育背景及工作经历:

  • 1916年4月30日出生于美国密歇根州。
  • 1936年毕业于密歇根大学并获得数学和电子工程学士学位。
  • 1940年分别获得麻省理工学院电子工程硕士学位和数学博士学位。
  • 1941年加入贝尔实验室数学部。
  • 1956年成为麻省理工学院(MIT)客座教授,并于1958年成为终生教授,1978年成为名誉教授。
  • 2001年2月26日去世,享年84岁。

学术成就:

  1. 创立信息论

1940年在普林斯顿高级研究所(The Institute for Advanced Study at Princeton)期间开始思考信息论与有效通信系统的问题。经过8年的努力,香农在1948年6月和10月在《贝尔系统技术杂志》(Bell System Technical Journal)上连载发表了具有深远影响的论文《通讯的数学原理》。1949年,香农又在该杂志上发表了另一著名论文《噪声下的通信》。在这两篇论文中,香农阐明了通信的基本问题,给出了通信系统的模型,提出了信息量的数学表达式,并解决了信道容量、信源统计特性、信源编码、信道编码等一系列基本技术问题。两篇论文成为了信息论的奠基性著作。

  1. 奠定了数字电路的理论基础

香农硕士论文题目是《A Symbolic Analysis of Relay and Switching Circuits》(继电器与开关电路的符号分析)。在该论文中,香农注意到电话交换电路与布尔代数之间的类似性,即把布尔代数的“真”与“假”和电路系统的“开”与“关”对应起来,并用1和0表示。于是他用布尔代数分析并优化开关电路,这就奠定了数字电路的理论基础。哈佛大学的Howard Gardner教授说,“这可能是本世纪最重要、最著名的一篇硕士论文。”1940年香农在MIT获得数学博士学位,而他的博士论文却是关于人类遗传学的,题目是《An Algebra for Theoretical Genetics》(理论遗传学的代数学)。这说明香农的科学兴趣十分广泛,后来他在不同的学科方面发表过许多有影响的文章。

信源模型

信源种类:

  • 汉语:字构成词,词构成句子
  • 字母:字母构成单词,单词构成句子
  • 图像:像素点构成
  • 声音:采样值构成

符号(基本消息):字,字母,像素点,采样值。如 aabbcc\cdots

符号集合(基本消息集合):如 {a,b,c,,z}\{a,b,c,\cdots,z\}

符号串(消息):如 Hello World

信源发出消息的过程 \Leftrightarrow 从一个基本消息集合取出基本消息的过程

对认识主体来说,信源在某一时刻输出的符号是随机的,可以用概率统计的数学方法描述

其余分类:

  • 单符号信源:信源每次输出一个符号,用离散随机变量描述
  • 多符号信源:信源每次输出多个符号(符号序列),用离散随机矢量描述
  • 离散信源:信源符号取值离散,包括单符号和多符号信源
  • 连续信源:信源符号取值连续,用随机过程描述

单符号离散信源的数学模型

[XP(X)]={x1,x2,,xi,,xnp(x1),p(x2),,p(xi),,p(xn)}p(xi)0i=1np(xi)=1\left[\begin{array}{c} \boldsymbol{X} \\ P(\boldsymbol{X}) \end{array}\right]=\left\{\begin{array}{cccccc} &x_1, \quad &x_2, &\cdots, \quad &x_i, &\cdots, \quad &x_n \\ &p\left(x_1\right), &p\left(x_2\right), &\cdots, &p\left(x_i\right), &\cdots, &p\left(x_n\right) \end{array}\right\} \quad p\left(x_i\right) \geq 0 \quad \sum_{i=1}^n p\left(x_i\right)=1

X\boldsymbol{X}:随机变量

p(xi)p(x_i):先验概率

单符号连续信源的数据模型

[XP(X)]=[x(a,b)p(x)]abp(x)dx=1p(x)0\left[\begin{array}{c} \boldsymbol{X} \\ P(\boldsymbol{X}) \end{array}\right]=\left[\begin{array}{c} x \in(a, b) \\ p(x) \end{array}\right] \quad \int_a^b p(x) d x=1 \quad p(x) \geq 0

X\boldsymbol{X}:随机变量

p(x)p(x):概率密度

多符号离散信源的数学模型

[XNP(αi)]=[α1,α2,αqNP(α1),P(α2),,P(αqN)]\left[\begin{array}{l} X^N \\ P\left(\alpha_i\right) \end{array}\right]=\left[\begin{array}{cccc} &\alpha_1, &\alpha_2, &\cdots &\alpha_{q^N} \\ & P\left(\alpha_1\right), &P\left(\alpha_2\right), &\cdots, &P\left(\alpha_{q^N}\right) \end{array}\right]

αi=(ai1ai2aiN)(i1,i2,iN=1,2,q)\alpha_i=\left(a_{i_1} a_{i_2} \cdots a_{i_N}\right) \quad\left(i_1, i_2, \cdots i_N=1,2, \cdots q\right)

P(αi)=P(ai1ai2aiN)=ik=1NP(aik)P\left(\alpha_i\right)=P\left(a_{i_1} a_{i_2} \cdots a_{i_N}\right)=\prod_{i_k=1}^N P\left(a_{i_k}\right)

i=1qNP(αi)=ik=1NP(aik)=1\sum_{i=1}^{q^N} P\left(\alpha_i\right)=\sum \prod_{i_k=1}^N P\left(a_{i_k}\right)=1

自信息

定义:若一个随机事件的概率为 p(xi)p(x_i),数学定义

I(xi)=f(p(xi))=logp(xi)I\left(x_i\right)=f\left(p\left(x_i\right)\right)=-\log p\left(x_i\right)

f()f(\cdot):信息测度。

性质

  • 非负

  • 单调递减

  • p(xi)=0p\left(x_i\right)=0 时,I(xi)I\left(x_i\right) \rightarrow \infty,不可能事件

  • p(xi)=1p\left(x_i\right)=1 时,I(xi)=0I\left(x_i\right)=0,确定性事件

含义

当事件 xix_i 发生以前,表示事件 xix_i 发生的不确定性。

当事件 xix_i 发生以后,表示事件 xix_i 提供的信息量。

单位:取决于对数的底

单位 表示
22 比特 bit\text {bit}
ee 奈特 nat\text {nat}
1010 哈特 hat\text {hat}

1 nat =1.44 bit ,1 hat =3.32bit1 \text { nat }=1.44 \text { bit }, 1 \text { hat }=3.32 \mathrm{bit}

联合自信息

二维联合集 XY\boldsymbol{XY} 上的元素(xiyjx_iy_j)的联合自信息:

I(xiyj)=logp(xiyj)I\left(x_i y_j\right)=-\log p\left(x_i y_j\right)

其中

[XYP(XY)]={x1y1,,x1ym,x2y1,,x2ym,,xny1,,xnymp(x1y1),,p(x1ym),p(x2y1),,p(x2ym),,p(xny1),p(xnym)}0p(xiyj)1i=1nj=1mp(xiyj)=1\left[\begin{array}{c} X Y \\ P(X Y) \end{array}\right]=\left\{\begin{array}{cccccccccc} &x_1 y_1, &\cdots, &x_1 y_m, &x_2 y_1, &\cdots, &x_2 y_m, &\cdots, &x_n y_1, &\cdots, &x_n y_m \\ &p\left(x_1 y_1\right), &\cdots, &p\left(x_1 y_m\right), &p\left(x_2 y_1\right), &\cdots, &p\left(x_2 y_m\right), &\cdots, &p\left(x_n y_1\right), &\cdots &p\left(x_n y_m\right) \end{array}\right\} \quad 0 \leq p\left(x_i y_j\right) \leq 1 \quad \sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right)=1

条件自信息

二维联合集 XY\boldsymbol{XY} 上的元素(xiyjx_iy_j)的条件自信息:

I(xi/yj)=log2p(xi/yj)I\left(x_i / y_j\right)=-\log _2 p\left(x_i / y_j\right)

关系计算

I(xiyj)=log2p(xi)p(yj/xi)=I(xi)+I(yj/xi)=log2p(yj)p(xi/yj)=I(yj)+I(xi/yj)\begin{aligned} I\left(x_i y_j\right) & =-\log _2 p\left(x_i\right) p\left(y_j / x_i\right)=I\left(x_i\right)+I\left(y_j / x_i\right) \\ & =-\log _2 p\left(y_j\right) p\left(x_i / y_j\right)=I\left(y_j\right)+I\left(x_i / y_j\right) \end{aligned}

XXYY 相互独立时

I(xiyj)=log2p(xi)log2p(yj)=I(xi)+I(yj)I\left(x_i y_j\right)=-\log _2 p\left(x_i\right)-\log _2 p\left(y_j\right)=I\left(x_i\right)+I\left(y_j\right)

Ex6:某地男青年中有 25%25\% 是大学生,他们其中 75%75\% 有驾照,而当地所有男青年中有驾照的比例为一半。问:当得知“某地有驾照的某位男青年是大学生”的消息时,我们获得多少信息量?

解:

设随机事件有

  • xix_i:”男青年是大学生“
  • yiy_i:”男青年有驾照“

xi/yix_i/y_i:”男青年有驾照是大学生“

yi/xiy_i/x_i:”男青年是大学生有驾照“

其中 p(xi)=14p(x_i)=\frac{1}{4}p(yi/xi)=34p(y_i/x_i)=\frac{3}{4}p(yi)=12p(y_i)=\frac{1}{2}

p(xi/yj)=p(xiyj)p(yj)=p(yj/xi)p(xi)p(yj)=38p\left(x_i / y_j\right)=\frac{p\left(x_i y_j\right)}{p\left(y_j\right)}=\frac{p\left(y_j / x_i\right) p\left(x_i\right)}{p\left(y_j\right)}=\frac{3}{8}

I(xi/yj)=log2p(xi/yj)=1.415bitI\left(x_i / y_j\right)=-\log _2 p\left(x_i / y_j\right)=1.415 b i t

互信息

[XP]={x1,x2,,xi,,xnp(x1),p(x2),,p(xi),,p(xn)}ip(xi)=1\left[\begin{array}{c} X \\ P \end{array}\right]=\left\{\begin{array}{cccccc} &x_1, & x_2, &\cdots, & x_i, &\cdots, & x_n \\ &p\left(x_1\right), &p\left(x_2\right), &\cdots, &p\left(x_i\right), &\cdots, &p\left(x_n\right) \end{array}\right\} \quad \sum_i p\left(x_i\right)=1

[YQ]={y1,y2,,yj,,ymq(y1),q(y2),,q(yj),,q(ym)}jq(yj)=1\left[\begin{array}{l} Y \\ Q \end{array}\right]=\left\{\begin{array}{cccccc} &y_1, &y_2, &\cdots, &y_j, &\cdots, &y_m \\ &q\left(y_1\right), &q\left(y_2\right), &\cdots, &q\left(y_j\right), &\cdots, &q\left(y_m\right) \end{array}\right\} \quad \sum_j q\left(y_j\right)=1

先验概率:信源发出消息 xix_i 的概率 p(xi)p(x_i)

后验概率:信宿收到消息 yiy_i 后推测信源发出消息 xix_i 的概率 p(xiyi)p(x_i|y_i)

互信息xix_i 的后验概率 p(xi)p(x_i) 与先验概率比值的对数称为 yiy_ixix_i 的互信息,即

I(xi;yj)=log2p(xi/yj)p(xi)(i=1,2,,n;j=1,2,,m)I\left(x_i ; y_j\right)=\log _2 \frac{p\left(x_i / y_j\right)}{p\left(x_i\right)} \quad(i=1,2, \cdots, n ; j=1,2, \cdots, m)

I(xi;yj)=log2p(xi)+log2p(xi/yj)=I(xi)I(xi/yj)I\left(x_i ; y_j\right)=-\log _2 p\left(x_i\right)+\log _2 p\left(x_i / y_j\right)=I\left(x_i\right)-I\left(x_i / y_j\right)

I(xi;yj)=log2p(xi/yj)log2p(xi)=log2p(xiyi)p(yi)log2p(xi)=log2p(xiyi)log2p(yi)log2p(xi)=I(xi)+I(yi)I(xiyj)I\left(x_i ; y_j\right)=\log _2 p\left(x_i / y_j\right)-\log _2 p\left(x_i\right)=\log _2 \frac{p\left(x_i y_i\right)}{p\left(y_i\right)}-\log _2 p\left(x_i\right)=\log _2 p\left(x_i y_i\right)-\log _2 p\left(y_i\right)-\log _2 p\left(x_i\right)=I\left(x_i\right)+I\left(y_i\right)-I\left(x_i y_j\right)

物理意义

I(xi;yj)=I(xi)I(xi/yj)I\left(x_i ; y_j\right)=I\left(x_i\right)-I\left(x_i / y_j\right)

  • I(xi;yj)I\left(x_i ; y_j\right):互信息,收到 yiy_i 而得到关于 xix_i 不确定度的减少量
  • I(xi)I\left(x_i\right):自信息,信宿收到 yiy_i 前,对信源发 xix_i 的不确定度
  • I(xi/yj)I\left(x_i / y_j\right):条件自信息,信宿收到 yiy_i 后,对信源发 xix_i 的不确定度

性质

  • 互易性:I(xi;yj)=I(yj;xi)I\left(x_i ; y_j\right)=I\left(y_j ; x_i\right)

证明:

I(xi;yj)=log2p(xi/yj)p(xi)=log2p(xiyi)p(yj)p(xi)=log2p(xiyj)p(xi)p(yi)I\left(x_i ; y_j\right)=\log _2 \frac{p\left(x_i / y_j\right)}{p\left(x_i\right)}= \log _2 \frac{\frac{p\left(x_i y_i\right)}{p\left(y_j\right)}}{p\left(x_i\right)}=\log _2 \frac{p\left(x_i y_j\right)}{p\left(x_i\right)p\left(y_i\right)}

I(xi;yj)=log2p(yj/xi)p(yj)=log2p(xiyi)p(xj)p(yi)=log2p(xiyj)p(xi)p(yi)I\left(x_i ; y_j\right)=\log _2 \frac{p\left(y_j / x_i\right)}{p\left(y_j\right)}=\log _2 \frac{\frac{p\left(x_i y_i\right)}{p\left(x_j\right)}}{p\left(y_i\right)}=\log _2 \frac{p\left(x_i y_j\right)}{p\left(x_i\right)p\left(y_i\right)}

可知二者相等。

或由 I(xi;yj)=I(xi)+I(yi)I(xiyj)I\left(x_i ; y_j\right)=I\left(x_i\right)+I\left(y_i\right)-I\left(x_i y_j\right)

  • 当事件 xix_iyjy_j 统计独立,互信息为零
  • 任何两事件之间的互信息不可能大于其中任一事件的自信息

I(xi;yj)I(xi);I(xi;yj)I(yj)I\left(x_i ; y_j\right) \leq I\left(x_i\right) ; \quad I\left(x_i ; y_j\right) \leq I\left(y_j\right)

条件互信息

三维联合集 XYZ\boldsymbol{XYZ} 上的元素(xiyjzkx_iy_jz_k)的条件互信息:给定条件 zkz_k 下,xix_iyjy_j 之间的互信息。

I(xi;yj/zk)=I(x/z)I(x/yz)=logp(xi/yjzk)p(xi/zk)I\left(x_i ; y_j / z_k\right)=I(x / z)-I(x / y z)=\log \frac{p\left(x_i / y_j z_k\right)}{p\left(x_i / z_k\right)}

推论1:

I(x;yz)=I(x;z)+I(x;y/z)I\left(x ; y z\right)=I\left(x ; z\right)+I\left(x ; y / z\right)

推论2:

I(x;y/z)I(x;y)=I(y;z/x)I(y;z)=I(z;x/y)I(z;x)\begin{aligned} I(x ; y / z)-I(x ; y) & =I(y ; z / x)-I(y ; z) =I(z ; x / y)-I(z ; x) \end{aligned}

Ex7:理想信道模型如下图

UU X=Y=y1y2y3X=Y=y_1 y_2 y_3 p(ui)p\left(u_i\right)
u1u_1 000000 14\frac{1}{4}
u2u_2 001001 14\frac{1}{4}
u3u_3 010010 18\frac{1}{8}
u4u_4 011011 18\frac{1}{8}
u5u_5 100100 116\frac{1}{16}
u6u_6 101101 116\frac{1}{16}
u7u_7 110110 116\frac{1}{16}
u8u_8 111111 116\frac{1}{16}

分别求出 p(ui/y1=0)p\left(u_i / y_1=0\right)p(ui/y1=0,y2=1)p\left(u_i / y_1=0, y_2=1\right)p(ui/y1y2y3=010)p\left(u_i / y_1 y_2 y_3=010\right)

I(u3;y1=0)I\left(u_3 ; y_1=0\right)I(u3;y2=1/y1=0)I\left(u_3 ; y_2=1 / y_1=0\right)I(u3;y3=0/y1=0;y2=1)I\left(u_3 ; y_3=0 / y_1=0 ; y_2=1\right)I(u3;y1y2y3=010)I\left(u_3 ; y_1 y_2 y_3=010\right)

解:

对于无记忆无损信道

p(ui/y1=0)p\left(u_i / y_1=0\right)

p(u1/y1=0)=p(y1=0/u1)p(u1)p(y1=0)=1×142×14+2×18=13=p(u2/y1=0)\begin{aligned} p\left(u_1 / y_1=0\right) & =\frac{p\left(y_1=0 / u_1\right) p\left(u_1\right)}{p\left(y_1=0\right)}\\ &=\frac{1 \times \frac{1}{4}}{2 \times \frac{1}{4}+2 \times \frac{1}{8}} \\ & =\frac{1}{3}=p\left(u_2 / y_1=0\right) \end{aligned}

p(u3/y1=0)=p(y1=0/u3)p(u3)p(y1=0)=1×182×14+2×18=16=p(u4/y1=0)\begin{aligned} p\left(u_3 / y_1=0\right) & =\frac{p\left(y_1=0 / u_3\right) p\left(u_3\right)}{p\left(y_1=0\right)}\\ &=\frac{1 \times \frac{1}{8}}{2 \times \frac{1}{4}+2 \times \frac{1}{8}} \\ & =\frac{1}{6}=p\left(u_4 / y_1=0\right) \end{aligned}

p(ui/y1=0)=p(y1=0/ui)p(ui)p(y1=0)=0×1162×14+2×18=0(i=5,6,7,8)\begin{aligned} p\left(u_i / y_1=0\right) & =\frac{p\left(y_1=0 / u_i\right) p\left(u_i\right)}{p\left(y_1=0\right)}\\ &=\frac{0 \times \frac{1}{16}}{2 \times \frac{1}{4}+2 \times \frac{1}{8}} \\ & =0 \quad (i=5,6,7,8) \end{aligned}

p(ui/y1=0,y2=1)=p(y1=0,y2=1/ui)p(ui)p(y1=0,y2=1)={1×182×18=12(i=3,4)0×p(ui)2×18=0(i=1,2,5,6,7,8)\begin{aligned} p\left(u_i / y_1=0, y_2=1\right) & =\frac{p\left(y_1=0, y_2=1 / u_i\right) p\left(u_i\right)}{p\left(y_1=0, y_2=1\right)}=\left\{\begin{array}{c} \frac{1 \times \frac{1}{8}}{2 \times \frac{1}{8}}&=\frac{1}{2} \quad (i=3,4) \\ \frac{0 \times p\left(u_i\right)}{2 \times \frac{1}{8}}&=0 \quad (i=1,2,5,6,7,8) \end{array}\right. \end{aligned}

p(ui/y1y2y3=010)=p(y1y2y3=010/ui)p(ui)p(y1y2y3=010)={1×1818=1(i=3)0×p(ui)18=0(i=1,2,4,5,6,7,8)\begin{aligned} p\left(u_i / y_1 y_2 y_3=010\right) & =\frac{p\left(y_1 y_2 y_3=010 / u_i\right) p\left(u_i\right)}{p\left(y_1 y_2 y_3=010\right)}=\left\{\begin{array}{c} \frac{1 \times \frac{1}{8}}{\frac{1}{8}}&=1 \quad (i=3) \\ \frac{0 \times p\left(u_i\right)}{\frac{1}{8}}&=0 \quad (i=1,2,4,5,6,7,8) \end{array}\right. \end{aligned}

I(u3;y1=0)=logp(u3/y1=0)p(u3)=log243=0.415 bitI\left(u_3 ; y_1=0\right)=\log \frac{p\left(u_3 / y_1=0\right)}{p\left(u_3\right)}=\log _2 \frac{4}{3}=0.415 \ \text{bit}

I(u3;y2=1/y1=0)=logp(u3/y1=0,y2=1)p(u3/y1=0)=log1216=log23=1.585 bit \begin{gathered}I\left(u_3 ; y_2=1 / y_1=0\right)=\log \frac{p\left(u_3 / y_1=0, y_2=1\right)}{p\left(u_3 / y_1=0\right)} =\log \frac{\frac{1}{2}}{\frac{1}{6}}=\log _2 3=1.585 \text { bit }\end{gathered}

I(u3;y2=1/y1=0,y2=1)=logp(u3/y1y2y3=010)p(u3/y1=0,y2=1)=log112=log22=1 bit \begin{gathered}I\left(u_3 ; y_2=1 / y_1=0, y_2=1\right)=\log \frac{p\left(u_3 / y_1 y_2 y_3=010\right)}{p\left(u_3 / y_1=0, y_2=1\right)}=\log \frac{1}{\frac{1}{2}}=\log _2 2=1 \text { bit }\end{gathered}

I(u3;y1y2y3=010)=logp(u3/y1y2y3=010)p(u3)=log118=log28=3 bit \begin{aligned} I\left(u_3 ; y_1 y_2 y_3=010\right)=\log \frac{p\left(u_3 / y_1 y_2 y_3=010\right)}{p\left(u_3\right)}=\log \frac{1}{\frac{1}{8}}=\log _2 8=3 \text { bit }\end{aligned}

I(u3;y1y2y3=010)=I(u3;y1=0)+I(u3;y2=1/y1=0)+I(u3;y3=0/y1=0,y2=1)=0.415+1.585+1=3 bit \begin{aligned} I\left(u_3 ; y_1 y_2 y_3=010\right)=I\left(u_3 ; y_1=0\right) +I\left(u_3 ; y_2=1 / y_1=0\right)+I\left(u_3 ; y_3=0 / y_1=0, y_2=1\right)=0.415+1.585+1=3 \text { bit }\end{aligned}

Ex8:实际信道模型如上图,一个等概率信源有八种消息符号,用四比特码字序列编码,码字中每一个二进制符号经信道输出可得二元符号 yy,已知条件概率(信道特性)为 P00=P11=1εP_{00}=P_{11}=1-\varepsilonP01=P10=εP_{01}=P_{10}=\varepsilon。这里 P00P_{00} 定义为信道转移概率 P(y=0/x=0)P(y=0 / x=0),以此类推。当实验结果得 y=0000\vec{y}=0000 时,求:

(1)第一位码测定后所得的关于 x\vec{x} 的自信息

(2)第二,第三,第四位码测定后各得多少关于 x1\vec{x}_1 的自信息

(3)全部结果 y=0000\vec{y}=0000 关于 x1\vec{x}_1 的自信息

(4)讨论 ε=0\varepsilon=0ε=12\varepsilon=\frac{1}{2} 时上述各自信息的情况

解:

x1=0000,x2=0011,x3=0101,x4=0110x5=1001,x6=1010,x7=1100,x8=1111\begin{aligned} & \overrightarrow{x_1}=0000, \overrightarrow{x_2}=0011, \overrightarrow{x_3}=0101, \overrightarrow{x_4}=0110 \\ & \overrightarrow{x_5}=1001, \overrightarrow{x_6}=1010, \overrightarrow{x_7}=1100, \overrightarrow{x_8}=1111\end{aligned}

P(xi)=1nP\left(\overrightarrow{x_i}\right)=\frac{1}{n}

P(x=0)=P(x=1)=12P(x=0)=P(x=1)=\frac{1}{2}

I(xi)=log21n=log2n=log28=3bitI\left(\overrightarrow{x_i}\right)=-\log _2 \frac{1}{n}=\log _2 n=\log _2 8=3 \text{bit}

(1)第一位码测定后所得的关于 x\vec{x} 的自信息

解法一:I(x1;y1=0)=I(x1)I(x1/y1=0)I\left(\overrightarrow{x_1} ; y_1=0\right)=I\left(\overrightarrow{x_1}\right)-I\left(\overrightarrow{x_1} / y_1=0\right)

  • P(y1=0/x1)=P(y1=0/x1=0)=1εP\left(y_1=0 / \overrightarrow{x_1}\right)=P\left(y_1=0 / x_1=0\right)=1-\varepsilon

  • P(y1=0)=P(x=0)(1ε)+P(x=1)ε=12P\left(y_1=0\right)=P(x=0)(1-\varepsilon)+P(x=1) \varepsilon=\frac{1}{2}

由上两式,得到

P(x1/y1=0)=P(y1=0/x1)p(x1)p(y1=0)=(1ε)1812=1ε4P\left(\overrightarrow{x_1} / y_1=0\right)=\frac{P\left(y_1=0 / \overrightarrow{x_1}\right) p\left(\overrightarrow{x_1}\right)}{p\left(y_1=0\right)}=\frac{(1-\varepsilon) \cdot \frac{1}{8}}{\frac{1}{2}}=\frac{1-\varepsilon}{4}

I(x1/y1=0)=log2P(x1/y1=0)=log21ε4=log241εbitI\left(\overrightarrow{x_1} / y_1=0\right)=-\log _2 P\left(\overrightarrow{x_1} / y_1=0\right)=-\log _2 \frac{1-\varepsilon}{4}=\log _2 \frac{4}{1-\varepsilon} \text{bit}

即:

I(x1;y1=0)=I(x1)I(x1/y1=0)=3log241ε=log[2(1ε)]bitI\left(\overrightarrow{x_1} ; y_1=0\right)=I\left(\overrightarrow{x_1}\right)-I\left(\overrightarrow{x_1} / y_1=0\right)=3-\log _2 \frac{4}{1-\varepsilon}=\log [2(1-\varepsilon)] \text{bit}

解法二:I(x1;y1=0)=I(y1=0)I(y1=0/x1)I\left(\overrightarrow{x_1} ; y_1=0\right)=I\left(y_1=0\right)-I\left(y_1=0 / \overrightarrow{x_1}\right)

P(y1=0)=P(x=0)(1ε)+P(x=1)ε=12P\left(y_1=0\right)=P(x=0)(1-\varepsilon)+P(x=1) \varepsilon=\frac{1}{2}

I(y1=0)=log2P(y1=0)=log212=1bitI\left(y_1=0\right)=-\log _2 P\left(y_1=0\right)=-\log _2 \frac{1}{2}=1 \text{bit}

P(y1=0/x1)=P(y1=0/x1=0)=1εP\left(y_1=0 / \overrightarrow{x_1}\right)=P\left(y_1=0 / x_1=0\right)=1-\varepsilon

I(y1=0/x1)=log2P(y1=0/x1)=log21ε=log211εbitI\left(y_1=0 / \overrightarrow{x_1}\right)=-\log _2 P\left(y_1=0 / \overrightarrow{x_1}\right)=-\log _2 1-\varepsilon=\log _2 \frac{1}{1-\varepsilon} \text{bit}

即:

I(x1;y1=0)=I(y1=0)I(y1=0/x1)=1log211ε=log[2(1ε)]bitI\left(\overrightarrow{x_1} ; y_1=0\right)=I\left(y_1=0\right)-I\left(y_1=0 / \overrightarrow{x_1}\right)=1-\log _2 \frac{1}{1-\varepsilon}=\log [2(1-\varepsilon)] \text{bit}

(2)第二,第三,第四位码测定后各得多少关于 x1\vec{x}_1 的自信息

依次对应要求:

  • I(x1;y2=0/y1=0)I\left(\overrightarrow{x_1} ; y_2=0 / y_1=0\right)
  • I(x1;y3=0/y1=y2=0)I\left(\overrightarrow{x_1} ; y_3=0 / y_1=y_2=0\right)
  • I(x1;y4=0/y1=y2=y3=0)I\left(\overrightarrow{x_1} ; y_4=0 / y_1=y_2=y_3=0\right)

求条件互信息的四种方法:

I(x;y/z)=I(x/z)I(x/yz)=I(y/z)I(y/xz)=I(x;yz)I(x;z)=I(y;xz)I(y;z)\begin{aligned} I(x ; y / z) & =I(x / z)-I(x / y z) \\ & =I(y / z)-I(y / x z) \\ & =I(x ; y z)-I(x ; z) \\ & =I(y ; x z)-I(y ; z) \end{aligned}

(3)全部结果 y=0000\vec{y}=0000 关于 x1\vec{x}_1 的自信息

(4)讨论 ε=0\varepsilon=0ε=12\varepsilon=\frac{1}{2} 时上述各自信息的情况

信源熵—平均自信息

自信息是一个随机事件概率的函数, 不能用作整个信源的信息测度。为了研究整个事件集合或符号序列(如信源)的平均的信息量(总体特征),就需要引入新的概念。定义自信息的数学期望为信源的平均信息量,称为信源的平均自信息量(信源熵)

H(X)=E[I(ai)]=i=1npiI(ai)=i=1npilogp(ai)H(X)=E\left[I\left(a_i\right)\right]=\sum_{i=1}^n p_i I\left(a_i\right)=-\sum_{i=1}^n p_i \log p\left(a_i\right)

单位:  bit / 信源 \text { bit } / \text { 信源 }

物理意义

  • 信源输出前,信源的平均不确定性
  • 信源输出后,每个符号所携带的平均信息量
  • 反映了变量 XX 的随机性。

Ex9:天气预报,有两个信源,求 H(X)H(X)H(Y)H(Y)

[XP]={a1,a20.8,0.2}[YQ]={b1b20.5,0.5}\left[\begin{array}{c} X \\ P \end{array}\right]=\left\{\begin{array}{cc} a_1, & a_2 \\ 0.8, & 0.2 \end{array}\right\} \quad\left[\begin{array}{l} Y \\ Q \end{array}\right]=\left\{\begin{array}{cc} b_1 & b_2 \\ 0.5, & 0.5 \end{array}\right\}

解:

H(X)=0.8log20.80.2log20.2=0.702 bit / 符号 H(X)=-0.8 \log _2 0.8-0.2 \log _2 0.2=0.702 \text { bit } / \text { 符号 }

H(Y)=0.5log20.50.5log20.5=1bit/ 符号 H(Y)=-0.5 \log _2 0.5-0.5 \log _2 0.5=1 \mathrm{bit} / \text { 符号 }

显然 H(X)<H(Y)H(X)<H(Y),说明信源 YY 比信源 XX 的平均不确定性要大,信源 YY 提供信息发送能力比信源 XX 大,对于离散信源来说,等概率分布是发送信息能力最大的必要条件。

Ex10:设某信源输出四个符号,其符号集合的概率分布为,求则其熵。

S={s1s2s3s4p1p2p3p4}={s1s2s3s412141818}S=\left\{\begin{array}{cccc} s 1 & s 2 & s 3 & s 4 \\ p 1 & p 2 & p 3 & p 4 \end{array}\right\}=\left\{\begin{array}{cccc} s 1 & s 2 & s 3 & s 4 \\ \frac{1}{2} & \frac{1}{4} & \frac{1}{8} & \frac{1}{8} \end{array}\right\}

解:

H(S)=i=14pilogpi=12log2+14log4+28log8=1.75 bit / 符号 H(S)=-\sum_{i=1}^4 p_i \log p_i=\frac{1}{2} \log 2+\frac{1}{4} \log 4+\frac{2}{8} \log 8=1.75 \text { bit } / \text { 符号 }

Ex11:电视屏上约有 500×600=3×105500 \times 600=3 \times 10^5 个格点,按每点有 1010 个不同的灰
度等级考虑,则共能组成 n=103×105n=10^{3 \times 10^5} 个不同的画面。按等概率 1103×105\frac{1}{10^{3 \times 10^5}} 计算,求平均每个画面可提供的信息量。

解:

H(X)=i=1np(xi)log2p(xi)=log2103×105=3×105×3.32bit/ 画面 \begin{aligned} H(X) & =-\sum_{i=1}^n p\left(x_i\right) \log _2 p\left(x_i\right)=-\log _2 10^{-3 \times 10^5} =3 \times 10^5 \times 3.32 \mathrm{bit} / \text { 画面 } \end{aligned}

联合熵

联合熵:联合离散符号集合 XYXY 上的每个元素对 xiyjx_iy_j 的联合自信息量的数学期望。

H(XY)=i=1nj=1mp(xiyj)I(xiyj)=i=1nj=1mp(xiyj)log2p(xiyj)H(X Y)=\sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right) I\left(x_i y_j\right)=-\sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right) \log _2 p\left(x_i y_j\right)

条件熵

条件熵:联合离散符号集合 XYXY 上的条件自信息量的数学期望。

H(Y/X)=E[I(yj/xi)]=j=1mi=1np(xiyj)I(yj/xi)=j=1mi=1np(xiyj)log2p(yj/xi)\begin{aligned} H(Y / X)=E\left[I\left(y_j / x_i\right)\right] & =\sum_{j=1}^m \sum_{i=1}^n p\left(x_i y_j\right) I\left(y_j / x_i\right)=-\sum_{j=1}^m \sum_{i=1}^n p\left(x_i y_j\right) \log _2 p\left(y_j / x_i\right) \end{aligned}

思考:求条件熵时为什么要用联合概率加权?

p(yj/xi)=p(yjxi)p(xi)j=1np(yj/xi)=1p\left(y_j / x_i\right)=\frac{p\left(y_j x_i\right)}{p\left(x_i\right)} \quad \sum_{j=1}^n p\left(y_j / x_i\right)=1

已知特定事件 xix_i 发生,事件 yjy_j 发生的不确定度为:

I(yj/xi)=logp(yj/xi)I\left(y_j / x_i\right)=-\log p\left(y_j / x_i\right)

统计集合 YY 中的所有元素

H(Y/xi)=j=1mp(yj/xi)logp(yj/xi)H\left(Y / x_i\right)=-\sum_{j=1}^m p\left(y_j / x_i\right) \log p\left(y_j / x_i\right)

再统计集合 XX 中的所有元素

H(Y/X)=i=1np(xi)H(Y/xi)=i=1np(xi)j=1mp(yj/xi)logp(yj/xi)=i=1nj=1mp(xi)p(yj/xi)logp(yj/xi)=i=1nj=1mp(xiyj)logp(yj/xi)\begin{aligned} & H(Y / X)\\ &=\sum_{i=1}^n p\left(x_i\right) H\left(Y / x_i\right)\\ &=-\sum_{i=1}^n p\left(x_i\right) \sum_{j=1}^m p\left(y_j / x_i\right) \log p\left(y_j / x_i\right) \\ & =-\sum_{i=1}^n \sum_{j=1}^m p\left(x_i\right) p\left(y_j / x_i\right) \log p\left(y_j / x_i\right)\\ &=-\sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right) \log p\left(y_j / x_i\right) \end{aligned}

熵函数

H(X)=H(p1,p2pn)=i=1npilogpii=1npi=1pi0(i=1,2,,n)H(X)=H\left(p_1, p_2 \cdots p_n\right)=-\sum_{i=1}^n p_i \log p_i \quad \sum_{i=1}^n p_i=1 \quad p_i \geq 0 \quad(i=1,2, \ldots, n)

二元熵H(X)=H(p,1p)=H(p)H(X)=H(p, 1-p)=H(p)

数学特性

对称性

H(p1,p2,pn)=H(pn,p1,p2,pn1)H\left(p_1, p_2, \ldots p_n\right)=H\left(p_n, p_1, p_2, \ldots p_{n-1}\right)

非负性

H(X)0H(X) \geq 0

拓展性

Limε0Hn+1(p1,p2,,pnε,ε)=Hn(p1,p2,,pn)\operatorname{Lim}_{\varepsilon \rightarrow 0} H_{n+1}\left(p_1, p_2, \ldots, p_n-\varepsilon, \varepsilon\right)=H_n\left(p_1, p_2, \ldots, p_n\right)

确定性

H(1,0)=H(0,1)=H(1,0,0,0)=0H(1,0)=H(0,1)=H(1,0,0, \ldots 0)=0

可加性

H(XY)=H(X)+H(Y/X)H(XY)=H(Y)+H(X/Y)\begin{aligned} & H(X Y)=H(X)+H(Y / X) \\ \\ & H(X Y)=H(Y)+H(X / Y) \end{aligned}

证明:

定义

[XP]={a1,a2,,anp1,p2,,pn}[YQ]={b1,b2,,bmq1,q2,,qm}[XYR]={c11,c12,,cnmr11,r12,,rnm}\left[\begin{array}{c} X \\ P \end{array}\right]=\left\{\begin{array}{c} a_1, a_2, \ldots, a_n \\ p_1, p_2, \ldots, p_n \end{array}\right\} \quad\left[\begin{array}{c} Y \\ Q \end{array}\right]=\left\{\begin{array}{c} b_1, b_2, \ldots, b_m \\ q_1, q_2, \ldots, q_m \end{array}\right\} \quad\left[\begin{array}{c} X Y \\ R \end{array}\right]=\left\{\begin{array}{c} c_{11}, c_{12}, \ldots, c_{n m} \\ r_{11}, r_{12}, \ldots, r_{n m} \end{array}\right\}

这里有

  • cij=aibjc_{i j}=a_i b_j
  • pji=p(y=bj/x=ai)p_{j i}=p\left(y=b_j / x=a_i\right)
  • qij=q(x=ai/y=bj)q_{i j}=q\left(x=a_i / y=b_j\right)
  • rij=pipji=qiqijr_{i j}=p_i p_{j i}=q_i q_{i j}

H(X)=H(p1,p2,,pn)=i=1npilogpiH(X)=H\left(p_1, p_2, \ldots, p_n\right)=-\sum_{i=1}^n p_i \log p_i

H(Y)=H(q1,q2,,qm)=j=1mqjlogqjH(Y)=H\left(q_1, q_2, \ldots, q_m\right)=-\sum_{j=1}^m q_j \log q_j

H(XY)=H(r11,r12,,rnm)=i=1nj=1mrijlogrijH(X Y)=H\left(r_{11}, r_{12}, \ldots, r_{n m}\right)=-\sum_{i=1}^n \sum_{j=1}^m r_{i j} \log r_{i j}

H(Y/X)=E[I(bj/ai)]=i=1nj=1mrijlogpji=i=1nj=1mpipjilogpjiH(Y / X)=E\left[I\left(b_j / a_i\right)\right]=-\sum_{i=1}^n \sum_{j=1}^m r_{i j} \log p_{j i}=-\sum_{i=1}^n \sum_{j=1}^m p_i p_{j i} \log p_{j i}

H(X/Y)=E[I(ai/bj)]=i=1nj=1mrijlogqij=i=1nj=1mpipjilogqijH(X / Y)=E\left[I\left(a_i / b_j\right)\right]=-\sum_{i=1}^n \sum_{j=1}^m r_{i j} \log q_{i j}=-\sum_{i=1}^n \sum_{j=1}^m p_i p_{j i} \log q_{i j}

H(XY)=i=1nj=1mrijlogrij=i=1nj=1mpipjilogpipji=i=1npij=1mpjilogpii=1nj=1mpipjilogpji (注: j=1mpji=1)=i=1npilogpii=1npij=1mpjilogpji=H(X)+i=1npi[j=1mpjilogpji]=H(X)i=1nj=1mpipjilogpji=H(X)+H(Y/X)\begin{aligned} H(X Y)&=-\sum_{i=1}^n \sum_{j=1}^m r_{i j} \log r_{i j}\\ &=-\sum_{i=1}^n \sum_{j=1}^m p_i p_{j i} \log p_i p_{j i} \\ & \left.=-\sum_{i=1}^n p_i \sum_{j=1}^m p_{j i} \log p_i-\sum_{i=1}^n \sum_{j=1}^m p_i p_{j i} \log p_{j i} \quad \text { (注: } \sum_{j=1}^m p_{j i}=1\right) \\ & =-\sum_{i=1}^n p_i \log p_i-\sum_{i=1}^n p_i \sum_{j=1}^m p_{j i} \log p_{j i} \\ & =H(X)+\sum_{i=1}^n p_i\left[-\sum_{j=1}^m p_{j i} \log p_{j i}\right]\\ &=H(X)-\sum_{i=1}^n \sum_{j=1}^m p_i p_{j i} \log p_{j i}\\ &=H(X)+H(Y / X) \end{aligned}

H(XY)=H(X)+H(Y/X)H(X Y)=H(X)+H(Y / X)

同理可证 H(XY)=H(Y)+H(X/Y)H(X Y)=H(Y)+H(X / Y)

推广:

H(X1,X2,X3,,XN)=H(X1)+H(X2/X1)+H(X3/X1X2)++H(XN/X1X2XN1)=i=1NH(X1/X1X2Xi1)\begin{aligned} H\left(X_1, X_2, X_3, \ldots, X_N\right) & =H\left(X_1\right)+H\left(X_2 / X_1\right)+H\left(X_3 / X_1 X_2\right)+\cdots+H\left(X_N / X_1 X_2 \cdots X_N 1\right)\\ &=\sum_{i=1}^N H\left(X_1 / X_1 X_2 \cdots X_{i-1}\right) \end{aligned}

H(X1X2XN)H(X1)+H(X2)++H(XN)H\left(X_1 X_2 \cdots X_N\right) \leq H\left(X_1\right)+H\left(X_2\right)+\cdots+H\left(X_N\right)

极值性

Hn(p1,p2,,pn)lognH_n\left(p_1, p_2, \ldots, p_n\right) \leq \log n

证明:

H(X)=H(p1,p2pn)=i=1npilogpii=1npi=1H(X)=H\left(p_1, p_2 \cdots p_n\right)=-\sum_{i=1}^n p_i \log p_i\quad \sum_{i=1}^n p_i=1

利用拉格朗日乘子构造函数:

G(p1,p2pn,λ)=i=1npilnpi+λ(i=1npi1)G\left(p_1, p_2 \cdots p_n, \lambda\right)=-\sum_{i=1}^n p_i \ln p_i+\lambda\left(\sum_{i=1}^n p_i-1\right)

分别对 pip_iλ\lambda,并令导数为 00,得

{Gpi=lnpi1+λ=0i=1npi1=0\left\{\begin{array}{l} \frac{\partial G}{\partial p_i}=-\ln p_i-1+\lambda=0 \\ \\ \sum_{i=1}^n p_i-1=0 \end{array}\right.

p1=p1==pn=1/np_1=p_1=\cdots=p_n=1 / n

推论1:任何概率分布下的信息熵一定不会大于它对其他概率分布下自信息的数学期望。

Hn(p1,p2,,pn)i=1npilogqiH_n\left(p_1, p_2, \ldots, p_n\right) \leq-\sum_{i=1}^n p_i \log q_i

证明:

利用 lnxx1\ln x \leq x-1,易得

lnqipiqipi1\ln \frac{q_i}{p_i} \leq \frac{q_i}{p_i}-1

i=1npilnqipi=i=1npilnpi+i=1npilnqii=1npi(qipi1)=i=1nqii=1npi=0\sum_{i=1}^n p_i \ln \frac{q_i}{p_i}=-\sum_{i=1}^n p_i \ln p_i+\sum_{i=1}^n p_i \ln q_i \leq \sum_{i=1}^n p_i\left(\frac{q_i}{p_i}-1\right)=\sum_{i=1}^n q_i-\sum_{i=1}^n p_i=0

i=1npilnpi+i=1npilnqi0-\sum_{i=1}^n p_i \ln p_i+\sum_{i=1}^n p_i \ln q_i \leq 0

i=1npilnpii=1npilnqi-\sum_{i=1}^n p_i \ln p_i \leq -\sum_{i=1}^n p_i \ln q_i

换底 logax=lnxlna\log _a x=\frac{\ln x}{\ln a},得 lnx=lnalogax\ln x=\ln a \log _a x

lnai=1npilogpilnai=1npilogqi-\ln a \sum_{i=1}^n p_i \log p_i \leq -\ln a \sum_{i=1}^n p_i \log q_i

i=1npilogpii=1npilogqi-\sum_{i=1}^n p_i \log p_i \leq -\sum_{i=1}^n p_i \log q_i

Hn(p1,p2,,pn)i=1npilogqiH_n\left(p_1, p_2, \ldots, p_n\right) \leq-\sum_{i=1}^n p_i \log q_i

推论2:条件熵一定不会大于无条件熵。

H(Y)H(Y/X)H(Y) \geq H(Y / X)

利用推论1,pip_i 替换为 pjip_{ji}qiq_i 替换为 qjq_j

j=1mpjilogpjij=1mpjilogqj-\sum_{j=1}^m p_{j i} \log p_{j i} \leq-\sum_{j=1}^m p_{j i} \log q_j

i=1npi=1\sum_{i=1}^n p_i=1

i=1npij=1mpjilogpjii=1npij=1mpjilogqj=j=1m[i=1npipji]logqj=j=1mqjlogqj\begin{aligned} -\sum_{i=1}^n p_i \sum_{j=1}^m p_{j i} \log p_{j i} \leq-\sum_{i=1}^n p_i \sum_{j=1}^m p_{j i} \log q_j =-\sum_{j=1}^m\left[\sum_{i=1}^n p_i p_{j i}\right] \log q_j=-\sum_{j=1}^m q_j \log q_j \end{aligned}

H(Y)H(Y/X)H(Y) \geq H(Y / X)

上凸性:设 f(X)=f(x1,x2,,xn)f(X)=f\left(x_1, x_2, \ldots, x_n\right) ,对于任意一个小于 11 的正数 α(0<α<1)\alpha(0<\alpha<1) 以及 f(X)f(X) 定义域内的两个任意矢量 X1X_1X2X_2,满足

f[αX1+(1α)X2]αf(X1)+(1α)f(X2)f\left[\alpha X_1+(1-\alpha) X_2\right] \geq \alpha f\left(X_1\right)+(1-\alpha) f\left(X_2\right)

则称 f(X)f(X) 为定义域上的上凸函数。

f[αX1+(1α)X2]>αf(X1)+(1α)f(X2)f\left[\alpha X_1+(1-\alpha) X_2\right] > \alpha f\left(X_1\right)+(1-\alpha) f\left(X_2\right):严格上凸函数。

证明1:

H[αpi+(1α)qi]=i=1n[αpi+(1α)qi]log[αpi+(1α)qi]=i=1nαpilog[αpi+(1α)qi]i=1n(1α)qilog[αpi+(1α)qi]=αi=1npilog{pipi[αpi+(1α)qi]}(1α)i=1nqilog{qiqi[αpi+(1α)qi]}=αi=1npilogpi(1α)i=1nqilogqiαi=1npilogαpi+(1α)qipi(1α)i=1nqilogαpi+(1α)qiqi=αH(P)+(1α)H(Q)αi=1npilogαpi+(1α)qipi(1α)i=1nqilogαpi+(1α)qiqi\begin{aligned} & H\left[\alpha p_i+(1-\alpha) q_i\right]\\ & =-\sum_{i=1}^n\left[\alpha p_i+(1-\alpha) q_i\right] \log \left[\alpha p_i+(1-\alpha) q_i\right] \\ & =-\sum_{i=1}^n \alpha p_i \log \left[\alpha p_i+(1-\alpha) q_i\right]-\sum_{i=1}^n(1-\alpha) q_i \log \left[\alpha p_i+(1-\alpha) q_i\right] \\ & =-\alpha \sum_{i=1}^n p_i \log \left\{\frac{p_i}{p_i}\left[\alpha p_i+(1-\alpha) q_i\right]\right\}-(1-\alpha) \sum_{i=1}^n q_i \log \left\{\frac{q_i}{q_i}\left[\alpha p_i+(1-\alpha) q_i\right]\right\} \\ & =-\alpha \sum_{i=1}^n p_i \log p_i-(1-\alpha) \sum_{i=1}^n q_i \log q_i-\alpha \sum_{i=1}^n p_i \log \frac{\alpha p_i+(1-\alpha) q_i}{p_i} -(1-\alpha) \sum_{i=1}^n q_i \log \frac{\alpha p_i+(1-\alpha) q_i}{q_i}\\ & =\alpha H(P)+(1-\alpha) H(Q)-\alpha \sum_{i=1}^n p_i \log \frac{\alpha p_i+(1-\alpha) q_i}{p_i}-(1-\alpha) \sum_{i=1}^n q_i \log \frac{\alpha p_i+(1-\alpha) q_i}{q_i} \end{aligned}

利用推论1

H(p1,p2,,pn)i=1npilog[αpi+(1α)qi]H\left(p_1, p_2, \ldots, p_n\right) \leq-\sum_{i=1}^n p_i \log \left[\alpha p_i+(1-\alpha) q_i\right]

i=1npilog[αpi+(1α)qi]+i=1npilogpi0-\sum_{i=1}^n p_i \log \left[\alpha p_i+(1-\alpha) q_i\right]+\sum_{i=1}^n p_i \log p_i \geq 0

依次证明:

αi=1npilogαpi+(1α)qipi=αi=1npilog[αpi+(1α)qi]+αi=1npilogpi0-\alpha \sum_{i=1}^n p_i \log \frac{\alpha p_i+(1-\alpha) q_i}{p_i}=-\alpha \sum_{i=1}^n p_i \log \left[\alpha p_i+(1-\alpha) q_i\right]+\alpha \sum_{i=1}^n p_i \log p_i \geq 0

(1α)i=1nqilogαpi+(1α)qiqi=(1α)i=1nqilog[αpi+(1α)qi]+(1α)i=1nqilogqi-(1-\alpha) \sum_{i=1}^n q_i \log \frac{\alpha p_i+(1-\alpha) q_i}{q_i}=-(1-\alpha) \sum_{i=1}^n q_i \log \left[\alpha p_i+(1-\alpha) q_i \right] + (1-\alpha) \sum_{i=1}^n q_i \log q_i

所以

H[αP+(1α)Q]αH(P)+(1α)H(Q)H[\alpha P+(1-\alpha) Q] \geq \alpha H(P)+(1-\alpha) H(Q)

证明2:

Jensen\text {Jensen} 不等式:若 ff[a,b][a,b] 上凸函数,则对于任意 x[a,b]x \in [a,b]λi>0(i=1,2,,n)\lambda_i>0 \quad (i=1,2,\cdots,n)i=1nλi=1\sum_{i=1}^n \lambda_i=1,满足

f(i=1nλixi)i=1nλif(xi)f\left(\sum_{i=1}^n \lambda_i x_i\right) \geq \sum_{i=1}^n \lambda_i f\left(x_i\right)

αi=1npilogαpi+(1α)qipi(1α)i=1nqilogαpi+(1α)qiqi-\alpha \sum_{i=1}^n p_i \log \frac{\alpha p_i+(1-\alpha) q_i}{p_i}-(1-\alpha) \sum_{i=1}^n q_i \log \frac{\alpha p_i+(1-\alpha) q_i}{q_i}

讨论上面的式子,简记 αpi+(1α)qipi=piα\frac{\alpha p_i+(1-\alpha) q_i}{p_i}=p_i^\alpha

上式简写为

αi=1npilogpiαpi(1α)i=1nqilogqiαqi-\alpha \sum_{i=1}^n p_i \log \frac{p_i^\alpha}{p_i}-(1-\alpha) \sum_{i=1}^n q_i \log \frac{q_i^\alpha}{q_i}

i=1npilogpiαpilogi=1npipiαpi\sum_{i=1}^n p_i \log \frac{p_i^\alpha}{p_i} \leq \log \sum_{i=1}^n p_i \frac{p_i^\alpha}{p_i}i=1nqilogqiαqilogi=1nqiqiαqi\sum_{i=1}^n q_i \log \frac{q_i^\alpha}{q_i} \leq \log \sum_{i=1}^n q_i \frac{q_i^\alpha}{q_i}

上式αlogi=1npipiαpi(1α)logi=1nqiqiαqi=αlogi=1npiα(1α)logi=1nqiα=0\begin{aligned} \text{上式} \geq & -\alpha \log \sum_{i=1}^n p_i \frac{p_i^\alpha}{p_i}-(1-\alpha) \log \sum_{i=1}^n q_i \frac{q_i^\alpha}{q_i} \\ = & -\alpha \log \sum_{i=1}^n p_i^\alpha-(1-\alpha) \log \sum_{i=1}^n q_i^\alpha \\ = & 0 \end{aligned}

各种熵函数的互换关系

联合熵:H(XY)H(X Y)

信息熵:H(X)H(X)

条件熵:H(X/Y)H(X/Y)

三者关系

H(XY)=H(X)+H(Y/X)=H(Y)+H(X/Y)H(X Y)=H(X)+H(Y / X)=H(Y)+H(X / Y)

H(XY)H(X)+H(Y)H(X Y) \leq H(X)+H(Y)

H(X/Y)H(X),H(Y/X)H(Y)H(X / Y) \leq H(X), H(Y / X) \leq H(Y)

等号成立的关系的条件是 XXYY 相互独立。

从信息传输系统角度看熵的意义

  • H(X)H(X)XX 的信息熵,表示信源边每个符号的平均信息量(信源熵)。
  • H(Y)H(Y)YY 的信息熵,表示信宿边每个符号的平均信息量(信宿熵)。
  • H(X/Y)H(X/Y):条件熵,表示在信宿接收到 YY 后,信源 XX 尚存的平均不确定性。
  • H(Y/X)H(Y/X):噪声熵,表示在已知信源发出 XX 后,信宿 YY 尚存的平均不确定性。
  • H(XY)H(XY):联合熵,表示整个信息传输系统的平均不确定性。

熵函数的应用

信息增益(Information Gain):通常在决策树中,用来表示采用某种属性分类后不确定性的减少量,依据这个减少量衡量是否采用该属性进行分类:

IG=H(X)H(XY)I G=H(X)-H(X \mid Y)

构造好的决策树关键在于如何选择属性。最常用的分类属性选择指标就是信息增益。从候选属性中选择属性信息增益大的属性,其划分子集的纯度高(平均熵值小),分类能力强。

平均互信息

平均互信息:回忆互信息 I(xi;yj)I\left(x_i;y_j\right),定量描述输入随机变量发出某个具体消息 xix_i,输出变量出现某个具体消息 yjy_j,流经信道的信息量。随 xix_iyjy_j 变化而变化的随机变量,不能从整体上作为信道中信息流通的测度。需要从整体的角度出发,在平均意义上度量每通过一个符号流经信道的平均信息量。

在联合概率空间 P(XY)P(XY) 中统计平均值为 YYXX 的平均互信息量为:确定通过信道的信息量的多少,也称为信道传输率

I(X;Y)=jip(xiyj)I(xi;yj)=jip(xiyj)logp(xi/yj)p(xi)I(X ; Y)=\sum_j \sum_i p\left(x_i y_j\right) I\left(x_i ; y_j\right)=\sum_j \sum_i p\left(x_i y_j\right) \log \frac{p\left(x_i / y_j\right)}{p\left(x_i\right)}

  • I(X;Y)=H(X)H(X/Y)I(X ; Y)=H(X)-H(X / Y)
  • I(X;Y)=H(Y)H(Y/X)I(X ; Y)=H(Y)-H(Y / X)
  • I(X;Y)=H(X)+H(Y)H(XY)I(X ; Y)=H(X)+H(Y)-H(X Y)

性质:

非负性:通过一个信道总能传递一些信息,最差的条件下,输入输出完全独立,不传递任何信息,平均互信息等于 00,但决不会失去已知的信息。

I(X;Y)0I(X ; Y) \geq 0

对称性:立足的观测点不同。

I(X;Y)=I(Y;X)I(X ; Y)=I(Y ; X)

极值性:一般来说,平均互信息总是小于信源的熵,只有当信道是无损信道时,平均互信息才等于信源的熵。

I(X;Y)H(X)I(X ; Y) \leq H(X)

凸状性:对于 I(X;Y)I(X ; Y),是 P(X)P(X) 的上凸函数,是 P(Y/X)P(Y/X) 的下凸函数。

信道容量的理论基础:对于固定的信道,平均互信息 I(X;Y)I(X ; Y) 是信源概率分布 P(X)P(X) 的上凸函数。

信息率失真理论的基础:对于固定的信源,平均互信息 I(X;Y)I(X ; Y) 信道传递概率分布 P(Y/X)P(Y/X) 的下凸函数。

互信息的应用

信息不可增性原理:

I(X;Z)I(X;Y)I(X ; Z) \leq I(X ; Y)

测试系统:

I(X;Y1)I(X;Y1Y2)I(X;Y1Ym)I\left(X ; Y_1\right) \leq I\left(X ; Y_1 Y_2\right) \leq \cdots \leq I\left(X ; Y_1 \cdots Y_m\right)

信源编码:

I(X;Y)={Imax(X;Y) 质量问题 Imin (X;Y) 效率问题 I(X ; Y)=\left\{\begin{array}{l} I_{\max }(X ; Y) \Rightarrow \text { 质量问题 } \\ I_{\text {min }}(X ; Y) \Rightarrow \text { 效率问题 } \end{array}\right.

信号处理系统:

I(X;Y)H(X/Y)I(X ; Y) \uparrow \Longrightarrow H(X / Y) \downarrow

测试估计系统:

# 信息散度

https://zh.wikipedia.org/wiki/相对熵

序列的信息熵

平均符号熵

HL(XL)=1LH(X1X2XL)H_L\left(\overrightarrow{X_L}\right)=\frac{1}{L} H\left(X_1 X_2 \ldots X_L\right)

平稳性:序列的统计规律与统计的起始点无关。

p[x1(t),x2(t),xL(t)]=p[x1(t+τ),x2(t+τ),xL(t+τ)]p\left[x_1\left(t\right), x_2\left(t\right), \cdots x_L\left(t\right)\right]=p\left[x_1\left(t+\tau\right), x_2\left(t+\tau\right), \cdots x_L\left(t+\tau\right)\right]

遍历性(各态历经性):序列的集合平均等于它的时间平均特性。

E[X]=i=1npixi=limN1Ni=1NxiE[X]=\sum_{i=1}^n p_i x_i=\lim _{N \rightarrow \infty} \frac{1}{N} \sum_{i=1}^N x_i

当信源具备上述特性,可以定义有记忆信源的极限熵概念。

极限熵:对于平均符号熵 HL(XL)H_L\left(\overrightarrow{X_L}\right),当 LL \rightarrow \infty,平均符号熵取极限值,即

H(XL)=limL1LH(XL)=limLH(XL/X1X2XL1)H_{\infty}\left(\overrightarrow{X_L}\right)=\lim _{L \rightarrow \infty} \frac{1}{L} H\left(\overrightarrow{X_L}\right)=\lim _{L \rightarrow \infty} H\left(X_L / X_1 X_2 \ldots X_{L-1}\right)

也称作熵率。极限熵提供了一种不用统计联合概率计算平均符号熵的方法。

马尔可夫信源

mm 阶马尔可夫信源:在任何时刻,信源发出消息符号的概率仅与它前面发出的 mm
符号有关,而与其更前面的符号无关。

P(xL/xL1xL2xLm)P\left(x_L / x_{L-1} x_{L-2} \cdots x_{L-m}\right)状态转移概率

mm 个符号序列 Xm\overrightarrow{X_m} 状态:si={xL1xL2xLm}s_i=\left\{x_{L-1} x_{L-2} \cdots x_{L-m}\right\}

发出一个新符号 xLx_L

mm 个符号序列 Xm\overrightarrow{X_m} 状态:sj={xLxL1xLm+1}s_j=\left\{x_L x_{L-1} \cdots x_{L-m+1}\right\}

则有:

P(xL/xL1xL2xLm+1xLmx2x1)=P(xL/xL1xL2xLm+1xLm)=P(xL/si)=P(sj/si)\begin{aligned} & P\left(x_L / x_{L-1} x_{L-2} \ldots x_{L-m+1} x_{L-m} \ldots x_2 x_1\right) \\ = & P\left(x_L / x_{L-1} x_{L-2} \ldots x_{L-m+1} x_{L-m}\right) \\ = & P\left(x_L / s_i\right) \\ = & P\left(s_j / s_i\right) \end{aligned}

如果马尔可夫信源符号集共有 nn 个符号,则信源共有 nmn^m 个不同状态,即马尔可夫信源的状态空间:

[s1s2snmp(s1)p(s2)p(snm)]\left[\begin{array}{cccc} s_1 & s_2 & \cdots & s_{n^m} \\ p\left(s_1\right) & p\left(s_2\right) & \cdots & p\left(s_{n^m}\right) \end{array}\right]

信源的冗余度

后续更新…

重点总结

自信息I(ai)=deflogpiI\left(\boldsymbol{a}_{\boldsymbol{i}}\right) \stackrel{\mathrm{def}}{=}-\log \boldsymbol{p}_{\boldsymbol{i}}

(自)互信息I(ai;bj)= def logP(ai/bj)p(ai)=I(ai)I(ai/bj)I\left(\boldsymbol{a}_i ; \boldsymbol{b}_j\right) \stackrel{\text { def }}{=} \log \frac{P\left(a_i / b_j\right)}{p\left(a_i\right)}=I\left(a_i\right)-I\left(a_i / b_j\right)

信息熵H(X)= def E[I(ai)]=i=1npilogpiH(X) \stackrel{\text { def }}{=} E\left[I\left(a_i\right)\right]=-\sum_{i=1}^n p_i \log p_i

联合熵H(XY)=i=1nj=1mp(xiyj)I(xiyj)=i=1nj=1mp(xiyj)log2p(xiyj)H(X Y)=\sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right) I\left(x_i y_j\right)=-\sum_{i=1}^n \sum_{j=1}^m p\left(x_i y_j\right) \log _2 p\left(x_i y_j\right)

条件熵H(Y/X)=E[I(yj/xi)]=j=1mi=1np(xiyj)I(yj/xi)=j=1mi=1np(xiyj)logp(yj/xi)\begin{aligned} H(Y / X)=E\left[I\left(y_j / x_i\right)\right] & =\sum_{j=1}^m \sum_{i=1}^n p\left(x_i y_j\right) I\left(y_j / x_i\right)=-\sum_{j=1}^m \sum_{i=1}^n p\left(x_i y_j\right) \log p\left(y_j / x_i\right)\end{aligned}

熵函数性质

信息熵/联合熵/条件熵关系

平均互信息I(X;Y)=jip(xiyj)I(xi;yj)=jip(xiyj)logp(xi/yj)p(xi)I(X ; Y)=\sum_j \sum_i p\left(x_i y_j\right) I\left(x_i ; y_j\right)=\sum_j \sum_i p\left(x_i y_j\right) \log \frac{p\left(x_i / y_j\right)}{p\left(x_i\right)}

I(X;Y)=H(X)H(X/Y)=H(Y)H(Y/X)=H(X)+H(Y)H(XY)I(X ; Y)=H(X)-H(X / Y)=H(Y)-H(Y / X)=H(X)+H(Y)-H(X Y)

信息散度D(P//Q)= def XP(x)logP(x)Q(x)D(P / / Q) \stackrel{\text { def }}{=} \sum_X P(x) \log \frac{P(x)}{Q(x)}

数据处理定理

I(X;Z)I(Y;Z),I(X;Z)I(X;Y)I(X;Y1)I(X;Y1Y2)I(X;Y1Ym)\begin{aligned} & I(X ; Z) \leq I(Y ; Z), I(X ; Z) \leq I(X ; Y) \\ & I\left(X ; Y_1\right) \leq I\left(X ; Y_1 Y_2\right) \leq \cdots \leq I\left(X ; Y_1 \cdots Y_m\right) \end{aligned}