首先引用百度百科的蒙特卡洛分析法介绍
蒙特卡罗分析法(Monte Carlo method)(统计模拟法),是一种采用随机抽样(Random Sampling)统计来估算结果的计算方法,可用于估算圆周率,由约翰·冯·诺伊曼提出。由于计算结果的精确度很大程度上取决于抽取样本的数量,一般需要大量的样本数据,因此在没有计算机的时代并没有受到重视。
蒙特卡洛分析法利用正方形的面积跟正方形体内最大圆的面积比值为4:π来求解π的值!
假设给正方形内部随机仍石子,那么落在圆里面的概率为四分之π,当石子数量足够大的时候就可以i精确的求出π的值了!现代利用计算机可以很容易的模拟出来。以下是代码展示.
import random #引入随机数库
N = 1000*1000 #定义循环次数
list=0.0 #落入圆中的石子数量
for i in range(1,1+N):
x = random.uniform(0,1) #模拟石子的随机坐标
y = random.uniform(0,1)
if (pow(x,2) + pow(y,2)) <= 1:
list = list+1
print(4*(list/N)) #打出圆周率