C语言实现测量数据误差处理.docVIP

误差处理程序（C语言） 电子信息工程学院 通信1009子豪 实验目的 实现对输入数据的误差处理： 剔除粗大误差。 判断累进性系统误差和周期性系统误差。 计算平均值，方差，不确定度。 程序代码 #includestdio.h #includemath.h double SUM(double x[],int n); double AVRG(double x[],int n); double SD(double x[],int n); int PauTa(double x[],int n); int Chauvenet(double x[],int n); int Grubbs_1(double x[],int n); int Grubbs_2(double x[],int n); static int n; static double a[500]; int main() { int i,choose,leap=1; double avg,sd,v[500],M=0,AH=0,vmax=0; double P,PX[]={12.706,4.303,3.182,2.776,2.571,2.447,2.365,2.306,2.262,2.228,2.131,2.086,2.060,2.042,2.021,2.000,1.980,1.960}; printf(请输入数据总个数：\n); scanf(%d,n); printf(请输入数据：\n); for(i=0;in;i++)scanf(%lf,a[i]); avg=AVRG(a,n); sd=SD(a,n); printf(\n输入数据的平均值为%lf，标准差为%lf\n,avg,sd); while(leap) { printf(请选择粗大误差的检验法：\n1.莱特检验法\n2.肖维纳检验法\n3.格拉布斯检验法(置信概率99%%)\n4.格拉布斯检验法(置信概率95%%)\n5.停止检验\n); scanf(%d,choose); if(choose==1n10)printf(数据总量小于10，不能使用莱特检验法。\n); else if(choose==2n5)printf(数据总量小于5，不能使用肖维纳检验法。\n); else switch(choose) { case 1 : leap=PauTa(a,n);break; case 2 : leap=Chauvenet(a,n);break; case 3 : leap=Grubbs_1(a,n);break; case 4 : leap=Grubbs_2(a,n);break; case 5 : leap=0;break; default: printf(输入错误，请重新选择\n);break; } } avg=AVRG(a,n); sd=SD(a,n); for(i=0;in;i++) { v[i]=a[i]-avg; if(vmax*vmax=v[i]*v[i])vmax=v[i]; } printf(粗大误差剔除完毕，剩余数据如下(括号内为残差)：\n); for(i=0;in;i++)printf(%f(%f) ,a[i],v[i]); printf(平均值为%lf，标准差为%lf\n,avg,sd); for(i=0;in/2;i++)M+=v[i]; for(i=(n+1)/2;in;i++)M-=v[i]; M=sqrt(M*M); vmax=sqrt(vmax*vmax); if(M=vmax)printf(\n根据马利科夫判据，存在累进性系统误差。\n); else printf(\n根据马利科夫判据，不存在累进性系统误差。\n); for(i=0;in-1;i++)AH+=v[i]*v[i+1]; if((AH*AH)=((n-1)*sd*sd*sd*sd))printf(\n根据阿卑-赫梅特准则，存在周期性系统误差。\n); else printf(\n根据阿卑-赫梅特准则，不存在周期性系统误差。\n); if(n120)P=PX[17]; else if(n60)P=PX[16]; else if(n40)P=PX[15]; else if(n30)P=PX[14]; else if(n25)P=PX[13]; else if(n20)P=PX[12]; else if(n15)P=PX[11]; else if(n10)P=PX[10]; else P=PX[n-1]; printf(\n该组数组的置信区间为[%lf,%l