matlab费诺编码程可实现

函数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=(t+1):jx(k,r)='1';endd=t;f1(i,d,p,r+1);%递归调用及相互调用f2(d+1,j,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)%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组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);elseendendendreturn;主函数clc;clearall;fprintf('请输入信源数组：');p=input('P=');%输入信源数组s=sum(p(1,:));if(s~=1)error('信源符号概率和不等1')end[M,N]=size(p);fori=1:N-1%按概率分布大小对信源排序forj=i+1:Nifp(i)p(j)m=p(j);p(j)=p(i);p(i)=m;endendendx=f1(1,N,p,1);fprintf('编码后所得码字:\n');disp(x)%显示按概率降序排列的码字