- 1、本文档共34页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
5.蒙特卡罗方法
计算物理 蒙特卡罗方法 蒲丰投针(1/5) 蒲丰投针(2/5) 蒲丰投针(3/5) 蒲丰投针(4/5) 蒲丰投针(5/5) 收敛性、误差和优缺点(1/4) 收敛性、误差和优缺点(2/4) 收敛性、误差和优缺点(3/4) 收敛性、误差和优缺点(4/4) 任意分布的随机数(1/11) 任意分布的随机数(2/11) 任意分布的随机数(3/11) 任意分布的随机数(4/11) 任意分布的随机数(5/11) 任意分布的随机数(6/11) 任意分布的随机数(7/11) 任意分布的随机数(8/11) 任意分布的随机数(9/11) 任意分布的随机数(10/11) 任意分布的随机数(11/11) 粒子输运问题(1/7) 粒子输运问题(2/7) 粒子输运问题(3/7) 粒子输运问题(4/7) 粒子输运问题(5/7) 粒子输运问题(6/7) 粒子输运问题(7/7) 随机过程模拟(1/3) 随机过程模拟(2/3) 随机过程模拟(3/3) 梅氏抽样(1/1) 作业 例:a 粒子衰变的蒙特卡罗模拟 随机现象的观测数据?概率分布函数 观测数据:每隔 Dt半衰期 测量一次放射的 a 粒子数,共测量 N=2608 次,测得 k 个粒子的次数为 nk 408 5 16 10 532 4 27 9 525 3 45 8 383 2 139 7 203 1 273 6 57 0 nk k nk k 分析现象,建立模型 Dt 内的衰变数 k 是随机事件的发生次数,是随机量 k 在有限的平均值 a 上下波动 k 可看作大量( N 足够大)独立(每个原子核的衰变不受其它原子核的影响)试验的结果 事件有相同的很小的概率(原子核衰变的机会相等) k 的上述特征表明: k 近似服从泊松分布 √ 蒙特卡罗模拟 产生 2608 个泊松分布的随机数,统计其中数值等于 k 的个数并与观测值 nk 比较 √ * * 3/lesson/ComputationalPhysics 蒙特卡罗方法 蒲丰投针 收敛性、误差和优缺点 任意分布的随机数 粒子输运问题 随机过程模拟 梅氏抽样 √ 蒙特卡罗方法 又称随机抽样技巧或统计试验方法 以概率统计理论为基础的 能够比较逼真地描述事物的特点及物理实验过程 解决一般数值方法难以解决的问题 随着电子计算机的发展而发展 首先在核武器的试验与研制中得到了应用 蒲丰投针 法国数学家蒲丰的1777年出版的著作:“在平面上画有一组间距为 d 的平行线,将一根长度为 l (ld) 的针任意掷在这个平面上,求此针与平行线中任一条相交的概率。” √ 步骤 在桌面上画出间距为 2d 的平行线 准备长度为 2l (ld) 的针 向桌面随机投针 如果针与平行线相交,则计数器 n 加 1 计算:计数器 n 与总投针数 N 的比例(视作相交概率 P ) 概率分析 P = ? 各条平行线地位等同,仅考虑某条平行线附近的情况 平行线方向的 x 坐标对概率没影响 针的中点的 y 坐标在线之间等概率落入(均匀分布在 [0, d]),仅当 y?l 才可能针-线相交 针-线的夹角 q 均匀分布在 [0, p],q 与 y 独立 x y √ x y 概率 P = 2l / (p d) ,可计算圆周率 3.1795 859 2,520 0.5419 1925 Reina 3.1415929 1,808 3,408 0.83 1901 Lezzerini 3.1595 489 1,030 0.75 1884 Fox 3.137 382 600 1.0 1860 De Morgan,C 3.1554 1,218 3,204 0.6 1855 Smith 3.1596 2,532 5,000 0.8 1850 Wolf p 值 相交数 总投数 针长 时间 实验者 √ 关于蒙特卡罗方法的分析和总结 基本思想 确定所求问题的解是某事件的概率(或某随机变量的数学期望、或与概率/数学期望有关的量,如 p ) 通过试验方法,得出事件发生的频率(或该随机变量若干个具体观察值的算术平均值,如 P ),求解 数学期望与概率:当随机变量的取值仅为 1 或 0 时,它的数学期望就是事件的概率;反之亦然 数学期望与算术平均值 用随机试验的方法计算积分,将积分看作服从分布密度函数 f (r) 的随机变量 g (r) 的数学期望 通过试验,得到 N 个观察值 r1, r2, …, rN (从 f (r) 中抽取 N 个子样 r1, r2, …, rN ),求 g (r) 的算术平均值 √ 试验方法和次数 试验方法不一定可行 精确的近似解需要巨量的试验次数,但人工方法作大量的试验相当困难,甚至是不可能的 用计算机模拟随机试验过程,完成巨量的次数 抽象 x 的分布密度函数: q 的分布密度函数: 产生任
您可能关注的文档
- 2016年黑龙江省高新技术企业认定申报指南-七台河科技局.ppt
- 2016年教学工作计划.doc-临沂大学.doc
- 2015年微光夜视技术重点实验室招聘.ppt
- 2016高考英语书面表达训练专题.ppt.ppt
- 2016级新生开学典礼院长讲话-浙江工商大学信息与电子工程学院.ppt
- 20偶氮苯的光化异构体-山西大学化学实验中心.ppt
- 2016年第二季度党风廉政教育课件.ppt
- 2017年度国家自然科学基金项目资金申请预算编报三.ppt
- 2、护理诊断.ppt
- 23寇福东庞培氏病发病原理与治疗.ppt
- 《西亚》(教案及反思)-2024-2025学年鲁教版(五四学制)(2024)地理六年级下册.docx
- 防溺水七不准课件.pptx
- 2024-2025学年人教PEP版小学英语五年级下册(全册)知识点归纳(附目录).docx
- 《印度》(教案及反思)-2024-2025学年鲁教版(五四学制)(2024)地理六年级下册.docx
- 2024-2025学年人教精通版(三起)小学英语四年级下册(全册)知识点归纳.docx
- 2024-2025学年人教精通版(三起)小学英语六年级下册(全册)知识点复习要点.docx
- 2024-2025学年人教精通版(三起)小学英语五年级下册(全册)知识点归纳.docx
- 2024-2025学年外研版(三起)小学英语四年级下册(全册)知识点归纳.docx
- 2024-2025学年外研版(三起)小学英语六年级下册(全册)知识点归纳.docx
- 防溺水与自救课件.pptx
最近下载
- Hi5622V100 (海思AX3000低成本路由器wifi芯片)用户指南.pdf
- 2023年保安员笔试考前题库500题(含答案).pdf VIP
- 2024年联通智家工程师(初级)认证理论备考试题库(附答案).pdf VIP
- 2012款陆风X8_汽车使用手册用户操作图解驾驶指南车主车辆说明书电子版.doc
- 2024年中考数学压轴题型(安徽专用)专题08 解答题压轴题(几何综合(一)(学生版).docx
- 三农村劳动力转移就业职业技能培训方案.doc VIP
- 2016年大学日语专业四级考试真题.doc VIP
- 三农村劳动力转移培训指导书.doc VIP
- (完整版)军事地形学.pptx
- 联通智家工程师初级认证理论备考试题及答案.doc VIP
文档评论(0)