资源描述:
《matlab费诺编码程序.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、%函数f1存放于f1.mfunctionx=f1(i,j,p,r)globalx;x=char(x);if(j<=i)return;elseq=0;fort=i:j%对于区间[i,j]自上而下求累加概率值q=p(t)+q;y(t)=q;endfort=i:j%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组v(t)=abs(y(t)-(q-y(t)));endfort=i:jif(v(t)==min(v))%求该数组中最小的一个值来确定分界点位置fork=i:t%赋值码字x(k,r)='0';endfork=
2、(t+1):jx(k,r)='1';endd=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);elseendendendreturn;%函数f2存放于f2.mfunctionx=f2(i,j,p,r)globalx;x=char(x);if(j<=i)return;elseq=0;fort=i:j%对于区间[i,j]自上而下求累加概率值q=p(t)+q;y(t-i+1)=q;endfort=1:j-(i-1)%把所有自上而下的累加概率值与该区间
3、总概率值减该累加概率值之差取绝对值存在一数组v(t)=abs(y(t)-(q-y(t)));endfort=1:j-(i-1)if(v(t)==min(v))%求该数组中最小的一个值来确定分界点位置d=t+i-1;fork=i:d%赋值码字x(k,r)='0';endfork=(d+1):jx(k,r)='1';endf2(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);elseendendendreturn;主函数clc;clearall;fprint
4、f('………………费诺编码程序………………');fprintf('小组成员:苏、杨、张');fprintf('请输入信源符号的个数:');N=input('N=');%输入信源符号的个数s=0;l=0;H=0;fori=1:Nfprintf('请输入第%d个符号的概率:',i);p(i)=input('p=');%输入信源符号概率分布矢量,0
5、
6、p(i)>=1error('请注意P的范围是0
7、(s~=1)error('信源符号概率和不等1')endtic;fori=1:N-1%按概率分布大小对信源排序forj=i+1:Nifp(i)
8、printf('信息熵:H(X)=');disp(H)%显示信息熵fprintf('编码效率:η=');disp(n)%显示编码效率%fprintf('计算耗时time=%f',toc);