实战宝典 原创大赛 仪器问答 仪友会 总帖：8043206今日发帖：436 我的社区 签到
快速导航 采购仪器 提醒
实验室认可/资质认定
版主:
入住本版
专家:
居民列表

主题：【第九届原创】蒙特卡洛模拟评定不确定实例

浏览0 回复10 电梯直达
普通用户63
头像
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
发表于：2016/08/20 21:32:27 楼主 管理 分享 倒序浏览 只看楼主 回复 私聊
维权声明：本文为zzg2002原创作品，本作者与仪器信息网是该作品合法使用者，该作品暂不对外授权转载。其他任何网站、组织、单位或个人等将该作品在本站以外的任何媒体任何形式出现均属侵权违法行为，我们将追究法律责任。
随着jjf1059.2的实施，大家越来越关心蒙特卡洛模拟法（MCM）在不确定度评定中的应用，现在举一个较简单的例子，大家一起学习下。
首先蒙特卡洛法法的提出，是为了克服jjf1059.1（GUM）存在的缺点。
MCM法的原理 :

        MCM法无论测量模型的函数满足何种分布，依据各个分量的分布，生成相应的随机数对分析结果的均值和标准差（不确定度）进行模拟；MCM法克服了GUM法对测量模型函数的分布都满足正态分布这一理论缺陷。
MCM法的步骤

  1、确定测试模型

  2、确定各个分量的分布

          通过测量工具的检定证书

          引用标物的标物证书

          根据以往的经验和惯例对某个分量分布的处理

    3、依据测试模型进行模拟，模拟次数一般取n>10^6

          随机数的生成可以借助matlab等工具软件；也可以自行编程解决，依托开发工具提供产生的均匀分布随机数来产生其他分布的随机数。

    4、模拟结果升序排列，计算模拟结果的均值和标准偏差，然后依据给定的概率划定不同的包含区间，以最短的包含区间为评定结果。
蒙特卡罗模拟法评定不确定度的简单例子
假设测试模型为z=x+y，x符合分布标准正态分布x~（0,1）；y符合（0,1）之间的均匀分布U（0,1），给出置信概率95%下的z的置信区间和标准差。
    1、确定测试函数

    2、确定模拟次数

    3、按要求生成各分量的随机数，依据测试模型计算模拟结果

    4、计算模拟结果平均值和标准差

    5、依据给定的置信概率，划分置信区间，以最短的置

          信区间作为评定结果



评定结果：
    模拟次数10000次，均值：0.505，标准差：1.04,95%置信概率下的置信区间（-1.55,2.47）。
这个模型均值的理论值0.5，标准差的理论值是1.04。可见，用MCM法模拟10000次，模拟值和理论值吻合较好。
这个模型用GUM法评定的结果为：置信概率95%（k=2），置信区间为：（-1.58，2.58）。两种评定方法结果接近，但是GUM法认为分布函数是对称的，和实际不符。
为您推荐
近期热榜
热门活动
您可能想找: 气相色谱仪(GC) 询底价
选参数 看心得 找厂商 查方案
专属顾问快速对接
立即提交
可能感兴趣
刘彦刚
pxsjlslyg
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 6:52:16 沙发 管理 分享 倒序浏览 只看楼主 回复 私聊
你使用matlab吗？这个很难哦！
赞贴
0
收藏
0
拍砖
0
举报
普通用户63
zzg2002
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 7:05:18 板凳 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 刘彦刚(pxsjlslyg) 发表:
你使用matlab吗？这个很难哦！
用VB编程，用matlab太繁琐了，而且不能进行多次模拟。
赞贴
0
收藏
0
拍砖
0
举报
刘彦刚
pxsjlslyg
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 7:09:48 马扎 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 普通用户63(zzg2002) 发表:
用VB编程，用matlab太繁琐了，而且不能进行多次模拟。
一般来说，你现在给出的是X与Y不相关的情况，你试想如果输入量X和Y强相关呢？你该怎样处理？能与不相关怎样去区别呢？
赞贴
0
收藏
0
拍砖
0
举报
刘彦刚
pxsjlslyg
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 7:16:03 地毯 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 普通用户63(zzg2002) 发表:
用VB编程，用matlab太繁琐了，而且不能进行多次模拟。


我不会VB，更不会matlab，而是一现成的专用用软件：
赞贴
0
收藏
0
拍砖
0
举报
zyl3367898
zyl3367898
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 8:18:29 地板 管理 分享 倒序浏览 只看楼主 回复 私聊
听了微课堂，讲的真好。
赞贴
0
收藏
0
拍砖
0
举报
普通用户63
zzg2002
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 10:58:59 6楼 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 刘彦刚(pxsjlslyg) 发表:
一般来说，你现在给出的是X与Y不相关的情况，你试想如果输入量X和Y强相关呢？你该怎样处理？能与不相关怎样去区别呢？
给出相关系数，然后生成有相关性的随机数就可以了。
赞贴
0
收藏
0
拍砖
0
举报
普通用户63
zzg2002
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/21 11:30:38 7楼 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 刘彦刚(pxsjlslyg) 发表:
我不会VB，更不会matlab，而是一现成的专用用软件：

这个软件是一款通用的mcm模拟软件，评定结果不满足jjf1059.2的要求：jjf1059.2需要最短置信区间；分量分布类型不全，没有常用的梯形分布和反正弦分布。

赞贴
0
收藏
0
拍砖
0
举报
刘彦刚
pxsjlslyg
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/22 9:26:20 8楼 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 普通用户63(zzg2002) 发表:
用VB编程，用matlab太繁琐了，而且不能进行多次模拟。
你是每次评定都用VB编程，还是先用VB编好一个程序，具体评定时再给出模型和对应的输入量呢？
赞贴
0
收藏
0
拍砖
0
举报
刘彦刚
pxsjlslyg
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/22 9:29:31 9楼 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 zyl3367898(zyl3367898) 发表:
听了微课堂，讲的真好。


我可惜没法入群没有听到讲课。我加你微信（我的手机号13979972566），你能分享给我吗？
赞贴
0
收藏
0
拍砖
0
举报
2016/8/22 9:56:14 Last edit by pxsjlslyg
普通用户63
zzg2002
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2016/8/22 9:30:14 10楼 管理 分享 倒序浏览 只看楼主 回复 私聊
原文由 刘彦刚(pxsjlslyg) 发表:
你是每次评定都用VB编程，还是先用VB编好一个程序，具体评定时再给出模型和对应的输入量呢？
当时是后者。
赞贴
0
收藏
0
拍砖
0
举报
手机版： 蒙特卡洛模拟评定不确定实例
猜你喜欢最新推荐热门推荐更多推荐
头像
高级回复 发表评论
品牌合作伙伴
岛津
日立科学仪器
珀金埃尔默仪器（上海）有限公司（PerkinElmer）
日本电子株式会社
丹纳赫
安捷伦
赛默飞世尔科技
普析通用
欧波同
天美
天瑞仪器
德国耶拿
海能技术
马尔文帕纳科
磐诺科技
上海仪电科仪
梅特勒托利多
聚光科技
莱伯泰科
盛瀚
多宁生物
丹东百特
科哲
卓立汉光
屹尧科技
华谱科仪
宝德仪器
优莱博
HORIBA
布鲁克核磁