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

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

ID:14270974

大小:41.50 KB

页数:9页

时间:2018-07-27

简单短序列的算术编码的matlab实现_第1页
简单短序列的算术编码的matlab实现_第2页
简单短序列的算术编码的matlab实现_第3页
简单短序列的算术编码的matlab实现_第4页
简单短序列的算术编码的matlab实现_第5页
资源描述:

《简单短序列的算术编码的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

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。