c语言随机误差分析.docVIP

误差理论与数据处理 姓名： 班级： 学号：随机误差分析 随机误差： 在实际相同条件下，多次测量同一量值时，其绝对值和符号无法预计的测量误差。 #include stdio.h #include stdlib.h #include math.h double d[100]; double x1=0.0,x2=0.0,x3=0.0,average=0.0; // t分布表 α=0.01 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+=datas[k]; } average/=datas_num; return average; } //求残差 double surplus(int num,double a[],double x) { int i,j;double c[100],w; for(i=0;inum;i++) { c[i]=a[i]-x; d[i]=c[i]; } printf(求得的残余误差为:\n ); for(j=0;jnum;j++) printf(%f ,d[j]); for(i=0;inum;i++) { w+=d[i]; } return(w); } //校核算术平均值 void check(int num,double w) {double A; printf(该算术平均值末位数的一个单位是：); scanf(%lf,A); if(num%2==0) { if(wnum/2.0*A) { printf(该算术平均值不正确!\n); } else printf(该算术平均值正确!\n); } else { if(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+=d[j]; for(j=t;jnum;j++) n+=d[j]; z=m-n; } else { t=(num+1)/2; for(j=0;jt;j++) m+=d[j]; for(j=t;jnum;j++) n+=d[j]; z=m-n;printf(%lf,z); } if(z=0.002) printf(无系统误差\n); else printf(有系统误差\n); } //求标准差 std(int num,double data[]) {double m=0.0;int k; for(k=0;knum;k++) { m+=data[k]*data[k]; } x1=sqrt(m/(num-1)); x2=x1/(sqrt(num)); } //判别粗大误差 void judge2(int w) {int i; for(i=0;iw;i++) { if(d[i]3*x1) printf(该数据有粗大误差\n); else printf(该数据无粗大误差\n); } } //求极限误差 limit(int num) {double t=0.0; t=T[num-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+=datas[i]; } k/=datas_num; y=k+x