粗大误差处理程序.docVIP

粗大误差处理程序 #include stdio.h #include math.h #include stdlib.h int i,N; double a[20],v[20],c,a1,a2,a3,k,g,r,r1,r2; void laiyite(); void roman(); void geluobusi(); void dikesong(); double average(double arr[],int m); double bzhch(double arr1[],int m); void kzhi(int n1); void gzhi(int n2); void rzhi(int n3,double arr[],int n); void main() { int n,w; printf(\n*****************************粗大误差处理****************************\n\n); printf(请输入数据个数: ); scanf(%d,N); printf(\n请?输?入?数簓据Y: ); for(i=0;iN;i++) scanf(%lf,a[i]); S:{ printf(\n\n***************************菜单**************************\n); printf(\t\t\t1.莱以特准则\n\t\t\t2.罗曼洛夫斯基准则ò\n\t\t\t3.格罗布斯准则\n\t\t\t4.狄克松准则\n); printf(**********************************************************\n); printf(\n请选择判别的准则： [ ]\b\b); scanf(%d,n); switch(n) { case 1: laiyite();break; case 2: roman();break; case 3: geluobusi();break; case 4: dikesong();break; } printf(\n\n要继续吗\t1.继续 2.退出? [ ]\b\b); scanf(%d,w); switch(w) { case 1: goto S;break; case 2:break; } } } void laiyite() //1.莱以特准则 { double b; L1:{ c=average(a,N); b=bzhch(a,N); for(i=0;iN;i++) v[i]=a[i]-c; for(i=0;iN;i++) { v[i]=abs(v[i]); if(v[i]3*b) { printf(\n剔除值残余误差： %8.3lf 剔除值 %8.3lf\n,v[i],a[i]); for(;iN;i++) a[i]=a[i+1]; a[N-1]=\0; N--; goto L1; } } } printf(\n标准差： %8.3lf \n,b); printf(\n输出最后数值： ); for(i=0;iN;i++) printf(%8.3lf ,a[i]); } void roman() //2.罗曼洛夫斯基准则 { double t=100,x[20],b,ave,b2,sum=0,s2; int d; printf(\n请?输?入?显?著?度è：阰 ); scanf(%lf,a1); for(i=0;iN;i++) x[i]=a[i]; L2:{ kzhi(N); for(i=0;iN;i++) { if(x[i]t) { t=x[i]; d=i; } } for(i=d;iN;i++) x[i]=x[i+1]; x[N-1]=\0; ave=average(x,N-1); b2=bzhch(x,N-2); if(abs(t-ave)k*b2) { pri