实验数据误差分析和数据处理13.pdf
-第二章 实验数据误差分析和数据处理 第一节 实验数据的误差分析 由于实验方法和实验设备的不完善，周围环境的影响，以及人的观察力，测量程序等限制，实验观测值和真值之间，总是存在一定的差异。人们常用绝对误差、相对误差或有效数字来说明一个近似值的准确程度。为了评定实验数据的准确性或误差，认清误差的来源及其影响，需要对实验的误差进展分析和讨论。由此可以判定哪些因素是影响实验准确度的主要方面，从而在以后实验中，进一步改良实验方案，缩小实验观测值和真值之间的差值，提高实验的准确性。一、误差的根本概念 测量是人类认识事物本质所不可缺少的手段。通过测量和实验能使人们对事物获得定量的概念和发现事物的规律性。科学上很多新的发现和突破都是以实验测量为根底的。测量就是用实验的方法，将被测物理量与所选用作为标准的同类量进展比拟，从而确定它的大小。1.真值与平均值 真值是待测物理量客观存在确实定值，也称理论值或定义值。通常真值是无法测得的。假设在实验中，测量的次数无限多时，根据误差的分布定律，正负误差的出现几率相等。再经过细致地消除系统误差，将测量值加以平均，可以获得非常接近于真值的数值。但是实际上实验测量的次数总是有限的。用有限测量值求得的平均值只能是近似真值，常用的平均值有以下几种:(1)算术平均值 算术平均值是最常见的一种平均值。设1x、2x、nx为各次测量值，n代表测量次数，则算术平均值为 nxnxxxxniin 121(2-1)(2)几何平均值 几何平均值是将一组n 个测量值连乘并开n次方求得的平均值。即 nnxxxx 21几(2-2)3均方根平均值 nxnxxxxniin 1222221均 (2-3)(4)对数平均值 在化学反响、热量和质量传递中，其分布曲线多具有对数的特性，在这种情况下表征平均值常用对数平均值。设两个量1x、2x，其对数平均值 21212121lnlnlnxxxxxxxxx对 (2-4)-应指出，变量的对数平均值总小于算术平均值。当1x/2x2 时，可以用算术平均值代替对数平均值。当1x/2x=2，对x=1.443,x1.50,(对x-x)对x=4.2%,即1x/2x2，引起的误差不超过4.2%。以上介绍各平均值的目的是要从一组测定值中找出最接近真值的那个值。在化工实验和科学研究中，数据的分布较多属于正态分布，所以通常采用算术平均值。2.误差的分类 根据误差的性质和产生的原因，一般分为三类：1系统误差 系统误差是指在测量和实验中未觉察或未确认的因素所引起的误差，而这些因素影响结果永远朝一个方向偏移，其大小及符号在同一组实验测定中完全一样，当实验条件一经确定，系统误差就获得一个客观上的恒定值。当改变实验条件时，就能发现系统误差的变化规律。系统误差产生的原因：测量仪器不良，如刻度不准，仪表零点未校正或标准表本身存在偏差等；周围环境的改变，如温度、压力、湿度等偏离校准值；实验人员的习惯和偏向，如读数偏高或偏低等引起的误差。针对仪器的缺点、外界条件变化影响的大小、个人的偏向，待分别加以校正后，系统误差是可以去除的。2偶然误差 在已消除系统误差的一切量值的观测中，所测数据仍在末一位或末两位数字上有差异，而且它们的绝对值和符号的变化，时而大时而小，时正时负，没有确定的规律，这类误差称为偶然误差或随机误差。偶然误差产生的原因不明，因而无法控制和补偿。但是，倘假设对*一量值作足够屡次的等精度测量后，就会发现偶然误差完全服从统计规律，误差的大小或正负的出现完全由概率决定。因此，随着测量次数的增加，随机误差的算术平均值趋近于零，所以屡次测量结果的算数平均值将更接近于真值。3过失误差 过失误差是一种显然与事实不符的误差，它往往是由于实验人员粗心大意、过度疲劳和操作不正确等原因引起的。此类误差无规则可寻，只要加强责任感、多方警觉、细心操作，过失误差是可以防止的。3、精细度、准确度和准确度 反映测量结果与真实值接近程度的量，称为精度亦称准确度。它与误差大小相对应，测量的精度越高，其测量误差就越小。精度应包括精细度和准确度两层含义。1精细度：测量中所测得数值重现性的程度，称为精细度。它反映偶然误差的影响程度，精细度高就表示偶然误差小。2准确度 测量值与真值的偏移程度，称为准确度。它反映系统误差的影响精度，准确度高就表示系统误差小。3准确度精度 它反映测量中所有系统误差和偶然误差综合的影响程度。在一组测量中，精细度高的准确度不一定高，准确度高的精细度也不一定高，但准确度高，则精细度和准确度都高。为了说明精细度与准确度的区别，可用下述打靶子例子来说明。如图2-1所示。图 2-1(a)中表示精细度和准确度都很好，则准确度高；图 2-1(b)表示精细度很好，但准确度却不高；图 2-1(c)表示精细度与准确度都不好。在实际测量中没有像靶心那样明确的真值，-而是设法去测定这个未知的真值。学生在实验过程中，往往满足于实验数据的重现性，而忽略了数据测量值的准确程度。绝对真值是不可知的，人们只能订出一些国际标准作为测量仪表准确性的参考标准。随着人类认识运动的推移和开展，可以逐步逼近绝对真值。a b c 图 2-1 精细度和准确度的关系 4、误差的表示方法 利用任何量具或仪器进展测量时，总存在误差，测量结果总不可能准确地等于被测量的真值，而只是它的近似值。测量的质量上下以测量准确度作指标，根据测量误差的大小来估计测量的准确度。测量结果的误差愈小，则认为测量就愈准确。1绝对误差 测量值X和真值0A之差为绝对误差，通常称为误差。记为：0AXD(2-5)由于真值0A一般无法求得，因而上式只有理论意义。常用高一级标准仪器的示值作为实际值A以代替真值0A。由于高一级标准仪器存在较小的误差，因而A不等于0A，但总比X更接近于0A。X与A之差称为仪器的示值绝对误差。记为 AXd(2-6)与d相反的数称为修正值，记为 XAdC(2-7)通过检定，可以由高一级标准仪器给出被检仪器的修正值C。利用修正值便可以求出该仪器的实际值A。即 CXA(2-8)2相对误差 衡量*一测量值的准确程度，一般用相对误差来表示。示值绝对误差d与被测量的实际值A的百分比值称为实际相对误差。记为%100AdA (2-9)以仪器的示值X代替实际值A的相对误差称为示值相对误差。记为%100XdX (2-10)一般来说，除了*些理论分析外，用示值相对误差较为适宜。3引用误差 为了计算和划分仪表准确度等级，提出引用误差概念。其定义为仪表示值的绝对误差与量程围之比。%100%100nAXd量程范围示值绝对误差 (2-11)d-示值绝对误差；nX-标尺上限值-标尺下限值。4算术平均误差 算术平均误差是各个测量点的误差的平均值。ndi平ni,2,1(2-12)n测量次数；-id为第 i次测量的误差。5标准误差 标准误差亦称为均方根误差。其定义为 ndi2(2-13)上式使用于无限测量的场合。实际测量工作中，测量次数是有限的，则改用下式 12ndi (2-14)标准误差不是一个具体的误差，的大小只说明在一定条件下等精度测量集合所属的每一个观测值对其算术平均值的分散程度，如果的值愈小则说明每一次测量值对其算术平均值分散度就小，测量的精度就高，反之精度就低。在化工原理实验中最常用的U形管压差计、转子流量计、秒表、量筒、电压等仪表原则上均取其最小刻度值为最大误差，而取其最小刻度值的一半作为绝对误差计算值。5、测量仪表准确度 测量仪表的准确等级是用最大引用误差又称允许误差来标明的。它等于仪表示值中的最大绝对误差与仪表的量程围之比的百分数。%100%100maxmaxnnXd量程范围最大示值绝对误差(2-15)式中：ma*仪表的最大测量引用误差；dma*仪表示值的最大绝对误差；*n标尺上限值标尺下限值。通常情况下是用标准仪表校验较低级的仪表。所以，最大示值绝对误差就是被校表与标准表之间的最大绝对误差。测量仪表的精度等级是国家统一规定的，把允许误差中的百分号去掉，剩下的数字就称为仪表的精度等级。仪表的精度等级常以圆圈的数字标明在仪表的面板上。例如*台压力计的允许误差为1.5%，这台压力计电工仪表的精度等级就是1.5，通常简称1.5 级仪表。仪表的精度等级为 a，它说明仪表在正常工作条件下，其最大引用误差的绝对值ma*不能超过的界限，即%100maxmaxaXdnn(2-16)由式(2-16)可知，在应用仪表进展测量时所能产生的最大绝对误差简称误差限为 nXad%max(2-17)而用仪表测量的最大值相对误差为 XXaXdnnn%maxmax(2-18)由上式可以看出，用只是仪表测量*一被测量所能产生的最大示值相对误差，不会超过仪表允许误差a%乘以仪表测量上限*n与测量值*的比。在实际测量中为可靠起见，可用下式对仪表的测量误差进展估计，即 XXanm%(2-19)-例 2-1 用量限为 5A，精度为 0.5 级的电流表，分别测量两个电流，I1=5A,I2=2.5A,试求测量I1和I2的相对误差为多少.由此可见，当仪表的精度等级选定时，所选仪表的测量上限越接近被测量的值，则测量的误差的绝对值越小。例2-2 欲测量约90V 的电压，实验室现有0.5 级0-300V 和1.0 级0-100V 的电压表。问选用哪一种电压表进展测量为好.用0.5级0-300V 的电压表测量90V的相对误差为 用1.0级0-100V 的电压表测量90V的相对误差为 上例说明，如果选择得当，用量程围适当的 1.0 级仪表进展测量，能得到比用量程围大的0.5 级仪表更准确的结果。因此，在选用仪表时，应根据被测量值的大小，在满足被测量数值围的前提下，尽可能选择量程小的仪表，并使测量值大于所选仪表满刻度的三分之二，即*2*n/3。这样就可以到达满足测量误差要求，又可以选择精度等级较低的测量仪表，从而降低仪表的本钱。二、有效数字及其运算规则 在科学与工程中，该用几位有效数字来表示测量或计算结果，总是以一定位数的数字来表示。不是说一个数值中小数点后面位数越多越准确。实验中从测量仪表上所读数值的位数是有限的，而取决于测量仪表的精度，其最后一位数字往往是仪表精度所决定的估计数字。即一般应读到测量仪表最小刻度的十分之一位。数值准确度大小由有效数字位数来决定。1、有效数字 一个数据，其中除了起定位作用的0外，其他数都是有效数字。如 0.0037 只有两位有效数字，而 370.0 则有四位有效数字。一般要求测试数据有效数字为 4 位。要注意有效数字不一定都是可靠数字。如测流体阻力所用的 U 形管压差计，最小刻度是 1mm，但我们可以读到0.1mm，如 342.4mmHg。又如二等标准温度计最小刻度为 0.1，我们可以读到 0.01，如 15.16。此时有效数字为 4 位，而可靠数字只有三位，最后一位是不可靠的，称为可疑数字。记录测量数值时只保存一位可疑数字。为了清楚地表示数值的精度，明确读出有效数字位数，常用指数的形式表示，即写成一个小数与相应10的整数幂的乘积。这种以10 的整数幂来记数的方法称为科学记数法。如 75200 有效数字为4位时，记为7.520*105 有效数字为3 位时，记为7.52*105 有效数字为2位时，记为7.5*105 0.00478 有效数字为4位时，记为4.780*10-3 有效数字为3位时，记为4.78*10-3 有效数字为2位时，记为4.7*10-3 2、有效数字运算规则 1记录测量数值时，只保存一位可疑数字。2当有效数字位数确定后，其余数字一律舍弃。舍弃方法是四舍六入，即末位有效数字后边第一位小于 5，则舍弃不计；大于 5 则在前一位数上增 1；等于 5 时，前一位为奇数，则进 1 为偶数，前一位为偶数，则舍弃不计。这种舍入原则可简述为：小则舍，大则入，正好等于奇变偶。如：保存4 位有效数字 3.717293.717;-5.142855.143 7.623567.624 9.376569.376 3在加减计算中，各数所保存的位数，应与各数中小数点后位数最少的一样。例如将24.65 0.0082 1.632三个数字相加时，应写为 24.65+0.01+1.63=26.29。4在乘除运算中，各数所保存的位数，以各数中有效数字位数最少的那个数为准；其结果的有效数字位数亦应与原来各数中有效数字最少的那个数一样。例如：0.012125.641.05782 应写成 0.012125.641.06=0.328。上例说明，虽然这三个数的乘积为0.3281823，但只应取其积为0.328。5在对数计算中，所取对数位数应与真数有效数字位数一样。三、误差的根本性质 在化工原理实验常直接测量或间接测量得到有关的参数数据，这些参数数据的可靠程度如何.如何提高其可靠性.因此，必须研究在给定条件下误差的根本性质和变化规律。1、误差的正态分布 如果测量数列中不包括系统误差和过失误差，从大量的实验中发现偶然误差的大小有如下几个特征：1绝对值小的误差比绝对值大的误差出现的时机多，即误差的概率与误差的大小有关。这是误差的单峰性。2绝对值相等的正误差或负误差出现的次数相当，即误差的概率一样。这是误差的对称性。3极大的正误差或负误差出现的概率都非常小，即大的误差一般不会出现。这是误差的有界性。4随着测量次数的增加，偶然误差的算术平均值趋近于零。这叫误差的低偿性。根据上述的误差特征，可疑的出误差出现的概率分布图，如图 2-2 所示。图中横坐标表示偶然误差，纵坐标表示个误差出现的概率，图中曲线称为误差分布曲线，以)(xfy 表示。其数学表达式有高斯提出，具体形式为：22221xey(2-20)或 22xhehy(2-21)上式称为高斯误差分布定律亦称为误差方程。式中 为标准误差，h 为准确度指数，和 h的关系为 21y(2-22)假设误差按函数关系分布，则称为正态分布。越小，测量精度越高，分布曲线的峰越高切窄；越大，分布曲线越平坦且越宽，如图 1-3所示。由此可知，越小，小误差占的比重越大，测量精度越高。反之，则大误差占的比重越大，测量精度越低。-2、测量集合的最正确值 在测量精度一样的情况下，测量一系列观测值1M,2M,3M,nM所组成的测量集合，假设 图 2-2 误差分布 其平均值为mM,则各次测量误差为 miiMMx,i=1、2n，当采用不同的方法计算平均值时，所得到误差值不同，误差出现的概率亦不同。假设选取适当的计算方法，使误差最小，而概率最大，由此计算的平均值为最正确值。根据高斯分布定律，只有各点误差平方和最小，才能实现概率最大。这就是最小乘法值。由此可见，对于一组精度一样的观测值，采用算术平均得到的值是该组观测值的最正确值。图2-3 不同 的误差分布曲线 3、有限测量次数中标准误差 的计算 由误差根本概念知，误差是观测值和真值之差。在没有系统误差存在的情况下，以无限屡次测量所得到的算术平均值为真值。当测量次数为有限时，所得到的算术平均值近似于真值，称最正确值。因此，观测值与真值之差不同于观测值与最正确值之差。令真值为A，计算平均值为a，观测值为M，并令d=M-a，D=M-A，则 因为 0naMinaMi 代入nAMDii中，即得 nDAai 223 将式223式代入di=Mi-a 中得 nDDnDAMdiiiii)(224 将式224两边各平方得 对i 求和2222)()(2nDnnDDdiiii 因在测量中正负误差出现的时机相等，故将 Di2展开后，D1D2、D1D3，为正为负的数目相等，彼此相消，故得 从上式可以看出，在有限测量次数中，自算数平均值计算的误差平方和永远小于自真值计算的误差平方和。根据标准误差的定义 式中 Di2代表观测次数为无限多时误差的平方和，故当观测次数有限时，12ndi225 4可疑观测值的舍弃 由概率积分知，随机误差正态分布曲线下的全部积分，相当于全部误差同时出现的概率，-即 121222dxepx226 假设误差*以标准误差 的倍数表示，即*=t，则在t 围出现的概率为 2(t)，超出这个围的概率为1-2(t)。(t)称为概率函数，表示为 ttdtet02221)(227 2(t)与 t 的对应值在数学手册或专著中均附有此类积分表，读者需要时可自行查取。在使用积分表时，需 t 值。由表 2-1 和图2-4给出几个典型及其相应的超出或不超出|*|的概率。由表 2-1 知，当 t=3,|*|=3 时，在 370 次观测中只有一次测量的误差超过 3 围。在有限次的观测中，一般测量次数不超过十次，可以认为误差大于 3，可能是由于过失误差或实验条件变化未被觉察等原因引起的。因此，但凡误差大于 3 的数据点予以舍弃。这种判断可疑实验数据的原则称为3 准则。5函数误差 上述讨论主要是直接测量的误差计算问题，但在许多场合下，往往涉及间接测量的变量，所谓间接测量是通过直接测量的量之间有一定的函数关系，并根据函数被测的量，如传热问题中的传热速率。因此，间接测量值就是直接测量得到的各个测量值的函数。其测量误差是各个测量值误差的函数。图 2-4 误差分布曲线的积分 表2-1 误差概率和出现次数 t|*|=t 不超出|*|的 概率2(t)超出|*|的概率 1-2(t)测量次数 n 超出|*|的测量次数 0.67 0.67 0.49714 0.50286 2 1 1 1 0.68269 0.31731 3 1 2 2 0.95450 0.04550 22 1 3 3 0.99730 0.00270 370 1 4 4 0.99991 0.00009 11111 1（1）函数误差的一般形式 在间接测量中，一般为多元函数，而多元函数可用下式表示：y=f(*1,*2,*n)228 式中 y间接测量值；*i直接测量值。由台劳级数展开得 nnxxfxxfxxfy2211229 或 iniixxfy1-它的最大绝对误差为iniixxfy1 230 式中 ixf 误差传递系数；*i直接测量值的误差；y 间接测量值的最大绝对误差。函数的相对误差 为 231 2*些函数误差的计算 函数y=*z 绝对误差和相对误差 由于误差传递系数1,1zfxf，则函数最大绝对误差 y=|*|+|z|232 相对误差 zxzxyyr 233 函数形式为wxzKy，*、z、w 为变量 误差传递系数为：wKzxy 函数的最大绝对误差为 wwKxzzwKxxwKzy2 234 函数的最大相对误差为 wwzzxxyyr 235 现将*些常用函数的最大绝对误差和相对误差列于表22 中。例2-3 用量热器测定固体比热容时采用的公式 OpHpCttmttMC2)()(2102 式中 M量热器水的质量 m被测物体的质量 t0 测量前水的温度 t1 放入量热器前物体的温度 t2 测量时水的温度 CpH2O水的热容，4.187Kj/(kg.K)测量结果如下：M=2500.2gm=62.310.02g-t0=13.520.01t1=99.320.04 t2=17.790.01 试求测量物的比热容之真值，并确定能否提高测量精度。解：根据题意，计算函数之真值，需计算各变量的绝对误差和误差传递系数。为了简化计算，令0=t2-t0=4.27,1=t1t2=81.53,.方程改写为OpHpCmMC210 表2-2 *些函数的误差传递公式 函数式 误 差 传 递 公 式 最大绝对误差y 最大相对误差r 321xxxy|)|(|321xxxy yyr/21xxy|)|(|21xxy yyr/21xxy|)|(|1221xxxxy)(2211xxxxr 321xxxy|)|(|132231321xxxxxxxxxy)(332211xxxxxxr nxy )(1xnxyn)(xxnr nxy )1(11xxnyn)1(xxnr 21/xxy )(222112xxxxxy)(2211xxxxr cxy xcy)(xxr xylg xx0.4343y yyr/xyln xxy yyr/各变量的绝对误差为 各变量的误差传递系数为 函数的绝对误差 =3.5210-30.21.4110-20.02+0.2060.021.0810-20.05 =0.70410-30.28210-3+4.1210-3-0.5410-3=4.0010-3 J/(gK)880.0187.453.8131.6227.4250pC J/(gK)故真值Cp=0.87980.0003 J/(gK)由有效数字位数考虑以上的测量结果清度已满足要求。假设不仅考虑有效数字位数，尚需从比拟各变量的测量精度，确定是否有可能提高测量精度。则本例可从分析比拟各变量的相对误差着手。-各变量的相对误差分别为 其中以0的相对误差为 0.468%，误差最大，是 M 的 5.85 倍，是 m 的 14.63 倍。为了提高 Cp的测量精度，可改善0的测量仪表的精度，即提高测量水温的温度计精度，如采用贝克曼温度计，分度值可达0.002，精度为0.001。则其相对误差为 由此可见，变量的精度根本相当。提高0精度后Cp的绝对误差为 Cp=3.5210-30.21.4110-20.02+0.2060.0021.0810-20.05=0.70410-30.28210-3+0.41210-3-0.5410-3=2.9410-4J/(gK)系统提高精度后，Cp的真值为 Cp=0.87980.0003 J/(gK)
