蒙特卡洛数值模拟-计算定积分的两种方法

原理

法1（撒米法）计算积分

1.他的逻辑原理其实很简单，对于函数图像围成的面积，我们总可以用一个矩形覆盖住，然后在上面撒下一把米（随机抽样），然后去数，那些米粒落在了函数图像围成的面积内，将得到的数据除以米粒的总数，得到一个积分区域面积与矩形面积的比例。这时，矩形的面积好算，将其乘上比例，就得到了积分区域面积。

法2（数学期望）计算积分

1.数学推导

代码

# 基于蒙特卡洛求定积分的两种方法
# 法1
decisionCondition = function(a){
  a[2]-log(a[1])
}
# 用于判断在面积内的函数

n = 100000
x = runif(n,5,9)
y = runif(n,0,log(9))
A = cbind(x,y)
A
b = apply(A, 1, decisionCondition)
# length(which(b<0))
mean(b<0)*4*log(9) # 注意这里，将算出来的面积乘上总面积得到积分区域面积

# 法2
n=10000
x = runif(n,5,9)
mean(log(x))*4 # 乘以（b-a）