不确定度与大数据处理

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

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

3、)(置信概率-68.3%)真值的最佳式平;总滁准差最佳估计值其中S(X) JZ(Xi -X)2平均值的标准偏差 k(k-1)(并不是只做一次测量)例1：某观察量的n次独立测量的结果是X, X2,，Xn。试用方差合成公式证明平均值的标准偏差是样本标准偏差S(X)S(X) =-。n、X解：X =一L由题知X相互独立，则根据方差合成公式有nu(X)u2(X1)u2(Xn)n利用样本标准偏差白定义，可知u(X)=S(Xlu(X) =S(X)=S2(X)S2(X)ni=1,2, nnS2(X)= S(X) n . n3 .系统误差与仪器误差(限)(1)系统误差：在同一被测量的多次测量过程中，保持恒定或以

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

5、估计；钢板尺、螺旋测微计的仪器误差限按其最小分度的1/2计算。 指针式仪表：A«= a%Nm式中M是电表的量程，a是准确度等级。数字仪表：仪=a%N<+b%|n 或仪=a%N：+n字式中a是数字式电表的准确度等级， NX是显示的读数，b是误差的绝对项系数， M是仪表的?t度值，n代表仪器固定项误 差，相当于最小量化单位的倍数。电阻箱：&i=Zai%R+R0式中R是残余电阻，R是第i个度盘的示值，a是相应电阻度盘的准确度级别。直流电位差计：仪=a% (Ux 业)10式中a是电位差计的准确度级别，U是标度盘示值，U0是有效量程的基准值，规定为该量程中最大的10的整数哥。直流

6、电桥:仪=a%(RX +R0)10式中R是电桥标度盘示值，a是电桥的准确度级别， R是有效量程的基准值，意义同上。3 3) B类不确定度的处理仪、灵敏度误差灵和估计误差限估。其中灵敏在物理实验中，B类不确定度的来源通常包括以下三种：仪器误差 020 2度误差可表不为灵=02 = 0.2SAn/AxB类不确定度与各种误差限之间的关系为ub=: 。.34 .不确定度的合成(1)直接测量 X： Ua(x) , Ub(x)u(x) Ru2(x)为2 (x)(称为合成不确定度)(2)间接测量y=f(x1, x2,，xn)其中x1, x2,，xn为相互独立的直接测量量u(y)=；Z(£ )2u2

7、(xi)或吗=|£(1)12(为)i 二xiy . i 为(3)最终结果表述形式：Nu(N)= (单位)结果有效数字的确定原则：不确定度u(N)只保留一位有效数字；测量结果N与不确定度u(N)小数位数对齐。所对应的=50 58例2:用分光计测棱镜材料的折射率公式为n=2。已测得A=60力'±2',黄光(汞灯光源)Asin 一23',则黄光所对应的折射率n±j(n)= 1.6479 a0.0007。A 、60 0 50 58sin- sin； A、. A解： n =2- =2：=1.6479lnn=lnsin-lnsin.A. 60 022s

8、in sin 22dnnA5 11cos(dAd、)222A 、sin2cosA1dA2 2A sin 21 ,二2(ctgA 1 A 、-ctg)dA ctgd、222u(n) 1 A 、 A 2 21 1 2 A 、2/.、-=4(ctg ctg 2) u (A) 4ctg 2u ()60 0 50 582-ctg60 0、2/21809)(一')260 二1 + 24ctg60 0 50 582180)2 n=0.000426u(n) =n 幽=1.6479 0.000426 =0.0007实用文档5.有效数字及其运算法则(1)有效数字：由若干位可靠数字加一位可疑数字构成。在不计

9、算不确定度的情况下，结果的有效数字由运算法则决定。(2)运算法则 加减法：以参加运算各量中有效数字最末一位位数最高的为准并与之取齐。N=A+BGD 则 u(N)= u1 2(A) u2(B) u2(C) u2(D)取决于u(A)、u(B)、u(C)、u(D)中位数最高者，最后结果与之对齐。 乘除法：以参加运算各量中有效数字最少的为准，结果的有效数字个数与该量相同。AB N -CD,则u(N)向(A) 1+a9)7+冰)7+中(口)1 Eg 1石后g取决于其中相对不确定度最大者，即有效数字个数最少者。混合四则运算按以上原则按部就班执行。例3：某物理量的计算公式为丫=k1 1.6d/H表示式中分母

10、的值具有 4位有效数字，正确测 H的方法是(d )。,其中k为常数，1.6为准确数，H= 16cmx d=0.1500cm。若使丫的(a)(c)解:用游标卡尺估读到 cm千分位用米尺只读到 mniu(d)1.6d 1.6 0.1500 八=0.015(b)16分母用米尺估读到cm百分位 用米尺只读到cm位1+16d/.015为4位有效数字H即H只需2位有效数字即可，故应选 (d) 特殊函数的有效数字：根据不确定度决定有效数字的原则，从不丢失有效位数的前提出发，通过微分关系传播 处理。例4：tg452' =1.00116423解:令 y=tg x ,其中 x=45 2'最多可取几

11、位有效数字？取.：x=1=0.00029(rad)601802cos 45 2M0.000294.00058即小数点后第四位产生误差tg45 2 =1.0012 ,有五位有效数字。例5:二x bb一， ,e , 一,对实验数据进行处理的计算结果如下表所示。S Sx=0.28144mmb=5.9325mmb' =0.7855mmS=27.65cmS' =75.90cmu(b)/b =0.025(b' )/b' =0.025(S) =0.5cm(S) =0.5cmx)=2.010 M10-4mm(b)=0.005mm(b' )=0.005mm(S) =0.0

12、5cm(S) =0.05cm双棱镜测波长的计算公式为O注：下标1代表来自方法误差，下标 2代表来自仪器误差。要求：(1)给出测量结果的正确表述(包括必要的计算公式)(2)定量讨论各不确定度的分量中，哪些是主要的，哪些是次要的，哪些是可以忽略的？如果略去次要因素 和可以忽略项的贡献，不确定度的计算将怎样简化？结果如何？解：(1)也亚二。281445.9325 0.7855 =5.86716 10Mmm(276.5 759.0)x 2b 2b S S S Sd' djx) db db dS=十-T -dS实用文档22222暇+陶+黑黑、藁1 =0.0111其中U(Ax)=2.010父10&

13、quot; =0.000714 ；x 0.28144"(b)4(b)/J314(b)0.025ch”。=X=- =0.00722_2_2_2 4 2b 2b _273 b2,/3t ,U(b2 _ 产 I + 1U2(b)IU2(b).：2(b)/,30.0050 00024a' 一名."|L2b . "iW = = =0.0002432b 2b 2 5.9325,3U1 (b) _a(b” J3-1 M &(b)_°.025 _0 007222 2b'2b'2召b，2v3|U2(b) _& (b ) / 於0.0

14、05一-0 00184、2b'2b'2 X0.7855V3'U1(S)& (S) / <30.5一 S十S- S +S，l u.uury(27.65 +75.90) J3|U2(S) _a2(S)/v3-0.05-0 000279|S +S，S +S'(27.65 +75.90)V3-. 2.2. 2u(b) U1(b)U2(b)IL2b,: , I 2b,: . " |L 2b-u(S) 2 二 U1(S)_ 2 U2(S) 2S ' S'：|_S ' S'：" g SU1(S) Ai(S)/V

15、3 S +S' - S +S' u2(S) .：2(S) / ,3 S S S S='= 0.00279(27.65 75.90),30.05=0.000279(27.65 75.90) , 3U(S ) 2 _ U1(Si 2 U2 2 s s,:|_s ' s,:" ILs-S_ U( )=587 _7nm于是得U( )= 3=5.86716 10 H 0.0111=6.53 10-6mm_4避)_和_U1(S2 ,可以忽略的因S S S S(2)由前面的计算可知，不确定度主要来自U1 (b)和 U1(b) ,次要因素是 Ubj2b 2b2b素是3

16、、山®、粤和2。x 2b S S S SJ2. -2U1U1(b)=：1(b) =0.0102IL 2b . IL 2b .6b贝 U 有u( )=6nm± u( )=587 +6nm 比严格计算的结果稍小但相差无几。二、数据处理方法1,列表法：按一定规律把数据列成表格。列表原则：(1)表格的标题栏中 注明物理量的名称、符号和单位 ；(2)记录原始数据(如记录刻度数，而不是记录长度)；(3)简单处理结果(如算出长度)或函数关系；(4)参数和说明(如表格名称、仪器规格、环境参数、常量以及公用单位等)2.作图法：把实验数据用自变量和因变量的关系作成曲线，以便反映它们之间的变化规

17、律或函数关系。作图要点：(1)原始数据列表表示 见列表法；(2)用坐标纸作图，图纸大小以不损失有效数字和能包括所有点为最低要求，因此至少应保证坐标纸的最小分格 (通常为1mm以下的估计位与实验数据中最后一位数字对应；(3)选好坐标轴并标明有关物理量的名称(或符号)、单位和坐标分度值。其中分度比例一般取1、2、5、10较好，以便于换算和描点；(4)实验数据点以+、X、 、 、等符号标出，一般不用细圆点“ ”标示实验点；光滑连接曲线 并使实验点匀称实用文档地分布于曲线两侧(起平均的作用)；(5)图解法求直线斜率和截距时，应：在线上取点(不能使用实验点)；所取两点要相距足够远 (以提高精度)；在图上

18、要注明所取点的坐标 。例6:拉伸法测弹性模量的载荷一一伸长曲线如图所示，图上至少有 不规范。它们是 坐标轴应标注物理量和单位 ，轴上缺少分度值 目标记标出 ， 曲线应光滑连接 ，计算点坐标标注不规范 。3.最小二乘法与一元线性回归法(1)最小二乘法：对等精密度测量若存在一条最佳的拟合曲线，那么5处绘制错误或 实验点应以醒各测量值与这条曲线上对应点之差的平方和应取极小值。例7:试用最小二乘原理推导直线方程y=kx中回归系数k的计算公式。n解：根据最小二乘原理应有£ (yi -欧)=mini 1nnnn即 v (yi -kxi )2 =0 .二 2(yi -kxi)(-x ) =07.

19、Xi yi -kv x =0.k i 1i 1i Ji 1于是得、Xi yixy(2) 一元线性回归法:由最小二乘原理，应有二 k 2x yi -(a bxi) Saia.二 kx yi -(a bx。2I ; b i m解之得2 xi设直线方程y=a+bx,其中自变量_k. _.一 2.' yi 一(a bxj =mini 1x的误差可略=0k三 2yi (a+bxi)(1)=0k£2匹一(a +bxi)(为)=0kkak +bZ xi =£ yiiTiT：kkkax x bx x2 八 xiyi,i 1idi =1" xyi -ky xiyi，xj2

20、-kx x2“ xiyi ' xiyi yx2a =；272(“为)-Q x：xy - xy-22x - x=y -bxr物理意义r| =1|r| 两 1r >0r <0r : 0(3)相关系数r：用于检验x和y之间是否存在线性关系。y = a +bx通过全部实验点xi、yi之间线性相关强烈yi随xi增加而增加yi随xi增加而减小xyi之间无线性关系 t拟合直线为与x轴平行的直线例8:根据所给相关系数r作出实验点分布草图:实用文档 r=0.9993 r=0.015(4)回归法使用要点： 自变量x测量误差可略，即应选择测量精度较高的物理量作自变量； 因变量y为等精度测量 或近

21、似等精度测量，即 u(y。近似相等；作线性关系的检验：利用物理规律或作图等其它方法确认线性关系的存在；或检验相关系数是否满足| r | .1。例9:实验线路及测量数据如下，用一元线性回归法计算电压表内阻R (写出计算公式即可),E RR ()20.050.0100.0200.0300.0400.0V (V)2.802.722.602.382.202.04解：根据线路图可彳导E = R + RV =里R RVRVV VulR1 |_°11, -°1L = b.005,0.00025 】计算R 1精度：R-2°.° 400.0可知R的精度较高VU(1/V)

22、u(V) 0.01 0.011"A” j-=- I, = 0.0036,0.005、1/V V 12.80 2.041111故将公式变形为-=R+- 令 -=y R三x并设直线方程y=a+bxVERVEV贝 U有a = b =-=-RV =亘EERVRVb4.逐差法(1)测量次数为偶数的逐差法r 一设自变量和因变量之间存在线性关系，并有一组实验数据：3X1' ,xn,xn" x2n¥，yk,yn+，y2n隔n项逐差，可得到b1 =yn、-y1，bn = y2n yn取平均值xn 1 -七x2n - *n1 .n对于自变重x等间隔分布的情况，有 xn+ -x

23、i =Anx 于是 b=乙(yn+-yi)ni nx i =1求得 后，可由公式E yi =a+b£xi求出a yi -bZ xi(不要求计算不确定度)t (C)85.080.075.070.065.060.055.050.0R ()0.36220.35650.34990.34370.33800.33240.32700.3215例10:已知R=R(1+t),实验数据如下，用逐差法求电阻温度系数解：R=R+ R0t 并设 y=a+bx则有a= R0b= R0=a而利用逐差法可得:ba =一 a it =t i +4-1 i(c)R= Ri+4-R120.00.0242220.00.02

24、41320.00.0229420.0平均20.00.02220.0233 5于是有R 0.02335 =0.0011675.t20.011=1(% Rt) =1(2.7312 -0.0011675 540.0) =0.2626n80.0011675与=4.45 父10，(1/C)a 0.2626i12345平均N= n i +5- n500500500500500500d500=X+5-X (mm)0.161600.160500.15970 :0.158850.157850.159702d5002 0.15970500= 638.8(nm)u( ) u(d 500) 2 u(N) 2一二,d5

25、00_ N其中U(d500 ) = ：, Ua (d 500 ) ub (d500 )ua (d 500 ) 二% (di -d 500 )2 i 500=0.000648(mm)ub (d 500 )5 4 0.00005 :0.0000289(mm) 、311:迈克尔逊干涉仪实验数据处理条纹吞吐n0100200300400M镜位置X (mm)34.4830534.5158534.5483034.58060P 34.61300条纹吞吐n500600700800900M镜位置X (mm)34.6446534.6763534.7080034.73945P 34.77085解法由逐差法可得u()=

26、 638.8.0.159700.000648 20.5 u(N)=ub(N) =- =0.289 -J30.2892，一 二2.6(nm) 5002u( )=639 二3(nm)解法二： 由逐差法可得i12345平均N'=(ni+5-ni )/5100100100100100100d100=(X+5-X) /5 (mm)0.0323200.0321000.0319400.0317700.0315700.031940=2"031940 =638.83)u(Z)Ru(d200)T + Lu(Nj fN100九 H d100 1! N'其中u(d100) ua(d100)

27、ub (d100 )ua (d100 ) =“ (di -d10。)25 40.00005 一= 0.000130(mm)ub(d100 )0.00000577(mm)5 3.1u(d100) T - F(N ) I2，d100_ N= 638.8(0.5u(N ) =5 (N ) = =0.0577、5V30.00013020.0577 22 2 ： 2.6(nm)0.03194021002u( )=639 :3(nm)(2)测量次数为奇数的逐差法 处理原则：去掉中间的数据。例12：重新处理迈克尔逊干涉仪实验数据条纹吞吐n0100200300400M镜位置X (mm)34.4830534.5158534.54830P 34.5806034.61300条纹吞吐n5006007008009001000M镜位置X (mm)34.6446534.6763534.7080034.7394534.7708534.80280解：去掉中间的数据后为条纹吞吐n0100200300400M镜位置X (mm)34.4830534.5158534.5483034.5806034.61300条纹吞吐n6007008009001000M镜位置X (mm)34.6763534.7080034.73945P 34.7708534.80280由逐差法可得i12345平均N=ni+6- ni