蒙特卡罗方法是一种计算方法。原理是通过大量随机样本,去了解一个系统,进而得到所要计算的值。最经典的莫过于圆周率π的计算了。
圆周率π的计算
单位圆面积为π,并且内切与边长为2的正方形,用随机数,10000个点往下面投,落到圆里面的概率乘以4就是π的值
说明:等式左边,是圆与正方形面积之比;等式右边,是落在单位圆内的点与落在正方形内的点的比值。二者相等。
分析
1、以圆心为原点,建立二维坐标系;
2、随便投一个点,这个记为(x, y),落在圆内,则到圆心的距离 x2 + y2 <= 1;
3、投点10000次,落入圆内次数n,π值为n / 10000 * 4;
Python代码
分析好了,代码就很容易实现了,