数据压缩 第四次作业

5.给定如表4-9所示的概率模型，求出序列a₁a₁a₃a₂a₃a₁的实值标签。

算术编码代码

clc,clear all;  
symbol=['abc'];          %输入一串序列 
pr=[0.2 0.3 0.5];        %各字符出现的概率  
temp=[0.0 0.2 0.5 1.0];  %各字符的累积概率
orignal=temp;  
in=input('input a string of abc:');  
n=length(in);  
%编码  
for i=1:n  
    width=temp(4)-temp(1);  
    w=temp(1);  
    switch in(i)  
        case 'a'  
            m=1;  
        case 'b'  
            m=2;  
        case 'c'  
            m=3;  
        otherwise  
            error('do not input other character');  
    end  
    temp(1)=w+orignal(m)*width;  
    temp(4)=w+orignal(m+1)*width;  
    left=temp(1);  
    right=temp(4);  
    fprintf('left=%.6f',left);  
    fprintf('    ');  
    fprintf('right=%.6f\n',right);  
end  
encode=(temp(1)+temp(4))/2 

结果：

input a string of abc:'aacbca'
left=0.000000    right=0.200000
left=0.000000    right=0.040000
left=0.020000    right=0.040000
left=0.024000    right=0.030000
left=0.027000    right=0.030000
left=0.027000    right=0.027600

encode =

    0.0273

结论：  a1a1a3a2a3a1的实值标签为 0.0273

 

6.对于表4-9给出的概率模型，对于一个标签为0.63215699的长度为10的序列进行解码。

算术解码代码：

encode=0.63215699;     %序列标签
n=10;    %序列长度
pr=[0.2 0.3 0.5];        %各字符出现的概率  
temp=[0.0 0.2 0.5 1.0];  %各字符的累积概率
orignal=temp; 
decode=['0'];  
for i=1:n  
    fprintf('tmp=%.6f\n',encode);  
    if(encode>=orignal(1)& encode<orignal(2))  
        decode(i)='a';  
        t=1;  
    elseif(encode>=orignal(2)& encode<orignal(3))  
        decode(i)='b';  
        t=2;  
    else  
        decode(i)='c';  
        t=3;  
    end  
    encode=(encode-orignal(t));  
    encode=encode/pr(t);  
end  
decode  

结果：

tmp=0.632157
tmp=0.264314
tmp=0.214380
tmp=0.047933
tmp=0.239666
tmp=0.132219
tmp=0.661093
tmp=0.322185
tmp=0.407284
tmp=0.690947

decode =

cbbabacbbc

结论：  标签为0.63215699的长度为10的序列为  a3a2a2a1a2a1a3a2a2a3