mdmacode.m

自己用matlab写的算术编码的程序

%算术编码
%信源为二元信源,码长为20
clear;
close all;
s=input('s=');
P=[1/4 3/4];%概率分布p(0)=1/4,p(1)=3/4
F=[0 1/4];%累积分布函数F(0)=0，F(1)=1/4
Fs=0;As=1;%初始化
for k=1:1:20%信源序列长度
    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([50,49]);
c=num2bin(q,Fs);%将累积分布概率转化为二进制
c_B=c(2:L+1);%取小数点后长度为l的码字
%判断l位以后是否有尾数，若有尾数就进位到第l位
c_D=bin2dec(c_B);
c2=c(L+2:50);
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