蒙特卡罗法用于分析测试不确定度评定研讨.pptVIP

;内容结构;1、蒙特卡洛模拟法介绍;a)蒙特卡洛模拟法计算圆周率 利用单位圆与边长为1的正方形面积之比来计算π的近似值。 具体思想如下: 如图所示,单位圆的1/4为一个扇形,它是边长为1的正方形的一部分。考虑扇形面积在正方形面积中所占的比例k,得出其结果为π/4,然后乘以4就可以得到π的值。 另一方面，假设生成符合均匀 分布（U（0,1）)的随机数对（x，y）， 共生成n个，落入扇形的区域的数有m个， 那么： ;两式结合： 算法： 0、计数器置0，m=0；设置总模拟次数n 1、生成均匀分布的随机数对（x，y） 2、判断（x，y）是否落在扇形内，如果落在扇形内 m=m+1 3、判断模拟次数是否为n：小于或者等于n，回到1，继续 ；大于n，终止，计算4*m/n，输出结果 ;模拟结果;b)蒙特卡洛模拟计算定积分;举例：;;;二、蒙特卡洛法用于分析方法不确定度评定;;;;;二、与不确定度评定相关的计量学名词;;举例：;f)真值、测量结果以及不确定度间的关系;3、不确定度评定方法;GUM法评定不确定度的原理 GUM法评定不确定度时认为：分析结果符合正态分布，依据各个分量通过模型得值来计算得到分析结果；利用方差传播定理和测试模型来计算分析结果的标准差s；通过各分量不确定度和自由度计算分析结果分布的自由度n，然后在给定的置信概率下，确定包含因子K，最后按照要求把评定结果表示出。 注释：由于正态分布的方差未知，测试结果（x）和真值（μ)的关系为： ;b)GUM法不确定度评定的步骤 1、确定测试模型：把间接测量量表示为直接测量量的函数 2、因果图：进行各分量不确定度来源分析 3、各分量不确定度评定：A类和B类 4、各分量自由度确定 5、各分量灵敏度计算（用测试模型对各分量求一阶偏导） 6、依据各分量测试值，计算测试结果x；计算 合成不确定 度U 7、合成自由度计算n 8、在给定的置信概率下，结合自由度n差t分布分位数表， 确定包含因子k 9、评定结果表示：（x-kU，x+kU），同时报告置信概率 ;c) GUM法存在的缺点 GUM法做出的分析结果符合正态分布的假设，在对某些测试模型的处理中与实际情况偏差较大。有的测试模型明显不符合正态分布。 测试模型函数的分布理论上可以依据各个分量的分布通过卷积计算出来，但是在计算过程中需要用到积分，有的函数“不可积”，因而测试模型函数的分布不能准确获得。 ;;;f)MCM法的步骤 1、确定测试模型 2、确定各个分量的分布 通过测量工具的检定证书 引用标物的标物证书 根据以往的经验和惯例对某个分量分布的处理 3、依据测试模型进行模拟，模拟次数一般取n106 随机数的生成可以借助matlab等工具软件； 也可以自行编程解决，依托开发工具提供产生的均匀 分布随机数来产生其他分布的随机数。 4、模拟结果升序排列，计算模拟结果的均值和标准偏 差，然后依据给定的概率划定不同的包含区间，以 最短的包含区间为评定结果。 ;蒙特卡罗模拟法评定不确定度的简单例子;;蒙特卡洛模拟法评定车内空气中甲苯检测的不确定度评定;