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 3446x系列前面板“Statistics”功能可显示平均值、标准差、峰峰值和样本数;通过SCPI的CALC:AVER:AVER?、CALC:AVER:SDEV?可读取。注意区分总体标准差(分母n)与样本标准差(分母n−1),仪器菜单中通常可切换。

二、概率分布:不确定度评定的“形状”假设

标准差只描述离散程度,要给出置信概率还必须知道分布形状。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,正态分布按置信概率查包含因子。

Keysight测量数据 → 统计量 → 不确定度分量的映射关系 原始测量数据 34465A 重复读数 x₁, x₂, …, xₙ CSV / SCPI 读取 分辨力 1 µV NPLC = 100 n = 100 次 统计量计算 x̄ = Σxᵢ / n s² = Σ(xᵢ−x̄)²/(n−1) s = √s² u_A = s / √n r = cov(xᵢ,xⱼ)/(sᵢsⱼ) NumPy: ddof=1 不确定度分量 A类:u_A = s/√n(正态) B类:u_B = a/√3(均匀) B类:u_B = a/√6(三角) 协方差项:2·cᵢ·cⱼ·u(xᵢ,xⱼ) 合成:u_c² = Σ(cᵢuᵢ)² + 2ΣΣ… k=2 → 约95.45% 概率分布形状与包含因子 k 的对应关系 正态分布 k=2 → 95.45% k=3 → 99.73% 均匀(矩形)分布 u = a / √3 分辨力、年准确度 三角分布 u = a / √6 两均匀量之和 相关输入量必须计入协方差项;忽略正相关会低估 u_c,忽略负相关会高估 u_c

五、一个可复算的示例

设用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不确定度评定中不可回避的原因。

常见错误:把仪器显示的标准差直接当作平均值的标准不确定度;把均匀分布的半宽a直接当作标准不确定度;对明显相关的输入量按独立处理。三者都会导致合成不确定度偏离合理值。