2. 测量不确定度的数学基础
在Keysight测试系统中,不确定度评定最终要落到数字上:一个测量结果、一个合成标准不确定度、一个包含因子和一段置信区间。支撑这些数字的数学骨架,是标准差、方差、概率分布、协方差与相关系数。本节不重复概率论教材的推导,而是把这些概念放回Keysight测量数据的语境中,说明它们如何被表达、如何被计算、以及如何影响最终的不确定度结果。
一、从测量数据到统计量:均值、方差与标准差
对同一被测件在重复性条件下做n次测量,得到一组观测值x₁, x₂, …, xₙ。Keysight仪器(如34465A数字万用表、53230A计数器)在多次触发采集模式下会直接给出这组读数,或通过SCPI命令读取。最常用的统计量是算术平均值:
x̄ = (1/n) · Σ xᵢ
它作为被测量的最佳估计值。方差s²衡量数据围绕均值的离散程度,样本方差用n−1作分母(贝塞尔校正),因为均值本身由同一组数据估计,损失了一个自由度:
s² = Σ(xᵢ − x̄)² / (n − 1)
s = √s²
标准差s与测量值同量纲,是A类标准不确定度u_A的直接来源:当以n次测量的平均值作为结果时,u_A = s/√n。这一点在Keysight数据中很关键——仪器显示的“标准差”通常是单次测量的s,而平均值的不确定度要再除以√n。例如用34465A在10 V量程、100 NPLC下重复测量100次,读数标准差s = 2.5 µV,则平均值的标准不确定度为0.25 µV。若误把2.5 µV当作平均值的不确定度,会高估10倍。
二、概率分布:不确定度评定的“形状”假设
标准差只描述离散程度,要给出置信概率还必须知道分布形状。Keysight不确定度评定中最常用三类分布:
- 正态分布:重复性测量、随机噪声主导的量,其平均值近似正态。包含因子k=2对应约95.45%置信概率,k=3对应约99.73%。
- 均匀(矩形)分布:仅知某误差限±a而不知其分布时采用,标准不确定度u = a/√3。例如仪器年准确度指标、分辨力引入的不确定度。
- 三角分布:两个均匀分布量之和,或中心值最可能的情形,u = a/√6。
分辨力是最典型的均匀分布来源。若34465A在1 V量程显示分辨力为1 µV,则半宽a = 0.5 µV,分辨力标准不确定度u = 0.5/√3 ≈ 0.289 µV。这个值虽小,但在低电平测量中常与噪声不确定度同量级,不可忽略。
三、协方差与相关系数:输入量之间的“联动”
当输出量y = f(x₁, x₂, …, x_N)的多个输入量彼此不独立时,合成方差中必须加入协方差项:
u_c²(y) = Σ (∂f/∂xᵢ)² u²(xᵢ) + 2 ΣΣ (∂f/∂xᵢ)(∂f/∂xⱼ) u(xᵢ,xⱼ)
相关系数r(xᵢ,xⱼ) = u(xᵢ,xⱼ) / [u(xᵢ)·u(xⱼ)],取值−1到+1。r=0表示不相关,r=±1表示完全线性相关。在Keysight系统中,相关性常来自:同一台仪器校准证书给出的多个修正值(共享同一校准源)、同一环境温度同时影响多个通道、同一参考标准溯源链。忽略正相关会低估合成不确定度,忽略负相关则可能高估。
处理相关的实用做法有两种:一是把相关输入量合并为一个共同分量重新建模;二是显式估计协方差。例如用同一台3458A校准两台被测表,两台表的修正值误差高度相关,此时不能简单按独立处理。
四、Keysight数据中的表达与计算工具
Keysight仪器与软件对上述统计量的表达方式各有约定。34465A/34470A的统计功能给出平均值、标准差、峰峰值;Truevolt系列还提供直方图和趋势图,可直观判断分布是否近似正态。BenchVue、PathWave BenchVue和Keysight 89600 VSA软件可导出CSV,用Python/NumPy计算:
import numpy as np
x = np.loadtxt("dmm_readings.csv", delimiter=",")
mean = np.mean(x)
s = np.std(x, ddof=1) # 样本标准差
uA = s / np.sqrt(len(x)) # 平均值的A类不确定度
cov = np.cov(x1, x2) # 协方差矩阵
r = np.corrcoef(x1, x2) # 相关系数矩阵
注意NumPy默认ddof=0为总体标准差,做不确定度评定时应设ddof=1。对于B类分量,则按分布假设换算:均匀分布除以√3,三角分布除以√6,正态分布按置信概率查包含因子。
五、一个可复算的示例
设用34465A测量某10 V参考源,重复测量n=25次,得s = 3.0 µV,则u_A = 3.0/√25 = 0.6 µV。仪器1 V量程分辨力1 µV,均匀分布u_res = 0.5/√3 ≈ 0.289 µV。校准证书给出修正值不确定度U = 4 µV(k=2),则u_cal = 2 µV。若这三项互不相关,合成标准不确定度:
u_c = √(0.6² + 0.289² + 2.0²) ≈ √(0.36 + 0.0835 + 4.0) ≈ 2.11 µV
取k=2,扩展不确定度U ≈ 4.22 µV。若校准修正值与另一台同源仪器的修正值存在r=0.8的相关,则需在合成中增加协方差项,结果会明显不同。这正是协方差与相关系数在Keysight不确定度评定中不可回避的原因。