一种蒙特卡罗方法在测量不确定度评定中的新算法.pdfVIP

t计i与测试技拳》2D11年第38卷第10期 一种蒙特卡罗方法在测量不确定度评定中的新算法 ANew MomeCanoMethodto Co玎啦llt口tion以Method对Using El，口胁口绍朋名倪s髓彤mP咒珞【砌cP，纪觑秒 崔伟群 (中国计量科学研究院北京，100013) 摘要：通过对现有利用蒙特卡罗方法进行不确定度评定过程的分析，指出了且前评定过程中存在理论上的不衔接，并给出了一种全新的利用蒙特卡 罗方法进行不确定度评定的算法．弥补了现有方法的缺憾。 关键词：蒙特卡罗；不确定度；不确定度评定 O引言 的矩形分布； 蒙特卡罗方法(MCM)的基本思想是依据大数定律， (3)选择蒙特卡洛试验样本量的大小肘。 当求解某种随机事件出现的概率，或者是某个随机变量 (4)从输入量x的PDF中抽取肘个样本值戈，，．『=l， 的期望值时，可以通过某种“实验”的方法，以这种事件出 2，…，M 现的频率估计这一随机事件的概率或期望，并将其作为 (5)对每个样本巧，计算相应y的模型值力=写，．『= 问题的解。 l，2，…，肘 近年来，蒙特卡罗方法被大量应用在测量不确定度 评定中，核心是求取模拟样本的标准差。其依据是分布 (6)由肘个，7计算y的样本均值为歹=玄聋勺及样 传播，即能够通过输入量的概率分布确定输出量的概率 本方差为s2()，)=面气蚤(勺一歹)2。 分布。在这一点上一般没有疑义，但是蒙特卡罗方法应 假设实验员通过上述方法进行一次评定，评出其标 用于不确定度传播的前提条件是输出量的样本均值的方 准不确定度为s，()，)。但是实验员为了更准确，又利用 差就应该是实际中要求的方差，然而蒙特卡洛方法未给 上述方法进行了一次评定，评出其标准不确定度为s2 出在模拟“实验”的过程中可以模拟方差的结论，这是因 (，)。那么如何处理这两次评定的结果呢?取平均值 为不同的“实验”过程会引入不同的误差，从而产生不同 吗?显然不是，这是因为在概率统计中，有明确的方法处 的样本方差。这从蒙特卡罗方法中给出的减小方差的各 理这两次评定的结果。对于第一次模拟有一组随机变量 种算法中可以看出。 工l，对与第二次模拟也有一组随机变量x2，则两次模拟 目前使用的蒙特卡罗方法进行不确定度评定一般按 下述步骤进行： 的结果为，，：垒≮土，其方差为： (1)建立y和Xl’．··粕之间的模型y=，(工l，…， 粕)； (2)利用可获信息，为置设定概率分布(肋F)等；笪[茎!二墨[墨!](茎2二墨[茎2321 ，n 、工7 2 (3)选择蒙特卡洛试验样本量的大小肘。 考虑到实验员两次评定的独立性，可以认为两组模 (4)从输入量玉的只凹厶．(五)中抽取肘个样本值， 拟数据不相关，则其标准不确定度为： i=1，2，…，Jv，_『=1，2，…，肘； (5)对每个样本矢量(茗l’石2，，…，髫舾)，计算相应y的 s2(，，)：煎学垃 (2) 模型值力=以茗v，锄，…