蒙特卡洛方法求定积分及python实现(转)

Python017

蒙特卡洛方法求定积分及python实现(转),第1张

蒙特卡洛方法计算定积分

计算定积分

利用蒙特卡洛计算方法,核心步骤是求取随机的 g(X1),………,g(Xn),n∈[a,b],由数学期望和大数定理可以近似计算定积分,公式为

原函数:

导函数:

计算导函数在[10,15]上的定积分;

Python

用蒙特卡洛方法计算的定积分:

直接用原函数计算的定积分:

偏差程度为:

此算法诞生的背景是:

蒙特卡洛算法理论基础是概率论,实际就是暴力计算逼近理想结果。正是在以上两个背景下,它刚好得到了极大的应用和发展。

上面叙述的是定义,我来描述一个例子:蒙特卡罗法计算圆周率