Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

jingang159 发表于 2006-5-16 11:30

诚心求教——算术编码解码

我调用下面3个子程序后没有结果，请大虾帮我看看下面3个子程序有没有问题，如果没问题，如何用一个主程序来调用它们来完成算术编码和解码。不是不会调用，是调用后，MATLAB马上死掉，动不了了。 程序如下： %****************************************************************************************** % 算术编码（Arithmetic Coding） %****************************************************************************************** function varargout=ArithCode(x,p) % 参考文献： % 数据压缩原理与应用（第二版）[美]David Salomon 著 pp:71-82 Len=length(x); CodeLen=0; % 初试码长设为0 Code=0; % 初试化编码结果 Low=0+(1-0)*p(1).low; High=0+(1-0)*p(1).high; % 对Low和High重定标 Low=floor(Low*10^4); High=floor(High*10^4); % 逐个进行编码 for i=2:Len for j=1:length(p) if x(i)==p(j).value Num=p(j).num; end end % 还原为小数 Low=Low/10^4; if High==9999 High=1; else High=High/10^4; end % LowT,HighT分别用语暂存Low和High LowT=Low; HighT=High; % Low和High更新 Low=LowT+(HighT-LowT)*p(Num).low; High=LowT+(HighT-LowT)*p(Num).high; Low=floor(Low*10^4); High=floor(High*10^4-1); % 一旦Low和High最左边的数字相等，则将最左边的数字写入输出流 if floor(Low/1000)==floor(High/1000) CodeLen=CodeLen+1;% 码长加1 Code=Code*10+floor(Low/1000); % Low的右端移入0，High的右端移入9 Low=(Low-floor(Low/1000)*1000)*10+0; High=(High-floor(High/1000)*1000)*10+9+1; end end CodeLen=CodeLen+4; % 加上最终四位码子作为总的码长 Code=Code*10^4+Low; % 编码最后4位 varargout{1}=CodeLen; if nargout>=2 varargout{2}=Code; end %****************************************************************************************** % 计算符号的概率，累计频率以及频率范围 %****************************************************************************************** function varargout=CumFreq(x) xt=x; % 暂存原始数据 maxx=99; % 定义最大值，避免重复扫描 count1=1;count2=1; Len=length(x); % 定义结构体数组用于存储序号，符号，该符号的概率以及累计频率的上下限 p=struct('num',[],'value',[],'probability',[],'low',[],'high',[]); for i=1:Len if xt(i)==maxx i=i+1; else for j=i+1:Len if xt(j)==x(i) count1=count1+1; % 频率加1 xt(j)=maxx; end end p(count2).num=count2; p(count2).value=x(i); p(count2).probability=count1/Len; count2=count2+1; count1=1; end end % 对应的频率范围 p(1).high=1.0 p(1).low=p(1).probability; for i=2:length(p) p(i).high=p(i-1).low; p(i).low=p(i-1).low-p(i).probability; end varargout{1}=p; %****************************************************************************************** % 算术解码（Arithmetic Decoding） %****************************************************************************************** function varargout=ArithDecode(p,CodeLen,Code) % 参考文献： % 数据压缩原理与应用（第二版）[美]David Salomon 著 pp:71-82 Len=0; Low=0000;High=9999; CodeLenT=CodeLen-4; CodeT=floor(Code/10^CodeLenT); % 码字初始化 %逐个进行解码 while CodeT~=Low index=((CodeT-Low+1)*10-1)/(High-Low+1); index=floor(index)/10; for i=1:length(p) if index>=p(i).low & index<p(i).high Len=Len+1; x(Len)=p(i).value; Low=Low+(High-Low+1)*p(i).low; High=Low+(High-Low+1)*p(i).high; end end % 一旦Low和High最左边的数字相等，则将最左边的数字写入输出流 if floor(Low/1000)==floor(High/1000) % Low的右端移入0，High的右端移入9 Low=(Low-floor(Low/1000)*1000)*10+0; High=(High-floor(High/1000)*1000)*10+9+1; CodeLen=CodeLen-1; Code=Code-floor(Low/1000)*10^CodeLen; CodeLenT=CodeLenT-1; CodeT=floor(Code/10^CodeLenT); end end varargout{1}=Len; if nargout>=2 varargout{2}=x; end
[此贴子已经被作者于2006-5-21 14:42:57编辑过]

hxjiangsha 发表于 2006-5-16 15:46

一种改进了的算术编码！ 不错， 很值得大家研究研究一下！
[此贴子已经被作者于2006-5-16 15:47:49编辑过]

jingang159 发表于 2006-5-17 22:47

happy教授最近很忙吗？还有哪些高手能帮我看看程序！！

happy 发表于 2006-5-20 14:02

回复：（jingang159）求教happy教授——算术编码解码问...

代码是否正确没功夫看 你可以找一个算术编码的代码比较验证一下 前面我好像给过一个，你搜索一下 用主程序调用是什么意思啊？调用函数你不会？晕 另外问问题的时候不要指名道姓谁来回答，这样即使有知道的人看到了也不一定回答你

jingang159 发表于 2006-5-21 14:51

不是不会调用，可能是我的子程序有问题吧，调用后MATLAB死掉了，动不了！！~~~

zhenyulu18 发表于 2009-1-20 14:57

楼主的程序中ArithCode(x,p),中x,p分别代表什么？

ch_j1985 发表于 2009-1-20 22:45

回复 7楼 zhenyulu18 的帖子

这个估计得自己看书了，LZ好久都没有来了！
《数据压缩原理与应用（第二版）》[美]David Salomon 著 pp:71-82

shazi 发表于 2009-5-8 17:06

请Happy教授给我一个算术编码的代码，好不好？

页: [1]

声振论坛's Archiver

诚心求教——算术编码解码

回复：（jingang159）求教happy教授——算术编码解码问...

回复 7楼 zhenyulu18 的帖子