4. 测量不确定度评定框架

在R&S测试系统中,无论是接收机参考电平的校准、频谱仪幅度平坦度的验证,还是信号源输出功率的溯源,最终给出的校准结果都不能只是一个孤立的数值。按照GUM(《测量不确定度表示指南》)与JJF 1059《测量不确定度评定与表示》的要求,完整的测量结果应表述为“被测量的最佳估计值 ± 扩展不确定度”,并说明包含因子与置信概率。本节围绕这一框架,讲清A类评定、B类评定、合成标准不确定度与扩展不确定度的计算逻辑与落地方法。

一、从测量模型出发

不确定度评定的第一步不是套公式,而是写出测量模型。以“用功率计校准信号源在1 GHz处的输出功率”为例,设信号源设定值为 Pset,功率计读数为 Pread,则被校信号源实际输出功率可表示为:

P = P_read + δP_meter + δP_mismatch + δP_drift + δP_repeat

其中 δP_meter 为功率计校准因子引入的修正量,δP_mismatch 为失配修正,δP_drift 为功率计漂移,δP_repeat 为重复性修正。每一个输入量都是不确定度来源,评定工作就是为它们分别赋一个标准不确定度 u(xi),再按传播律合成。

二、A类评定:用重复观测数据说话

A类评定指对在规定条件下重复测量得到的一系列观测值进行统计分析。对同一被测量独立重复测量 n 次,得到 x1, x2, …, xn,则单次测量的实验标准差为:

s(x) = sqrt( Σ(x_i - x̄)² / (n - 1) )

若以 n 次测量的算术平均值作为结果,则平均值的标准不确定度为:

u_A = s(x̄) = s(x) / sqrt(n)

例如对某信号源1 GHz输出功率重复测量10次,得到读数(单位dBm):10.02、10.05、9.98、10.01、10.03、9.99、10.04、10.00、10.02、10.01。计算得平均值约10.015 dBm,单次标准差约0.022 dB,则平均值标准不确定度 uA ≈ 0.022/√10 ≈ 0.007 dB。注意:A类评定得到的是标准不确定度,不是扩展不确定度,且自由度通常取 n−1。

三、B类评定:基于证书、手册与经验

B类评定针对无法重复观测、但可依据校准证书、技术手册、历史数据或经验判断的不确定度来源。核心是把已知的“半宽”a 转换为标准不确定度 uB = a/k,其中 k 取决于分布假设:

  • 正态分布且给出置信概率(如95%),取 k = 1.96;
  • 均匀(矩形)分布,取 k = √3 ≈ 1.732;
  • 三角分布,取 k = √6 ≈ 2.449;
  • 反正弦(U形)分布,取 k = √2 ≈ 1.414。

仍以上例为例:功率计校准证书给出扩展不确定度 U = 0.05 dB(k = 2),则 uB1 = 0.05/2 = 0.025 dB;失配引入的反射系数不确定度按均匀分布估计半宽0.03 dB,则 uB2 = 0.03/√3 ≈ 0.017 dB;功率计年漂移按均匀分布估计半宽0.02 dB,则 uB3 = 0.02/√3 ≈ 0.012 dB。B类评定同样输出标准不确定度,其自由度往往按经验估计,必要时可用JJF 1059给出的公式反算。

四、合成标准不确定度

当各输入量相互独立时,合成标准不确定度为各标准不确定度分量的方和根:

u_c = sqrt( u_A² + u_B1² + u_B2² + u_B3² )

代入上例:uc = √(0.007² + 0.025² + 0.017² + 0.012²) ≈ √(0.000049 + 0.000625 + 0.000289 + 0.000144) ≈ √0.001107 ≈ 0.033 dB。若输入量之间存在相关性,则需引入灵敏系数与协方差项,R&S测试系统中同一台功率计同时提供多个修正量时尤其要注意这一点。

五、扩展不确定度与结果表述

扩展不确定度 U = k · uc,k 通常取2(对应约95%置信概率)或按有效自由度查 t 分布表。上例取 k = 2,则 U = 2 × 0.033 ≈ 0.066 dB。最终校准结果可表述为:信号源1 GHz输出功率为10.015 dBm,扩展不确定度 U = 0.066 dB(k = 2,约95%置信概率)。

JJF 1059强调:不确定度评定报告应包含测量模型、各分量来源与评定方法、uc、U、k 及置信概率,缺一不可。只给 U 而不给 k 是不完整的。
GUM / JJF 1059 不确定度评定框架 测量模型:P = P_read + δP_meter + δP_mismatch + δP_drift + δP_repeat 识别全部不确定度来源,明确输入量与输出量的函数关系 A类评定(统计方法) 重复测量 n 次 → 算术平均值 x̄ 单次标准差 s(x) = √[Σ(xᵢ−x̄)²/(n−1)] 平均值标准不确定度 u_A = s(x)/√n 示例:u_A ≈ 0.007 dB(n=10) B类评定(非统计方法) 依据证书/手册/经验 → 半宽 a u_B = a / k(k 由分布决定) 正态 k=1.96|均匀 k=√3|三角 k=√6 示例:u_B1=0.025,u_B2=0.017,u_B3=0.012 dB 合成标准不确定度 u_c = √(u_A² + u_B1² + u_B2² + u_B3²) 示例:u_c ≈ 0.033 dB(各分量独立) 扩展不确定度 U = k·u_c,k=2 → U ≈ 0.066 dB(约95%)

六、评定过程中的常见判断

实际执行时,有几个判断点直接影响结果。第一,A类与B类不是按“重要程度”划分,而是按“是否用统计方法获得”划分,同一来源在不同场合可能归入不同类别。第二,B类评定的分布假设要有依据:校准证书通常按正态处理,仪器分辨力、失配等按均匀分布处理更稳妥。第三,当某分量明显小于其他分量的1/3时,可评估其是否可忽略,但应在报告中说明。第四,合成前应检查各分量是否已统一到同一单位与同一量纲,R&S测试系统中dB与线性值混用是常见错误来源。

不要把A类评定得到的标准差直接当作扩展不确定度使用,也不要在未说明 k 的情况下给出 U。JJF 1059要求不确定度表述具备可追溯性与可复现性。

七、可直接套用的评定记录表

为便于在R&S测试系统校准中落地,下表给出一个可填写的不确定度分量汇总示例,字段与本节示例对应,实际使用时按被测对象替换来源与数值即可。

序号不确定度来源类别半宽 a 或 s(x)分布 / k标准不确定度 uᵢ备注
1重复性As(x)=0.022 dBn=100.007 dB平均值
2功率计校准因子BU=0.05 dB正态 k=20.025 dB证书
3失配B0.03 dB均匀 k=√30.017 dB反射系数估算
4功率计漂移B0.02 dB均匀 k=√30.012 dB年漂移
合成标准不确定度 u_c0.033 dB
扩展不确定度 U(k=2)0.066 dB

把这张表与前面的测量模型、A/B类评定方法、合成与扩展公式连起来看,就是GUM与JJF 1059在R&S测试系统校准中的完整评定链条:先建模,再分源,后合成,最后按约定包含因子给出结果。掌握这条链条,才能让校准报告中的每一个数值都站得住脚。