3. 测量模型建立实战

测量不确定度评定的第一步,也是最容易被轻视的一步,就是建立测量模型。所谓测量模型,就是把一个真实的物理测量问题,翻译成可计算的数学函数关系:输出量 Y 是被测量,输入量 X1、X2、…、XN 是影响 Y 的所有量,模型写成 Y = f(X1, X2, …, XN)。后续所有 A 类、B 类评定、合成标准不确定度、扩展不确定度,全部建立在这个函数关系之上。模型错了,后面的计算再精确也没有意义。

一、从物理问题到函数关系:先写“测量原理”,再写“修正项”

建立模型不要一上来就写公式,而应先回答三个问题:被测量到底是什么?它是通过什么物理原理得到的?还有哪些量会系统性地影响这个结果?

以“用游标卡尺测量圆柱体直径”为例。被测量是直径 d,测量原理是卡尺两爪直接接触圆柱面读数。若只考虑读数,模型可写成 d = d读。但真实测量中,卡尺存在零位误差 δ零,测量时温度偏离参考温度 20 ℃ 会带来热膨胀影响,圆柱体与卡尺线膨胀系数不同,修正量约为 d·(α件 − α尺)·(t − 20)。于是更完整的模型为:

d = d读 + δ零 + d·(α件 − α尺)·(t − 20)

注意这里 d 同时出现在等式两边,严格说这是一个隐式模型。工程上常把修正项中的 d 用 d读 近似,得到显式模型:

d = d读 + δ零 + d读·(α件 − α尺)·(t − 20)

这一步“近似”本身要记录在模型说明中,因为它决定了后续灵敏度系数怎么求。

二、输入量识别:把“影响量”与“被测量”分开

输入量识别最容易犯的错误,是把“被测量的重复读数”当成多个输入量。实际上,重复测量得到的一组读数,属于同一个输入量 d读 的多次观测,其分散性用 A 类评定处理,而不是拆成 d读1、d读2。

识别输入量时,建议按“直接测量量—修正量—影响量—常数”四类逐一排查:

  • 直接测量量:卡尺读数 d读,由多次测量取平均得到;
  • 修正量:零位误差 δ零,由校准证书给出,属于 B 类;
  • 影响量:温度 t、线膨胀系数 α件、α尺,它们不直接等于被测量,但通过物理规律影响结果;
  • 常数:参考温度 20 ℃,按定义取固定值,不引入不确定度。

对于本例,若温度控制在 (20 ± 2) ℃,α件 = 11.5×10−6 ℃−1,α尺 = 16.5×10−6 ℃−1,直径约 25 mm,则温度修正项最大约为 25 × (−5×10−6) × 2 = −0.00025 mm,即 0.25 μm。这个量级是否可忽略,取决于被测直径的公差要求。若公差为 ±0.01 mm,则 0.25 μm 约占 2.5%,通常不能忽略;若公差为 ±0.1 mm,则可考虑忽略并说明理由。

判断某个影响量能否忽略,不要凭感觉,而应估算其最大可能贡献,与目标不确定度或公差比较。一般当某项贡献小于合成标准不确定度的 1/3 时,可考虑忽略,但需在报告中说明。

三、输出量定义:明确“测的是什么”和“结果给谁用”

输出量定义要落到具体物理量和具体位置。例如“测量圆柱体直径”应明确:是测量某一横截面的直径,还是沿轴向多个截面直径的平均值?如果是后者,输出量应定义为 d = (1/n)Σdi,输入量包括每个截面的直径,模型结构随之改变。

再如“测量电阻”时,若用伏安法,输出量是 R = U/I;若用四线制开尔文测量,输出量仍是 R,但输入量中引线电阻的影响被大幅削弱,模型可简化为 R = U/I 加上接触电势修正。输出量定义不同,输入量清单和模型形式都会不同。

四、模型结构可视化:输入量如何汇入输出量

下面这张图把本例的测量模型画成“输入量—函数关系—输出量”的结构关系图。图中每个输入量都标注了类型和典型来源,便于后续逐项评定。

测量模型结构:圆柱体直径 d 的输入量汇入关系 输入量 Xᵢ d读 直接测量量 · A类 δ零 零位修正 · B类 t 温度 · B类 α件 线膨胀系数 · B类 α尺 线膨胀系数 · B类 测量函数 f d = d读 + δ零 + d读·(α件 − α尺) ·(t − 20) 参考温度 20 ℃ 为定义常数,不评不确定度 修正项中 d 用 d读 近似 输出量 Y 圆柱体直径 d 被测量 · 最终报告量 单位:mm 输入量识别原则:直接测量量、修正量、影响量、常数四类逐一排查;重复读数属于同一输入量的多次观测,不拆分为多个输入量。

五、实战检查清单:模型是否可用

模型写完后,用下面几个问题自检:

检查项 判断标准
输出量是否唯一、明确 能说清测的是什么物理量、在什么位置、什么条件下
输入量是否覆盖主要影响 每个输入量都能对应到具体测量步骤或物理效应
函数关系是否可计算 给定输入量数值,能算出输出量估计值
修正项是否显式表达 系统误差修正已写入模型,而非事后“加一点”
忽略项是否有依据 被忽略的影响量已估算最大贡献并说明理由

以本例为例,若最终模型为 d = d读 + δ零 + d读·(α件 − α尺)·(t − 20),输入量共 5 个:d读、δ零、t、α件、α尺。其中 d读 用 10 次重复测量平均,A 类标准不确定度按实验标准差除以 √10 计算;其余 4 个按 B 类评定,依据校准证书、温度计精度和材料手册数据。这样,一个物理测量问题就完整地转化成了可继续评定的数学模型。

常见错误:把“卡尺分辨力”直接当成输入量。分辨力是读数 d读 的不确定度来源之一,应通过 B 类评定计入 d读 的标准不确定度,而不是单独列一个输入量。输入量是“量”,不是“误差源”。

建立测量模型不是一次性的形式工作,而是贯穿整个评定过程的“骨架”。后续每引入一项新的不确定度来源,都应先问:它对应哪个输入量?是否已经出现在模型中?如果不在,是补充进模型,还是证明它可以忽略。保持模型与评定过程一致,才能保证最终给出的扩展不确定度真实可信。