公开/公告号CN1912861A
专利类型发明专利
公开/公告日2007-02-14
原文格式PDF
申请/专利权人 上海市计量测试技术研究院;
申请/专利号CN200510028587.8
申请日2005-08-08
分类号G06F17/00(20060101);
代理机构31128 上海世贸专利代理有限责任公司;
代理人李浩东
地址 200233 上海市宜山路716号
入库时间 2023-12-17 18:16:49
法律状态公告日
法律状态信息
法律状态
2011-10-12
未缴年费专利权终止 IPC(主分类):G06F17/00 授权公告日:20081210 终止日期:20100808 申请日:20050808
专利权的终止
2008-12-10
授权
授权
2007-04-11
实质审查的生效
实质审查的生效
2007-02-14
公开
公开
[技术领域]
本发明涉及计量测试技术领域,尤其是涉及一种用蒙特卡罗统计模拟计算合成不确定度的方法。
[背景技术]
不确定度的概念在计量测试技术领域占有着重要的地位。任一有效的测量必须带有不确定度的评定,以便于人们可以对这一测量过程的优劣进行评估、比较和复现这一测量过程。它在科技领域成为不可缺少的部分,所以国际上于1993年由国际标准化组织(ISO)等7个国际组织的名义联合发布了《测量不确定度表示指南》(Guide to the Expression of Uncertainty in Measurement简称GUM)。我国也于1999年发布了JJF1059-1999《测量不确定度评定与表示》,以便于与国际同步,可对重要的测量过程进行统一的评估,也保证了计量领域的量值传递的可行性和统一性。
但是,不确定度伴随着任一测量过程。许多领域,如工业、商业以及卫生和安全等领域,需要伴随着这个测量值的不确定度有一个确切的置信区间,以期望将其作为依据,预测项目的可行性。人们期望得到的数据(不确定度)既不夸大可预知性,即在概率要求确定的情况下,缩小不确定度;也不缩小预知性(如果发生相反的情况),而期望获得正确的预测。由于大多数测量结果的概率分布是很难用数学方法得出的,所以上述标准根据不确定度传播率导出的合成标准不确定度uc乘以复盖因子kp得到的扩展不确定度Up,并不确切地与置信水平的概率相关联(如p=95%,,表示大约应有95%的测量结果落入此区间)。
[发明内容]
本发明所要解决的技术问题是提供一种可使最终的合成不确定度与置信水平精确关联的用蒙特卡罗统计模拟计算合成不确定度的方法,它包括如下步骤:a、引入测量输入量和输出量函数关系的数学模型Y=f(x1,x2,…xn),上式中Y为输出量、x1,x2,…xn为输入量;b、引入输入量x1,x2,…xn不确定度的概率分布和参数;c、根据输入量x1,x2,…xn不确定度的概率分布和参数,选择模拟量ξ1,ξ2,…ξn;d、选择对模拟量ξ1,ξ2,…ξn的检验方法;e、在计算机上进行模拟计算输出量Y、输出量Y的标准偏差和与置信水平精确相关的合成不确定度。本发明所提出的任务还可进一步通过如下技术方案加以实现:在计算机上进行模拟计算是将各随机数ξ1,ξ2,…ξn代入Y=f(x1,x2,…xn)得到yi=f(ξ1i,ξ2i,…ξni),各随机数的取值范围在各不确定度分量分布的全范围,相应得到的yi即为输出量的可能值,i=1,2,…n,在计算机上产生n个模拟值y1,y2,…yn,输出量
本发明可使最终的合成不确定度与置信水平精确关联,特别是高置信水平的合成不确定度(置信水平95%或99%等)。如果根据不确定度传播率导出的合成标准不确定度uc乘以复盖因子kp得到的扩展不确定度Up,由于输出量的分布往往不可知,而不确切地与置信水平的概率相关联。本发明将在许多领域,如工业、商业、医药卫生等领域发挥很大的作用。
[具体实施方式]
表1为模拟计算的计算机程序流程。
表2为模拟计算Y=f(D0,Mx,Dx)的分布。
表3为模拟计算的10次重复计算值(不同模拟次数)。
实施例:本实施例是一个校准扫描电子显微镜记录图像上显微标尺的实例,其中测量数学模型即为计算校准值的计算公式:
M=D0×Mx/Dx
式中:
M——显微标尺的校准值;
D0——标准物质的长度值;
Mx——记录图像上显微标尺的长度测量值;
Dx——标准物质的长度测量值。
各输入量不确定度的确定和产生随机数的方法
在测量的数学模型f(D0,Mx,Dx)=D0×Mx/Dx已建立的情况下,其中各输入量D0,Mx,Dx的不确定度假设确定如下:
标准物质的长度值D0=4.6μm±0.05μm,不确定度服从正态分布N(a,σ),这里a=4.6μm;σ=0.05μm。标准物质的长度值用模拟数ξ1i表示,选择通常蒙特卡罗方法介绍的正态分布N(a,σ)的随机数产生的方法,可以得到:
记录图像上显微标尺的长度测量值Mx=15mm±0.5mm,不确定度服从(a,b)区间的均匀分布,这里a=15mm-0.5mm=14.5mm;b=15mm+0.5mm=15.5mm。记录图像上显微标尺的长度测量值用模拟数ξ2i表示,选择通常蒙特卡罗方法介绍的对任意区域(a,b)上均匀分布的随机数产生的方法产生,可以得到:
ξ2i=a+(b-a)r2i=14.5+r2i(其中r2i为(0,1)上均匀分布的随机数)
标准物质的长度测量值Dx=36mm±0.5mm,不确定度服从在(a,b)区间的均匀分布,这里a=36mm-0.5mm=35.5mm;b=36mm+0.5mm=36.5mm。标准物质长度测量值用模拟数ξ3i表示,选择通常蒙特卡罗方法介绍的对任意区域(a,b)上均匀分布的随机数产生的方法产生,可以得到:
ξ3i=a+(b-a)r3i=35.5+r3i(其中r3i为(0,1)上均匀分布的随机数)模拟数的检验
由于本实例是采用美国Microsoft公司的产品Visual Basic 6.0软件编制的计算机程序来进行模拟计算,随机数的产生是该软件所附的随机数发生器产生的。可以认为这样获得的在(0,1)上均匀分布的随机数系列是通过各项参数检验、均匀性检验和独立性检验的,所以在计算程序中,我们仅对在该随机数基础上构造起来的在任意区域(a,b)上均匀分布的随机数和正态分布N(a,σ)的随机数的平均值
对于任意区域(a,b)上均匀分布随机数的数学期望和方差表示为:
E(y)=(a+b)/2
σ2=E(y2)-[E(y)]2=(b3-a3)/3-(a+b)2/4
由中心极限定理知统计量:
当N充分大时,渐近地服从N(0,1)分布。
对于正态分布N(a,σ)的随机数的数学期望和方差直接就是a和σ,所以由中心极限定理得到统计量:
当N充分大时,渐近地服从N(0,1)分布。
取显著性水平α=0.1,则当|u|>1.645时,计算机程序将重新产生随机数,一直到产生的随机数通过参数检验然后继续进行模拟计算。
模拟计算的计算机程序流程参见表1。
模拟计算结果:
模拟计算Y=f(D0,Mx,Dx)的分布参见表2。
Y的模拟计算平均值:E(y)=1.91670410155901μm;模拟计算次数:n=100000;Y的模拟计算标准偏差:S=4.62415096429435E-02μm;模拟计算值的置信概率为95%的合成不确定度:U95=0.086317873244385μm。
本发明还可以采用多次重复计算取得平均值,以便获得更高的精度。模拟计算的10次重复计算值(不同模拟次数)参见表3。
如果取n=100000的10次重复计算平均值为最终结果,可得校准值及不确定度为:
图像上显微标尺校准值M=1.917μm,标准偏差s=0.046μm,U95(M)=0.085μm。
根据《测量不确定度表示指南》(Guide to the Expression of Uncertainty inMeasurement简称GUM)和中国国家标准JJF1059-1999《测量不确定度评定与表示》中规定的不确定度的传播率
表1
表2
(续下页)(接上页)
表3
机译: 蒙特卡罗器件仿真方法和用于执行相同方法的蒙特卡罗器件仿真程序
机译: 蒙特卡罗模拟的统计改进系统
机译: 蒙特卡罗模拟的统计改进系统