蒙特卡洛算法详讲.pdf资源-CSDN文库

版权申诉

30 浏览量 2022-06-17 23:48:08 上传评论收藏 773KB PDF 举报

资源推荐

资源详情

资源评论

Monte Carlo 法

§8.1 概述

Monte Carlo 法不同于前面几章所介绍的确定性数值方法，它是用来解决数

学和物理问题的非确定性的（概率统计的或随机的）数值方法。Monte Carlo 方

法（MCM），也称为

统计试验方法

，是理论物理学两大主要学科的合并：即随机

过程的概率统计理论（用于处理布朗运动或随机游动实验）和位势理论，主要是

研究均匀介质的稳定状态

[1]

。它是用一系列随机数来近似解决问题的一种方法，

是通过寻找一个概率统计的相似体并用实验取样过程来获得该相似体的近似解

的处理数学问题的一种手段。运用该近似方法所获得的问题的解 in spirit 更接近

于物理实验结果，而不是经典数值计算结果。

普遍认为我们当前所应用的 MC 技术，其发展约可追溯至 1944 年，尽管在

早些时候仍有许多未解决的实例。MCM 的发展归功于核武器早期工作期间 Los

Alamos（美国国家实验室中子散射研究中心）的一批科学家。Los Alamos 小组

的基础工作刺激了一次巨大的学科文化的迸发，并鼓励了 MCM 在各种问题中的

应用

[2]-[4]

。“Monte Carlo”的名称取自于 Monaco（摩纳哥）内以赌博娱乐而闻名

的一座城市。

Monte Carlo 方法的应用有两种途径：仿真和取样。仿真是指提供实际随机

现象的数学上的模仿的方法。一个典型的例子就是对中子进入反应堆屏障的运动

进行仿真，用随机游动来模仿中子的锯齿形路径。取样是指通过研究少量的随机

的子集来演绎大量元素的特性的方法。例如，

f (x)

在

a  x  b

上的平均值可以

通过间歇性随机选取的有限个数的点的平均值来进行估计。这就是数值积分的

Monte Carlo 方法。MCM 已被成功地用于求解微分方程和积分方程，求解本征

值，矩阵转置，以及尤其用于计算多重积分。

任何本质上属随机组员的过程或系统的仿真都需要一种产生或获得随机数

的方法。这种仿真的例子在中子随机碰撞，数值统计，队列模型，战略游戏，以

及其它竞赛活动中都会出现。Monte Carlo 计算方法需要有可得的、服从特定概

率分布的、随机选取的数值序列。

§8.2 随机数和随机变量的产生

[5]－[10]全面的论述了产生随机数的各类方法。其中较为普遍应用的产生随

机数的方法是选取一个函数

g(x)

，使其将整数变换为随机数。以某种方法选取

，并按照

k 1

 g(x

)

产生下一个随机数。最一般的方程

g(x)

具有如下形式：

（8.1）

其中



初始值或种子（

 0

）

g(x)  (ax  c) mod m

a 

乘法器（

a  0

）

c 

增值（

c  0

）

m 

模数

对于

数位的二进制整数，其模数通常为

。例如，对于 31 位的计算机

即可

取

311

。这里

, a

和

都是整数，且具有相同的取值范围

m  a, m  c, m  x

。

所需的随机数序





便可由下式得

n1

 (ax

 c) mod m

（8.2）

该序列称为

线性同余序列

。例如，若

 a  c  7

且

m  10

，则该序列为

7 ， 6 ， 9 ， 0 ， 7 ， 6 ， 9 ， 0……

（8.3）

可以证明，同余序列总会进入一个循环套；也就是说，最终总会出现一个无休止

重复的数字的循环。（8.3）式中序列周期长度为 4。当然，一个有用的序列必是

具有相对较长周期的序列。许多作者都用术语

乘同余法

和

混合同余法

分别指代

c  0

和

c  0

时的线性同余法。选取

, a, c

和

的法则可参见[6,10]。

这里我们只关心在区间

(0,1)

内服从均匀分布的随机数的产生。用字符

来

表示这些数字，则由式（8.2）可得

U 

（8.4）

这样

仅在数组



0,1/ m,2 / m,......, (m 1) / m



中取值。（对于区间（0,1）内的随机

数，一种快速检测其随机性的方法是看其均值是否为 0.5。其它检测方法可参见

[3,6]。）产生区间

(a,b)

内均匀分布的随机数

，可用下式

X  a  (b  a)U

（8.5）

用计算机编码产生的随机数（利用式（8.2）和（8.4））并不是完全随机的；

事实上，给定序列种子，序列的所有数字

都是完全可预测的。一些作者为强调

这一点，将这种计算机产生的序列称为

伪随机数

。但如果适当选取

a, c

和

，序

列

的随机性便足以通过一系列的统计检测。它们相对于真随机数具有可快速产

生、需要时可再生的优点，尤其对于程序调试。

Monte Carlo 程序中通常需要产生服从给定概率分布

F(x)

的随机变量

。

该步可用[6]，[13]－[15]中的几种方法加以实现，其中包括

直接法

和

舍去法

。

直接法（也称反演法或变换法），需要转换与随机变量

蒙特卡洛算法详讲.pdf

蒙特卡洛算法详讲.docx

蒙特卡洛算法-基于光的波粒二象性猜想的数学仿真.pdf

MATLAB算法-马尔可夫链蒙特卡洛算法详解，附代码.pdf

基于禁忌搜索的车联网蒙特卡洛定位算法.pdf

大数据-算法-波函数蒙特卡洛方法在量子开放系统中的应用.pdf

数学建模-蒙特卡洛算法资料合集.rar

WSN定位蒙特卡洛盒MCB算法MATLAB源码.pdf

基于Matlab_Simulink的杀爆弹蒙特卡洛飞行落点仿真.pdf

神经网络的反过程——蒙特卡洛实现想象力.pdf

基于蒙特卡洛算法的电动汽车有序充电控制策略研究.pdf

基于蒙特卡洛的医学图像重建体积计算算法GPU加速研究.pdf

基于蒙特卡洛算法的电动汽车充电需求负荷研究.pdf

基于蒙特卡洛仿真和并行粒子群优化算法的携行备件优化.pdf

缺失数据处理的期望-极大化算法与马尔可夫蒙特卡洛方法.pdf

中国新能源汽车动力电池报废量预测与对策建议——基于蒙特卡洛模拟的测算.pdf

Matlab常用算法大集合.zip

大作业基于python蒙特卡洛算法实现黑白棋MiniAlphaGo.zip

蒙特卡洛方法及应用.pdf

算法大全常用数值算法MATLAB数学建模算法灰色算法蒙特卡洛神经网络图论算法遗传算法资料大集合.zip

使用卡尔曼滤波修正蒙特卡洛算法的电动汽车充电负荷预测研究.pdf

云计算-直接模拟蒙特卡洛计算的并行算法研究.pdf

蒙特卡洛算法模拟随机数代码.txt

围棋AI，基于蒙特卡洛树搜索(MCTS)算法 Java版.zip

基于Matlab的弹道蒙特卡洛仿真研究.pdf

第十七届全国大学生智能车竞赛技术报告

32个uniapp项目源码 涵盖商城团购等

反编译脚本解包 wxapkg

video downloadhelper 合作应用2.0.10版本.

CSDN积分的获取方式

最新资源

32个uniapp项目源码涵盖商城团购等