bianma1.m

来自「信息论与编码中」· M 代码 · 共 27 行

clc;clear;
a=[0.25,0.25,0.20,0.15,0.10,0.05];%初始概率，也可为输入
k=length(a);y=0;
s=zeros (k,1);b=zeros(k,1);p=zeros(8,1);
for m=1:k
   s(m)=y;
   y=y+a(m);                     %求和，
   b(m)=ceil(-log2(a(m)));       %取整，求位数
   z=zeros(b(m),1);
   x=s(m);
   p=b2d10 (x);                  %调用子函数
   
for r=1:b(m)
    z(r)=p(r);
   end
  
  disp('输出结果为:')            %结果显示
  disp('初始概率'),  disp(a(m))  
  disp('求和结果'),  disp(s(m))
  disp('编码位数'),  disp(b(m)) 
  disp('最终编码'),  disp(z') 

end