第三次课误差理论月

主要内容：1.系统误差2.误差的合成3.最小二乘原理4.曲线的拟合5.补充：不确定度1一、随机误差的合成按方和根法得到它们的标准误差:

二、系统误差的合成1、恒定系差的合成可按代数和法合成：当误差项数较多时，一般情况下按方和根法合成较好。第九节误差的合成22.变值系差的合成第j个系差的误差区间为[aj,bj]系统不确定度为：ej=1/2*(bj-aj)标准误差为：σj=ej/kj（系统不确定度或极限误差与置信系数之比）合成方法：（1）线性相加法：

e=e1+e2+e3+••••••+en3（2）方和根法：

e=√e12+e22+e32+••••••+en2（3）广义方和根法：将各变值系差的系统不确定度转换成相应的标准误差，用方和根法合成后，得到总的标准误差，再转化为总的系统不确定度。4三、随机误差与系统误差的合成1.线性相加法:g=e+Δ

线性相加的结果，综合不确定度g偏大。2.方和根法：g=√e2+Δ23.广义方和根法5四、最后结果的表示（1）随机不确定度（又称A类不确定度）与系统不确定度（又称B类不确定度）在结果中分别标明。最后结果可表示为：M=（±Δ,±e）Where:M为被测量的测量值或计算结果；e及Δ分别是相应的系统、随机不确定度。（2）用随机不确定度与系统不确定度合成后的综合不确定度表示之。最后结果可写为：M±g

。6例：标准活塞式压力计实验测得各种误差因数引起的压力的极限误差值如下。求总的不确定度（压力P=ma/S，单位均略）。（1）恒定系差：ε=+0.2，由系统安装误差引起。（2）系统不确定度：e1=10.3，是由活塞有效面积S引起的：e2=3.2，来自砝码及活塞质量（m）；e3=0.5，是由重力加速度a的误差引起的。（3）随机不确定度：Δ1=11.6，是由活塞有效面积引起的；Δ2=4.8，是由砝码及活塞质量（m）引起的。7解：设引起误差的各个因数是相互独立的，按照方和根法合成之。总的系统不确定度为：=e=√e12+e22+e32=

√10.32+3.22+0.52=10.8总的随机不确定度为：

Δ=√Δ12+Δ22=√11.62+4.82=12.6故活塞压力计总的不确定度g及修正量c为：g=√e2+Δ2=√10.82+12.62≈16.6≈17C=-ε=-0.28例1：某弹簧管压力表的测量范围为0~1.6MPa，准确度等级为2.5级，校验时出现的最大绝对误差为0.05MPa，问这块仪表是否合格？为什么？9例2：用标准压力表来校验工业压力表时，应如何选用标准压力表精度等级？可否用一台精度等级为0.2，量程为0~25MPa的标准表来校验精度等级为1.5级，量程为0~2.5MPa的压力表，为什么？10第十节最小二乘原理最小二乘原理：欲得真值的最佳估计值，应使各测量值xi的残差vi的平方之和为最小。真值x0的最佳估计值即算术平均值x，具有残差平方和最小值的特性。由于残差均是实数，各个残差的平方必为正数，故残差的平方和为最小值就保证了相应的标准偏差及方差为最小值，同时也说明了测量数据的离散度也是最小的，精度较高。11第十一节曲线拟和的概念

实际工作中由于各随机因素的干扰，所得到的数据往往不同程度存在着误差。因此，插值方法只能适用那些误差可以忽略不记的情况，当误差较大而不能忽略时，又如何通过这些观测数据确定其内在的变化规律呢？本节所介绍的曲线拟合就是解决这一问题的主要方法之一。12如图所示，常常需要从一组获得的数据点中，寻找变量与变量之间的变化规律．用几何方法来解释，就是用已知平面内的一组点，来确定一条曲线，使该曲线能在整体上刻画这组点的变化趋势而不需通过每个点，我们称这种方法为曲线拟合，所求出的曲线称为拟合曲线。1314测量不确定度

重要性与政策：1、测量的质量以不确定度表示2、国家标准及ISO都强调不确定度的重要性3、中国实验室国家认可委员会的测量不确定政策指出：实验室应有能力进行不确定度的评定15误差（Error）准确度（Accuracy）不确定度（Uncertainty）测量结果测得值减真值测量结果与真值的一致程度有定义批量生产仪器没有定义示值最大允许误差没有定义计量标准装置有定义（示值误差）没有定义没有定义但需要给出16定义不确定度定义：与测量结果相联系的参数，表征合理地赋予被测量之值的分散性。此参数可以是用标准差或其倍数，也可以是给定概率下置信区间的半宽。不确定度恒为正值，由方差得出时取其正平方根。17标准不确定度：用标准偏差表示的测量结果的不确定度。（不确定度的）A类评定：对观测列进行统计分析以评定不确定度的方法。（不确定度的）B类评定：评定标准不确定度的非统计分析方法。影响量：不是被测量但对测量有影响的量。

18合成标准不确定度：当结果由若干其它量得来时，该测量结果的标准不确定度等于这些量的方差和协方差加权的正平方根，权的大小取决于这些量的变化及测量结果影响的程度。扩展不确定度：确定测量结果区间的量，期望测量结果以合理地赋予的较高置信水平包含在此区间内。包含因子：为获得扩展不确定度，作为合成不确定度乘数的数字因子。19

测量不确定度的结构:

20◆A类标准不确定度：

（由于随机效应导致的不确定度）

实验标准（偏）差

s是单次观测值Xi的实验标准（偏）差，是m次测量所得算术平均值的实验标准（偏）差，它是分布的标准（偏）差的估计值。为易于区别前者用表示，后者用表示，故：21◆B类标准不确定度：

（由于系统效应导致的不确定度）

不同于A类对观测列进行统计分析的方法来评定标准不确定度，称为不确定度B类的评定，有时也称B类不确定度评定。B类不确定度评定是根据经验和资料及假设的概率分布估计的标准（偏）差表征，也就是说其原始数据并非来自观测列的数据处理，而是基于实验或其他信息来估计，含有主观鉴别的成分。B类不确定度的信息来源一般有：221.以前的观测数据；2.对有关技术资料的测量仪器特性的了解和经验；3.生产企业提供的技术说明文件；4.校准证书（检定证书）或其他文件提供的数据、准确度的等别或级别，包括目前仍在使用的极限误差、最大允许误差等；5.手册或某些资料给出的参考数据及其不确定度；6.规定试验方法的国家标准或类似技术文件中给出的重复性限或复现性。237.测量仪器的示值不够准确；8.标准物质的标准值不够准确；9.引用的数据或其它参量的不够准确；10.取样的代表性不够，即被测样本不能完全代表所定义的被测量；11.化学分析中的基体效应，分析空白，干扰影响，回收率及反映效率等系统影响；12.测量方法和测量程序的近似和假设；13.其他因素……24B类不确定度的评定方法:B类评定的标准不确定度可以用下述四种方法得到：（Ⅰ）若有关资料（如：计量校准/检定证书，仪器说明书等）给出估计值的扩展不确定度为其估计值标准不确定度的倍，则标准不确定度：

◆B类标准不确定度：

（由于系统效应导致的不确定度）

25（Ⅱ）如估计值的扩展不确定度不是按标准不确定度的倍给出，而是给出了置信概率P为90%、95%、99%的置信区间的半宽、、，除非另有说明，一般按照正态分布考虑评定其标准不确定度。对应于上述三种置信概率的包含因子分别为1.64、1.96或2.58，即：标准不确定度为26

表-1:正态分布情况下置信概率与包含因子之间的关系27（Ⅲ）如已知信息表明之值接近正态分布，并以0.68的概率落于的对称范围之内，按表-1，，则。（Ⅳ）若已知估计值分散区间的半宽为，且落在至范围内的概率p为100%，通过对分布的估计，可以得出标准不确定度为：28式中：为测量值概率分布区间半宽度,因为与分布状态有关，见右表：29

合成标准不确定度

式中：为：二个输入量的相关系数为：灵敏系数（当被测量Y为相互独立的输入量的线性函数时，灵敏系数为+1或-1）

30在下列情况下可以不考虑各输入量之间的相关性：(1)可以确认各输入量之间相互独立无关；(2)虽然某些输入量相互之间可能存在相关性，但相关性较弱；（3）虽然某些输入量之间存在强相关性，但这些输入量所对应的不确定度分量对合成不确定度的贡献不大。31◆合成标准不确定度当:;且或-1时,

32结束选定包含因子一般为2~333◆合成标准不确定度用Up

表示扩展不确定度,即合成标准不确定度的自由度称为有效自由度，如果是两个以上多个估计方差分量的合成，即，则即使当每个均服从正态分布的输入量的估计值，变量可以近似为t分布，其有效自由度可以由韦尔奇-萨特思韦特（Welch-Satterthwaite）公式计算：34◆合成标准不确定度上式也可用与相对标准不确定度的合成：35评定不确定度的基本程序

可用下述框图表示

分析引起不确定度的各种因素评定各个不确定度分量A类标准不确定度B类标准不确定度评定合成标准不确定度扩展不确定度报告测量不确定度建立数学模型36建立数学模型

37测量不确定度的评定程序

第一步，按A类评定方法评出A类标准不确定度；第二步，按B类评定方法评出各B类标准不确定度；第三步，按照合成不确定度的合成方法合成标准不确定度；38测量不确定度的评定程序第四步，扩展不确定度的评定方法：（1）按照包含因子k=2~3给出扩展不确定度（2）按照标准不确定度乘以给定概率p的包含因子，从而得到扩展不确定度。与的分布有关。当估计接近正态分布时，采用t分布临界值，，一般采用p值为99%和95%。在多数情况下，采用p=95%。39测量不确定度的评定程序对某项计量标准的检定或校准，根据规定可采用p=99%。当充分大时，可以近似认为、，从而分别得出、。（3）当只给出扩展不确定度U时，不必评定各分量及合成标准不确定度的自由度及有效自由度。除另有说明,对Y可能值的分布作正态分布的估计，虽未计算，但可估计40测量不确定度的评定程序值并不太小时，则大约置信概率近似为95%的区间的半宽，而大约置信概率近似为99%的区间的半宽。（4）如果可以确定Y可能值的分布不是正态分布，而是接近于其他某种分布，则决不能按或计算U或。（5）当Y可能值近似为矩形分布，则包含因子与之间的关系应为：41概率分布情况的估计（参考件）

►正态分布：1)重复条件或复现条件下多次测量的算术平均值的分布；2)被测量量Y用扩展不确定度Up给出，而对其分布又没有特殊指明时，估计值Y的分布；3)被测量Y的合成标准不确定度中，相互独立的分量较多，它们之间的大小也比较接近时，Y的分布；424)被测量Y的合成标准不确定度中，相互独立的分量中，存在两个界限值接近的三角分布，或4个界限值接近的均匀分布时；5)被测量Y的合成标准不确定度的相互独立的分量中，量值较大的分量（起着决定作用的分布）接近正态分布时。43►矩型（均匀）分布：1)数据修约导致的不确定度；2)数字式测量仪器对示值量化（分辨率）导致的不确定度；3)测量仪器由于滞后、摩擦效应导致的不确定