python 蒙特卡罗法求π值

论搬砖的艺术

已于 2022-03-18 17:08:08 修改

阅读量3.6k

点赞数 3

分类专栏： Python 文章标签： python

于 2022-03-18 17:05:15 首次发布

本文链接：https://blog.csdn.net/qq_53817374/article/details/123579469

版权

Python 专栏收录该内容

19 篇文章 1 订阅

订阅专栏

这篇博客介绍了两种不同的Python代码实现，通过生成随机点并计算落入单位圆内的比例来估算π的值。第一种方法直接比较点的坐标是否在圆内，第二种方法则检查y坐标是否小于等于x的平方。这种方法展示了概率和几何在计算数学常数上的应用。

摘要由CSDN通过智能技术生成

题目如下：

这里使用函数random.random()随机产生10000个点，范围为0-1的浮点数。

求π的方法：通过求出四分之一的圆的面积，然后反解出π。

面积的求法：通过（落在圆中的随机点的数量/总的随机点的数量）乘正方形的面积）

代码如下

import random
random.seed(180)  # 设置随机数种子，也可以不设置
hit = 0
n = 10000  # 总的点数
for i in range(10000):
    x,y = random.random(),random.random()
    if (x**2 + y**2) < 1:         # 检验是否落在半径为1的圆内
       hit += 1
       
print("{:.2f}".format((hit/10000)*4))    # 反向解出π的值

结果：

例二：

代码如下：

import random
hit = 0
random.seed(10)
for x in range(100000):
    x,y = random.random(),random.random()
    if y<=x**2:
        hit += 1
print("{:.3f}".format(hit/100000))