xinxilun.rar_信息论与编码_算术编码

%主程序，可以更改信源字母表、信源序列和概率值，实现任意元的编码 clear all close all format long e; symbol=['abcd']; ps=[0.5 0.25 0.125 0.125]; inseq=('baacdabaa'); %输出十进制小数编码结果 codeword=arithencode(symbol,ps,inseq) %十进制小数转换为二进制数 for i=1:length(inseq) temp(i)=find(inseq(i)==symbol); As(i)=ps(1,temp(i)); end pro=prod(As); innum=codeword; N=ceil(abs(log2(1/pro))) %二进制小数应取的长度N if (innum>1)|(N == 0) %判断输入的有效性 disp('error!'); return; end count=0; tempnum=innum; record=zeros(1,N); while(N) count=count+1; %长度小于N if(count>N) N=0; end tempnum=tempnum*2; %小数转换为二进制,乘2取整 if tempnum>1 record(count)=1; tempnum=tempnum-1; elseif(tempnum==1) record(count)=1; N=0; %stop loop else record(count)=0; end end %输出二进制编码结果 codebinary=record %输出解码结果 outseq=arithdecode(symbol,ps,codeword,length(inseq))