香农的matlab编码

%clearall;clc;p=input('Pleaseinputp(x):');%只能对单符号二进制编码p=sort(p,'descend');n=length(p);Hx=0;fori=1:n%求解熵Hx=Hx-(p(i))*log2(p(i));endhp=0;xn=[0];fori=2:n%求解p(aj)xn(i)=hp+p(i-1);hp=xn(i);endzn=[];fori=1:n%求解位数Kzn(i)=ceil(-log2(p(i)));endaa=0;fori=1:n%求解均值码Kaa=aa+p(i)*zn(i);endnn=Hx/aa*100;fori=1:ncn=0;tn=zeros(1,zn(i));bn=0;%tnis单个编码forj=1:zn(i)bn=2^j;cn=cn+1/bn;%bn是不同阶的判断条件ifxn(i)cncn=cn-1/bn;continue;endifxn(i)=cntn(j)=1;endendtn%显示编码enddisp(['信息熵(H(X))=',num2str(Hx)]);disp(['平均码长K=',num2str(aa)]);disp(['编码效率n=',num2str(nn)]);