实验六 信息编码.docVIP

实验六 香农—费诺—埃利斯码的ＭＡＴＬＡＢ语言实现 实验学时：3 实验类型：（演示、验证、综合、√设计、研究） 实验要求：（√必修、选修） 一、实验目的 通过本次试验了解香农—费诺—埃利斯码的编码方法和原理并用matlab实现。 二、实验内容 编程实现香农—费诺—埃利斯码算法 三、实验原理、方法和手段 1. 编码原理：香农—费诺—埃利斯码不是分组码，它根据信源符号的积累概率分配码字，不是最佳码，但它的编码和译码效率都很高。 步骤1、求修正累计概率。 2、转换成二进制小数形式。 3、求码长。 4、截取码字。 2. 编码如下： p=input(please input a number:); %提示输入界面 n=length(p); for i=1:n if p(i)0 fprintf(\n The probabilities in huffman can not less than 0!\n); p=input(please input a number:) %如果输入的概率数组中有小于0的值，则重新输入概率数组 end end if abs(sum(p)-1)0 fprintf(\n The sum of the probabilities in huffman can more than 1!\n); p=input(please input a number:) %如果输入的概率数组总和大于1，则重新输入概率数组 end A=p; A=fliplr(sort(A));%降序排列 [m,n]=size(A); for i=1:n B(i,1)=A(i);%生成B的第1列 end %生成B第2列的元素 a=sum(B(:,1))/2; for k=1:n-1 if abs(sum(B(1:k,1))-a)=abs(sum(B(1:k+1,1))-a) break; end end for i=1:n %生成B第2列的元素 if i=k B(i,2)=0; else B(i,2)=1; end end %生成第一次编码的结果 END=B(:,2); END=sym(END); %生成第3列及以后几列的各元素 j=3; while (j~=0) p=1; while(p=n) x=B(p,j-1); for q=p:n if x==-1 break; else if B(q,j-1)==x y=1; continue; else y=0; break; end end end if y==1 q=q+1; end if q==p|q-p==1 B(p,j)=-1; else if q-p==2 B(p,j)=0; END(p)=[char(END(p)),0]; B(q-1,j)=1; END(q-1)=[char(END(q-1)),1]; else a=sum(B(p:q-1,1))/2; for k=p:q-2 if abs(sum(B(p:k,1))-a)=abs(sum(B(p:k+1,1))-a); break; end end for i=p:q-1 if i=k B(i,j)=0; END(i)=[char(END(i)),0]