思路
正方形的面积为:
圆的面积为:
圆的面积比上正方形的面积为:π / 4
使用蒙特卡洛法在正方形内随机撒点,落在圆内的点 / 落在正方形内的点(全部的点),就约等于圆的面积 / 正方形的面积 = π / 4
代码实现以第一象限的1/4圆为例。
代码
import randomdef count_pi(n):#这里用1/4圆写代码逻辑i = 0count = 0# n 为传入的总点数量while i < n:# 随机产生x,y坐标# random.random()用于生成一个0到1的随机符点数: 0 <= n < 1.0x = random.random()y = random.random()# 如果x平方 + y平方 < 1,说明在圆内if (pow(x, 2) + pow(y, 2)) < 1:count += 1i += 1# π的值为:4 * (落在圆内的点/总的点)return 4 * (count / n)pi = count_pi(100000)print(pi)
运行结果
运行100次,取平均值。
蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。