《倪育才 不确定度》由会员分享，可在线阅读，更多相关《倪育才 不确定度（22页珍藏版）》请在金锄头文库上搜索。
1、如何建立数学模型讲授人：中国计量科学研究院研究员 倪育才应中国计量编辑部的邀请，根据近几年来在各种“测量不确定度评定”学习班上讲课的经验 和对初学者经常感到困惑的问题或不容易理解的问题，现写若干短篇并在中国计量上连载。为使这 些短篇相互独立，每一篇讲述一个问题，合起来看又有一定的系统性。经与编辑部商讨，初步确定下述 几个问题:1 如何建立数学模型2 极差法和贝塞尔法之间的比较3 被测量 Y 可能值分布的判定4 包含因子 k 的选择5 测量不确定度评定在不同应用中的差别6 测量误差的基本概念7 测量不确定度的基本概念8 测量误差和测量不确定度的差别在测量不确定度评定中，建立数学模型也称为测量模型
2、化，目的是要建立满足测量不确定度评定所 要求的数学模型，即建立被测量Y和所有各影响量X间的函数关系，其一般形式可写为：Y=f（X1， X2， Xn）可以说，建立数学模型是进行测量不确定度评定最关键的第一步，也是许多初学者在进行测量不确 定度评定时遇到的第一个困难。测量不确定度表示指南（GUM）在摘要介绍测量不确定度评定步骤时，首先就提到要建立数学模 型，并说 ：“ The function f should contain every quantity, including all corrections and correction factors, that can contribute
3、a significant component of uncertainty to the result of measurement. ”。其意是数学模型f中应包含所有对测量结果的不确定度有影响的修正值和修正因子。 也就是说，数学模型中应包含所有应该考虑的影响量，而每一个影响量将对测量结果贡献一个值得考虑 的不确定度分量。因此一个好的数学模型，其中所包含的影响量和此后不确定度评定中所考虑的每一个 不确定度分量应该是一一对应的。这样建立起来的数学模型，既能用来计算测量结果，又能用来全面地 评定测量结果的不确定度。要找出每一个影响量与被测量之间的函数关系，往往是很困难的，有时简直不可能得到两者关
4、系的 解析表达式。于是许多初学者往往将测量中用来获得被测量的计算公式作为数学模型而列出。例如在各 种测量中，最经常采用的方法之一是比较测量。将被测量值y和参考标准所提供的标准量值s相比较， 通过测量两者之差可以计算出被测量y。于是在已经发表的各种测量不确定度评定的文章中，经常见 到将y=x+4作为数学模型的情况。但在进行不确定度评定时，则又往往脱离数学模型而重新考虑各个不 确定度分量。这样的数学模型对测量不确定度评定实际上毫无帮助。在某些特殊情况下（例如某些检测项目）将计算公式作为数学模型可能是允许的，但一般说来不要把 数学模型简单地理解为就是计算测量结果的公式，也不要理解为就是测量的基本原理
5、公式。两者之间经 常是有区别的。从原则上说，似乎所有对测量结果有影响的输入量都应该在计算公式中出现，但实际情况却不然。 有些输入量虽然对测量结果有影响，但由于信息量的缺乏，在具体测量时无法定量地计算它们对测量结 果的影响。也有些输入量由于对测量结果的影响很小而被忽略，故在测量结果的计算公式中也不出现， 但它们对测量结果的不确定度的影响却可能是必须考虑的。因此如果仅从计算公式出发来进行不确定度 评定，则上述这些不确定度分量就可能被遗漏。当然，在某些特殊情况下如果所有其他不确定度贡献因 素的影响都可以忽略不计时，数学模型也可能与计算公式相同。对于不同的被测量和不同的测量方法，数学模型的具体形式可能
6、差别很大，但实际上都可以用一种 比较系统的方式来给出数学模型，或者说可以给出数学模型的通式。根据测量误差的定义:误差=测量结果-真值。同时误差又可以分为随机误差和系统误差两类，且三 者之间的关系为:误差=系统误差+随机误差。于是可以得到:真值=测量结果-误差=测量结果-系统误差-随机误差由于修正值等于负的误差，于是上面的关系式就成为: 真值=测量结果-系统误差-随机误差=测量结果+系统误差的修正值+随机误差的修正值实际上，真值就是想得到的被测量的测量结果，于是上式可写成被测量=测量结果+系统误差的修正值+随机误差的修正值例1:对于常见的量块比较测量，若l为标准量块的长度， l为测得的两量块的长
7、度差，于是被测s量块长度l的计算公式为：Xl =1 + lxs由于测量时量块的温度通常会偏离标准参考温度20C，考虑到温度和线膨胀系数对测量结果的影 响，计算公式成为:l =l + l+l 6 a 0 +l a 6 Qx s s x s s式中a和Q分别表示线膨胀系数和对标准参考温度20C的偏差；脚标“”、“ ”分别表示标准量sx块和被测量块；以及6 Q =Q -Q 和6 a =a -a 。s x s x考虑到量块测量点可能偏离量块测量面中心点对测量结果的影响，数学模型成为:l =l + l+l 6 a Q +l a 6 Q +6 lx s s x s s将此数学模型和上面给出的通式相比较就可
8、以发现，等式右边的第一、二项l + l即是由测量得到s的未修正的测量结果。等式右边的第三、四项l 6aQ +l a 6 Q是对由温度偏差所引入的系统误差的s x s s修正值，在本例中这两项的数值十分小而可以忽略，但它们对测量结果不确定度的影响是必须考虑的。 等式右边的最后一项6 l,是表示由于测量点可能偏离量块中心对测量结果的影响。测量点的偏离对测 量结果引入随机误差，因此最后一项实际上是对该随机误差的修正值。由下图可见两者之间的对应关系。鮫测駅=测Rt结集+霍统谋歪的搐正慎+R&机谋產的修正值1 t:I4=+辭 + Si例 2:砝码校准，将被测砝码的质量与具有相同标称值的标准砝码相比较。若
9、被校准砝码和标准砝码的折算质量分别为m和m ,测得两者的质量差为卜m，于是被校准砝码折算质量 m的计算公式xsx为：m =m +A mxs考虑到标准砝码的质量自最近一次校准以来可能产生的漂移 m ,质量比较仪的偏心度和磁效应的d影响 m，以及空气浮力对测量结果的影响6 B后，其数学模型成为：cm =m +A m+6 m +6 m +6 Bx s d c模型中等式右边的第一、二项为未修正的测量结果。该测量不存在值得考虑的系统误差，也就是说， 在数学模型中不存在对系统误差的修正值。等式右边的第三、四、五项为对三项随机误差分量的修正量。铁灣二测结乐+垂统误爰的修正值+随机误莖的修正侑t :：:与数学
10、模型通式之间的对应关系为： J 2宀禺在建立数学模型时，未修正的测量结果和系统误差的修正值通常都能比较容易地得到解析形式的数 学表达式。惟有随机误差的修正值无法得到其解析形式的表达式。因此只能在数学模型中简单地加上一 项，表示对随机误差的修正值。根据随机误差的定义，无限多次测量结果的随机误差的平均值等于零， 因此这些项的数学期望为零。也就是说，增加这些修正值后不会对被测量的数值有影响。需要知道的是 这些修正值的可能取值范围，通常可以由测量者的经验或辅助的实验测量得到。再由假定的概率分布， 可以通过 B 类评定估算出它们的标准不确定度。有些测量，其计算公式中可能仅包含各影响量的积和商，即被测量可
11、以用下述函数形式表示:式中的系数c并非灵敏系数，而是比例常数，且指数Pi可以为正数或负数。在这种情况下，需要增 加的不是修正值，而是相乘的修正因子。此时，数学模型的通式可以表示为:被测量等于未修正测量结 果的计算公式乘以由于系统误差引入的修正因子（它们的数学期望值不等于 1），再乘以由于随机误差引 入的修正因子（它们的数学期望值等于1）。有些领域，例如化学分析领域，经常出现这种类型的数学模型。例 3:在用原子吸收光谱法测定陶瓷容器中镉的溶出量的实例中，被测量为被醋酸溶液浸泡的容器单位表面积镉的溶出量r,它可以表示为:式中：P0醋酸浸取液中镉的质量浓度；VL醋酸浸取液体积；av被醋酸溶液浸泡的容
12、器表 面积。考虑到还有三项随机误差在上述公式中未反映出来，它们分别是浸泡温度、浸泡时间和醋酸的体积分 数对测量结果的影响，于是最后采用的数学模型成为：_ 在该数学模型中，S是未修正的测量结果，t 、和f d分别是相对于三项随机误差的修正因子，temp time acid它们的数学期望均等于1。在本例中不存在值得考虑的系统误差。由此可见，写出符合要求的数学模型并不难，关键还是要找到所有能影响测量结果的误差来源。一般 先根据测量的最基本原理导出被测量的基本计算公式，然后考察该计算公式是否已经对所有的系统误差 进行了修正，否则就补充加入其余未考虑的系统误差分量的修正值（或乘以修正因子），最后再加上对
13、所 有随机误差分量的修正值（或乘以修正因子）。只要对测量工作有一定程度的了解，写出计算公式和系统 误差修正值的函数形式对大部分测量人员并不困难，因此要做的仅是简单地将所有需要考虑的随机误差 的修正值（或修正因子）补充进入数学模型。必须注意，即使对于相同的被测量和相同的测量方法，数学模型也不是一成不变的。随着所选择的 影响量的不同，对测量不确定度评定所要求的严密程度的不同，其数学模型也可能会有所不同。此外，对于测量仪器和量具的常规检定或校准来说，还必须注意两者在数学模型上可能存在的微小差 别。当被测对象是测量仪器时，由于仪器本身一般不提供标准量值，其量值需要用其他测量标准进行标 定。故在进行测量
14、不确定度评定时，被测量应该是测量仪器的示值误差E,因此其数学模型需写成示值X误差的形式，即“E =”。当被测对象是实物量具时，由于实物量具本身能提供一个标准量值，故在X进行测量不确定度评定时，被测量既可以是其相对于标称值的偏差（相当于示值误差），也可以是它所提 供的量值。也就是说，其数学模型既可以写成“y=”的形式，也可以写成“E =”。由于两者之 丿X间仅相差一个标称值，而标称值是一个规定值而不存在不确定度，因此两种数学模型在不确定度评定时 毫无差别。极差法和贝塞尔法之间的比较讲授人：中国计量科学研究院研究员 倪育才标准不确定度的 A 类评定定义为:“用对观测列进行统计分析的方法，来评定标准
15、不确定度”。国 家计量技术规范JJF1059-1999测量不确定度评定与表示中介绍了两种A类评定的方法，贝塞尔法和 极差法。1. 贝塞尔法当在重复性或复现性条件下，对被测量X进行n次独立观测。若得到的测量结果分别为气，x2,， x，n次测量的平均值为#。于是用贝塞尔公式可以求出单次测量结果x.的实验方差S2（x.）和实验标准 nii差 S（x） 。.2. 极差法当在重复性或复现性条件下，对被测量X进行n次独立观测。若n个测量结果中最大值和最小值之差为R（称为极差），在可以估计X接近正态分布的条件下，单次测量结果的实验标准差s（x.v）可近似.地表示为:s(x.)=R/C=u(x.).式中系数C为极差系数。极差系数之值与测量次数n的大小有关。表1给出极差法的极差系数和自由度与测量次数的关系。1及自由度裂- - 1n2356789C1.13L692 062. 332. 532.702.852.97V0. 91.82. 73, 64.5 5.36-06.8既然随机变量 X 的标准偏差可以用两种方法得到，就不可避免地会提出两种
