电子测量大作业要点分析.doc

电子测量大作业 姓 名 杨嘉明 学 号 班 级 _ 自动化 1106 指导教师 王化深 题目一C语言实现测量数据误差处理 题2-21：参考例2-2-6的解题过程，用c语言或matlab设计测量数据误差处理的通用程序，要求如下： （1）提供测试数据输入、粗大误差判别准则选择等的人机界面； （2）编写程序使用说明； （3）通过实例来验证程序的正确性。 程序使用说明及实例验证 程序如下: #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[