控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 蒙特卡罗计算机模拟应用一例.PDFVIP

第 16 卷 第 4 期 大 　学 　物 　理 　实 　验 　 Vol16 No4 2003 年 12 月出版 PHYSICAL EXPERIMENT OF COLLEGE Dec 2003 ( ) 文章编号 :1007 - 2934 2003 04 - 0059 - 02 蒙特卡罗计算机模拟应用一例 聂士忠 (石油大学 ,东营 ,257062) 摘 　要 　通过圆周率的计算 ,说明蒙特卡罗计算机模拟方法的应用过程。 关键词 　蒙特卡罗 ;计算机模拟 中图分类号 　O4 - 39 　　　　文献标识码 :A 1 　投针试验 著名的Buffon 投针实验是一种求 π近似值的方法 ,该方法 是在平坦桌面上划一组相距为 l 的平行线 ,然后向桌面随意地 α 投掷长度为 l 的细针 ,设细针与平行线的垂直方向的夹角为 , α l ·| cos | α αχ 则细针与平行线相交的概率为 = | cos | 。由于 是在 l 1 π π [ 0 , ] 间均匀分布的, 所以细针与平行线相交的概率等于 ∫| 0 π 图 1 　投针示意图 2 α α cos | d = π。设进行了 N 次投针试验, 有 M 次与平行线相交, 当 N 足够大时, 细针与平行线相交的频率就等于以上的概率, 于是得到计算 π的公式为 2 N π = M 2 　计算方法和结果 21 　模型 在二维平面上画三条相距为 05 、长度为 l 的平行线, 取细针长度为05 , 如图 1 所示, 因为本文所用随机数在[ 0 , 1] 之间, 所以平行线的有效长度为 l , 也就是说所有投掷试验 都等效于在上述边长为 l 的正方形区域内进行的。 22 　算法 ( ) 1 为计算作准备, 取总投针次数初始值 N = 0 , 相交次数 M = 0 , 设定总投针试验次 数 Nmax 。 ( ) 2 由蒙特卡罗方法产生两个[ 0 , 1] 间均匀分布的随机数, 并作为随机构造的细针的 收稿日期 :2003 - 01 - 20 —59 — © 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. ( ) 一个端点的坐标 x 1 、y 1 。 ( ) 3 再产生一个随机数, 作为细针另一端点的横坐标 x2