matlab费诺编码程序.docVIP

%函数f1存放于f1.m function x=f1(i,j,p,r) global x; x=char(x); if(j=i) return; else q=0; for t=i:j %对于区间[i,j]自上而下求累加概率值 q=p(t)+q;y(t)=q; end for t=i:j%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组 v(t)=abs(y(t)-(q-y(t))); end for t=i:j if(v(t)==min(v)) %求该数组中最小的一个值来确定分界点位置 for k=i:t %赋值码字 x(k,r)=0; end for k=(t+1):j x(k,r)=1; end d=t; f1(i,d,p,r+1); %递归调用及相互调用 f2(d+1,j,p,r+1); f1(d+1,j,p,r+1); f2(i,d,p,r+1); else end end end return; %函数f2存放于f2.m function x=f2(i,j,p,r) global x; x=char(x); if(j=i) return; else q=0; for t=i:j %对于区间[i,j]自上而下求累加概率值 q=p(t)+q;y(t-i+1)=q; end for t=1:j-(i-1)%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组 v(t)=abs(y(t)-(q-y(t))); end for t=1:j-(i-1) if(v(t)==min(v)) %求该数组中最小的一个值来确定分界点位置 d=t+i-1; for k=i:d %赋值码字 x(k,r)=0; end for k=(d+1):j x(k,r)=1; end f2(d+1,j,p,r+1);%递归调用及相互调用 f1(i,d,p,r+1); f2(i,d,p,r+1); f1(d+1,j,p,r+1); else end end end return; 主函数 clc;clear all; fprintf(………………费诺编码程序………………\n ); fprintf(小组成员：苏、杨、张\n \n ); fprintf(请输入信源符号的个数：); N=input(N=);%输入信源符号的个数 s=0;l=0;H=0; for i=1:N fprintf(请输入第%d个符号的概率：,i); p(i)=input(p=);%输入信源符号概率分布矢量,0p(i)1 if p(i)=0|| p(i)=1 error(请注意P的范围是0P1) end s=s+p(i) H=H+(- p(i)*log2(p(i)));%计算信源信息熵 end if (s~=1) error(信源符号概率和不等1) end tic; 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 x=f1(1,N,p,1); for i=1:N %计算平均码长 L(i)=length(find(x(i,:))); l=l+p(i)*L(i); end n=H/l; %计算编码效率 fprintf(编码后所得码字:\n); disp(x) %显示按概率降序排列的码字 fprintf(平均码长:K=\n); disp(l)% 显示平均码长 fprintf(信息熵:H(X)=\n); disp(H)%显示信息熵 fprintf(编码效率: η=\n); disp(n) %显示编码效率 %fprintf(计算耗时time= %f\n,toc);