费诺编码的matlab实现

...wd......wd......wd...多媒体技术实验报告学院：城南学院姓名：学号：指导教师：尹波时间：2015年11月25日教师评语：教师评语：成绩评阅教师日期实验一：费诺编码的matlab实现1．实验目的掌握费诺编码的思想和具体方法。用MATLAB语言实现费诺编码。2.实验原理及编码思想：费诺编码属于概率匹配编码，但不是最正确的编码方法。在编N进制码时首先将信源消息符号按其出现的概率依次由大到小排列开来，并将排列好的信源符号按概率值分N大组，使N组的概率之和近似一样，并对各组赋予一个N进制码元0、1……N-1。之后再针对每一大组内的信源符号做如上的处理，即再分为概率和一样的N组，赋予N进制码元。如此重复，直至每组只剩下一个信源符号为止。此时每个信源符号所对应的码字即为费诺码。具体过程如下：[1]将信源消息符号按其出现的概率大小依次排列：P1>=P2>=…>=Pn。[2]依次排列的信源符号按概率值分为两大组，使两个组的概率之和近似一样，并对各组赋予一个二进制码元“0”和“1[3]使划分后的两个组的概率之和近似一样，并对各组赋予一个二进制符号“0”和“1[4]如此重复，直至每个组只剩下一个信源符号为止。[5]信源符号所对应的码字即为费诺码。例：有一幅40个像素组成的灰度图像，灰度共有5级，分别用符号A，B，C，D和E表示。40个像素中出现灰度A的像素数有15个，出现灰度B的像素数有7个，出现灰度C的像素数有7个，其余情况见表。符号ABCDE出现的次数157765编码00011010110111费诺编码方法属于概率匹配编码，具有如下特点：1、概率大，则分解次数小；概率小则分解次数多。这符合最正确码原则。2、码字集合是唯一的。3、分解完了，码字出来了，码长也有了,即先有码字后有码长。因此，费诺编码方法又称为子集分解法。3.程序设计思路4.程序代码clc;clear;A=[0.19,0.18,0.17,0.16,0.13,0.10,0.06,0.01];A=fliplr(sort(A));%降序排列[m,n]=size(A);fori=1:nB(i,1)=A(i);%生成B的第1列end%生成B第2列的元素a=sum(B(:,1))/2;fork=1:n-1ifabs(sum(B(1:k,1))-a)<=abs(sum(B(1:k+1,1))-a)break;endendfori=1:n%生成B第2列的元素ifi<=kB(i,2)=0;elseB(i,2)=1;endend%生成第一次编码的结果END=B(:,2)';END=sym(END);%生成第3列及以后几列的各元素j=3;while(j~=0)p=1;while(p<=n)x=B(p,j-1);forq=p:nifx==-1break;elseifB(q,j-1)==xy=1;continue;elsey=0;break;endendendify==1q=q+1;endifq==p|q-p==1B(p,j)=-1;elseifq-p==2B(p,j)=0;END(p)=[char(END(p)),'0'];B(q-1,j)=1;END(q-1)=[char(END(q-1)),'1'];elsea=sum(B(p:q-1,1))/2;fork=p:q-2ifabs(sum(B(p:k,1))-a)<=abs(sum(B(p:k+1,1))-a);break;endendfori=p:q-1ifi<=kB(i,j)=0;END(i)=[char(END(i)),'0'];elseB(i,j)=1;END(i)=[char(END(i)),'1'];endendendendp=q;endC=B(:,j);D=find(C==-1);[e,f]=size(D);ife==nj=0;elsej=j+1;endendBAENDfori=1:n[u,v]=size(char(END(i)));L(i)=v;endavlen=sum(L.*A)H=-sum(A.*log2(A));disp(['ENtropyofthesource:',