信息论与编码实验程序与结果图(matlab)

t)=[char(ENDc1(y)),'0'];Рt=t+1;РENDc1=ENDc;РendРA%排序后的原概率序列РENDc%编码结果Рfor i=1:nР [a,b]=size(char(ENDc(i)));Р L(i)=b;РendРavlen=sum(L.*A)%平均码长Рselen=(L-avlen).^2%?Рmselen=sum((selen).*A)%码长均方差РH=-A*(log2(A'))%?РP=H/avlen%?Рsubplot(2,1,2)Рh=stem(1:n,selen);%Рaxis([0 n+1 0 max(selen)+0.1]);Рset(h,'MarkerFaceColor','blue','linewidth',2)Рxlabel('信源向下排');ylabel('方差值selen');Рhold onРplot(0:n+1,mselen*ones(1,n+2),'r','linewidth',2);Рhold offРlegend('每个码长与平均码长的方差','码长均方差');Р实验结果РA = 0.4000 0.2000 0.2000 0.1000 0.1000РENDc =[ 11, 1, 0, 101, 100]Рavlen = 1.8000Рselen = 0.0400 0.6400 0.6400 1.4400 1.4400Рmselen =0.5600РH =2.1219РP =1.1788РA =0.4000 0.2000 0.2000 0.1000 0.1000РENDc =[ 0, 10, 111, 1101, 1100]Рavlen =2.2000Рselen = 1.4400 0.0400 0.6400 3.2400 3.2400Рmselen =1.3600РH =2.1219РP =0.9645Р结果图