简单短序列的算术编码的matlab实现

《简单短序列的算术编码的matlab实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、简单短序列的算术编码的MATLAB实现正确实现的算术编码算法压缩能力Shannond定理描述的理论极限，是目前已知的压缩能力最强的无损压缩算法。   不过，由于算术编码算法的实现比较复杂，使用它作为默认压缩算法的应用程序还相当少。在Unix平台上非常流行的bzip2(这个工具有命令行模式的Windows版本)使用的就是经过修改的算术编码算法。   目前为止还没有使用算术编码作为默认压缩算法的Windows应用程序，WinRAR和WinIMP能够支持bzip2的解压。除此之外，在最新的JPEG标准中

2、也用到了经过修改的算术编码压缩算法，但JPEG所用的那种算法受专利保护，因此使用时必须获得授权。    在之后的文章会很好的研究这个算法的实现：现在给出一个简单的实例：%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%算术编码过程实例ssbm.m%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%about={...'本实例说明：''字符串不能太长，程序不加判断，请注意溢出；''本实例只限定少数字符串abcde；''实例只是说明一下算术编

3、码过程。'};disp(about);str=input('请输入编码的字符串(本程序仅仅是一个实例，请仅输入abcde）:');l=0;r=1;d=1;%初始间隔%程序限定字符为：a、b、c、d、ep=[0.20.30.10.150.25];%字符的概率分布，sum(p)=1n=length(str);disp('abcde')disp(num2str(p))fori=1:n  switchstr(i)    case'a'      m=1;    case'b'      m=2;    c

4、ase'c'      m=3;    case'd'      m=4;    case'e'      m=5;    otherwise      error('请不要输入其它字符！');    end    %判断字符    pl=0;pr=0;     forj=1:m-1        pl=pl+p(j);     end     forj=1:m        pr=pr+p(j);     end     %概率统计      l=l+d*pl;      r=l+d*(pr-p

5、l);      strl=strcat('输入第',int2str(i),'符号的间隔左右边界：');      disp(strl);      formatlong      disp(l);disp(r);      d=r-l;end运行过程如下：'本实例说明：''字符串不能太长，程序不加判断，请注意溢出；''本实例只限定少数字符串abcde；''实例只是说明一下算术编码过程。'请输入编码的字符串(本程序仅仅是一个实例，请仅输入abcde）:'aaabded'abcde0.20.30.1

6、0.150.25输入第1符号的间隔左右边界：00.200000000000000输入第2符号的间隔左右边界：00.040000000000000输入第3符号的间隔左右边界：00.008000000000000输入第4符号的间隔左右边界：0.0016000000000000.004000000000000输入第5符号的间隔左右边界：0.0030400000000000.003400000000000输入第6符号的间隔左右边界：0.0033100000000000.003400000000000输入第

7、7符号的间隔左右边界：0.0033640000000000.003377500000000 %I=imread('001.bmp')%imshow(I);clearI=[33113312;23313232;12333312];%I=[111100101110];[m,n]=size(I);%第一列为灰度值,第二列为个数,第三列为概率百分数，应该也可以用imhisttable=tabulate(I(); %注意的是，tabulate要求I的元素必须为非负整数%否则，以采用如下方法求解%如[123；1

8、22]，则统计出结果1是2个，2是3个，3是1个%sortM=sort(M();                                                             %uniqueM=([diff(sortM);1]>0);                                                  %count=[sortM(uniqueM)diff(find([1;uniqueM]))]   %即color,p如下所示c