Skip to content

16.4.1 测量误差及其分布 ​

16.4.1.1 测量误差的定性特征 ​

欲根据起因描述测量误差, 需要区分下述三类误差:

(1)粗误差由读数不准确或错误引起, 这类误差是可排除的.

(2) 系统测量误差由测量设备的不精确或测量方法引起, 包括读取数据的方法, 以及测量系统的测量误差, 这类误差并非都可以避免.

(3)统计误差或随机测量误差源于不易或不可能控制的测量条件的随机变化, 也可能由观察事件的某种随机性产生.

在测量误差理论中, 通常假定粗误差和系统测量误差是可排除的, 故只有统计性质和随机测量误差被纳入到误差计算中.

16.4.1.2 测量误差的密度函数 ​

1. 测量规约 ​

为计算不确定性特征, 必须假设把测量结果作为主要记法列入测量记录中, 从而可得到不确定数值的频率或密度函数 f(x) 或累计频率或分布函数 F(x) . 所考虑的随机变量 X 的实现用 x 表示.

2. 误差密度函数 ​

对于测量误差性质的特别假定, 导致了误差分布的密度函数具有一些特殊性质.

(1)连续密度函数 由于随机测量误差可在某区间内取任意值, 故可用连续密度函数 f(x) 描述.

(2) 偶密度函数 若绝对值相同但符号不同的测量误差是等可能的, 则密度函数是偶函数: f(−x)=f(x) .

(3) 单调递减密度函数 若具有较大绝对值的测量误差与具有较小绝对值的测量误差相比可能性小,则当 x>0 时,密度函数 f(x) 是单调递减的.

(4) 期望有限 误差绝对值的期望一定是有限数:

(16.189)E(|X|)=∫−∞∞|x|f(x)dx<∞.

误差的不同性质导致密度函数的类型各异.

3. 正态分布误差 ​

(1) 密度函数和分布函数 在大多数实际情况下, 可假设测量误差的分布是期望 μ=0 、方差为 σ2 的正态分布,即测量误差的密度函数 f(x) 和分布函数 F(x) 是

(16.190a)f(x)=1σ2πe−x22σ2

和

(16.190b)F(x)=1σ2π∫−∞xe−t22σ2dt=Φ(xσ),

其中 Φ(x) 是标准正态分布的分布函数 (参见第 1070 页 (16.73a) 和 1458 页表 21.17). 在 (16.190a, 16.190b) 情形下, 也称为正态误差

(2) 几何表示 密度函数 (16.190a) 如图 16.18(a) 所示, 图中标出了拐点和重心. 图 16.18(b) 给出了不同方差对应的密度函数. 拐点的横坐标为 ±σ ,半区域的重心在 ±η 处. 当 x=0 时,函数取得最大值 1/(σ2π) . 当 σ2 增大时,曲线变宽, 曲线下的面积恒等于 1. 分布图显示, 小误差经常发生, 大误差很少出现.

0193686d-91c3-7222-a100-f59d7e5e597d_105_364_1289_910_358_0.jpg

4. 描述正态分布误差的参数 ​

除了方差 σ2 ,或者也被称为均方误差或标准差的标准偏差 σ ,还有其他一些参数可用于描述正态分布误差,如精确度 h 、平均误差 η 和可能误差 γ .

(1) 精确度 除了方差 σ2 ,精确度

(16.191)h=1σ2

可用来描述正态分布的宽度. 高斯曲线越窄, 精确度越高 (图 16.18(b)). 利用试验值 σ~ 或由测量值得到的 σ~x 取代 σ ,精确度可用来说明测量方法的准确性.

(2)平均误差 误差绝对值的期望 η 定义为

(16.192)η=E(|X|)=2∫0∞xf(x)dx=2πσ.

(3)可能误差 误差绝对值的边界 γ 满足性质

(16.193a)P(|X|≤γ)=12

称为可能误差. 这意味着

(16.193b)∫−γγf(x)dx=2Φ(γσ)−1=2π∫0γ/σe−t2/2dt=12,

其中 Φ(x) 是标准正态分布的分布函数. 条件 (16.193b) 是非线性方程,可借助于计算机代数系统近似求解,以得到 γ/σ . 其结果是

(16.193c)γσ≈0.6745

(4) 给定的误差界 若已知误差的上界 a>0 ,则由 (16.193b) 可计算出误差落在区间 [−a,a] 内的概率是

(16.194)P(|X|≤a)=2Φ(aσ)−1.

(5) 标准差、平均误差、可能误差和精确度之间的关系 若误差服从正态分布, 则根据 (16.193c), 下述关系式成立:

(16.195a)η=2πσ(16.195b)γ=0.6745σ,(16.195c)h=12σ

16.4.1.3 测量误差的量化特征 ​

1. 真实值及其近似值 ​

一个可测量的真实值 xw 通常是未知的. 故通过测量值 xi(i=1,2,⋯,n) 实现的随机变量期望,可作为 xw 的估计值,从而可把下述均值视为 xw 的近似值.

(1) 算术平均

(16.196a)x¯=1n∑i=1n⁡xi.

如果测量值以绝对频率 hj 分成 k 组,组平均值为 x¯j(j=1,2,⋯,k) ,则

(16.196b)x¯=∑j=1k⁡hjx¯j.

(2) 加权平均

(16.197)x¯(g)=∑i=1n⁡gixi∑i=1n⁡gi.

此时,对单次测量值以权重因子 gi(gi>0) 进行加权 (参见第 1112 页 16.4.1.6,1.).

2. 测量序列中的单次测量误差 ​

(1) 测量序列中单次测量的真实误差 是真实值 xw 和测量值之差. 由于真实值通常是未知的,测量结果为 xi 的第 i 次测量之真实误差 εi 也是未知的:

(16.198)εi=xw−xi.

(2) 测量序列中单次测量的平均误差或表观误差 是算术平均值和测量值 xi 之差:

(16.199)vi=x¯−xi.

(3) 单次测量的均方误差或标准差 由于 n 次测量真实误差 εi 之和的期望值, 以及表观误差 vi 之和的期望值为 0 (与它们多大无关),故用误差平方和表示:

(16.200a)ε2=∑i=1n⁡εi2(16.200b)v2=∑i=1n⁡vi2.

从实用观点看,既然只有 vi 的值可由测量过程确定,我们仅对 (16.200b) 感兴趣. 因此, 测量序列中单次测量的均方误差可定义为

(16.200c)σ~=∑i=1n⁡vi2/(n−1).

σ~ 是误差分布标准差 σ 的近似值.

在正态误差分布情形下,可得到 σ~=σ :

(16.200d)P(|ε|≤σ~)=2Φ(1)−1=0.6826,

即真实值的绝对值不超过 σ 的概率约为 68% .

(4) 单次测量的可能误差 是数 γ ,由于

(16.201a)P(|ε|≤γ)=12,

即误差绝对值不超过 γ 的概率是 50% . 横坐标 ±γ 把密度函数下方左右两边的区域分成相等的两部分 (图 16.18(a)).

在正态误差分布情形下, 根据 (16.193c), 有

(16.201b)γ≈0.6745σ≈0.6745σ~=γ~.

(5) 单次测量的平均误差 是数 η ,它是误差绝对值的期望:

(16.202a)η=E(|ε|)=∫−∞∞|x|f(x)dx.

在正态误差分布情形下,有 η=2πσ 和

(16.202b)P(|ε|≤η)=2Φ(ησ)−1=2Φ(2π)−1≈0.5751.

误差不超过 η 值的概率约为 57.6% . 密度函数下方左右两侧区域重心的横坐标是 ±η (图 16.18(a)). 在正态误差分布情形下:

(16.202c)η=2πσ≈2πσ~=η~.

3. 测量序列算术平均值的误差 ​

测量序列算术平均值 x¯ 的误差由单次测量误差给出.

(1)均方误差或标准差

(16.203)σ~AM=∑i=1n⁡vi2/[n(n−1)]=σ~n.

(2) 可能误差

(16.204)γ~AM=0.6745∑i=1n⁡vi2/[n(n−1)]=0.6745⋅σ~n.

(3) 平均误差

(16.205)η~AM=0.7979∑i=1n⁡vi2/[n(n−1)]=0.7979⋅σ~n.

(4) 误差的可达水平 由上述 (16.203) (16.205) 定义的三类误差, 与对应的单次测量误差 (16.200c),(16.201b),(16.202c) 成比例,且与 n 的平方根的倒数成比例. 当达到一定次数后, 再增加测量次数则没有实际意义. 更有效的方法是提高测量方法的精确度 h(16.191) .

4. 绝对误差和相对误差 ​

(1) 绝对不确定性、绝对误差 测量结果的不确定性用误差 εi,vi,σi,γi,ηi 或 ε,v,σ,γ,η 来描述,它们度量了测量的可靠性. 由绝对误差给出的绝对不确定性概念, 对所有类型的误差和误差传播的计算都有意义 (参见第 1114 页 16.4.2). 作为测量数量, 它们的维度相同.

引入 “绝对” 误差这个词汇是为了避免与相对误差概念混淆,经常使用记号 Δxi 或 Δx 表示. "绝对" 与绝对值的概念完全不同: 它指可测量没有符号限制的数值 (比如长度、重量、能量).

(2)相对不确定性、相对误差 用相对误差给出的相对不确定性, 是关于可测量数值测量方法品质的度量. 由于相对误差是绝对误差和可测量数值之商, 故与绝对误差相反,相对误差没有维度. 若测量值未知,可用量 x 的平均值代替:

(16.206a)δxi=Δxix≈Δxix¯.

相对误差大多数用百分数给出, 也称为百分误差:

(16.206b)δxi/%=δxi⋅100%.

5. 最大绝对误差和最大相对误差 ​

(1) 最大绝对误差 欲确定量 z ,且 z 是可测量 x1,x2,⋯,xn 的函数,即 z= f(x1,x2,⋯,xn) ,则计算所产生的误差时必须考虑到函数 f . 有两种不同方法检测误差. 第一种方法是运用统计误差分析,即通过最小二乘法 (使 ∑(zi−z)2 最小) 平滑数据值. 第二种方法是确定绝对误差量的上界 Δzmax . 对于 n 个独立变量 xi , 有

(16.207)Δzmax=∑i=1n⁡|∂∂xif(x1,x2,⋯,xn)|Δxi,

其中平均值 x¯i 被 xi 取代.

(2)最大相对误差 最大相对误差是最大绝对误差除以测量值 (多数情况下用

z 的平均值):

(16.208)δzmax=Δzmaxz≈Δzmaxz¯.

16.4.1.4 使用误差界确定测量结果 ​

只有也给出了期望误差时, 才可能对测量结果进行实际解释, 误差估计值和误差界是测量结果的组成部分. 根据数据必须弄清楚误差的类型、置信区间和显著性水平.

(1)定义误差 单次测量值可用下述形式给出

(16.209a)x=xi±Δx≈xi±σ~,

且均值满足

(16.209b)x=x¯±ΔxAM≈x¯±σ~AM.

对于 Δx ,最常使用标准差 σ~ ,有时也用 γ~ 和 η~ .

(2) 任意置信限法则 若总体服从 N(μ,σ2) 分布,根据 (16.96),量 T=X−xwσ~AM 服从自由度为 f=n−1 的 t 分布 (见 (16.97b)). 对于预期的显著性水平 α 或接受概率 S=1−α ,未知量 xw=μ 在 t 分位数 tα/2;f 处的置信限是

(16.210)μ=x¯±tα/2;f⋅σ~AM,

即真实值 xw 以概率 S=1−α 落在上述置信限给出的区间内.

在大多数情况下,以尽可能低的水平保持测量序列的容量 n 很有意义. 当 1−α 值减小,测量次数 n 增大时,置信区间的长度 2tα/2;fσ~AM 减小. 由于 σ~AM 对应 1/n 成比例地减小,自由度 f=n−1 的分位数 tα/2;f 也对应 1/n 成比例地减小 (对于 5 到 10 之间的 n 值,可参见第 1463 页表 21.20),对于上述 n 值,置信区间的长度对应 1/n 成比例减小.

16.4.1.5 精确度相同时直接测量的误差估计 ​

若 n 次测量都有相同的标准差 σi ,则测量的精确度相同, h= 常数. 在这种情况下, (16.200c), (16.201b) 和 (16.202b) 给出了由最小二乘法产生的误差量.

  • 下表给出的测量序列中,包含了精确度相同的 n=10 次直接测量数据,试确定其最终结果.
x¯=1.580,σ~=∑i=1n⁡vi2/(n−1)=0.0131,σ~AM=σ~/n=0.004.

xi

1.592

1.581

1.574

1.566

1.603

1.580

1.591

1.583

1.571

1.559

vi⋅103

-12

-1

+6

+14

-23

0

-11

-3

+9

+21

vi2⋅106

144

1

36

196

529

0

121

9

81

441

最终结果: x=x¯±σ~AM=1.580±0.004 .

16.4.1.6 精确度不同时直接测量的误差估计 ​

1. 加权测量 ​

如果直接测量值 xi 根据不同测量方法得到,或表示具有相同均值 x¯ 、不同方差 σ~i2 的单次测量的均值,则计算加权平均

(16.211)x¯(g)=∑i=1n⁡gixi/∑i=1n⁡gi

其中 gi 定义为

(16.212)gi=σ~2σ~i2.

此时, σ~ 是任意正值,通常是最小的 σ~i . 它可作为偏差的加权单位,即对于 σ~i=σ~ , 有 gi=1 . 由 (16.210) 可推出,测量的权重较大时,会产生较小的偏差 σ~i .

2. 标准差 ​

加权单位的标准差可估计为

(16.213)σ~(g)=∑i=1n⁡givi2/(n−1).

必须确保 σ~(g)<σ~ . 相反,若 σ~(g)>σ~ ,则存在带有系统偏差的 xi 值.

单次测量的标准差是

(16.214)σ~i(g)=σ~(g)gi=σ~(g)σ~σ~i,

其中预期 σ~i(g)<σ~i .

加权平均的标准差是

(16.215)σ~AM(g)=σ(g)/∑i=1n⁡gi=∑i=1n⁡givi2/((n−1)∑i=1n⁡gi).

3. 误差描述 ​

可通过 1111 页 16.4.1.4 描述误差,或通过误差的定义,或根据自由度为 f 的 t 分位数.

x¯i

σ~AMi

σ~AMi2

gi

zi

gizi

zi2

gizi2

1.573

0.010

1.0⋅10−4

0.81

−1.2⋅10−2

−9.7⋅10−3

1.44⋅10−4

1.16⋅10−4

1.580

0.004

1.6⋅10−5

5.06

−5.0⋅10−3

−2.5⋅10−2

2.50⋅10−5

1.26⋅10−4

1.582

0.005

2.5⋅10−5

3.24

−3.0⋅10−3

−9.7⋅10−3

9.0⋅10−6

2.91⋅10−5

1.589

0.009

8.1⋅10−5

1.00

+4.0⋅10−3

4.0⋅10−3

1.6⋅10−5

1.6⋅10−5

1.591

0.011

1.21⋅10−4

0.66

+6.0⋅10−3

3.9⋅10−3

3.6⋅10−5

2.37⋅10−5

(x¯i)mσ~

∑i=1n⁡gi

∑i=1n⁡gizi

∑i=1n⁡gizi2

=1.583=0.009

=10.7

=3.6⋅10−2

=3.1⋅10−4

计算可得 (x¯i)m=1.5830 ,选取 x0=1.585 和 σ~=0.009 ,且 zi=x¯i− x0,gi=σ~2/σ~i2 ,得到 z¯=−0.0036 和 x¯=x0+z¯=1.582 . 标准差是 σ~(g)= ∑i=1n⁡givi2/(n−1)=0.0088<σ~ 和 σ~x=σ~AM=0.0027 . 最终结果是 x=x¯±σ~x=1.585±0.0027.

version 1.24.0