不确定度与数据处理

1、用标准偏差S代替标准差S(x)单次测量的标准偏差结果表述:xiS(x(置信概率 68.3%)真值的估计值单次测量标准差最佳估计值S(x)的物理意义：在有限次测量中，每个测量值平均所具有的标准偏差。 通常不严格区分标准差与标准偏差，统称为标准差。(2)平均值的标准差真值的最佳估计值是平均值，故结果应表述为：(并不是只做一次测量)x/真值的最佳估计值S( x )(置信概率 68.3%)准差最佳估计值其中S(x)平均值的标准偏差 , k(k 1)不确定度与数据处理误差与不确定度1 .误差与不确定度的关系(1)误差：测量结果与客观真值之差x=x-A公认值一如物理常数等其中A称为真值，一般不可能准确知道

2、，常用约定真值代替：标准值一更高精度仪器测量结果理论值一理论公式计算结果对一个测量过程，真值 A的最佳估计值是平均值 X。在上述误差公式中，由于 A不可知，显然 x也不可知，对误差的最佳估计值是不确定度u(x)o(2)不确定度：对误差情况的定量估计，反映对被测量值不能肯定的程度。通常所说“误差” 一般均为“不确定度”含义。不确定度分为A、B两个分量，其中A类分量是可用统计方法估计的分量，它的主要成分是随机误差。2 .随机误差：多数随机误差服从正态分布。定量描述随机误差的物理量叫标准差。(1)标准差与标准偏差标准差.真值A不可知,尸L且测量次数 k为有限次.,实际上也不可知，于是:例1：某观察量

3、的n次独立测量的结果是X1, X2, ,Xn。试用方差合成公式证明平均值的标准偏差是样本标准偏差的S(X)S(X).n X解： X -L由题知Xi相互独立，则根据方差合成公式有 nu(X)u2(X1)u2(Xn)n利用样本标准偏差白定义，可知u(Xi)=S(X)u(X) S(X)S2(X)S2(X)ni=1,2, ,n.nS2(X)nS(X).n3 .系统误差与仪器误差(限)(1)系统误差：在同一被测量的多次测量过程中，保持恒定或以可以预知方式变化的那一部分误差称为系统误差。 已被确切掌握了其大小和符号的系统误差，称为可定系统误差；对大小和符号不能确切掌握的系统误差称为未定系统误 差。前者一般

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

5、%Nx+n字式中a是数字式电表的准确度等级，Nx是显示的读数，b是误差的绝对项系数， Nm是仪表的满度值,误差，相当于最小量化单位的倍数。n代表仪器固定项电阻箱:仪=ai % Ri R0式中Ro是残余电阻，Ri是第i个度盘的示值，ai是相应电阻度盘的准确度级别。直流电位差计：仪=a% (Ux U-0-)10式中a是电位差计的准确度级别，Ux是标度盘示值，U0是有效量程的基准值，规定为该量程中最大的10的整数哥。直流电桥:入 R。仪=a%Rx )10式中Rx是电桥标度盘示值，a是电桥的准确度级别， R0是有效量程的基准值，意义同上。(3) B类不确定度的处理在物理实验中，B类不确定度的来源通常包

6、括以下三种：仪器误差仪、灵敏度误差灵和估计误差限估。其中灵敏度误差可表不为0.20.2灵-Z-Sn/ xB类不确定度与各种误差限之间的关系为Ub4.不确定度的合成直接测量u(x)间接测量x： ua(x) , ub(x) ,u2(x) u；(x) y=f(x1, x2, xn)(称为合成不确定度)(3)u(y),f、2()u i xi(Xi)其中x1, x2,- u(y)或 y,Xn为相互独立的直接测量量lnf 2 2 / 、 )u (为)Xi最终结果表述形式:结果有效数字的确定原则：N u(N)= 不确定度测量结果(单位)u(N)只保留一位有效数字；N与不确定度u(N)小数位数对齐。例2：用分

7、光计测棱镜材料的折射率公式为Asinn 2.A sin 2。已测得A=60 0' 2',黄光(汞灯光源)所对应的=50 58'3',则黄光所对应的折射率n u(n)= 1.6479 0.0007解:dnAsinn 2.A sin 2A J h cos(dA22Asin260 0 50 58 sin2.60 0 sin21d ) cosAldA2 2 2Asin 一 21.6479A lnn lnsin 2A lnsin 2,A1 . Actg )dA ctg 22d 2u(n) 1 AA 2 21 . 2 A 24 (函2 ctg 2) u (A) 4 ctg

8、2u ()160 050 5860 0 2 2180 2(ctg ctg)2( 一 )2.42260u(n) n u(n 1.6479 0.000426 0.000712 60 0 50 58ctg 42180 2一)2 0.000426n5.有效数字及其运算法则(1)有效数字：由若干位可靠数字加一位可疑数字构成。n u(n)=1.6479 0.0007在不计算不确定度的情况下，结果的有效数字由运算法则决定。(2)运算法则 加减法：以参加运算各量中有效数字最末一位位数最高的为准并与之取齐。N=A+B-C-D,则 u(N) vu2(A) u2(B) u2(C) U2(D)取决于u(A)、u(B)

9、、u(C)、u(D)中位数最高者，最后结果与之对齐。 乘除法：以参加运算各量中有效数字最少的为准，结果的有效数字个数与该量相同。N 组 则 u(N u(A) 2 u(B) 2 u(C) 2 u(DPCDNA B C D取决于其中相对不确定度最大者，即有效数字个数最少者。混合四则运算按以上原则按部就班执行。例3：某物理量的计算公式为Y ,其中k为常数，1.6为准确数，H = 16cm, d=0.1500cm。若使丫的1 1.6d/H表示式中分母的值具有 4位有效数字，正确测 H的方法是(d )。(a)用游标卡尺估读到 cm千分位(b)用米尺估读到cm百分位(c)用米尺只读到 mm位(d)用米尺只

10、读到cm位1.6d 1.6 0.15001.6d解：0.015 分母 1 1.015为4位有效数子H 16H即H只需2位有效数字即可，故应选 (d)。 特殊函数的有效数字：根据不确定度决定有效数字的原则，从不丢失有效位数的前提出发，通过微分关系传播 处理。例4:tg45 2' =1.00116423最多可取几位有效数字？1解： 令 y=tgx ,其中 x=45 2'取 x 1 0.00029(rad)6018011,去、_»则 y x 2 0.00029 0.00058即小数点后第四位产生误差cos x cos 45 2tg45 2' =1.0012 ,有五位

11、有效数字。例5：双棱镜测波长的计算公式为xVbb ,对实验数据进行处理的计算结果如下表所示。S Sx=0.28144mmb=5.9325mmb'=0.7855mmS=27.65cmS'=75.90cmu( x)=2.010 10-4m m(b)/b=0.025(b')/b'=0.025(S) =0.5cm(S') =0.5cm(b)=0.005mm(b')=0.005mm(S) =0.05cm(S') =0.05cm注：下标1代表来自方法误差，下标 2代表来自仪器误差 要求：(1)给出测量结果的正确表述(包括必要的计算公式)(2)定量讨论

12、各不确定度的分量中，哪些是主要的，哪些是次要的，哪些是可以忽略的？如果略去次要因素和可以忽略项的贡献，不确定度的计算将怎样简化？结果如何?解：(1)0.28144 -5.9325。迺 5.86716 10 4mmS S(276.5 759.0)其中InIn xu()222u(b )2u(S)u( x)u(b)x2b2bS Su( x)2.01010 40.000714 ;1 , ,1d d( x)-Inb -Inb In(S S) - 22xx0.28144db db2b 2b2u(S)S SdSS S0.0111dSS SW (b)2b1(b)/、3132,3 b0.0252“30.0072

13、22bu2(b)2(b)/、.30.0050 0002432b2b2 5.9325,3u1(b )1(b)/、.311(b)0.025'0.007222b2b2、3 b2,3u2(b )2(b)/ .30.005n2b2b2 0.7855.3u1(S)1 (S) / , 30.50.00279S SS S(27.65 75.90).3u2(S)2(S)/ .30.050.000279S SS S(27.65 75.90).3u1(S)i(S)/ , 30.50.00279S SS S(27.65 75.90) . 3u2(S)2(S)/ .30.050.000279S SS S(27.

14、65 75.9)0) , 32u(b)2u1(b)2U2(b)2b2b2b2u(b )5(b) 2必9) 22b2b2b2u(S)2Ui(S)2U2(S)S SS SS S2u(S)2Ui(S)2U2(S)S SS SS S于是得 u( )=un 5.86716 10 4 0.0111=6.53 10-6mmu( )=587 7nmU1(S)和U1(S),可以忽略的因S S S S(2)由前面的计算可知，不确定度主要来自u1(b)和 u1(b) ,次要因素是U2(b )2b 2b2b素是一、处L吧和里Stx 2b S S S S22若只考虑主要项的贡献：山 u1(b) u1(b)1(b) 0.

15、01022b 2b6b则有 u( )=6nmu( )=587 6nm比严格计算的结果稍小但相差无几。二、数据处理方法1 .列表法：按一定规律把数据列成表格。列表原则：(1)表格的标题栏中 注明物理量的名称、符号和单位；(2)记录原始数据(如记录刻度数，而不是记录长度)；(3)简单处理结果(如算出长度)或函数关系；(4)参数和说明(如表格名称、仪器规格、环境参数、常量以及公用单位等)。2 .作图法：把实验数据用自变量和因变量的关系作成曲线，以便反映它们之间的变化规律或函数关系。作图要点：(1)原始数据列表表示 一一见列表法；(2)用坐标纸作图，图纸大小以不损失有效数字和能包括所有点为最低要求，因

16、此至少应保证坐标纸的最小分格(通常为1mm)以下的估计位与实验数据中最后一位数字对应；(3)选好坐标轴并标明有关物理量的名称(或符号)、单位和坐标分度值。其中分度比例一般取1、2、5、10较好，以便于换算和描点；(4)实验数据点以+、X、 、 、等符号标出，一般不用细圆点“ ”标示实验点；光滑连接曲线 并使实验点匀称 地分布于曲线两侧(起平均的作用)；(5)图解法求直线斜率和截距时，应：在线上取点(不能使用实验点)；所取两点要相距足够远 (以提高精度)；在图上要注明所取点的坐标 。例6:拉伸法测弹性模量的载荷一一伸长曲线如图所示，图上至少有 不规。它们是 坐标轴应标注物理量和单位，轴上缺少分度

17、值标记标出， 曲线应光滑连接，计算点坐标标注不规 。3.最小二乘法与一元线性回归法5处绘制错误或 实验点应以醒目各测量值与这条曲线上对应点之差的平方和应(1)最小二乘法：对等精密度测量若存在一条最佳的拟合曲线，那么 取极小值。例7:试用最小二乘原理推导直线方程解：根据最小二乘原理应有n(yik i 1kx,)20n(yii 1n2(yii 1y=kx中回归系数2kxi)minkxi)( Xi) 0k的计算公式。nxiyii 1xi20于是得xyxi y-2-xi(2) 一元线性回归法:由最小二乘原理,应有设直线方程ky=a+bx,其中自变量x的误差可略yi (a12bxi) minkyi即 a

18、i 1(abx)2k2yii 1(a bx/( 1)akxi1kyii 1解之得kbi"(a(3)相关系数r物理意义bxi )2r：yik2yii 1(a bxi)(4)kxiyii 1(xj2xiyikxi(x)xiyi2xiyi vx2xixy xy-2 xbx用于检验x和y之间是否存在线性关系。xyxyx2x2)(y2y2)y a bx®过全部实验点 xryi之间线性相关强烈 yi随xj增加而增加yi随xj增加而减小xyi之间无线性关系拟合直线为与x轴平行的直线例8:根据所给相关系数 r=-1r作出实验点分布草图: r=0.9993 r=0.015|r| 1。解：根据

19、线路图可得R ()20.050.0100.0200.0300.0400.0V (V)2.802.722.602.382.202.04Rv (写出计算公式即可)EVRRvRRvRv1 -计算R、-精度:Vu(R) Ru(1 /V)1/V0.10.1200'400.0u(V) V故将公式变形为LrERv则有1_ErVaRV0.005,0.000250.01 0.01 ， 2.80 2.04可知r的精度较高0.0036,0.005并设直线方程y=a+bx(4)回归法使用要点： 自变量x测量误差可略，即应选择测量精度较高的物理量作自变量； 因变量y为等精度测量 或近似等精度测量，即 u(yi)

20、近似相等；作线性关系的检验：利用物理规律或作图等其它方法确认线性关系的存在；或检验相关系数是否满足例9:实验线路及测量数据如下，用一元线性回归法计算电压表阻E RERvV4.逐差法(1)测量次数为偶数的逐差法隔n项逐差，可得到b1 小一y1xn 1 x1y a bx,并件-组实验数据：x1, xn , xn 1, x2ny1，, yk, yn 1,y2nhy2nbn必取平均值b 1nbi 1n 绚x2nxnn i 1n i 1 xn i 为设自变量和因变量之间存在线性关系对于自变量x等间隔分布的情况，有 xn求得同后，可由公式yi a b X1 n ，xinx 于b(yn i yi )n nx

21、 i 1求出a yi bxi于是有R 0.02335t 20.00.0011675例10：已知R=Ro(i+ t),实验数据如下，用逐差法求电阻温度系数(不要求计算不确定度)t (C)85.080.075.070.065.060.055.050.0R ()0.36220.35650.34990.34370.33800.33240.32700.3215解：R=Ro+ Ro t 并设 y=a+bx 则有a= Rob= Ro =a即-而利用逐差法可得:ai1234平均t=ti+4-ti ( C )20.020.020.020.020.0R= Ri+4-Ri ()0.02420.02410.02290

22、.02220.0233511a -( R b t) -(2.7312 0.0011675 540.0) 0.2626 n80.00116750.2626一一 一 34.45 10(1/C)条纹吞吐n500600700800900M2镜位置X (mm)34.6446534.6763534.7080034.7394534.77085解法其中2d 500N2 0.15970638.8(nm)500U(d500 )./Ua (d 500 ).2 /Ub (d500 )u( ) u(d 500)d d 500(d id 500 )U a (d 500 ) ,542 u(N)N0.000648(mm)0.

23、00005ub(d500) 0.0000289( mm),3例11:迈克尔逊干涉仪实验数据处理条纹吞吐n0100200300400M2镜位置X (mm)34.4830534.5158534.5483034.5806034.61300i12345平均N=ni+5-ni500500500500500500d500=Xi+5-Xi (mm)0.161600.160500.159700.158850.157850.15970由逐差法可得0.5 u(N) ub(N)0.2893u(2u(d50。)d500u(N)N0.000648 2638.8、?0.1597020.289222.6(nm)5002u(

24、 )=639 3(nm)i12345平均N'=( ni+5-ni)/5100100100100100100d100=(Xi+5-Xi ) /5 (mm)0.0323200.0321000.0319400.0317700.0315700.031940解法由逐差法可得u()22 u(d100)2d 100N2 0.031940 638.8(nm)100d100u(N )N其中U(di00 ). Ua (d 100 ).2 /Ub (d100 )U a (di00 )U b ( d100 )(di d100)25 40.00005 0.000130(mm) 0.00000577(mm)5.3

25、u(N ) Ub(N )0.5p 0.05775.3u()2u(d100 )d 1002u(N )N638.80.00013020.031940220.0577210022.6(nm)u( )=639 3(nm)测量次数为奇数的逐差法处理原则：去掉中间的数据。例12:重新处理迈克尔逊干涉仪实验数据条纹吞吐n0100200300400M2镜位置X (mm)34.4830534.5158534.5483034.5806034.61300条纹吞吐n5006007008009001000M2镜位置X (mm)34.6446534.6763534.7080034.7394534.7708534.80280条纹吞吐n0100200300400M2镜位置X (mm)34.4830534.5158534.5483034.5806034.61300解:去掉中间的数据后为条纹吞吐n6007008009001000M2镜位置X (mm)34.6763534.7080034.7394534.7708534.80280i12345平均N=ni+6-ni600600600600600600D600=Xi+6-Xi (mm)0.193