[数学]蒙特卡罗积分方法
- 格式:ppt
- 大小:880.00 KB
- 文档页数:82
蒙特卡洛模拟法求积分1. 引言蒙特卡洛模拟法是一种基于随机采样的数值计算方法,被广泛应用于求解各种数学问题。
其中之一便是利用蒙特卡洛模拟法求解积分。
本文将介绍蒙特卡洛模拟法的基本原理、步骤以及在求解积分中的应用。
2. 蒙特卡洛模拟法基本原理蒙特卡洛模拟法以概率统计为基础,通过生成大量的随机样本来近似计算一个问题的解。
其基本原理可以概括为以下几个步骤:•随机生成样本:根据问题的要求,生成符合一定概率分布的随机样本。
•计算函数值:将每个随机样本代入目标函数中进行计算,得到对应的函数值。
•统计平均:对所有函数值进行求和并取平均,得到近似解。
3. 求解积分的蒙特卡洛模拟法步骤在使用蒙特卡洛模拟法求解积分时,需要按照以下步骤进行操作:步骤1:确定积分范围需要明确要求解的积分范围。
假设要求解的积分为∫f(x)dx,其中x的范围从a到b。
步骤2:确定随机样本生成规则根据积分范围确定随机样本生成规则。
可以使用均匀分布或其他概率分布来生成随机样本,确保样本覆盖整个积分区间。
步骤3:生成随机样本使用确定的随机样本生成规则,生成足够数量的随机样本。
通常情况下,生成的样本数越多,计算结果越接近真实值。
步骤4:计算函数值将每个随机样本代入目标函数f(x)中进行计算,得到对应的函数值。
这相当于在积分区间上进行采样,并计算采样点处的函数值。
步骤5:统计平均对所有函数值进行求和并取平均,得到近似解。
根据大数定律,当样本数量充足时,平均值将趋近于真实解。
4. 蒙特卡洛模拟法求解积分示例以下是一个使用蒙特卡洛模拟法求解积分的示例:假设要求解的积分为∫x^2dx,积分范围为0到1。
步骤1:确定积分范围。
积分范围为0到1。
步骤2:确定随机样本生成规则。
使用均匀分布生成随机样本。
步骤3:生成随机样本。
生成足够数量的随机样本,例如10000个。
步骤4:计算函数值。
将每个随机样本代入目标函数f(x)=x^2中进行计算,得到对应的函数值。
步骤5:统计平均。
蒙特卡洛类方法
蒙特卡洛方法是一类随机化的计算方法,主要应用于求出高维度空间中的定积分或概率分布的特性。
该方法以随机样本为基础,通过大量生成且符合某种分布律的随机数,从中抽取样本,利用样本的统计性质来计算近似解。
常见的蒙特卡洛方法包括:
1.随机模拟法
在数学建模、广告投放、经济预测等领域,随机模拟(也称蒙特卡罗方法)已经成为了一个重要的工具。
其基本思想是,系统表现出的某些规律和性质可以用随机过程进行模拟和预测。
2.随机游走算法
随机游走是一种基于随机过程的数值计算算法,通过简单的偏随机移动来解决复杂问题,被广泛应用于物理、化学、生物学、金融等领域。
随机游走算法的核心思想是通过随机漫步遍历所有可能的状态,找到最终解。
3.马尔可夫链蒙特卡罗方法
马尔可夫链蒙特卡罗方法(MCMC)是一种近似随机模拟算法,用于计算高维空间中的积分和概率分布。
这种方法通过构造一个马尔可夫链来模拟复杂的概率
分布,并通过观察链的过程来获得所求的统计量。
4.重要性采样
重要性采样是一种通过迭代抽样来估算积分值或概率分布的方法。
它的基本思想是利用不同的概率分布来采样目标分布中的样本,从而增加目标分布中采样到重要样本的概率,从而提高采样的效率。
总之,蒙特卡洛方法在物理学、统计学、金融学、计算机科学、生物科学等众多领域都有广泛的应用,是一种很实用的工具。
蒙特卡洛法的基本原理蒙特卡洛法(Monte Carlo method)是一种基于随机抽样的数值计算方法,用于解决难以通过解析方法或传统数学模型求解的问题。
它在物理学、化学、工程学、计算机科学、金融学、生物学等领域都有广泛应用。
本文将介绍蒙特卡洛法的基本原理,包括随机数生成、统计抽样、蒙特卡洛积分、随机漫步等方面。
一、随机数生成随机数是蒙特卡洛法中的基本元素,其质量直接影响着计算结果的准确性。
随机数的生成必须具有一定的随机性和均匀性。
常见的随机数生成方法有:线性同余法、拉斯维加斯法、梅森旋转算法、反序列化等。
梅森旋转算法是一种广泛使用的准随机数生成方法,其随机数序列的周期性长、随机性好,可以满足大多数应用的需要。
二、统计抽样蒙特卡洛法利用抽样的思想,通过对输入参数进行随机取样,来模拟整个系统的行为,并推断出某个问题的答案。
统计抽样是蒙特卡洛方法中最核心的部分,是通过对概率分布进行样本抽取来模拟随机事件的发生,从而得到数值计算的结果。
常用的统计抽样方法有:均匀分布抽样、正态分布抽样、指数分布抽样、泊松分布抽样等。
通过对这些概率分布进行抽样,可以在大量随机取样后得到一个概率分布近似于输入分布的“抽样分布”,进而求出所需的数值计算结果。
三、蒙特卡洛积分蒙特卡洛积分是蒙特卡洛法的重要应用之一。
它利用统计抽样的思想,通过对输入函数进行随机抽样,计算其随机取样后的平均值,来估算积分的值。
蒙特卡洛积分的计算精度与随机取样的数量、抽样分布的质量等因素有关。
蒙特卡洛积分的计算公式如下:$I=\frac{1}{N}\sum_{i=1}^{N}f(X_{i})\frac{V}{p(X_{i})}$$N$为随机取样的数量,$f(X_{i})$为输入函数在点$X_{i}$的取值,$V$为积分区域的体积,$p(X_{i})$为在点$X_{i}$出现的抽样分布的概率密度函数。
通过大量的样本拟合,可以估算出$I$的值接近于真实积分的值。
用蒙特卡洛方法计算积分简介蒙特卡洛方法是一种通过随机抽样来计算数学问题的方法。
在计算积分时,蒙特卡洛方法可以提供一种简单而有效的解决方案。
方法步骤1. 确定积分范围:首先确定要计算的积分范围,并将其表示为一个多维的定积分。
2. 创建随机点:生成一组随机点,这些随机点需要在积分范围内均匀分布。
3. 判断点的位置:对于每个随机点,判断它是否在被积函数的曲线下方。
4. 计算积分值:计算在被积函数下方的点数与总随机点数的比例,并乘以积分范围的体积,得到积分的近似值。
优势和注意事项蒙特卡洛方法的优势在于其简单性和适用性广泛性。
然而,在使用蒙特卡洛方法进行积分计算时,需要注意以下几点:- 随机点的数量:随机点的数量越多,计算结果越精确,但计算时间也会增加。
- 积分范围的选择:选择合适的积分范围可以提高计算效率和准确性。
- 随机点的生成:生成随机点需要遵循均匀分布原则,以确保计算结果的准确性。
示例以下是使用蒙特卡洛方法计算积分的示例代码:import randomdef monte_carlo_integration(f, a, b, n):count = 0for _ in range(n):x = random.uniform(a, b)y = random.uniform(min(f(a), f(b)), max(f(a), f(b)))if 0 < y <= f(x):count += 1return count / n * (b - a) * (max(f(a), f(b)) - min(f(a), f(b)))def f(x):被积函数定义,根据实际情况修改return x**2a = 0 # 积分下限b = 1 # 积分上限n = # 随机点数量result = monte_carlo_integration(f, a, b, n)print("Approximate integral value:", result)注意:上述代码仅为示例,实际运行时请根据需要修改被积函数和参数。
蒙特·卡罗(Monte Carlo)法是一种统计模拟方法,通常是利用随机数来解决一些数值计算问题,本文要讲的就是利用蒙特·卡罗方法来求解数值积分。
基本思路首先我们知道定积分其实就是一个面积,将其设为I,现在我们就是要求出这个I。
我们的想法是通过在包含定积分的面积为S的区域(通常为矩形)内随机产生一些随机数,其数量为N,再统计在积分区域内的随机数,其数量为i,则产生的随机数在积分区域内的概率为iN,这与积分区域与总区域面积的比值IS应该是近似相等的,我们利用的就是这个关系,即IS≈iN最后即得所求定积分算式为:I=iNS代码部分有了上面的铺垫,我们就可以来写MATLAB代码了。
我们要求的定积分为∫0πsinxdx.对于上述积分我们很容易可以得到其解析解为2,下面我们来看用蒙特·卡罗方法得到的结果,输入代码% Monte Carlo% 蒙特卡洛法求定积分clearN = 1e4;x_min = 0; x_max = pi;f = @(x) sin(x);xx =x_min:0.01:x_max;x = x_min + (x_max-x_min)*rand(N,1);y_min = min(f(xx)); y_max = max(f(xx));y = y_min +(y_max-y_min)*rand(N,1);i = y < f(x);I = sum(i)/N*(x_max-x_min)*(y_max-y_min);% 画图plot(x,y,'go',x(i),y(i),'bo')axis([x_min x_max y_min y_max])hold onplot(xx,f(xx),'r-','LineWidth',2)。
蒙特卡洛积分法
蒙特卡洛积分法,是一种通过随机抽样计算数学积分的方法。
该方法可以用于解决各种类型的数学积分问题。
下面是详细的蒙特卡洛积分法的讲解。
一、定义
蒙特卡洛积分法是一种基于随机样本来估算积分的方法。
其基本思路是选取一组随机数(通常服从均匀分布),将这些数代入要求积分的函数中,并计算函数值总和,再将总和除以样本数量,得到积分的近似值。
由于蒙特卡洛积分法具有极高的灵活性,且无需求解复杂的数学方程,因此被广泛应用于各个领域的数学积分求解。
二、原理
蒙特卡洛积分法是通过概率统计的方法来近似计算积分。
具体方法是根据随机数的生成进行采样,将采样得到的样本代入积分式中,然后求和、取平均值得到近似积分值。
三、流程
1. 确定计算积分的区域和积分函数;
2. 通过随机抽样生成样本点;
3. 计算每个样本点对应的函数值;
4. 将所有函数值求和,再除以样本点数量,得到积分近似值。
四、优点
1. 适用性强:蒙特卡洛积分法对待求解的函数类型没有限制,而且求解的范围可以是多维的;
2. 精度高:积分结果的误差随着样本数的增加而减小;
3. 随机性强:由于样本点的随机分布,可以避免因为选取的样本点不具有代表性而导致的误差;
4. 可并行性强:蒙特卡洛法具有明显的并行性,可以充分利用现代计算机的多核处理器。
五、缺点
1. 运行速度慢:由于需要生成随机数,蒙特卡洛积分法的计算速度相对于传统数值积分较慢;
2. 精度依赖于样本数目:蒙特卡洛积分法的精度与样本数量有直接关系,因此需要大量样本数目;
3. 对函数的光滑度要求高:蒙特卡洛积分法在求解非光滑函数时精度会受到影响。
蒙特卡洛方法蒙特卡洛方法是一种以随机数代替确定性答案的方法,用来解决难以用传统数学方法求解的计算问题。
它的名字来自于摩纳哥的蒙特卡洛市,因为在二战时期,美国的原子弹计划曾在那里进行过试验。
现在,蒙特卡洛方法已经广泛应用于各种领域,包括统计学、计算机科学、物理学、金融等。
我们来举一个简单的例子来解释蒙特卡洛方法的基本原理。
假设我们要计算正方形中圆的面积,这个问题可以用传统的数学方法求解,而且结果是$π/4$。
但是,如果我们用蒙特卡洛方法求解这个问题,我们可以在正方形中随机生成很多点,并统计其中多少点在圆内。
如果我们生成的点足够多,那么圆内点的数量与总点数的比例就可以近似表示圆的面积与正方形面积之比,也就是$π/4$。
这种方法的优点在于,我们不需要事先知道圆的半径或面积,只需要用随机数模拟出圆内外的点,就可以得到一个近似的答案。
当然,随机生成的点的数量越多,计算结果就越精确。
蒙特卡洛方法的应用非常广泛,下面介绍几个例子:1. 在金融领域,蒙特卡洛方法被用来计算复杂的金融衍生品的价格。
金融衍生品是一种金融工具,其价值的变化受到其他金融资产的价格波动的影响。
这些衍生品的价格无法用传统的数学方法精确计算,因为它们涉及到多种不确定因素,如未来市场价格的波动、利率和货币汇率的变化等。
利用蒙特卡洛方法,可以在一个随机模拟的框架下,通过很多次模拟(通常是几千次)来计算期权的价格和各种可能结果出现的概率。
这些结果可以用来帮助投资者评估一种衍生品的实际价值。
2. 在科学计算中,蒙特卡洛方法可以用来求解很多复杂的数学问题,如高维积分、求解微分方程、求解偏微分方程等。
一个著名的例子就是蒙特卡洛积分法,它可以用来求解高维积分。
在这种方法中,我们首先确定积分范围(即多维空间中的一个区域),然后在这个区域中随机生成很多点,最后根据这些点的分布来估计积分的大小。
蒙特卡洛积分法的优点在于,它适用于复杂的积分问题,且收敛速度比传统的数值积分方法要快得多。
第五章蒙特卡洛方法在机器学习和强化学习中,蒙特卡洛方法是一类基于随机抽样的方法,用于估计未知概率分布的特征或求解复杂的问题。
在本章中,我们将介绍蒙特卡洛方法的基本原理和应用领域。
1.蒙特卡洛方法的原理蒙特卡洛方法是通过利用随机抽样的规律来估计未知概率分布的特征。
其基本原理如下:(1)随机抽样:根据已知概率分布进行随机抽样,得到一系列样本。
(2)样本推断:利用得到的样本进行统计推断,从而估计未知概率分布的特征。
(3)结果评估:通过对估计结果进行评估,得到对未知概率分布的特征的估计值。
2.蒙特卡洛方法的应用领域蒙特卡洛方法广泛应用于估计数学问题、求解优化问题以及模拟高维空间中的复杂系统。
以下是一些蒙特卡洛方法的应用领域的示例:(1)数值计算:蒙特卡洛方法可以用于计算复杂的数学问题,如计算积分、求解微分方程等。
通过随机抽样和统计推断,可以得到对问题的近似解。
(2)优化问题:蒙特卡洛方法可以用于求解优化问题,如最大化或最小化函数的值。
通过随机抽样和统计推断,可以找到函数的全局最优解或局部最优解。
(3)统计推断:蒙特卡洛方法可以用于估计未知概率分布的特征,如均值、方差、分位数等。
通过随机抽样和统计推断,可以得到这些特征的近似值。
(4)模拟与优化:蒙特卡洛方法可以用于模拟高维空间中的复杂系统,如金融市场、交通网络等。
通过随机抽样和统计推断,可以对系统的行为进行建模和优化。
3.蒙特卡洛方法的算法步骤蒙特卡洛方法的算法步骤如下:(1)随机抽样:根据已知概率分布进行随机抽样,得到一系列样本。
(2)样本推断:利用得到的样本进行统计推断,从而估计未知概率分布的特征。
常见的推断方法有样本平均法、样本方差法等。
(3)结果评估:通过对估计结果进行评估,得到对未知概率分布的特征的估计值。
常见的评估方法有置信区间估计、假设检验等。
4.蒙特卡洛方法的优缺点蒙特卡洛方法具有以下优点:(1)简单易实现:随机抽样和统计推断是蒙特卡洛方法的基本步骤,易于理解和实现。
蒙特卡罗方法计算定积分蒙特卡罗方法(Monte Carlo Method)是一种通过运用随机数抽样,利用计算机模拟实验来求解复杂问题的方法。
主要是利用概率统计中大数定律来求解,它可以帮助我们在数学上求解积分,快速准确的进行估计,被用在金融学,计算物理学,经济学,数学统计,技术分析等诸多领域。
一、蒙特卡罗方法计算定积分1、估计问题定义:计算一个未知函数积分。
2、随机抽样:由蒙特卡罗定理,我们在总体样本里进行随机抽样,随机抽取的样本可用与能够准确估计积分的概率多样性。
首先,创建一个具有期望值的随机变量,即函数中的随机变量,有且只有一个。
3、抽样的选择:抽样的选择非常重要,随机抽样的样本数量要远大于正常的定积分计算过程中要求的样本数量。
4、统计估计:通过蒙特卡罗方法,估计积分就是所抽样本的函数值的平均值乘以定积分范围。
二、蒙特卡罗方法容易出现的问题1、抽样样本量不够:如果抽样样本量不够,结果会出现较大误差,蒙特卡罗方法所估计值将可能与实际结果存在较大偏差。
2、估计值不够稳定:蒙特卡罗方法产生的结果一般存在很大的变动,估计值的结果可能会出现很大的波动,如果这种情况发生,就要调整抽样数量来达到稳定的结果。
3、结果不精确:由于蒙特卡罗方法依赖于随机样本,对精确度的要求很高,如果抽样数量不够,很可能出现精度较差的情况。
三、蒙特卡罗方法计算定积分的优点1、随机变量可定义:由于蒙特卡罗方法是基于随机变量的,所以可以通过定义方法来求出任意函数的定积分。
2、结果准确:在合理的的抽样数量下,蒙特卡罗方法的估计结果都基本准确。
3、实用性强:蒙特卡罗方法不仅实用于算法应用,还可以用于复杂估计。
四、总结蒙特卡罗方法是一种基于随机变量的,主要用于求解数学和经济类问题的方法。
它具有计算定积分快速准确、估计值结果可靠、实用性强等优点,是复杂问题求解的重要工具。
但同时也存在诸如抽样数量不足、估计值不稳定和精度较低等问题,因此在使用时要醉倒,确保估计结果的准确性。
蒙特卡罗(MonteCarlo)方法算积分❝蒙特卡罗(Monte Carlo)是摩纳哥最著名的一区,以豪华的赌场闻名于世,用它作为名字大概是因为随机性,就像赌博场里面的扔骰子的过程。
最早的「蒙特卡罗方法」是为了解决一些难求解的积分问题。
❞•「问题」•「蒙特卡洛方法」如果可以选择在的概率分布函数,则有:若在之间是均匀分布时,即,那么:这就是之前讲解的平均值法(点击跳转),另外随机投点法(点击跳转)也是「蒙特卡洛方法」. 一般均匀分布并不是好选择,因为如果在有不少点使得,那么这些点对的近似计算贡献很小,所以应尽可能少用这些点. 此时就需要采用「重要采样方法」选择合适的,从而提高精度,这部分内容我们后续会详细阐述,这次我们先分析「随机投点法」和「平均值法」的随机误差.•「误差分析」(1)「随机投点法」令且,则 iid . 由中心极限定理知:从而所以因此的随机误差为:.(2)「平均值法」由中心极限定理知:其中因此的随机误差为:,但其渐近方差更小.类似的,计算高维定积分的蒙特卡罗方法的随机误差也为,所以蒙特卡罗方法计算积分和维数关系不大,但数值积分则存在「维数诅咒」问题,这也是蒙特卡罗方法的「优势」.•「高维积分算例」「以下为Python代码」import numpy as npfrom scipy import integrate## (x1)^2(x2)^2(x3)^2 在 [0,1] 的积分a1,b1 = 0,1a2,b2 = 0,1a3,b3 = 0,1# 三重积分计算def f(x1,x2,x3):return x1**2 * x2**2 * x3**2I_exact, Error = integrate.tplquad(f,a1,b1,a2,b2,a3,b3)# 平均值法N = 10000x1_sample = a1 + (b1-a1)*np.random.rand(N)x2_sample = a2 + (b2-a2)*np.random.rand(N)x3_sample = a3 + (b3-a3)*np.random.rand(N)np.random.seed(1)h_x = f(x1_sample,x2_sample,x3_sample)I_approx_stat = (b3-a3)*(b2-a2)*(b1-a1)/N*np.sum(h_x)# 数值积分M = 200h1 = (b1-a1)/(M-1)h2 = (b2-a2)/(M-1)h3 = (b3-a3)/(M-1)x1 = np.linspace(a1,b1,M)x2 = np.linspace(a2,b2,M)x3 = np.linspace(a3,b3,M)x1_mesh, x2_mesh, x3_mesh = np.meshgrid(x1,x2,x3)I_approx_rec = np.sum( f(x1_mesh, x2_mesh, x3_mesh)*h1*h 2*h3 )print( '多重积分值:', I_exact )print( '\n平均值法结果:', I_approx_stat )print( '\n数值积分结果:', I_approx_rec )❝多重积分值:0.037037037037037035平均值法结果:0.03737256369148107数值积分结果:0.03788231093787493(大家可尝试画出:不同数量采样点对应的结果和真实值之间的关系图)❞。
蒙特卡洛方法是一种基于随机抽样的数值计算方法,广泛应用于金融学、物理学、工程学和计算机科学等领域。
它的原理是通过随机抽样来估计数学模型的结果,通过大量重复实验来逼近真实值。
在本文中,我们将探讨蒙特卡洛方法的原理、应用和局限,并共享个人对这一方法的理解和观点。
1. 蒙特卡洛方法的原理蒙特卡洛方法的核心思想是利用随机数来处理问题。
它通过生成大量的随机数,利用这些随机数的统计特性来近似求解问题。
在金融衍生品定价中,我们可以使用蒙特卡洛方法来模拟股票价格的随机漫步,从而估计期权合约的价格。
通过不断模拟股票价格的变化,并计算期权合约的价值,最终得到一个接近真实值的结果。
2. 蒙特卡洛方法的应用蒙特卡洛方法在金融领域被广泛应用于期权定价、风险管理和投资组合优化等问题。
在物理学中,蒙特卡洛方法可以用于模拟粒子的运动,求解无法用解析方法求解的复杂系统。
在工程学和计算机科学中,蒙特卡洛方法可以用于求解概率分布、优化问题和模拟系统行为。
3. 蒙特卡洛方法的局限虽然蒙特卡洛方法有着广泛的应用,但也存在一些局限性。
蒙特卡洛方法通常需要大量的随机抽样,计算成本较高。
随机性导致了结果的不确定性,需要进行大量的实验才能得到可靠的结果。
蒙特卡洛方法在高维问题和高精度要求下计算效率低下,需要借助其他数值方法进行辅助。
4. 个人观点和理解个人认为蒙特卡洛方法是一种非常强大的数值计算方法,能够解决复杂问题和高维问题。
它的随机性使得结果更加贴近真实情况,有利于处理实际情况中的不确定性和风险。
但是在实际应用中,需要注意随机抽样的方法和计算成本,并且需要结合其他数值方法进行验证和辅助,以确保结果的准确性和可靠性。
总结回顾蒙特卡洛方法是一种基于随机抽样的数值计算方法,通过大量重复实验来逼近真实值。
它在金融学、物理学、工程学和计算机科学等领域有着广泛的应用。
然而,蒙特卡洛方法也存在一些局限性,需要结合其他数值方法来弥补其不足。
个人认为蒙特卡洛方法是一种强大的数值计算方法,能够处理复杂和高维问题,但在实际应用中需要注意其随机性和计算成本。
蒙特卡洛方法求定积分一蒙特卡洛方法求定积分一蒙特卡洛方法是一种通过随机采样的方式来估算数值积分的方法。
它的基本思想是利用大量随机点的分布来逼近函数曲线下的面积,从而得到定积分的近似值。
本文将介绍蒙特卡洛方法的原理和实现步骤,并以一个具体的例子来演示如何使用该方法求解定积分。
具体地说,我们可以将定积分的计算问题转化为对一个随机变量的期望值的计算问题。
设要求解的定积分为:其中,f(x)是在积分区间上的一些函数。
我们可以将其写成一个概率密度函数的形式:其中,p(x)是一个在积分区间上的非负函数,满足:而f(x)可以被视为一个可测函数,满足其在随机变量分布上的期望存在。
根据随机变量的期望值性质,我们可以将上述的积分表达式写为:其中,X是一个在积分区间上的随机变量,其概率密度函数为p(x)。
利用蒙特卡洛方法求解定积分的步骤:1.选择一个适当的概率分布p(x)作为随机变量的密度函数。
2.从概率分布p(x)中生成N个随机样本x1,x2,...,xN。
3.计算函数的值f(x1),f(x2),...,f(xN)。
4.计算随机变量的期望值:通过大量的随机样本的估计值可以逼近积分的精确值。
下面我们以一个具体的例子来演示如何使用蒙特卡洛方法求解定积分。
我们要计算的积分为:我们选择一个均匀分布作为随机变量的概率分布,其密度函数为:最后,我们可以得到定积分的估计值为:这个结果可以视为原定积分的一个近似值,而且随着样本数量的增加,估计值会更加接近积分的精确值。
总结起来,蒙特卡洛方法通过随机采样的方式来求解定积分,其核心思想是利用随机样本的分布来逼近函数曲线下的面积,从而得到定积分的近似值。
虽然蒙特卡洛方法在计算过程中存在一定的误差,但是通过增加样本数量可以提高估计结果的精度。
因此,蒙特卡洛方法广泛应用于估计各种复杂的定积分问题,尤其是对于高维积分问题,蒙特卡洛方法具有更大的优势。
蒙特卡洛(MonteCarlo)法求定积分蒙特卡洛(Monte Carlo)法是⼀类随机算法的统称。
随着⼆⼗世纪电⼦计算机的出现,蒙特卡洛法已经在诸多领域展现出了超强的能⼒。
在机器学习和⾃然语⾔处理技术中,常常被⽤到的MCMC也是由此发展⽽来。
本⽂通过蒙特卡洛法最为常见的⼀种应⽤——求解定积分,来演⽰这类算法的核⼼思想。
⽆意识统计学家法则(Law of the unconscious statistician)这是本⽂后续会⽤到的⼀个定理。
作为⼀个预备知识,我们⾸先来介绍⼀下它。
先来看⼀下维基百科上给出的解释。
In probability theory and statistics, the law of the unconscious statistician (sometimes abbreviated LOTUS) is a theorem used to calculate the 期望值 of a function g(X)" role="presentation" style="position: relative;">g(X)g(X) of a 随机变量 X"role="presentation" style="position: relative;">XX when one knows the probability distribution of X" role="presentation" style="position: relative;">XX but one does not explicitly know the distribution of g(X)" role="presentation" style="position: relative;">g(X)g(X). The form of the law can depend on the form in which one states the probability distribution of the 随机变量 X"role="presentation" style="position: relative;">XX.If it is a discrete distribution and one knows its PMF function ƒX" role="presentation"style="position: relative;">ƒXƒX (but not ƒg(X)" role="presentation" style="position: relative;">ƒg(X)ƒg(X)), then the 期望值 of g(X)" role="presentation" style="position: relative;">g(X)g(X) isE[g(X)]=∑xg(x)fX(x)" role="presentation" style="position:relative;">E[g(X)]=∑xg(x)fX(x)E[g(X)]=∑xg(x)fX(x)where the sum is over all possible values x" role="presentation" style="position: relative;">xx of X" role="presentation" style="position: relative;">XX.If it is a continuous distribution and one knows its PDF function ƒX" role="presentation">ƒXƒX (but not ƒg(X)" role="presentation">ƒg(X)ƒg(X)), then the 期望值 of g(X)"role="presentation">g(X)g(X) isE[g(X)]=∫−∞∞g(x)fX(x)dx"role="presentation">E[g(X)]=∫∞−∞g(x)fX(x)dxE[g(X)]=∫−∞∞g(x)fX(x)dxLOTUS到底表达了⼀件什么事呢?它的意思是:已知随机变量X" role="presentation"style="position: relative;">XX的概率分布,但不知道g(X)" role="presentation" style="position: relative;">g(X)g(X)的分布,此时⽤LOTUS公式能计算出函数g(X)" role="presentation"style="position: relative;">g(X)g(X)的数学期望。