用蒙特卡洛法求圆周率pi

思路

正方形的面积为：
圆的面积为：
圆的面积比上正方形的面积为：π / 4
使用蒙特卡洛法在正方形内随机撒点，落在圆内的点 / 落在正方形内的点（全部的点），就约等于圆的面积 / 正方形的面积 = π / 4

代码实现以第一象限的1/4圆为例。

代码

import random
 
def count_pi(n):
    #这里用1/4圆写代码逻辑
    i = 0
    count = 0
    # n 为传入的总点数量
    while i < n:
        # 随机产生x,y坐标
        # random.random()用于生成一个0到1的随机符点数: 0 <= n < 1.0
        x = random.random()
        y = random.random()
        # 如果x平方 + y平方 < 1，说明在圆内
        if (pow(x, 2) + pow(y, 2)) < 1:
            count += 1
        i += 1
    # π的值为：4 * （落在圆内的点/总的点）
    return 4 * (count / n)
 
 
pi = count_pi(100000)
print(pi)

运行结果

运行100次，取平均值。

蒙特·卡罗方法（Monte Carlo method），也称统计模拟方法，一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。