误差试验报告_第1页
误差试验报告_第2页
误差试验报告_第3页
误差试验报告_第4页
误差试验报告_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、误差实验报告实验一 误差的基本概念一、实验目的通过实验熟悉MATLAB 的基本操作,了解误差的定义及表示法、熟悉误差的来源、误差分类以及有效数字与数据运算。二、实验原理1、误差的基本概念:所谓误差就是测量值与真实值之间的差,可以用下式表示误差=测得值-真值绝对误差:某量值的测得值和真值之差为绝对误差,通常简称为误差。绝对误差=测得值-真值相对误差:绝对误差与被测量的真值之比称为相对误差,因测得值与真值接近,故也可以近似用绝对误差与测得值之比值作为相对误差。相对误差=绝对误差/真值绝对误差/测得值2、精度反映测量结果与真值接近程度的量,称为精度,它与误差大小相对应,因此可以用误差大小来表示精度的

2、高低,误差小则精度高,误差大则精度低。3、有效数字与数据运算含有误差的任何近似数,如果其绝对误差界是最末位数的半个单位,那么从这个近似数左方起的第一个非零的数字,称为第一位有效数字。从第一位有效数字起到最末一位数字止的所有数字,不论是零或非零的数字,都叫有效数字。数字舍入规则如下:若舍入部分的数值,大于保留部分的末位的半个单位,则末位加1。若舍去部分的数值,小于保留部分的末位的半个单位,则末位加1。若舍去部分的数值,等于保留部分的末位的半个单位,则末位凑成偶数。即当末位为偶数时则末位不变,当末位为奇数时则末位加1。三、实验内容1、用自己熟悉的语言编程实现对绝对误差和相对误差的求解。实验程序:1

3、 -加=轧 371;2 -xl=6. 376:3 -fprintfC绝对误差为;)4 -d=s-x5 -fprintf/相对误差为;')6 -l=(d/x)实验结果:绝对误差为二 d -0i0050相对误差为:1 ;7,3481e-0042、按照数字舍入规则,用自己熟悉的语言编程实现对下面数据保留四位有效数字进行凑整。(保留四位有效数字可使用matlab控制运算精度函数 vpa)原有数据3.141592.717294.510503.215516.378501舍入后数据实验程序:1 -a=3. 14159, 2. 71729, 4,51050,3,21551, 6. 376S012 b=

4、vpa(a, 4)实验结果:3,142, 1.717* 4.E1L 3.216, 6,379实验二误差的基本性质与处理一、实验目的了解误差的基本性质以及处理方法二、实验原理(1)算术平均值对某一量进行一系列等精度测量,由于存在随机误差, 其测得值皆不相同,应以全部测得值的算术平均值作为最后的测量结果。1、算术平均值的意义:在系列测量中,被测量所得的值的代数和除以n而得的值成为算术平均值。nJi-ll 12 .ln i=11c 1x设11, 12,,1n为n次测量所得的值,则算术平均值nn算术平均值与真值最为接近,由概率论大数定律可知, 若测量次数无限增加, 则算术平均值x必然趋近于真值 L0。

5、Vi =1i x-1i 第i个测量值,i = 12,n;vi h的残余误差(简称残差)2、算术平均值的计算校核算术平均值及其残余误差的计算是否正确,可用求得的残余误差代数和性质来校核。nn_ vi ='、' 1i -nx残余误差代数和为:T Tn-V vi V当X为未经凑整的准确数时,则有: im 01)残余误差代数和应符合:nn-1i -' vi当if =nx ,求得的x为非凑整的准确数时,if 为零;nn% 1 i -、vi>nx ,求得的x为凑整的非准确数时,为正;其大小为求x时的余数。nnli -v Vi当i <nx ,求得的x为凑整的非准确数时,i

6、 为负;其大小为求 x时的亏数。2)残余误差代数和绝对值应符合:n£ M n当n为偶数时,y < 2 a;Z Vi'n -0.5 IA当n为奇数时,1g I2 J式中A为实际求得的算术平均值 x末位数的一个单位。(2)测量的标准差测量的标准偏差称为标准差,也可以称之为均方根误差。1、测量列中单次测量的标准差ni 1 n式中 n测量次数(应充分大)a 一测得值与被测量值的真值之差a-x2、测量列算术平均值的标准差:三、实验内容:对某一轴径等精度测量 8次,得到下表数据,求测量结果。序号li /mmvi /mm2 /2vi / mm124.674224.675324.673

7、424.676524.671624.678724.672824.674假定该测量列不存在固定的系统误差,则可按下列步骤求测量结果。1、算术平均值2、求残余误差3、校核算术平均值及其残余误差4、判断系统误差5、求测量列单次测量的标准差6、判别粗大误差7、求算术平均值的标准差8、求算术平均值的极限误差9、写出最后测量结果实验程序:1 - fprintfT 初始数据:)2 -1= 24. 674, 24. 675, 24, 673, 24. 676, 24, 671* 2d. 67 24. 672, 24, 6743 -珏nntfC L算数平均值4 -xl=itean(l)5 -"门足近&

8、quot;残余误差16 -v=l-xI7 - fprintfC 3.校核算术平均值及其残余误差,)8 a=sun(v);9 ah=abs (a)10 bh= (8/2)*0.0001It-if ah<=bh)12 -fprirrtfT核核结果;校核正确n.rf )13 -else14 -fprirrtfC校核结果;校核错误n措)15 end16 -fprrntf(J4.判断系统误差)17 xt=Eun(v( 1: 4) )-swi (v(5:8)13 -xf (xt<=0. 02)19 -fprintf。结果:不存在系统误差, rf n')20 -else21 -fpEin

9、tf ('结果:存在系统误差n n 322 end23-fprintf(J5.刎里列单次刎里的标准差)24-bz=sqrt (sum(v.2)/7)25 - fprintf (J 6<判别系统是否存在粗大误差rf )25p=sort (1):27-gO=2- 03;28-g 1 = (x 1 -p (1) )/bz29-g8=(p(8)-xl)/bz30-if .尊数平均值 -24. 67412.残余误差gO>=gl4&gO>-g8)31 -fprintf (结果;系统不存在粗大误差,rin?)32-else33 -fprmtf。结果:系统存在粗大误差n n

10、)34-end35 - fprintfC 7,则里列单次测量的标准差)36-sc=b2/ (s4rt (8)37 - fprmtf C 3.算术平均值的极B艮误差)33-1=2. 36: %t (7, 0. 05)39-j=t *sc40-fprintf C 9.最后的则里结果一)41-ll=X14-jx42-12=zl-js实验结果:初始数据:1 -24.674024.675024,673024, 676024.671024.678024,672024.6740-0. 00010.0009-0.00110. 0019-0.00310. 0039-0. 0021-0. 00013 .校核算术平

11、均值及其残余误差 ah =3.5527e-015bh =4.0000e-004校核结果:校核正确4 .判断系统误差xt =0. 0030结果:不存在系统误差5 .别里列单次测里的标准差 bz =0. 00226 .判别系统是否存在粗大误差g0 =2. 0300gl =1.4000g8 =1.7361结果:系统不存在粗大误差7 .测里列单次则里的标准差 sc =7.8?16e-0048 .算术平均值的极限误差 jx =0.00199 .最后的测里结果:11 =24.676012 =24.6723实验三线性参数的最小二乘法处理一、实验目的最小二乘法原理是一种在多学科领域中获得广泛应用的数据处理方法

12、。通过实验要求掌握最小二乘法基本原理、正规方程以及组合测量的最小二乘法处理办法。二、实验原理(1)测量结果的最可信赖值应在残余误差平方和为最小的条件下求出,这就是最小二 乘法原理。即$+/+.+/=最小(2)正规方程最小二乘法可以将误差方程转化为有确定解的代数方程组(其方程式的数目正好等于未知数的个数),从而可求解出这些未知参数。这个有确定解的代数方程组称为最小二乘法估 计的正规方程。(3)精度估计为了确定最小二乘估计量 x1,x2,-,xt的精度,首先需要给出直接测量所得测量数据的精 度。测量数据的精度也以标准差 仃来表示。因为无法求得 仃的真值,只能依据有限次的测 量结果给出。的估计值,所

13、谓精度估计,实际上是求出估计值。(4)组合测量是通过直接测量待测参数的各种组合量,然后对这些测量数据进行处理,从而求得待测参数的估计量,并给出其精度估计。三、实验内容如下图所示已知直接测量刻线的各种组合量,要求检定刻线A、B、C、D间距离x1x2、 x3 ,测量数据的标准差以及估计量的标准差。(1)x1x2x3l1 =2.018mm l2=1.986mml3=2.020mml4 = 4.020mm l5=3.984mml6=6.030mm实验程序:1 -11-2*018;2 -12=1.986;3 -13=2,020;4 -14=4.020;5 -15=3.984:6 -16二6. 030 :7 -1=2. 018:1. 936 :2. 020; 4. 020 : 3. 984 : 6. 030:8 -A=l 0 0:0 1 0;0 0 1; 1 1 0;0 1 1;1 1 1;9 -B=AJ ;10 - invC=inv(AJ *A):11 - *1;12 -zl=Z(l, 1):13 -i2=X(2, I);14 -x3=X(3, 1);15 -L=xl;x2;x3;xl+i2;x2-hc3;xl-hc2+x3;16 -V=LL :17 -fprintf ('则里结果的标准差:PIS -bzc=sqrt 1 (sldh(V.

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论