不确定度与数据处理

1、不确定度与数据处理、误差与不确定度x=x-A1 误差与不确定度的关系(1) 误差：测量结果与客观真值之差公认值一如物理常数等其中A称为真值，一般不可能准确知道，常用约定真值代替:标准值一更高精度仪器测量结果理论值一理论公式计算结果对一个测量过程，真值 A的最佳估计值是平均值 X。在上述误差公式中，由于 A不可知，显然 x也不可知，对误差的最佳估计值是不确定度u( x)。(2) 不确定度：对误差情况的定量估计，反映对被测量值不能肯定的程度。通常所说“误差”一般均为“不确定度”含义。不确定度分为 A B两个分量，其中A类分量是可用统计方法估计的分量，它的主要成分是随机误差。2.随机误差：多数随机误

2、差服从正态分布。定量描述随机误差的物理量叫标准差。(1) 标准差与标准偏差.(Xi A)2limkkT真值A不可知，且测量次数 k为有限次 二(Xi X)2k 1(置信概率标准差实际上也不可知，于是:用标准偏差S代替标准差结果表述:Xi/真值的估计值其中S(x) JS：x)单次测量的标准偏差68.3%)单次测量标准差最佳估计值S(x)的物理意义：在有限次测量中，每个测量值平均所具有的标准偏差。 通常不严格区分标准差与标准偏差，统称为标准差。(2 )平均值的标准差真值的最佳估计值是平均值，故结果应表述为：(置信概率 68.3%)(并不是只做一次测量)X / 真值的最佳估计值S( X)平均值的标准

3、差最佳估计值S(X)(Xi X)2k(k 1)平均值的标准偏差例1:的1 ，即某观察量的n次独立测量的结果是 X, X2, , Xno试用方差合成公式证明平均值的标准偏差是样本标准偏差S(X)S(X)、 n解:0由题知X相互独立，则根据方差合成公式有Xn利用样本标准偏差的定义，可知u(X)=S(X>v'S2(X)u(X) S(X)S2(X)i=1,2,nS2(X)u(X)S(X)nU2(XJU2(Xn)n3.系统误差与仪器误差(限)(1 )系统误差：在同一被测量的多次测量过程中，保持恒定或以可以预知方式变化的那一部分误差称为系统误差。 已被确切掌握了其大小和符号的系统误差，称为可

4、定系统误差；对大小和符号不能确切掌握的系统误差称为未定系统误 差。前者一般可以在测量过程中采取措施予以消除或在测量结果中进行修正；而后者一般难以作出修正，只能估计出它 的取值围。在物理实验中，对未定系统误差的估计常常利用仪器误差限来进行简化处理。(2 )仪器误差(限)：由国家技术标准或检定规程规定的计量器具的允许误差或允许基本误差，经过适当简化称为 仪器误差限，用以代表常规使用中仪器示值和(作用在仪器上的)被测真值之间可能产生的最大误差。长度测量仪器:游标卡尺的仪器误差限按其分度值估计；钢板尺、螺旋测微计的仪器误差限按其最小分度的1/2计算。指针式仪表：仪=a% Nm式中N.是电表的量程，a是

5、准确度等级。数字仪表：仪=a%N<+b%Nn或仪=a%N+n字式中a是数字式电表的准确度等级，NX是显示的读数，b是误差的绝对项系数，2是仪表的满度值，n代表仪器固定项误差，相当于最小量化单位的倍数。 电阻箱：仪=ai%Ri Roi式中F0是残余电阻，R是第i个度盘的示值，a是相应电阻度盘的准确度级别。 直流电位差计：仪=a% (Ux10式中a是电位差计的准确度级别，U是标度盘示值，U0是有效量程的基准值，规定为该量程中最大的10的整数幕。直流电桥：仪=a%(RX)10式中FX是电桥标度盘示值，a是电桥的准确度级别， R)是有效量程的基准值，意义同上。(3) B类不确定度的处理仪、灵敏度

6、误差灵和估计误差限估。其中灵敏在物理实验中，B类不确定度的来源通常包括以下三种：仪器误差度误差可表示为0.20.2n/ xB类不确定度与各种误差限之间的关系为Ub4.不确定度的合成(1)直接测量x :Ua(x) , Ub(x)u(x)间接测量、u；(x) u(x)y=f(X1, X2, xn)(称为合成不确定度)其中X1, X2, Xn为相互独立的直接测量量J / f 、2 2 ()Ui Xi最终结果表述形式：u(y)(Xi)(3)结果有效数字的确定原则：u(y) / ln f、2 2/ . ()u (Xi)y , i Xi (单位)用分光计测棱镜材料的折射率公式为=5058'解:dn

7、u(n)nN u(N=不确定度u(N只保留一位有效数字； 测量结果N与不确定度u(N小数位数对齐。Asin =nA-。已测得 A=60Asin2u( n)= 1.64790.0007。n60 050 58sin3'，则黄光所对应的折射率Asinn匚.Asin2A ,1, Acos (亍d AAsin 22d2.60 0sin2、 A 1 )cos dA 2 2Asin -21.64792(ctgA0'2'，黄光(汞灯光源)所对应的AInn lnsin 2A Insin2ctg -)dA ctg 2 20.000426u(n) n1.6479 0.0004260.0007

8、n u(n)=1.64790.0007n5有效数字及其运算法则(1) 有效数字：由若干位可靠数字加一位可疑数字构成。在不计算不确定度的情况下，结果的有效数字由运算法则决定。(2 )运算法则 加减法：以参加运算各量中有效数字最末一位位数最高的为准并与之取齐。N=A+BGD 则 u(N) u2(A) u2(B) u2(C) u2(D)取决于u(A)、u(B)、u(C)、u(D)中位数最高者，最后结果与之对齐。 乘除法：以参加运算各量中有效数字最少的为准，结果的有效数字个数与该量相同。N AB，则 u(N) u(A)CDA2u(D)D2 2u(B)u(C)BC取决于其中相对不确定度最大者，即有效数字

9、个数最少者。 混合四则运算按以上原则按部就班执行。例3：某物理量的计算公式为k1 1.6d/H，其中k为常数,1.6为准确数,H 16cm, d=0.1500cm。若使 Y 的表示式中分母的值具有 4位有效数字，正确测 H的方法是(d )。(a)用游标卡尺估读到 cm千分位(c)用米尺只读到mm位(d)解:型6 O'1500 0.015H16(b)用米尺估读到cm百分位用米尺只读到cm位分母11.015为4位有效数字H即H只需2位有效数字即可，故应选 (d) 特殊函数的有效数字：根据不确定度决定有效数字的原则，从不丢失有效位数的前提出发，通过微分关系传播 处理。例4：tg452'

10、; =1.00116423最多可取几位有效数字？解：令 y=tg x，其中 x=452'1取 x 10.00029(rad)60180则y1 12 x 2,-0.00029 0.00058即小数点后第四位产生误差cos2 xcos2 45 2tg452' =1.0012 ，有五位有效数字。例 5：双棱镜测波长的计算公式为汗，对实验数据进行处理的计算结果如下表所示。x=0.28144mmb=5.9325mmb' =0.7855mmS=27.65cmS' =75.90cm-4u(x)=2.01010mm(b)/b =0.025(b') /b' =0.

11、025(S =0.5cm(S')=0.5cm(b)=0.005m m(b' )=0.005mm(S)=0.05cm(S') =0.05cm注：下标1代表来自方法误差，下标 2代表来自仪器误差。要求：(1)给出测量结果的正确表述(包括必要的计算公式)。(2 )定量讨论各不确定度的分量中，哪些是主要的，哪些是次要的，哪些是可以忽略的？如果略去次要因素 和可以忽略项的贡献，不确定度的计算将怎样简化？结果如何？x . bb0.28144.5.93250.78554解：(1)5.8671610mmS S(276.5 759.0)11dd(x)dbdbdSdSInIn xIn bI

12、n bIn (SS)22x2b2bS S S Sr22 222u()u(x)u(b)u(b)u(S)u(S)0.0111Yx2b2bS SSS其中u( x)2.01010 40.000714 ;x0.281442b2b2、3b2 3U2(b)2(b)/ 30.0050 00022b2bU1 (b )1(b)/ 311(b)0.025Ui (b)12bU2(b )1(b)/ . 31(b)2 . 32b2(b)/ .3-2b2、3 b0.0052 0.0250.00722u(b)莎Ui(b)2b0.00722u(b)2b5(b )2b2U2(b)2b(b) 22bU1(S)1(S)/ 30.50

13、.002792u(S)2U1(S)U2(S)S SS S(27.6575.90) .3U2(S)2(S)/ 30.050.000279S SS SS SS SS S(27.6575.90). 3U1(S)1(S)/ 30.50.00279_ 2u(S)_ 25(S)S SS S(27.6575.90) . 3U2(S)U2(S)2(S)/ 30.050.000279S SS SS SS SS S(27.6575.90) . 30.001842b22于是得0.0111=6.53u()5.86716 10 4u()=巴-610 mm)=587 7nm则有由前面的计算可知，不确定度主要来自业和2b讐

14、，次要因素是讐U1(S)和U1(S)，可以忽略的因S S S SU2(b)2bu2 (S)和 U2 (S )若只考虑主要项的贡献:u()2U1(b)2U1(b )1(b)2bu()=6 nmu()=587 6nm0.01022b,6b比严格计算的结果稍小但相差无几。二、数据处理方法1. 列表法：按一定规律把数据列成表格。 列表原则：(1 )表格的标题栏中 注明物理量的名称、符号和单位 ；(2) 记录原始数据(如记录刻度数，而不是记录长度)；(3) 简单处理结果(如算出长度)或函数关系；(4) 参数和说明(如表格名称、仪器规格、环境参数、常量以及公用单位等)。例6:拉伸法测弹性模量的载荷一一伸长

15、曲线如图所示，图上至少有 不规。它们是 坐标轴应标注物理量和单位，轴上缺少分度值标记标出，曲线应光滑连接，计算点坐标标注不规。3 最小二乘法与一元线性回归法5处绘制错误或 实验点应以醒目各测量值与这条曲线上对应点之差的平方和应例7:试用最小二乘原理推导直线方程y=kx中回归系数k的计算公式。解：根据最小二乘原理应有n(yik i 1kxi )20n(yii 1n2(yii 1kxi )2minkxi)(为)0nXiyii 1nkxi20i 1于是得xi yiXixyx2(2) 元线性回归法:由最小二乘原理,应有设直线方程kyi (a1y=a+bx，其中自变量bxi )2 minx的误差可略ky

16、i即 a i 1k匚yib i 1(a(a解之得(3 )相关系数r物理意义bXi)2bXi)2Xik2yii 1k2yii 1Xiyi22(Xi) k XiXi yiXi yi yx2a 2 2(X )2k Xi2yi(a bxj( 1)(a bXj)( k )Xy xyX2x2bXr:用于检验x和y之间是否存在线性关系。y a bx通过全部实验点 之间线性相关强烈yiyi随xi增加而增加yi随xi增加而减小xi、yi之间无线性关系拟合直线为与akXi1xykyii 12Xikxiyii 1xyx2 X2)(y2y2)x轴平行的直线例8:根据所给相关系数 r=-1r作出实验点分布草图: r =

17、0.9993 r=0.015(1)最小二乘法：对等精密度测量若存在一条最佳的拟合曲线，那么 取极小值。(4) 回归法使用要点： 自变量x测量误差可略，即应选择测量精度较高的物理量作自变量； 因变量y为等精度测量 或近似等精度测量，即 u(y)近似相等； 作线性关系的检验：利用物理规律或作图等其它方法确认线性关系的存在；或检验相关系数是否满足|r| 1。解：根据线路图可得R ()20.050.0100.0200.0300.0400.0V (V厂2.802.722.602.382.202.04例9：实验线路及测量数据如下，用一元线性回归法计算电压表阻RV (写出计算公式即可)EVRVR RvERv

18、Vu(R).0.10.1计算R1丄精度:R20.0，400.0Vu(1/V)u(V)0.011/VV2.80故将公式:变形为1 R1令V ERvE则有a 1b1aEERvRvR Rv4.逐差法(1)测量次数为偶数的逐差法0.005,0.00025可知R的精度较高0010.0036,0.0052.041y R x并设直线方程y=a+bxV-Rvab设自变量和因变量之间存在线性关系y a bx ,并有组实验数据:yn 1y1隔n项逐差，可得到b1Xn 1X1,bny2n yn取平均值X2nXnX1, xn , xn 1, x2n .； y1, ,yk, yn 1,y2nb 1 n b 1 n yn

19、 i yi b bi 一ni 1 n i 1Xn i Xi对于自变量x等间隔分布的情况，有 xnXinX求得bl后，可由公式yia b求出n nX i 1yi bXi(yn i yi)例10：已知R=R)(1+t)，实验数据如下，用逐差法求电阻温度系数(不要求计算不确定度)t (C)85.080.075.070.065.060.055.050.0R ()0.36220.35650.34990.34370.33800.33240.32700.3215解：R=R)+ Rot 并设 y=a+bx则有a= Rob= Ro=a即而利用逐差法可得:ai1234平均t =t i +4- t i(C)20.0

20、20.020.020.020.0R= Ri +4- R()0.02420.02410.02290.02220.0233 5于是有0.0233520.00.0011675R b t)b 0.00116750.26264.45 10 3 (1/C)条纹吞吐n0100200300400M镜位置X (mm)34.4830534.5158534.5483034.5806034.61300例11:迈克尔逊干涉仪实验数据处理条纹吞吐n500600700800900M镜位置X (mm)34.6446534.6763534.7080034.7394534.77085解法一：由逐差法可得i12345平均N=ni+

21、旷 ni500500500500500500d5oo=X +5-X (mm)P 0.161600.160500.159700.158850.15785 :0.15970u()222d500N2 08(nm)U(d 500 )d500u(N)N其中I 22U(d500)Ua (d 500 ) Ub(d500)U a (d 500 )U b (d 500 ) 2(d id 500 )540.000648(mm)于是u(N) Ub(N)0.530.2890.000050.0000289( mm)u()2U(d500 )d5002u(N)N638.820.000648 0.1

22、5970220.28950022.6( nm)u()=6393( nm)i12345平均N'=( ni +5- n )/5100100100100100100d1oo=( Xi +5- Xi) /5 (mm):0.0323200.0321000.0319400.0317700.0315700.031940解法二：由逐差法可得u()22d100N2 0.031940100638.8( nm)u(d100)d100u(N )Nu(d100 ).Ua (d 100 ) Ub (d 100 )其中U a (d100 )U b ( d 100 )(didp。)25 40.000130(mm)于是

23、u()u(N ) Ub(N )2u(d100 )d 1000.000055、30.00000577(mm)0.55i30.05772u(N )N638.8 °.000130：O'0572V 0.031940210022.6 (nm)u()=6393( nm)(2)测量次数为奇数的逐差法 处理原则：去掉中间的数据。条纹吞吐n0100200300400例12：重新处理迈克尔逊干涉仪实验数据条纹吞吐n5006007008009001000M镜位置X (mm)34.6446534.6763534.7080034.7394534.7708534.80280M2 镜位置 X (mm) 3

24、4.4830534.5158534.5483034.5806034.61300M镜位置X (mm)34.4830534.5158534.5483034.5806034.61300i12345平均N=ni+6 ni600600600600600600D6oo=X +6-X (mm)P 0.193300.192150.191150.190250.18980 ：0.1913322d600N2 °19133637.8(nm)600u()u(d600 )d600u(N)N其中u()U( d600 ), U a (d 600 )Ub (d 600 )2(d i d 600 )5 40.00005

25、ub(d 600 )0.0000289(mm)Ua (d 600 )0.000636(mm)0 5u(N) ub(N)0.289v'3I22u(d600)U(N),d 600Nu()=6382(nm)" 20.0006366378,2 0.1913320.289 22.1( nm) 600解：去掉中间的数据后为条纹吞吐n0100200300400M镜位置X (mm)34.4830534.51585 :34.54830 :34.58060 :34.61300条纹吞吐n6007008009001000M镜位置X (mm)34.6763534.7080034.7394534.7708534.80280由逐差法可得(3)逐差法说明 逐差法多用在自变量等间隔测量且其测量误差可略去的情况，这样可简化计算。 使用逐差法