基于不确定度分析的计算结果有效数字确定方法

+士—~/1I1echnical
不/丄Column
基于不确定度分析的计算结果有效数字确定方法
张馨睿I,吴飒I,陈仁军食周学亮3
(1.北京航空航天大学，北京100083； 2.陆军装备部航空军事代表局，上海200233；
3.苏州天航长鹰科技发展有限公司，苏州215155)
摘要：本文提出一种基于不确定度传播原理对确定计算结果有效数字的方法。

当需要根据若干个变量间接计算得到时，每个自变量由于切尾修约产生的有效数字所带来的不确定度都会传递到计算结果中，通过不确定度传播规律对其进行分析，在此基础上进行计算结果有效位的确定。

本文对几种常用的计算公式给出了确定有效数字的具体说明，并通过实例证明根据常用计算法则与基于不确定度分析计算法则得出的有效位数的不同。

相对于常用的计算结果有效数字确定方法[1],本文提出的方法可以有效降低对计算结果错误使用的风险。

关键词：运算规则；不确定度；准确度有效数字
中图分类号：T文献标识码：A文章编号：1004-7204(2020)05-0161-06
Method for Determining Effective Number of Calculation Result Based on
Uncertainty Analysis
ZHANG Xin-rui1,WU Sa1,CHEN Ren-jun2,ZHOU Xue-liang3
(1.Beihang University,Beijing10083;
2.Aviation and Military Representative Bureau of the Ministry of Army Equipment,Shanghai200233;
3.Suzhou Tianhang Changying Technology Development Co.,Ltd.>Suzhou215155)
Abstract:This paper proposes a method for determining the effective number of calculation results based on the uncertainty propagation principle.When it needs to be calculated indirectly based on several variables,the uncertainty caused by the effective number of each independent variable due to trimming and rounding will be transferred to the calculation resuIt,and it will be analyzed by the uncertainty propagation law.On this basis,the effective bit of the calculation resuIt is determined.
This article gives specific instruetions for determining the effective number of several commonly used calculation formulas,and proves the difference between the effective number of digits according to the commonly used calculation rules and the calculation rules based on uncertainty analysis through pared with the commonly used method for determining the effective number of calculation resuIts,the method proposed in this paper can effectively reduce the risk of incorrect use of calculation resuIts.
Key words:algorithm；uncertainty；accuracy significant number
背景
在科学研究与工程实践中，不可避免地要进行大量的数据计算，也就是把一些数值作为已知量，通过规定的计算方法以数值的方式给出运算结果。

无论是作为已知量还是作为运算结果的数值都是通过有效数字的方式给出，也就是说这些数值在本质上具有一定的不确定性。

数据计算不仅应当按照计算方法给出理论结果，也应当把来自于数据源和计算方法的不确定度信息正确地传递给计算结果的使用者，一般是以有效数字位的方式完成数据结果和不确定度信息的传递。

因此，如何确定计算结果的有效位数是十分必要的，不正确的有效数字位会造成计算结果中所传递的不确定度被错误地使用。

正是
161环境技术/Environmental Technology
T echnical/~
Column/玖不二
因为有效数字位与不确定度存在着本质上的关联，因而有必要通过不确定度分析的方法进行有效数字位的确定。

1概念
1.1有效数字的含义
数据的来源一般有两种途径，一种是通过宜接测量得到的测量量，另一种是通过计算的得到的计算结果。

对于第一种情况，有效数字是反映被测量大小的全部数字，包含准确数字以及一位带有不确定度的估读数字。

有效数字的最末一位是估读的，因此只有最末一位数字具有不确定度，有效数字的不确定度代表了估读数字的不确定度，可以通过有效数字的不确定度大小来确定哪位数字是估读数字。

有效数字也在一定程度上表达了测量结果的不确定度叫其最末数字位数离小数点越近，被测值的不确定度越大。

对于第二种情况，通过对计算结果的修约得到其有效数字。

无论是直接测量值还是通过间接计算得到的数值，当以有效数字表示时，都要遵循只有最末位是估计值，其他位都是准确数字的要求。

只有这样才能保证正确反映出分析量的多少，才能使结果被正确使用。

1.2不确定度的含义
不确定度的概念首先应用在测量结果的表述中，用以合理表示被测值分散程度。

无论是直接测量量还是通过计算得到的测量值，其结果的完整表示应该是被测量量加上相应的不确定度。

不确定度反映了测量的质量，不确定度小，表明测量结果距真值的偏离范围就越小，质量就越高，其使用价值就越高，相反，不确定度越大，则偏离范围就越大，质量就越低，其使用价值就越低。

2测量不确定度评定
进行测量不确定度评定时，比较常用的方法是应用测量不确定度传播率的方法(GUM法)，具体过程包括建模、标准不确定度评定、计算合成标准不确定度和计算扩展不确定度。

2.1建模
被测量的值的不确定度取决于参与计算的其他分
量，确定这种取决关系即为建模。

被测量Y的最佳值
y由X”X2,•••,X n的最佳值x”x2,…，x…得到，即
y=f(xi,x2,x…)。

y的标准不确定度u(y)取决于x”
X2,…，九的标准不确定度。

2.2标准不确定度评定
A类评定方法是为了得到测量数据的标准差的最佳
估计值，通过一系列测量数据使用统计方法进行评定。

通常使用方法有贝塞尔法、最大残差法、极差法、最小
二乘法。

B类评定方法是通过非统计方法进行的，B类评定方
法中最常见的是使用来自于技术资料、证书或手册中提
供的数据。

这些数据往往提供了可能的取值范围，但其
概率分布要进行假设，根据已知信息可假设为正态分布、
均匀分布、反正弦分布、双均匀分布、三角分布，如果
信息缺乏，一般情况下可以假设为均匀分布。

2.3合成不确定度计算
合成不确定度是在不确定度传播公式下，自变量的
标准不确定度经函数关系传播到因变量所产生的不确定
度。

不确定度传播公式为：
(1)
假设自变量之间互不相关，则呛“对=0合成不确
定度公式为：
U C(JV)=U c=若(拾)%(冷)2(2)
2.4扩展不确定度计算
扩展不确定度表示测量值大部分合理地包含在此区
间内，扩展不确定度由合成不确定度与包含因子相乘求
得：
包含因子常用求取方法如下叫
2020年10月/October2020
162
+士—4^.~/1I1echnical 玖不PF二/1Column
1）当参与运算分量的自由度已知时，根据下式（4）计算合成不确定度的有效自由度
V：—（x;）的自由度。

由t分布求出包含因子P为置信水平，可取0.95,0.99, 0.997,可根据t分布表，通过p与v查表得tp（v）,即k。

2）当参与运算分量的自由度未知时，可通过赋值法求取k。

实际测量工作中，大多数测量量的分布可以近似为正态分布，对于正态分布的情况，包含因子取决于期望的置信概率。

在工程上一般的置信概率取值为90%、95%、99%,对应的包含因子分别为2、2.56、3。

其他分布情况，也有类似的幅值方法，具体可以参阅测量不确定度评定和数理统计相关书籍。

3计算结果的不确定度
计算结果的不确定度评定与测量的不确定度评定非常类似。

3.1建模
计算结果的来源公式即为模型。

3.2标准不确定度评定
计算所面临的情况不存在一系列测量数据，不使用统计方法进行评定,因此不使用A类评定方法。

计算过程的不确定度应使用B类评定方法进行评定，数据可来源于资料或前期计算结果，因此在进行标准不确定度评估时，参与计算的每个数值都是不确定度评定时的来源。

如果数值给出了不确定度信息，则可以利用该信息计算其标准不确定度。

如果该数值仅以有效数字给出，由于没有给出其他信息，故假设其为均匀分布，根据修约规则，其量值变化半宽度为a,标准不确定度等于寻如
3.3合成标准不确定度计算
计算结果的合成标准不确定度计算与测量不确定度计算类似。

合成标准不确定度是在不确定度传播公式下，参与计算分量的标准不确定度经灵敏度系数传播到计算结果所产生的不确定度。

不确定度传播公式为：
J n—I N
若诒）2心2+2a g云函唤宀）“（咖（初
（5）假设参与计算的分量之间互不相关，则r（x“Xj）=0,合成标准不确定度公式为：
3.4扩展不确定度计算
计算结果的合成标准不确定度计算与测量不确定度计算类似，扩展不确定度由合成不确定度与包含因子相乘求得：
U(N)=kx Uc(N)（7）k通过赋值法求得，取决于预期的置信概率，在工程上一般的置信概率取值为90%、95%、99%,对应的包含因子分别为2、2.56、3O
4计算结果的有效数字确定
4.1常规方法
数值计算是科研工作的核心，计算公式可以是常规的算术计算，也可以是非常复杂的计算方法。

下面就几种常见的计算方法，给出目前采用的有效数字确定原则[1]o
加减法：最终结果的有效数字的末位与参与运算变量中有效数字末位最高位者相同。

乘除法：最终结果的有效数字的位数与参与运算变量中有效数字位数最少者相同。

乘方运算：最终结果的有效数字位数与底数有效位数相同。

对数运算：最终结果的有效数字位数与真数的有效位数相同。

163环境技术/Environmental
Technology
echnical/-JA
丄Column/孑二
4.2基于不确定度分析的方法
有效数字的不确定度与最末一位修约数字有直接联系，无论是直接测量还是间接计算得到的结果，其不确定度都表现了计算结果在哪个数位上是不确定的，确定有效数字位数即确定修约数字在哪位，可以通过计算结果的大小与其不确定度来确定。

进行最终结果报告时，结果的末位应与不确定度的末位保持一致（如GJB 3756）o在进行计算时，根据修约规则，有效数字中只包含一位不确定数字就是最末一位修约数字，因而对于不确定度的关注只有一位数字即首位数字。

根据预期的置信概率对应的包含因子，得到扩展不确定度，然后根据通用的修约规则，得到只有一个数字位的数值，如该位数值大于等于5,则有效数字位进一位，否则有效数字位保持不变。

例如0.04对应的有效数字位是小数点后第二位，0.05对应的有效数字位是小数点后1位[5]°
4.3在常用计算方法中的应用
1）力口法：N=f（x”X2,±X2土…土 Xn
合成不确定度：
J（务）2则+（爲）2吨+••••••+（知2“玮
:___________________（8）
=y/uxl+U%2+............+
这是最简单的计算过程，当只有一个数值具有最少有效数字时,它将占主导地位,其他数值的影响可以忽略。

通过公式就可以明显看出。

这种情况下与常规方法是一致的。

只有当存在多个数值具有相同的最小有效数字时，才需要根据不确定度传播公式计算其合成标准不确定度。

2）乘法：N=f（x”x2,•••x n）=x1Xx2x•••xx…
合成不确定度：
-x n)2uxl+•••+(x1x2--x n_1)2izx2
（9）
3）除法：N=f（xi，X2,•••Xn）=Xi+X2+…+召
合成不确定度：
昭N）=J（签）2咗+…+（知2曲
彳gx…吩曲+…+(Xl X*X…X占)2歸
(10)
4）乘方（以平方运算为例）：N=f（x）=x2
合成不确定度：
%(Np u(x)2=|2%u(%)|(ID
5）乘方（以开方运算为例）：N=/（%）=五
合成不确定度：
（鴛）2吩）2=
6）对数:N=f（x）=log a x
合成不确定度：
如（2）打（野吩）勺佥吩）
7）指数:N=f（x）=a x
(12)
(13)
合成不确定度：
UcW=l(^)2u(x)2=la x In a u(x)|(i4)
5示例与分析
5.1加减法运算
在加减法运算中，如若参与运算的分量中最小有效数值只有一个时，其产生的不确定度占主导地位，其他数值产生的不确定度可以忽略，因此这种情况下，常用的确定方法与本文提出的确定方法得出的结果相同。

例1：N=xi+x2+x3+x4+x5+x6，X]=1.325,x2=1.364,
x3=1.102,x4=1.32,x5=1.364,%=1.182,计算结果N=7.657。

1）使用本文提出的确定方法：
X]、X2、X3、X4、X5、X6服从均匀分布，其标准不确定度为u1=u2=u3=u5=u6="署5=0.000289,u4==0.002887O u c（N）=yjuxl+uxl+uxj+uxl+uxj+uxl=0.002961,
由于只有一个数据占主导，且假设其为均匀分布，因而 认为计算结果也基本服从均匀分布,其扩展不确定度为，0.002961x2〜0.006故有效位进一位即0.01,按照修约规则N的最终表示结果是N=7.66。

2020年10月 /October2020
164
2)使用常用确定方法：
参与运算的有效数字中小数位数最少者为2,因此最终计算结果中保留的小数数位也是2位，最终计算结果为7.66O两种方法给出了同样的结果，原因是在加减法运算中，灵敏度系数是常数1,只有一个数据具有最小的有效位，在本文提出的方法中不确定度取决于最小有效位；常规方法中也取决于最小有效位。

例2:N=x1+x2+x3h-----x100,x1=x2=x1oo=O.O1325,
N=1.325
1)使用本文提出的方法：
Xi、X2、…、X100服从均匀分布，其标准不确定度为
iz1=U2=-,-=Uioo=5=0.000289,
73
Uc(N)=y/uxl+u%2+••-+u%i OO=10u%1=0.00289o
由于每个数据的不确定度都相同，且假设其为均匀分布，因而认为计算结果也基本服从均匀分布，其扩展不确定度为0.00289x2«0.005,故有效位即0.01,按照修约规则N的最终表示结果是N=1.33。

2)使用常用确定方法：
参与运算的最小有效数位为4,因此运算结果的有效数位也是4位,最终计算结果为1.325。

当参与加减法运算的变量修约数字所在位数相同且变量个数增加到一定程度时，采用两种方法得到的结果是不同的，原因是由于变量个数的增加使得扩展不确定度扩大了十倍，有效位也随之发生改变，此时常用的方法仍根据变量中小数位最小者确定有效位是不合理的。

5.2乘法运算
V=axbxh,a=5.375,b=4.230,h=4.110,结果为V=93.4459875O修约截尾产生的不确定度属于均匀分布，标准不确定度u*=警=0.000289、％=詈= 0.000289、u h=0.000289□
贝!luc(V)=』(必)2邂+(ah)2u舟+(ab)2u{=0.010436。

按正态分布95%置信概率的包含因子为2,则扩展不确定度为0.02,确定有效位是与不确定度首位，最终表示结果是V=93.45。

而采用常规方法，最终表示结果的有效位是所有数据中最小有效位数，表示结果是V=93.446O
采用两种方法得到的结果是不同的，原因在于采用常规方法时将不确定度传播中的灵敏度系数全部看作常数1,采用本文提出的方法，在进行不确定度传播计算时,真实的考虑了灵敏度系数，参与运算变量的灵敏度系数大于常数1,因此得到的结果会偏大。

5.3除法运算
v=l，s=40.265,t=1.62,计算结果v=25.07771605。

u==0.000289,u,==0.002887O
则九3)=J(护谐+(_訝诟=0.044326,按正态分布95%置信概率的包含因子为2,则扩展不确定度为0.044326x2«0.1,确定有效位是与不确定度首位，最终表示结果是v=25.1o而采用常规方法，最终表示结果的有效位是所有数据中最小有效位数，表示结果是v=25.08o o
5.4平方运算
E=v2,v=6.23,计算结果E=38.8129。

u,=号齐= 0.002887,u c(E v)=賈踌鬲=0.035972则。

按正态分布95%置信概率的包含因子为2,则扩展不确定度为0.035972x2〜0.1,确定有效位是与不确定度首位，最终表示结果是v=3&8。

而采用常规方法，最终表示结果的有效位是底数有效位数，表示结果是v=38.8o
5.5开方运算
T=VI,L=5.485m,计算结果T=2.342008, u L=嘗s=0.000289则％(T)=J(扣为嗨=o.ooo062,按
正态分布95%置信概率的包含因子为2,则扩展不确定度为0.000062x2〜0.0001,确定有效位是与不确定度首位，最终表示结果是T=2.3420。

而采用常规方法，最终表示结果的有效位是底数有效位数，表示结果是T=2.342O
5.6对数运算
y=log10x,x=20.11计算结果y=1.303412,标准不
(下转171页)
165环境技术/Environmental Technology
［3］杨琪莉，朱兰艳,李海涛•改进FAST特征点支持下的实时影像地标匹配算法卩］•计算机应用，2016,36⑸：1404-1409.
［4］徐志军，耿则勋,魏小峰，等•一种改进的透视不变匹配算法［J］.测绘科学技术学报，2016,33⑶:279-284,
［5］付偲，邓丽，卢根，等•基于快速视网膜关键点算法改进的图像匹配方法卩］•计算机工程与应用,2016,52(19):208-212.
【6］杨欢，沈晓军，李杰，等•基于张量黎曼度量的序列图像匹配光流场计算方法［J］.北京理工大学学报，2016,36(8):862-867,
［7］王峰，尤红建,傅兴玉，等•应用于多源SAR图像匹配的级联SIFT 算法［J］.电子学报，2016,44⑶:548_554.
［8］朱奇光，王梓魏，陈颖•基于全局特征与局部特征的图像分级匹配算法研究及应用卩］.中国机械工程，2016,27(16):2211-2217.
［刃果涛,熊文莉,侯培国,等•基于极曲线几何和支持邻域的鱼眼图像立体匹配［J］-光学精密工程，2016,24(8):2050-2058.
［10］高冠东，王晶，刘菲，等•一种基于极坐标变换的点模式匹配算法【J］.计算机工程与科学，2016,38⑵:331-337.
作者简介:
郭剑(1990-)，男，硕士，助教，研究方向：无人机航空摄影测量;
吴迪(1980-),男，本科，副教授，研究方向：精密工程测量;
贺国伟(1990-),男，硕士，助教，研究方向：开采沉陷。

(上接165页)
确定度吐=器=0.002887,则Uc(N)=J(盒尸诚= 0.0000623,按正态分布95%置信概率的包含因子为2,则扩展不确定度为0.0000623x2«0.0001,确定有效位是与不确定度首位，最终表示结果是y=1.3034。

而采用常规方法，最终表示结果的有效位数与真数的有效位数相同,表示结果是y=1.303°
6结论
本文基于不确定度传播原理对确定计算结果有效数字的方法，计算由自变量修约截尾产生的标准不确定度传播到测量值的合成不确定度，计算结果有效数字的确定与扩展不确定度相关，计算结果的有效位与扩展不确定度的首位相同。

采用常用方法与本文提出的方法得出的结果是有差异的，原因在于常规的方法由于没有考虑计算中各数据对不确定度的贡献是与计算方法中该变量的偏导数有关的，因而会造成一定的偏差，本文提出的方法考虑了在不确定度传播过程中灵敏度系数的影响，根据不确定度更准确地确定了有效位，文中举例说明两种确定法则的差异。

使用本文提出的方法进行有效位的确定仍然会造成不确定度的失真，因此在最终计算结果的表示形式上以有效位加不确定度共同表示较为合理。

参考文献：
［1］蒋志东•有效数字四则运算法则的研究卩］.科技信息,2011,
15(9):85.
［2］唐文强.大学物理实验教学中有效数字运算的几点讨论［J］.科技视
界，2014,(26):178+266.
⑶刘志敏•实验室认可中的不确定度和统计分析［M］.北京：中国标
准出版社，2007:68-76.
［4］陈红霞.化学有效数字的运算法则和修约规则［J］.科技创新导报，
2013,(5):181.
［5］张丽,潘华锦,齐建英•大学物理实验中测量结果不确定度的有效
数字位数【」］•测绘技术装备，2014,25(3):55.
作者简介：
张馨睿(1998-),女，北京航空航天大学可靠性与系统工程学院在读硕士，研究方向：环境试验方法与设备。

171环境技术/Environmental Technology。