实验三香农编码

《实验三香农编码》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、信息论与编码实验报告院系：哈尔滨理工大学荣成校区专业：电子信息工程学号：姓名：日期：2015年6月16日精品学习资料可选择pdf第1页，共5页-----------------------香农编码信息论与编码第三次实验报告一、实验目的和任务1、理解信源编码的意义；2、熟悉MATLAB程序设计；3、掌握香农编码的方法及计算机实现；4、对给定信源进行香农编码，并计算编码效率；二、实验原理介绍给定某个信源符号的概率分布，通过以下的步骤进行香农编码1、信源符号按概率从大到小排列；p1p2......pn2、确定满足下列不等式的整数码长Ki为-lb(pi)Kilb(pi)13、

2、为了编成唯一可译码，计算第i个消息的累加概率：l1pip(aK)K14、将累加概率Pi变换成二进制数；5、取Pi二进制数的小数点后Ki位即为该消息符号的二进制码字。三、实验设备1、计算机2、Matlab2014四、编码程序N=input('N=');%输入信源符号的个数s=0;%对行的概率进行初始化l=0;%对列的概率初始化精品学习资料可选择pdf第2页，共5页-----------------------H=0;%对信息熵初始化fori=1:Np(i)=input('p=');%输入信源符号概率分布矢量p(i)<1;%输入概率小于1s=s+p(i);%行概率求和H=

3、H+(-p(i)*log2(p(i)));%计算信源信息熵endifabs(s-1)>0%若行概率和>1则报错error('不符合概率分布')endfori=1:N-1%for循环按概率分布大小对信源排序forj=i+1:Nifp(i)

4、=l;%香农编码的信息率（这里L=1,m=2）n=H/r;%计算编码效率P(1)=0;%初始1列为0；fori=2:NP(i)=0;%初始1列为0；累加forj=1:i-1P(i)=P(i)+p(j);%从2列开始累加endend%计算累加概率fori=1:Nforj=1:wW(i,j)=fix(P(i)*2);%对每次的P（i）取码字,长度与a相同P(i)=P(i)*2-fix(P(i)*2);%将累加概率转化为L(i)位二进制码字endend精品学习资料可选择pdf第3页，共5页-----------------------W%显示码字r%显示平均码长n%显示编码

5、效率H%显示信息熵H五、实验结果>>xiangnongN=6p=0.25p=0.25p=0.20p=0.15p=0.10p=0.05W=000000100010000101101101111110r=2.7000n=0.8975H=2.4232精品学习资料可选择pdf第4页，共5页-----------------------六、实验心得此次实验让我明白了做程序需要先进行需求分析，然后画出流程图，根据流程图把每部分模块完成，然后再进行调试和编辑，提高了我做编程的能力精品学习资料可选择pdf第5页，共5页-----------------------