蒙特卡罗模拟8.pptVIP

蒙特卡罗模拟 一. 蒙特卡罗法计算定积分 二. 蒙特卡罗模拟试验次数的确定 例7.3.2 核反应堆屏蔽层设计问题 * * 蒙特卡罗(Monte－Carlo)模拟，又称蒙特卡罗方法、统计试验法等. M－C模拟是静态模拟，描述特定时间点上的系统行为. 模拟过程中不出现时间参数。 基本思想:把随机事件 （变量）的概率特征与 数学分析的解联系起来. 概率特征：随机事件的概率和随机变量的 数学期望等. 用试验方法确定 例7.3.1 用M－C 模拟求圆周率π的估计值. 1 1 0 设二维随机变量 (X, Y)在正方形内 服从均匀分布. (X, Y)落在圆内的概率为: 计算机上做n 次掷点试验： 产生n 对二维随机点(xi，yi) ，i＝1 ,2, …, n . xi 和yi 是RND 随机数对. 检查每对随机数是否满足: 相当于第i个随机点落在1/4圆内. 若有k 个点落在l／4圆内 随机事件“点落入1/4圆内”的 频率为 k/n 根据概率论中的大数定律, 事件发生的频率 依概率收敛于事件发生的概率p，即有 得圆周率π的估计值为 且当试验次数足够大时, 其精度也随之提高. 分析：实际上概率值为 恰为1/4圆的面积 频率法： 利用随机变量落进指定区域内的频率来计算定积分. 平均值