香浓编码实验.docVIP

实验一香农编码的实验报告

一、实验目的

1.了解香农编码的根本原理及其特点；

2.熟悉掌握香农编码的方法和步骤；

3.掌握C语言或者Matlab编写香农编码的程序。

二、实验要求

对于给定的信源的概率分布,按照香农编码的方法进行计算机实现.

三、实验原理

给定某个信源符号的概率分布，通过以下的步骤进行香农编码

1．信源符号按概率从大到小排列

2.对信源符号求累加概率，表达式：Gi=Gi-1+p(xi)

3.求自信息量，确定码字长度。自信息量I(xi)=-log(p(xi));码字长度取大于等于自信息量的最小整数。

4.将累加概率用二进制表示，并取小数点后码字的长度的码。

四、实验内容

离散无记忆信源符号S的概率分布：

S

P(S)=0.200.190.180.170.150.100.01

对离散无记忆信源分布S进行香农编码

1.画出程序设计的流程图，

开始计算平均码长按概率分布大小对信源排序判断概率和S是否等于1输入对应符号的概率P〔i〕输入信源符号的个数

开始

计算平均码长

按概率分布大小对信源排序

判断概率和S是否等于1

输入对应符号的概率P〔i〕

输入信源符号的个数

计算编码效率

计算编码效率

计算累加概率

计算累加概率

计算自信息量

计算自信息量

将累加概率转化为二进制码字

将累加概率转化为二进制码字

输出累加概率，码长，码字，自信息量

输出累加概率，码长，码字，自信息量，平均码长，编码效率

结束

结束

2.写出程序代码，

N=input(N=);%输入信源符号的个数

s=0;

l=0;

H=0;

fori=1:N

p(i)=input(p=);%输入信源符号概率分布矢量,p(i)1

s=s+p(i)

H=H+(-p(i)*log2(p(i)));I(i)=-log2(p(i));%计算信源信息熵

end

ifabs(s-1)0,

error(不符合概率分布)

end

fori=1N-1

forj=i+1:N

ifp(i)p(j)

m=p(j);

p(j)=p(i);

p(i)=m;

end

end

end%按概率分布大小对信源排序

fori=1:N

a=-log2(p(i));

ifmod(a1)==0

w=a;

else

w=fix(a+1);

end%计算各信源符号的码长

l=l+p(i)*w;%计算平均码长

end

l=l;

n=H/l;%计算编码效率

P(1)=0

fori=2:N

P(i)=0;

forj=1:i-1

P(i)=P(i)+p(j);

end

end%计算累加概率

fori=1:N

forj=1:w

W(i,j)=fix(P(i)*2);

P(i)=P(i)2-fix(P(i)*2);

end

end%将累加概率转化为L(i)位二进制码字

disp(W)%显示码字

disp(l)%显示平均码长

disp(n)%显示编码效率

disp(I)%显示自信息量

3.写出在调试过程中出现的问题，

问题1：自信量程序不会编写

问题2：累加概率时注意P(1)=0

问题3：程序运行时要依次输入各个符号概率

4.对实验的结果进行分析

由程序运行结果，得

2.32192.39592.47392.55642.73703.32196.6439

所以我们得到每个信源符号的自信息量为

根据公式,我们得到每个信源符号的码长为

由程序运行结果，

0000000

0011001

0110001

1