计算机模拟课件.pptVIP

若有k 个点落在l／4圆内 随机事件“点落入1/4圆内”的 频率为 k/n 根据概率论中的大数定律, 事件发生的频率 依概率收敛于事件发生的概率p，即有 得圆周率π的估计值为 且当试验次数足够大时, 其精度也随之提高. 计算机模拟全文共53页，当前为第32页。 分析：实际上概率值为 恰为1/4圆的面积 频率法： 利用随机变量落进指定区域内的频率来计算定积分. 平均值法： 利用随机变量的平均值(数学期望) 来计算定积分. 计算机模拟全文共53页，当前为第33页。 平均值法的算法如下： 产生RND 随机数：r1，r2，…，rn； （2）令 ui=a＋(b－a)ri，i=1，2，…，n； （3）计算 作为I 的估计值. 原理分析： 设随机变量ζ1，ζ2，…，ζn相互独立， 且ζi～U(0,1) {f(ξi)}，i=1，2，…，n 相互独立同分布 计算机模拟全文共53页，当前为第34页。 由(强)大数定律知 以概率为1 成立 当n 足够大时，得近似公式： 计算机模拟全文共53页，当前为第35页。 注： 平均值法本质上是用样本平均值作为 总体数学期望的估计。 二. 蒙特卡罗模拟试验次数的确定 M－C 模拟是一种试验近似方法 , 试验次数