编程实现算术编码算法

六、参考书
1.《信息论——基础理论及应用》傅祖芸，电子工业出版社
七、源代码
#include<math.h>
#include<string.h>
#include<stdio.h>
const double proc[]={0.10,0.10,0.10,0.1,0.1,0.1,0.1,0.1,0.15,0.05};
编程实现算术编码算法
中国地质大学计算机学院信息安全专业
信息论实验报告
实验三算术编码
一、实验内容
编程实现算术编码算法
二、实验环境
1.计算机
2.Windows 2000或以上
3.DEVC++
三、实验目的
1.进一步熟悉算术编码算法；
2.掌握C语言编程（尤其是数值的进制转换，数值与字符串之间的转换等）
四、实验要求
wei+=proc[temp++];
temp--;
areaEnd=areaBegin+wei*len;//计算新的空间
areaBegin=areaBegin+(wei-proc[temp]*len);
printf("%d",temp);
}
printf("\n");
int main()
{
if(readdat())
1.提前预习实验，认真阅读实验原理。
2.认真高效的完成实验，实验过程中服从实验室管理人员以及实验指导老师的管理。
3.认真填写实验报告。
五、实验原理
算术编码是把一个信源表示为实轴上0和1之间的一个区间，信源集合中的每一个元素都用来缩短这个区间。
1.算法流程
（1）输入信源符号个数，信源概率分布，还有需要编码的符号序列，（2）根据概率可以算出初始编码间隔，
double temp1=result;
int temp2;
for(int j=0;j<cordLength;j++)//是进制转换成二进制
{
temp1*=2;
temp2=int(temp1);
temp1-=temp2;
cord[j]=temp2;
printf("%d",temp2);
}
printf("\n");
double result,areaBegin,areaEnd;
int cord[1000],cordLength;
char str[1000];
int strLength=0;
bool readdat()
{
printf("***********固定模式***********\n");
printf("请输入字符串(0-9):\n");
double w=0.0,len;
areaBegin=0.0,areaEnd=1.0;
for(int i=0;i<strLength;i++)
{
int n=str[i]-'0',k;
w=0.0;
for(k=0;k<n;k++)
w+=proc[k];
len=areaEnd-areaBegin;
areaEnd=areaBegin+len*(w+proc[k]);
printf("字符输入错误!!!\n"
decord();
}}
return 0;
High——当前编码的上限，
Low——当前编码的下限，
high——中间变量，用来计算下一个编码符号的当前间隔的上限，low——中间变量，用来计算下一个编码符号的当前间隔的下限，d——当前间隔之间的距离。
（3）扫描需编码的符号序列，确定编码空间
第1个编码符号的当前间隔为其初始的编码间隔，
第i个编码符号的当前间隔为第i-1个编码后的[Low，High），第i+1个编码符号的当前间隔算法如下：high=Low+d*第i+1个初始编码符号对应的上限，low=Low+d*第i+1个编码符号对应的下限，然后High=high，Low=low，d=d*第i个编码符号的概率。
scanf("%s",str);
while(str[strLength]!='\0')
strLength++;
for(int i=0;i<strLength;i++)
if(str[i]>'9'||str[i]<'0')
return 1;
return 0;
}
void encord()
{
printf("编码:");
areaBegin=0.0,areaEnd=1.0;
wei=0.0;
int temp;
double len;
for(int j=0;j<strLength;j++)
{
temp=0;wei=0.0;len=areaEnd-areaBegin;
while(result-areaBegin>wei*len)
}
void decord()
{
printf("译码:\n");
result=0.0;
double wei=0.5;
for(int i=0;i<cordLength;i++,wei*+0.5)//二进制转换成十进制
result+=wei*cord[i];
printf("译码选取的数:%f\n",result);
areaBegin+=len*w;
}
result=areaBegin*0.01+areaEnd*0.99;//选择适合的点
cordLength=(int(-log(areaEnd-areaBegin)/log(2)))+1;
printf("编码位数:%d\n",cordLength);
printf("编码结果:");