matlab費諾編碼程序

1、%S數(shù)fl存放于fl.mfunctionx=f1(i,j,p,r)globalx;x=char(x);if(j<=i)return;elseq=0；fort=i:j%對于區(qū)間i,j自上而下求累加概率值q=p(t)+q;y(t尸q;endfort=i:j%把所有自上而下的累加概率值與該區(qū)間總概率值減該累加概率值之差取絕對值存在一數(shù)組v(t)=abs(y(t)-(q-y(t);endfort=i:jif(v(t)=min(v)%求該數(shù)組中最小的一個值來確定分界點位置fork=i:t%賦值碼字x(k,r)='0'endfork=(t+1):jx(k,r)='1'

2、endd=t;f1(i,d,p,r+1);%遞歸調(diào)用及相互調(diào)用f2(d+1,j,p,r+1);f1(d+1,j,p,r+1);f2(i,d,p,r+1);elseendendendreturn;%函數(shù)f2存放于f2.mfunctionx=f2(i,j,p,r)globalx;x=char(x);if(j<=i)return;elseq=0;fort=i:j%對于區(qū)間i,j自上而下求累加概率值q=p(t)+q;y(t-i+1)=q;endfort=1:j-(i-1)%把所有自上而下的累加概率值與該區(qū)間總概率值減該累加概率值之差取絕對值存在一數(shù)組v(t)=abs(y(t)-(q-y(t);e

3、ndfort=1:j-(i-1)if(v(t)=min(v)%求該數(shù)組中最小的一個值來確定分界點位置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);%遞歸調(diào)用及相互調(diào)用f1(i,d,p,r+1);f2(i,d,p,r+1);f1(d+1,j,p,r+1);elseendendendreturn;主函數(shù)clc;clearall;fprintf('費諾編碼程序n');fprintf('小組成員：蘇、楊、張nn');fprintf(&#

4、39;請輸入信源符號的個數(shù)：');N=input('N=');%輸入信源符號的個數(shù)s=0;l=0;H=0;fori=1:Nfprintf('請輸入第符號的概率：',i);p(i)=input('p=');%輸入信源符號卞S率分布矢量，0<p(i)<1ifp(i)<=0|p(i)>=1error('請注意P的范圍是0<P<1')ends=s+p(i)H=H+(-p(i)*log2(p(i);%計算信源信息嫡endif(s=1)error('信源符號概率和不等1')endti

5、c;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);fori=1:N%計算平均碼長L(i)=length(find(x(i,:);l=l+p(i)*L(i);endn=H/l;%計算編碼效率fprintf('編碼后所得碼字:n');disp(x)%顯示按概率降序排列的碼字fprintf('平均碼長:K=n');disp(l)%顯示平均碼長fprintf('信息嫡:H(X)=n');disp(H)%顯示信息嫡f