(完整word版)实验三香农编码_第1页
(完整word版)实验三香农编码_第2页
(完整word版)实验三香农编码_第3页
(完整word版)实验三香农编码_第4页
(完整word版)实验三香农编码_第5页
已阅读5页，还剩2页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、信息论与编码实验报告院系：哈尔滨理工大学荣成校区专业：电子信息工程学号：1姓名：日期：2015年6月16日香农编码信息论与编码第三次实验报告一、实验目的和任务1、理解信源编码的意义；2、熟悉MATLAB程序设计；3、掌握香农编码的方法及计算机实现；4、对给定信源进行香农编码，并计算编码效率；二、实验原理介绍给定某个信源符号的概率分布，通过以下的步骤进行香农编码1、信源符号按概率从大到小排列；Pl - P2 - -Pn2、确定满足下列不等式的整数码长Ki为-lb(Pi) < Ki < -lb(Pi) 13、为了编成唯一可译码，计算第i个消息的累加概率：l 1PiP(aK)K 14、将

2、累加概率P变换成二进制数；5、取P二进制数的小数点后 Ki位即为该消息符号的二进制码字。三、实验设备1、计算机2、Matlab2014四、编码程序N=input( 'N='); 瀚入信源符号的个数s=0;%对行的概率进行初始化1=0;%对列的概率初始化H=0;%寸信息嫡初始化for i=1:Np(i尸input( 'p=');瀚入信源符号概率分布矢量 p(i)<1;淅入概率小于1s=s+p(i);%亍概率求和H=H+(-p(i)*10g2(p(i);%计算信源信息嫡end if abs(s-1)>0*行概率和>1则报错error('不符

3、合概率分布')endfor i=1:N-1for j=i+1:Nif p(i)<p(j) m=p(j);p(j)=p(i);p(i)=m;endend%for循环按概率分布大小对信源排序%冒泡排序endfor i=1:Na=-log2(p(i);if mod(a,1)=0%确定满足不等式的整数 a,并领a为第i个码字的长度%a对1取余后赋值给aw=a;elsew=fix(a+1);end%计算各信源符号的码长l=l+p(i)*w;%计算平均码长endr=l;n=H/r;P(1)=0;for i=2:NP(i)=0;for j=1:i-1P(i)=P(i)+p(j);endendf

4、or i=1:Nfor j=1:wW(i,j)=fix(P(i)*2);P(i)=P(i)*2-fix(P(i)*2);end%香农编码的信息率(这里L=1,m=2 )% 计算编码效率%初始1列为0;%初始1列为0 ;累加%从2列开始累加% 十算累加概率%寸每次的P (i )取码字，长度与a相同%等累加概率转化为 L(i)位二进制码字endW%显示码字r%显示平均码长n%显示编码效率H%显示信息嫡五、实验结果>> xiangnongN=6p=0.25p=0.25p=0.20p=0.15p=0.10p=0.05W =000000100010000101101101111110r =2.7000 n =0.89752.4232六、实验心得此

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论