C语言算术编解码(附注释)

#/6#include<iostream.h>#include"math.h"//定义所需要用到的变量及数组charS[100],A[10];floatP[10],f[10],gFs;voidbianma(inta,inth)//编码程序{inti,j;floatfr;floatps=1;floatFs=0;floatSp[100],b[100],F[100];//以待编码的个数和字符个数为循环周期，将待编码的字符所对应的概率存入到Fs中for(i=0;i<h;i++){for(j=0;j<a;j++){if(S[i]==A[j]){Sp[i]=P[j];fr=f[j]；//将划分好的01)区间的对应点赋值给fr}}Fs=Fs+ps*fr;//从选择的子区间中继续进行下一轮的分割。不断的进行这个过程，直到所有符号编码完毕。ps*=Sp[i];/^Ps}cout<<"Fs="<<Fs<<endl;//显示最终的算术编码gFs=Fs;floatl=log(1/ps)/log(2);//计算算术编码的码字长度lif(l>(int)l)l=(int)l+1;elsel=int(l);〃将Fs转换成二进制的形式intd[20];intm=0;while(l>m){Fs=2*Fs;if(Fs>1){Fs=Fs-1;d[m]=1;}elseif(Fs<1)d[m]=0;else{d[m]=1;break;}m++;}intz=m;//解决有关算术编码的进位问题,给二进制数加1if(m>=l){while(1){d[m-1]=(d[m-1]+1)%2;//最后位力口1if(d[m-1]==1)break;elsem--;}}cout<<"s=";for(inte=0;e<z;e++)cout<<d[e];cout<<endl;}voidjiema(inta,inth)//解码程序{inti,j;floatFt,Pt;floatFs=0,Ps=1;for(i=0;i<h;i++)//以编码个数和符号个数为循环周期，对其进行解码{for(intj=a-1;j>-1;j--){Ft=Fs;Pt=Ps;Ft+=Pt*f[j]；//对进行逆编码Pt*=P[j];if(gFs>=Ft)//对其进行判断，并且将值存入到数组A中{Fs=Ft;Ps=Pt;cout<<A[j];break;}}}cout<<endl;}voidmain()2{8a<<"输入所要编码的符号的个数，并按回车跳转："<<endl;inta,i,h=0;cin>>a;cout<<"请输入符号及其相对应的概率值，并按回车跳转：<<endl;for(i=0;i<a;i++){charx;floaty;cin>>x;A[i]=x；//将字符依次存入数组A中cin>>y;P[i]=y；//将字符所对应的概率依次存入到数组P中}for(i=1;i<a;i++){f[0]=0;f[i]=f[i-1]+P[i-1]；//将要编码的数据映射到一个位于[0,1)的实数区间中}cout<<"请输入所要编码的符号序列，并以*结尾："<<endl;while(1)//这个while语句的作用是将要编码的字符存入到数组S中{charss;cin>>ss;if(ss=='*')break;//在以"*"为结尾的时候结束存储S[h++]=ss;}8a<