mathencode.m

霍夫曼编码、算术编码和行程编码的matlab实现

%算术编码
%信源为二元信源,信源序列、码长均可自定义
clear;
close all;
n=input('信源序列长度n=');
s=input('信源符号序列s=');
P=[input('p(0)=') input('p(1)=')]; 
F=[0 P(1,1)];%累积分布函数F(0)=0，F(1)=P(0)
Fs=0;As=1;%初始化
for k=1:1:n%信源序列长度
    Fs=Fs+As*F(1,s(1,k)+1);%信源序列累积分布函数
    As=As*P(1,s(1,k)+1);%信源序列对应的区间宽度
end
L=ceil(abs(log2(1/As)));%编码后码长
q=quantizer([3*n,3*n-1]);
c=num2bin(q,Fs);%将累积分布概率转化为二进制
c_B=c(2:L+1);%取小数点后长度为L的码字
%判断L位以后是否有尾数，若有尾数就进位到第L位
c_D=bin2dec(c_B);
c2=c(L+2:3*n);
c2_D=bin2dec(c2);
if c2_D~=0
    c_D=c_D+1;
    mc_B=dec2bin(c_D,L);
else
    mc_B=c_B;
end
disp(mc_B)