Monte Carlo模拟 引言(introduction) 均匀随机数的产生(Random number generation)

Presentation on theme: "Monte Carlo模拟 引言(introduction) 均匀随机数的产生(Random number generation)"— Presentation transcript:

1 Monte Carlo模拟 引言(introduction) 均匀随机数的产生(Random number generation)
任意分布的随机变量的抽样 Monte Carlo积分法 常用Monte Carlo模拟软件的使用 2019/7/13 Monte Carlo模拟

2 Monte Carlo模拟 第一章 引言 （Introduction) Monte Carlo方法 Monte Carlo方法简史
2019/7/13 Monte Carlo模拟

3 1.Monte Carlo方法 Monte Carlo方法： 亦称统计模拟方法（statistical simulation method）
一种采用统计抽样理论近似地求解物理或数学问题的方法 2019/7/13 Monte Carlo模拟

4 1.Monte Carlo方法 基本思想： 建立与待解问题相似或相关联的概率模型或概率过程，利用这种相似性把概率模型的某些特征量与待求解问题的解联系起来： 概率模型： 待求解问题的解： 随机事件的概率 定积分的值 随机变量的数学期望值 微分方程的解 对概率模型进行随机模拟或统计抽样，用所得样本得到这些特征量的估计值问题的近似解 2019/7/13 Monte Carlo模拟

5 Monte Carlo模拟 第一章 引言 （Introduction) Monte Carlo方法 Monte Carlo方法简史
2019/7/13 Monte Carlo模拟

6 2.Monte Carlo方法简史 Buffon投针实验 1768年，法国数学家Comte de Buffon利用投针实验估计的值 L d
2019/7/13 Monte Carlo模拟

7 2.Monte Carlo方法简史 Problem of Buffon’s needle:
If a needle of length l is dropped at random on the middle of a horizontal surface ruled with parallel lines a distance d > l apart, what is the probability that the needle will cross one of the lines? 2019/7/13 Monte Carlo模拟

8 2.Monte Carlo方法简史 Solution:
The positioning of the needle relative to nearby lines can be described with a random vector which has components: The random vector is uniformly distributed on the region [0,d)×[0,). Accordingly, it has probability density function 1/d. The probability that the needle will cross one of the lines is given by the integral 2019/7/13 Monte Carlo模拟

9 2.Monte Carlo方法简史 Enrico Fermi 1930年,利用Monte Carlo方法研究中子的扩散
并设计了一个Monte Carlo机械装置，Fermiac,用于计算核反应堆的临界状态 2019/7/13 Monte Carlo模拟

10 2.Monte Carlo方法简史 Stanislaw Ulam (1909-1984)
S. Ulam is credited as the inventor of Monte Carlo method in 1940s, which solves mathematical problems using statistical sampling. 2019/7/13 Monte Carlo模拟

11 Nicholas Metropolis (1915-1999)
2.Monte Carlo方法简史 Nicholas Metropolis ( ) The algorithm by Metropolis (and A Rosenbluth, M Rosenbluth, A Teller and E Teller, 1953) has been cited as among the top 10 algorithms having the "greatest influence on the development and practice of science and engineering in the 20th century." 2019/7/13 Monte Carlo模拟

12 2.Monte Carlo方法简史 The Name of the Game
Metropolis coined the name “Monte Carlo”, from its gambling Casino. Monte-Carlo, Monaco 2019/7/13 Monte Carlo模拟

13 Monte Carlo模拟 第一章 引言 （Introduction) Monte Carlo方法 Monte Carlo方法简史
2019/7/13 Monte Carlo模拟

14 3.Monte Carlo模拟的应用 Monte Carlo模拟的应用： 自然现象的模拟： 宇宙射线在地球大气中的传输过程；
高能物理实验中的核相互作用过程； 实验探测器的模拟 数值分析： 利用Monte Carlo方法求积分 2019/7/13 Monte Carlo模拟

15 3.Monte Carlo模拟的应用 Monte Carlo模拟在物理研究中的作用 2019/7/13 Monte Carlo模拟

16 3.Monte Carlo模拟的应用 注意以下两点： Monte Carlo方法与数值解法的不同:
数值解法:从一个物理系统的数学模型出发,通过求解一系列的微分方程来的导出系统的未知状态; Monte Carlo方法并非只能用来解决包含随机的过程的问题: 许多利用Monte Carlo方法进行求解的问题中并不包含随机过程 例如:用Monte Carlo方法计算定积分. 对这样的问题可将其转换成相关的随机过程, 然后用Monte Carlo方法进行求解 2019/7/13 Monte Carlo模拟

17 Monte Carlo模拟 第一章 引言 （Introduction) Monte Carlo方法 Monte Carlo方法简史
2019/7/13 Monte Carlo模拟

18 4.Monte Carlo算法的主要组成部分 Monte Carlo模拟的步骤：
根据欲研究的物理系统的性质，建立能够描述该系统特性的理论模型，导出该模型的某些特征量的概率密度函数； 从概率密度函数出发进行随机抽样，得到特征量的一些模拟结果； 对模拟结果进行分析总结，预言物理系统的某些特性。 2019/7/13 Monte Carlo模拟

19 4.Monte Carlo算法的主要组成部分 Monte Carlo算法的主要组成部分 概率密度函数(pdf)
必须给出描述一个物理系统的一组概率密度函数; 随机数产生器 能够产生在区间[0,1]上均匀分布的随机数 抽样规则 如何从在区间[0,1]上均匀分布的随机数出发,随机抽取服从给定的pdf的随机变量; 模拟结果记录 记录一些感兴趣的量的模拟结果 2019/7/13 Monte Carlo模拟

20 4.Monte Carlo算法的主要组成部分 误差估计 必须确定统计误差（或方差）随模拟次数以及其它一些量的变化； 减少方差的技术
利用该技术可减少模拟过程中计算的次数； 并行和矢量化 可以在先进的并行计算机上运行的有效算法 2019/7/13 Monte Carlo模拟