用蒙特卡洛算法编程求解Π的值

阿米修修

已于 2022-09-24 19:27:29 修改

阅读量1.1k

点赞数 3

文章标签：算法 python

于 2022-09-24 19:23:23 首次发布

本文链接：https://blog.csdn.net/weixin_53843696/article/details/127029409

版权

蒙特卡洛问题求PI值的原理

正方形的面积为4r^2
正方形内四分之一圆的面积为:Πr^2
正方形内四分之一圆的面积与正方形的面积的比值为Π/4
所以Π的值为4 * (四分之一圆的面积/正方形的面积)

在这里插入图片描述

import random
N = int(input("请输入要随机产生的点数："))
X=0
pi=0
for i in range(N):
    a = random.uniform(0,1) #随机产生0-1内的数
    b = random.uniform(0,1)
    if a*a + b*b <= 1: #随机产生的坐标点坐落在以原点为中心半径为1的四分之一圆内的个数(X^2+Y^2=1 表示以原点为中心半径为1的圆）
        X=X+1
pi=(X/N)*4
print('pi近似等于：{0}'.format(pi))