粗大误差C语言程序#(精选.)

误差分析的C语言实现学院：电子信息工程学院专业班级：通信1004学生姓名：童博学号：102840432012 年12 月26 日一、编程分析1、误差的基本概念所谓误差就是测量值与真实值之间的差，可以用下式表示误差=测得值-真值（1）绝对误差：某量值的测得值和真值之差为绝对误差，通常简称为误差。

绝对误差=测得值-真值（2）相对误差绝对误差与被测量的真值之比称为相对误差，因测得值与真值接近，故也可以近似用绝对误差与测得值之比值作为相对误差。

相对误差=绝对误差/真值≈绝对误差/测得值（3）引用误差所谓引用误差指的是一种简化和使用方便的仪器仪表表示值的相对误差，它以仪器仪表某一刻度点的示值误差为分子，以测量范围上限值或全量程为分母，所得的比值称为引用误差。

引用误差=示值误差/测量范围上限2、精度反映测量结果与真值接近程度的量，称为精度，它与误差大小相对应，因此可以用误差大小来表示精度的高低，误差小则精度高，误差大则精度低。

精度可分ⅰ准确度它反映测量结果中系统误差的影响程度ⅱ精密度它反映测量结果中随机误差的影响程度ⅲ精确度它反映测量结果中系统误差和随机误差综合的影响程度，其定量特征可以用测量的不确定度来表示。

3、有效数字与数据运算含有误差的任何近似数，如果其绝对误差界是最末位数的半个单位，那么从这个近似数左方起的第一个非零的数字，称为第一位有效数字。

从第一位有效数字起到最末一位数字止的所有数字，不论是零或非零的数字，都叫有效数字。

数字舍入规则如下：①若舍入部分的数值，大于保留部分的末位的半个单位，则末位加1。

②若舍去部分的数值，小于保留部分的末位的半个单位，则末位不变。

③若舍去部分的数值，等于保留部分的末位的半个单位，则末位凑成偶数。

即当末位为偶数时则末位不变，当末位为奇数时则末位加1。

（2）算术平均值对某一量进行一系列等精度测量，由于存在随机误差，其测得值皆不相同，应以全部测得值的算术平均值作为最后的测量结果。

1、算术平均值的意义在系列测量中，被测量所得的值的代数和除以n 而得的值成为算术平均值。

设 1l ，2l ，…,n l 为n 次测量所得的值，则算术平均值121...nin i l l l l x n n=++==∑ 算术平均值与真值最为接近，由概率论大数定律可知，若测量次数无限增加，则算术平均值x 必然趋近于真值0L 。

i v = i l -xi l ——第i 个测量值，i =1,2,...,;n i v ——i l 的残余误差（简称残差）2、算术平均值的计算校核算术平均值及其残余误差的计算是否正确，可用求得的残余误差代数和性质来校核。

残余误差代数和为：11n niii i v l nx ===-∑∑当x 为未经凑整的准确数时，则有1nii v==∑01）残余误差代数和应符合： 当1n ii l =∑=nx ，求得的x 为非凑整的准确数时，1nii v =∑为零；当1n ii l =∑>nx ，求得的x 为凑整的非准确数时，1nii v =∑为正；其大小为求x 时的余数。

当1nii l =∑<nx ，求得的x 为凑整的非准确数时，1nii v =∑为负；其大小为求x 时的亏数。

2）残余误差代数和绝对值应符合： 当n 为偶数时，1ni i v =∑≤2n A; 当n 为奇数时，1nii v =∑≤0.52n A ⎛⎫- ⎪⎝⎭ 式中A 为实际求得的算术平均值x 末位数的一个单位。

（3）测量的标准差测量的标准偏差称为标准差，也可以称之为均方根误差。

1、测量列中单次测量的标准差σ==式中 n —测量次数（应充分大）i δ——测得值与被测量值的真值之差σ=2、测量列算术平均值的标准差x σ=3、标准差的其他计算法1．别捷尔斯法nivσ=∑二、程序流程1、算术平均值2、求残余误差3、校核算术平均值及其残余误差4、判断系统误差5、求测量列单次测量的标准差6、判别粗大误差7、求算术平均值的标准差8、求算术平均值的极限误差9、写出最后测量结果三、编程实现源程序：#include "stdafx.h"#include "stdio.h"#include "stdlib.h"#include "math.h"double gebruce[]={1.15,1.46,1.67,1.82,1.94,2.03,2.11,2.18,2.23,2.28,2.33,2.37,2.41,2.44,2.48,2.50,2.53,2.56,2.58,2.60,2.62,2.64,2.66,2.74,2.81,2.87,2.96,3.17};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);//求平均数void canyu_error_fun(double datas[],double canyu_error_data[],double average,int datas_num);//求残余误差void revise_average_fun(double datas[],double canyu_error_data[],double average,int datas_num);void judge_system_error_fun(double canyu_error_data[],int datas_num);//判断系统误差//马利准则//阿卑准则void abei(double canyu_error_data[],double std_error_average,int datas_num);//求测量列单次测量的标准差double std_row_fun(double canyu_error_data[],int datas_num);//判别粗大误差void crassitude_error_fun(double datas[],double average,double std_besia,int datas_num);int liyiter(double canyu_error_data[],double std_row_besia,int datas_num);//求算数平均值的标准差double std_error_average_fun(double std_bersia,int datas_num);//求算数平均值的极限误差double limit_error_average_fun(double std_error_average,int datas_num);double myround(double val,int digits);int main(int argc, char* argv[]){//开辟数组的指针double *datas,*canyu_error_data;//保存的数据double average,canyu_error=0,std_besia=0,std_error_average=0,limit_error_average=0; //测试数据//double datas[]={24.774,24.778,24.771,24.780,24.772,24.777,24.773,24.775,24.774}; //数组长度以及一些常量int datas_num,j,flags;datas_num=9;printf("输入你要创建的数据容量：");scanf("%d",&datas_num);//创建数组，数据，残余误差datas=(double *)calloc(datas_num,sizeof(double));canyu_error_data=(double *)calloc(datas_num,sizeof(double));printf("请输入你要测试的数据：\n");//接收数据for(j=0;j<datas_num;j++){printf("%d:",j+1);scanf("%lf",datas+j);}// printf("接收数据成功!\n");average=myround(average_fun(datas,datas_num),3);canyu_error_fun(datas,canyu_error_data,average,datas_num);//平均值printf("****************************************************\n");printf("平均值为：%lf\n",average);printf("****************************************************\n");//残余误差printf("残余误差为：\n");for(j=0;j<datas_num;j++){printf("%d:",j+1);printf("%lf\n",canyu_error_data[j]);}printf("****************************************************\n");//校核算术平均值revise_average_fun(datas,canyu_error_data,average,datas_num);printf("****************************************************\n");//判别系统误差printf("马利可夫准则判定：\n");judge_system_error_fun(canyu_error_data,datas_num);printf("****************************************************\n");//测量列单次测量的标准差std_besia=std_row_fun(canyu_error_data,datas_num);printf("测量列单次测量的标准差为：%lf\n",std_besia);printf("****************************************************\n");//判别粗大误差crassitude_error_fun(datas,average,std_besia,datas_num);printf("****************************************************\n");flags=liyiter(canyu_error_data,std_besia,datas_num);printf("莱以特准则判断粗大误差：\n");if(flags==0){printf("不存在粗大误差\n");}else printf("存在粗大误差，应剔除\n");printf("****************************************************\n");//算数平均值的标准差std_error_average=std_error_average_fun(std_besia,datas_num);printf("算数平均值的标准差为：%lf\n",std_error_average);printf("****************************************************\n");//算数平均值的极限误差limit_error_average=limit_error_average_fun(std_error_average,datas_num);printf("算数平均值的极限误差为：%lf\n",limit_error_average);printf("****************************************************\n");//阿卑printf("阿卑准则判定:\n");abei(canyu_error_data,std_error_average,datas_num);printf("****************************************************\n");return 0;}//main//求平均数double average_fun(double datas[],int datas_num){int k=0;double average=0;for(k=0;k<datas_num;k++){average+=datas[k];}average/=datas_num;return average;}//求残余误差void canyu_error_fun(double datas[],double canyu_error_data[],double average,int datas_num) {int k=0;for(k=0;k<datas_num;k++){canyu_error_data[k]=datas[k]-average;}}//校核算术平均值void revise_average_fun(double datas[],double canyu_error_data[],double average,int datas_num){int k=0;double sum_canyu=0,sum_datas=0,sum_canyu_test,A;for(k=0;k<datas_num;k++){sum_canyu+=canyu_error_data[k];sum_datas+=datas[k];}sum_canyu_test=sum_datas-datas_num*average;printf("该算术平均值为:%lf\n",average);printf("该算术平均值末位数的一个单位是：");scanf("%lf",&A);if(datas_num%2==0){if(sum_canyu_test>datas_num/2.0*A){printf("该算术平均值不正确!\n");}else printf("该算术平均值正确!\n");}else{if(sum_canyu_test>(datas_num/2.0-0.5)*A){printf("该算术平均值不正确!\n");}else printf("该算术平均值正确!\n");}}//残余误差法判断系统误差void judge_system_error_fun(double canyu_error_data[],int datas_num){int k,j;double system_errors=0,system_error=0,data_error;j=(datas_num+1)/2;for(k=0;k<j;k++)system_errors+=canyu_error_data[k];for(k=0;k<datas_num;k++)system_error+=canyu_error_data[k];data_error=system_errors-system_error;printf("该系统误差为:%lf,",data_error);if(data_error<0.01)printf("系统误差较小，故不存在系统误差\n");else printf("存在系统误差\n");}//求测量列单次测量的标准差，方法二double std_row_fun(double canyu_error_data[],int datas_num){double std_row_besia=0,std_row_bejacks=0,u;int k=0;for(k=0;k<datas_num;k++){std_row_besia+=canyu_error_data[k]*canyu_error_data[k];std_row_bejacks+=double(fabs(canyu_error_data[k]));}std_row_besia=myround(sqrt(std_row_besia/((float)(datas_num-1))),4);std_row_bejacks=myround(1.253*std_row_bejacks/(sqrt(datas_num*(float)(datas_num-1))),4); u=myround(std_row_bejacks/std_row_besia-1,4);printf("贝赛尔-别捷尔斯判断系统误差:\n");if(double(fabs(u)<2/sqrt((double)datas_num-1))){printf("因为|u|=%lf<2/sqrt(%d-1)=%lf,",u,datas_num,2/sqrt((double)datas_num-1));printf("不存在系统误差\n");}else{printf("因为|u|=%lf>2/sqrt(%d-1)=%lf\n",u,datas_num,2/sqrt((double)datas_num-1));printf("存在系统误差\n");}return std_row_besia;}//判别粗大误差void crassitude_error_fun(double datas[],double average,double std_besia,int datas_num){double *datas_list,temp;double x1,xi,g0,g1,gi;//flag为α,flag=0时α=0.01，flag=1时α=0.05int k=0,j=0,flag=0;datas_list=(double *)calloc(datas_num,sizeof(double));for(k=0;k<datas_num;k++)datas_list[k]=datas[k];//冒泡排序法...由小到大排序for(k=0;k<datas_num;k++){for(j=0;j<datas_num-k;j++){if(datas_list[j]>datas_list[j+1]){temp=datas_list[j];datas_list[j]=datas_list[j+1];datas_list[j+1]=temp;}else continue;}}x1=datas_list[1];xi=datas_list[datas_num];x1=average-x1;xi=xi-average;gi=xi/std_besia;g0=gebruce[datas_num-3];g1=x1/std_besia;printf("使用格罗布斯判别准则判定粗大误差为：\n");printf("x1=%lf,\tx%d=%lf\n",datas_list[1],datas_num,datas_list[datas_num]);printf("x1=%lf,\tx%d=%lf\n",x1,datas_num,xi);printf("g1=%lf,\tg%d=%lf\tg0=%lf\n",gebruce[datas_num],datas_num,gi,g0);if(g1<gi&&gi<g0)printf("不存在粗大误差\n");else printf("存在粗大误差\n");free(datas_list);}//莱以特准则判别int liyiter(double canyu_error_data[],double std_row_besia,int datas_num) {int k=0,m=0;while(k<datas_num){if(fabs(canyu_error_data[k])>3*std_row_besia){m=1;break;}else k++;}return m;}//求算数平均值的标准差double std_error_average_fun(double std_bersia,int datas_num){double std_error_average;std_error_average=std_bersia/(sqrt((double)datas_num));return myround(std_error_average,4);}//求算数平均值的极限误差double limit_error_average_fun(double std_error_average,int datas_num) {double ta,limit_error_average;ta=T[datas_num-1];limit_error_average=ta*std_error_average;return myround(limit_error_average,4);}void abei(double canyu_error_data[],double std_error_average,int datas_num) {int k;double abei_sum=0;for(k=0;k<datas_num-1;k++){abei_sum+=canyu_error_data[k]*canyu_error_data[k+1];}abei_sum=fabs(abei_sum);if(abei_sum>sqrt((double)datas_num-1)*std_error_average*std_error_average) {printf("因为u=%lf>sqrt(%d-1)*δ*δ=%lf\n",abei_sum,datas_num, sqrt((double)datas_num-1)*std_error_average*std_error_average);printf("存在周期性误差\n");}else printf("不存在周期性误差\n");}double myround(double val,int digits){double d = pow(10.0, digits);return (int)( val * d + 0.5 ) / d;}四、程序使用说明根据界面提示输入所需测试数据个数再逐一输入测试数据即可判断出该组数据是否存在粗大误差五、程序实例1.输入数据个数为8个输入数据分别为：12.11312.11412.11712.11912.11512.11612.11812.301通过程序可以发现最后一位存在粗大误差2.输入数据个数为9个输入数据分别为：24.77424.77824.77124.78024.77224.77724.77324.77524.774通过程序验证该组数据无粗大误差最新文件仅供参考已改成word文本。