信息论课程设计

《信息论与编码》实验报告

班级：

学号：

姓名：sweety

信息与通信工程学院信息工程系

目录

实验一绘制信源熵函数曲线错误!未定义书签。

实验二哈夫曼编解码错误!未定义书签。

实验三离散信道容量错误!未定义书签。

实验一绘制信源熵函数曲线

一、实验目的

掌握离散信源熵的原理和计算方法。

熟悉matlab软件的基本操作，练习应用matlab软件进行信源熵函数曲线的绘制。理解信源熵的物理意义，并能从信源熵函数曲线图上进行解释其物理意义。

二、实验原理

离散信源相关的基本概念、原理和计算公式

产生离散信息的信源称为离散信源。离散信源只能产生有限种符号。

假定X是一个离散随机变量，即它的取值范围R={x1，x2，x3，…}是有限或可数的。设第i个变量x i发生的概率为p i=P{X=x i}。则：

定义一个随机事件的自信息量I（x i）为其对应的随机变量x i出现概率对数的负值。即：

I（x i）= -log2p(x i)

定义随机事件X的平均不确定度H（X）为离散随机变量x i出现概率的数学期望，即：

∑∑

-

=

=

i i

i i

i i

x

p

x

p

x

I

x

p

X

H)

(

log

)

(

)

(

)

(

)

(

单位为比特/符号或比特/符号序列。

平均不确定度H（X）的定义公式与热力学中熵的表示形式相同，所以又把平均不确定度H（X）称为信源X的信源熵。

必须注意一下几点：

某一信源，不管它是否输出符号，只有这些符号具有某些概率特性，必有信源的熵值；这熵值是在总体平均上才有意义，因而是个确定值，一般写成H（X），X 是指随机变量的整体（包括概率分布）。

信息量则只有当信源输出符号而被接收者收到后，才有意义，这就是给与信息者的信息度量，这值本身也可以是随机量，也可以与接收者的情况有关。

熵是在平均意义上来表征信源的总体特征的，信源熵是表征信源的平均不确定度，平均自信息量是消除信源不确定度时所需要的信息的量度，即收到一个信源符号，全部解除了这个符号的不确定度。或者说获得这么大的信息量后，信源不确定度就被消除了。信源熵和平均自信息量两者在数值上相等，但含义不同。 当某一符号x i 的概率p (x i )为零时，p (x i )log p (x i ) 在熵公式中无意义，为此规定这时的 p (x i )log p (x i ) 也为零。当信源X 中只含有一个符号x 时，必有p (x )=1，此时信源熵H （X ）为零。

例1-1，设信源符号集X ={0，1}，每个符号发生的概率分别为p (0)=p ，p (1)=q ，p+ q =1，即信源的概率空间为

??

????=?????? 1 0q p P X 则该二元信源的信源熵为：

H (X ) = - p log p – q log q = - p log p – (1- p ) log (1- p)

即：H (p) = - p log p – (1- p ) log (1- p)

其中0 ≤ p ≤1

P =0时，H (0) = 0

P =1时，H (1) = 0

MATLAB 二维绘图

例对函数y = f (x )进行绘图，则用matlab 中的命令plot(x , y )就可以自动绘制出二维图来。如果打开过图形窗口，则在最近打开的图形窗口上绘制此图；如果未打开图形窗口，则开一个新的图形窗口绘图。

例1-2，在matlab 上绘制余弦曲线图，y = cos x ，其中0 ≤ x ≤ 2。 >>x =0::2*pi ； %生成横坐标向量，使其为0，，，…，

>>y =cos(x )； %计算余弦向量

>>plot(x ,y ) %绘制图形

三、实验内容

用matlab 软件绘制二源信源熵函数曲线。根据曲线说明信源熵的物理意义。

四、实验要求

提前预习实验，认真阅读实验原理以及相应的参考书。

认真高效的完成实验，实验中服从实验室管理人员以及实验指导老师的管理。 认真填写实验报告。

五、实验代码：

p=::1;

h=-(p).*log2(p)-(1-p).*log2(1-p);

plot(p,h)

hold on ;

stem,1,'--');

xlabel('p')

title('二元信源熵')

六、实验结果：

七、物理意义：

信源熵为信息的不确定度性，而概率的大小决定了信息量的大小。

由上图可知概率为1时，信息量越小，不确定性最低。概率等于时熵最大。

实验二哈夫曼编码

一、实验目的

掌握哈夫曼编码的原理及编码步骤

练习matlab中哈夫曼编码函数的调用及通信工具箱的使用

二、实验原理

通信的根本问题是如何将信源输出的信息在接收端的信宿精确或近似的复制出来。为了有效地复制信号，就通过对信源进行编码，使通信系统与信源的统计特性相匹配。

若接收端要求无失真地精确地复制信源输出的信息，这样的信源编码即为无失真编码。即使对于一个小的时间段内，连续信源输出的信息量也可以是无限大的，

所以对其是无法实现无失真编码的；而离散信源输出的信息量却可以看成是有限的，所以只有离散信源才可能实现无失真编码。

凡是能载荷一定的信息量，且码字的平均长度最短，可分离的变长码的码字集合都可以称为最佳码。为此必须将概率大的信息符号编以短的码字，概率小的符号编以长的码字，使得平均码字长度最短。

变字长编码的最佳编码定理：在变字长码中，对于概率大的信息符号编以短字长的码；对于概率小的信息符号编以长字长的码。如果码字长度严格按照符号概率的大小顺序排列，则平均码字长度一定小于俺任何顺序排列方式得到的码字长度。

哈夫曼编码就是利用了这个定理，讲等长分组的信源符号，根据其概率分布采用不等长编码。概率大的分组，使用短的码字编码；概率小的分组，使用长的码字编码。哈夫曼编码把信源按概率大小顺序排列，并设法按逆次序分配码字的长度。在分配码字的长度时，首先将出现概率最小的两个符号相加，合成一个概率；第二步把这个合成的概率看成是一个新组合符号的概率，重复上述做法，直到最后只剩下两个符号的概率为止。完成以上概率相加顺序排列后，再反过来逐步向前进行编码。每一步有两个分支，各赋予一个二进制码，可以对概率大的编为0码，概率小的编为1码。反之亦然。

哈夫曼编码的具体步骤归纳如下：

统计n 个信源消息符号，得到n 个不同概率的信息符号。

将这n 个信源信息符号按其概率大小依次排序：

p (x 1) ≥ p (x 2)≥ …≥ p (x n )

取两个概率最小的信息符号分别配以0和1两个码元，并将这两个概率相加作为一个新的信息符号的概率，和未分配的信息符号构成新的信息符号序列。 将剩余的信息符号，按概率大小重新进行排序。

重复步骤3，将排序后的最后两个小概论相加，相加和与其他概率再排序。 如此反复重复n-2次，最后只剩下两个概率。

从最后一级开始，向前返回得到各个信源符号所对应的码元序列，即相应的码字，构成霍夫曼编码字。编码结束。

哈夫曼编码产生最佳整数前缀码，即没有一个码字是另一个码字的前缀，因此哈夫曼编码是唯一码。

编码之后，哈夫曼编码的平均码长为：

1()n i

i i K p x K ==∑

哈夫曼编码的效率为：

()=H x K

η=信源熵平均码长 例

该哈夫曼码的平均码长为

编码效率为： 三、实验内容

为某一信源进行哈夫曼编码。该信源的字符集为X={x 1, x 2, … x 6 }，相应的概率矢量为：P =,,,,,，即X ，P 的概率空间为：

123456 0.300.250.210.100.090.05X x x x x x x P ????=????????

根据哈夫曼编码算法对该信源进行哈夫曼编码。并计算其平均码长和编码效率。 调用matlab 哈夫曼编码函数进行哈夫曼编码，与人工编码结果做比较。 huffmandict 函数： 为已知概率分布的信源模型生成哈夫曼编解码索引表。 调用方法如下：

[dict ，avglen] = huffmandict (symbols, p)

[dict ，avglen] = huffmandict (symbols, p, N)

[dict ，avglen] = huffmandict (symbols, p, N, variance)

四、实验要求

提前预习实验，认真阅读实验原理以及相应的参考书。

认真高效的完成实验，实验中服从实验室管理人员以及实验指导老师的管理。 认真填写实验报告。

五、实验代码：

symbols = [1:6];

p = [.30 .25 .21 .10 .09 .05];

71

() 2.72i i i K p x K ===∑码元/符号() 2.610.9596/2.72H X K η===比特码元

[dict,avglen] = huffmandict(symbols,p)

temp=dict;

for i = 1:length(temp)

temp{i,2}= num2str(temp{i,2});

end

h=.*log2.*log2.*log2.*log2.*log2.*log2

n=h/avglen

temp

六、实验结果：

symbols=

1 2 3 4 5 6

P=

dict=

[1] [1x2 double]

[2] [1x2 double]

[3] [1x2 double]

[4] [1x3 double]

[5] [1x4 double]

[6] [1x4 double]

avglen=

n=

temp=

[1] ‘0 0’

[2] ‘0 1’

[3] ‘1 1’

[4] ‘1 0 1’

[5] ‘1 0 0 0’

[6] ‘1 0 0 1’

人工编码：其平均码长和编码效率分别为和%，与调用MATLAB哈夫曼树编码函数进行哈夫曼编码结果一致，但调用函数处理速度更快，效率更高，在进行大量数据运算时。优势明显。

实验三离散信道容量

一、实验目的

掌握离散信道容量的计算。

理解离散信道容量的物理意义。

练习应用matlab软件进行二元对称离散信道容量的函数曲线的绘制，并从曲线上理解其物理意义。

信道是传送信息的载体—信号所通过的通道。

信息是抽象的，而信道则是具体的。比如二人对话，二人间的空气就是信道；打电话，电话线就是信道；看电视，听收音机，收、发间的空间就是信道。 研究信道的目的：在通信系统中研究信道，主要是为了描述、度量、分析不同类型信道，计算其容量，即极限传输能力，并分析其特性。

二元对称信道BSC （Binary Symmetric Channel ）

二进制离散信道模型有一个允许输入值的集合X={0，1}和可能输出值的集合Y={0，1}，以及一组表示输入和输出关系的条件概率（转移概率）组成。如果信道噪声和其他干扰导致传输的二进序列发生统计独立的差错，且条件概率对称，即

(0/1)(1/0)(1/1)(0/0)1p Y X p Y X p p Y X p Y X p

======??======-? 这种对称的二进制输入、二进制输出信道称做二元对称信道（或二进制对称信道，简称BSC 信道），如下图所示：

信道容量公式：

{()}

max p x C I(X,Y)=

三、实验内容

BSC 信道是DMC 信道对称信道的特例，对于转移概率为P(0/1)=P(1/0)=p ，P(0/0)=P(1/01)=1-p ，求出其信道容量公式，并在matlab 上绘制信道容量C 与p 的曲线。

根据曲线说明其物理意义。

四、实验要求

提前预习实验，认真阅读实验原理以及相应的参考书。

认真高效的完成实验，实验中服从实验室管理人员以及实验指导老师的管理。 认真填写实验报告

五、实验代码：

p=linspace,,50);

c=log2(2)+p.*log2(p)+(1-p).*log2(1-p);

plot(p,c);

title('BSC 信道')

由图可知：p=0时信道容量是1bit/符号；当P=时，从输出得不到任何关于输入的信息，互信息为0即信道容量为0；信道容量以P=点为中心对称。

物理意义：P=0时，无误码概率，从接受的Y可以完全确定发送的X，信道的介入没有产生任何的损伤。P=时，说明损失的信息达到与输入符号信息熵相等的程度。可得C=0时，即信道上没能传送任何信息。