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.027600encode =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.690947decode =cbbabacbbc

結論：? 標簽為0.63215699的長度為10的序列為? a3a2a2a1a2a1a3a2a2a3

轉載于:https://www.cnblogs.com/rwg-xs/p/6033579.html

總結

以上是生活随笔為你收集整理的数据压缩 第四次作业的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。