蒙特卡洛方法是一种用野蛮粗暴的蛮力对抗精致数学的一种计算思维,能够将复杂数学问题转化为简单粗暴的重复步骤,在工程上有很多应用。我还用蒙特卡洛方法计算了圆周率,请看我另一篇博客。
张子豪 同济大学
蒙特卡洛方法计算自然常数e
蒙特卡洛方法计算自然常数e
原理
曲面四边形面积即为积分之后的值
原理:用蒙特卡洛方法随机在左图矩形方格中撒点,统计y=1/x内外点的个数,
由几何概型,估算出曲线下曲面四边形的面积。
同时,由定积分可算出这部分面积为ln2,即e**(估算出的面积)== 2,即可求出e。
撒点越多,e的计算值也越来越趋近2.71828的真实值。
源代码
# 张子豪 2019-3-14
import random
import matplotlib.pyplot as plt
import numpy as np
DARTS = 1024*1024 # 总共撒点的