香浓编码实验报告.docVIP

　香农编码实验报告 姓名：徐以刚 学号 专业班级：信计09.1 学院：理信学院 一 、实验目的 1. 了解香农编码的基本原理及其特点； 2. 熟悉掌握香农编码的方法和步骤； 3. 掌握C语言或者Matlab编写香农编码的程序。 二、实验要求 对于给定的信源的概率分布,按照香农编码的方法进行计算机实现. 三、实验原理 给定某个信源符号的概率分布，通过以下的步骤进行香农编码 1．信源符号按概率从大到小排列 2. 对信源符号求累加概率，表达式： Gi=Gi-1+p(xi) 3. 求自信息量，确定码字长度。自信息量I(xi)=-log(p(xi));码字长度取大于等于自信息量的最小整数。 4. 将累加概率用二进制表示，并取小数点后码字的长度的码 。 四、实验内容 离散无记忆信源符号S的概率分布： S P(S) = 0.20 0.19 0.18 0.17 0.15 0.10 0.01 对离散无记忆信源分布S进行香农编码 1.画出程序设计的流程图 2.写出程序代码， N=input(N=); %输入信源符号的个数 s=0; l=0; H=0; for i=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 if abs(s-1)0, error(不符合概率分布) end for i=1:N-1 for j=i+1:N if p(i)p(j) m=p(j); p(j)=p(i); p(i)=m; end end end %按概率分布大小对信源排序 for i=1:N a=-log2(p(i)); if mod(a,1)==0 w=a; else w=fix(a+1); end %计算各信源符号的码长 l=l+p(i)*w; %计算平均码长 end l=l; n=H/l; %计算编码效率 P(1)=0 for i=2:N P(i)=0; for j=1:i-1 P(i)=P(i)+p(j); end end %计算累加概率 for i=1:N for j=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.3219 2.3959 2.4739 2.5564 2.7370 3.3219 6.6439 所以我们得到每个信源符号的自信息量为 根据公式,我们得到每个信源符号的码长为 由程序运行结果， 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 1 1 0 0 0 1 1 0 0 1 0 0 0 1 0 1 1 1 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 0 我们得到每个信源符号的为对应的二进制数为