C语言误差处理(平均数、残差).doc

实验二 不等进度测量及评定#include#include#include#includedouble d100;double x1=0.0,x2=0.0,x3=0.0,average=0.0;double T=63.66,9.92,5.84,4.60,4.03,3.71,3.50,3.36,3.25,3.17,3.11,3.05,3.01,2.98,2.95,2.92,2.90,2.88,2.86,2.85,2.83, 2.82,2.81,2.80,2.79,2.78,2.77,2.76,2.76,2.75,2.70,2.68,2.66,2.65,2.64,2.63,2.63,2.58;double average_fun(double datas,int datas_num)int k=0;for(k=0;kdatas_num;k+)average+=datask;average/=datas_num;return average;double surplus(int num,double a,double x)int i,j;double c100,w;for(i=0;inum;i+)ci=ai-x;di=ci;printf(求得的残余误差为:n );for(j=0;jnum;j+)printf(%f ,dj);for(i=0;inum/2.0*A)printf(该算术平均值不正确!n);else printf(该算术平均值正确!n);elseif(w(num/2.0-0.5)*A)printf(该算术平均值不正确!n);else printf(该算术平均值正确!n);judge1(int num)int t,j;double m=0.0,n=0.0,z=0.0;if(num%2=0)t=num/2;for(j=0;jt;j+)m+=dj;for(j=t;jnum;j+)n+=dj;z=m-n;elset=(num+1)/2;for(j=0;jt;j+)m+=dj;for(j=t;jnum;j+)n+=dj;z=m-n;printf(%lf,z);if(z=0.002)printf(无根据怀疑此组存在系统误差n);elseprintf(怀疑此组存在系统误差n);std(int num,double data)double m=0.0;int k;for(k=0;knum;k+)m+=datak*datak;x1=sqrt(m/(num-1);x2=x1/(sqrt(num);void judge2(int w)int i;for(i=0;iw;i+)if(di3*x1)printf(该数据存在粗大误差n);elseprintf(该数据不存在粗大误差n);limit(int num)double t=0.0;t=Tnum-1;x3=t*x2;bear(int datas_num,double datas)double y=0.0,k=0.0;int i;for(i=0;idatas_num;i+)k+=datasi;k/=datas_num;y=k+x3;printf(结果为:);printf(%lf,y);void main()int num;int i;double a100;double x=0.0,w=0.0;printf(输入实验组数：);scanf(%d,&num);printf(输入测得的数据：);for(i=0;inum;i+)scanf(%lf,&ai);for(i=0;inum;i+)printf(%lf,ai);printf(n);x=average_fun(a,num);printf(平均数为：n);printf(%lf,x);printf(残余误差为：n);w=surplus(