利用Mathcad进行测量不确定度的蒙特卡洛方法评定.pdfVIP

利用 Mathcad 进行测量不确定度的蒙特卡洛方法评定 崔伟群 中国计量科学研究院, 北京 100013 摘要：简要介绍了利用 Mathcad 进行测量不确定度的蒙特卡洛方法评定过程中所涉及到的一 个关键函数，并给出了使用 Mathcad 进行不确定度蒙特卡洛方法评定的过程。该过程 编程简单，复用性高，使同类测量的不确定度蒙特卡洛方法评定具有了规范的可能。 关键词：蒙特卡洛方法；Mathcad ；不确定度 0. 引言 蒙特卡洛方法（MCM ）的基本思想是依据大数定律，当求解某种随机事件出现的概率， 或者是某个随机变量的期望值时，可以通过某种“实验”的方法，以这种事件出现的频率估 计这一随机事件的概率或期望，并将其作为问题的解。 近年来，蒙特卡洛方法被大量应用在测量不确定度评定[1-10] 中，核心是求取模拟样本的 标准差。其依据是分布传播，即能够通过输入量的概率分布确定输出量的概率分布。一般使 用蒙特卡罗方法进行不确定度评定的步骤如下： (1) 建立 Y 和X , X 之间的模型Y f (X ，，X )； 1 N 1 N (2) 利用可获信息，为X i 设定概率分布（PDF）等； (3) 选择蒙特卡洛试验样本量的大小M。 j 1,2,,M (4) 从输入量X 的PDF f (X ) 中抽取M 个样本值x ，i=1,2,...,N， ； i X i i ij (5) 对每个样本矢量 (x1j , x 2j , , xNj ) ，计算相应 Y 的模型值y j f (x1j ,x 2j , ,xNj ) ， j 1,2,,M 。 M 个 计算 Y 的估计值y 及 y 的标准不确定度u( y) 。 （6）由 y j 由于蒙特卡洛方法涉及大量的数值模拟和计算，所以必须借助计算机来实现。但是目 前多数编程工具和计算软件都要求使用者具有较高的程序开发水平，导致使用蒙特卡洛方法 进行不确定度评定难以大面积推广。 因而使用一款功能简单、界面友好、操作便捷、对内存要求低、计算速度快、可复用 性好的计算软件成为摆在测量不确定度评定人员面临的问题之一。 1. Mathcad 软件[11] Mathcad 是目前在科技和工程界上著名的四大数学软件之一，该软件集成了数值计算、 符号计算、图形界面、高级语言编程等强大功能。 由于该工具入门简单，并且不需要有编程基础，目前以作为数学、物理、力学和机械 工程等学科方面强有力的二次开发工具被大量应用。 Mathcad 系统具有丰富的内建函数与运算符，并提供外部接口以调用第三方软件，实现 与其他软件的交互和共享。 2. Mathcad 的蒙特卡洛模拟函数 Mathcad 自带一个蒙特卡洛模拟函数，专门用于数值模拟，其格式为 montecarlo(F, n, Rvals, [Limits, dist]) ，其中montecarlo 为蒙特卡洛模拟函数名，该函数有 3 个必选参数和 2 个可选参数。 其中参数 F 为一个实值函数，在不确定度评定中为测量模型的函数名；n 为要生成的蒙 特卡洛样本数；Rvals 为实值函数 F 中自变量的统计信息矩阵，第一列为每个自变量名称的 字符串，第二列为其期望值，第三列为其标准差 σ；Limits (可选) 为长度与 Rvals 相同的 矩阵，该矩阵的第一列和第二列分别为下限和上限，低于或高于此极限时，将裁减掉相应的 随机变量。如果未定义上下限，则必须使用 NaN 填充 Limits 的空元素；dist (可选) 是一 个长度与