使用蒙特卡罗方法求圆周率的值

方法如下 

from random import random
from time import perf_counter
DARTS = 1000*1000
hits = 0.0
start = perf_counter()
for i in range(1,DARTS + 1):
    x,y = random(),random()
    dist = pow(x**2+y**2,0.5)
    if dist <= 1.0:
        hits = hits + 1
pi = 4*(hits/DARTS)
print("圆周率的值是:{:.9f}".format(pi))
print("运行所用的时间是{:.6f}s".format(perf_counter()-start))
使用Python蒙特卡罗方法计算圆周率的近似可以通过以下步骤实现: 1. 首先,我们需要了解蒙特卡罗方法的原理。蒙特卡罗方法是一种基于随机抽样的统计模拟方法,通过生成大量的随机点来估计某个问题的解。在计算圆周率的问题中,我们可以通过在一个正方形内随机生成点,并统计落入圆内的点的数量来估计圆周率。 2. 创建一个正方形,并在其中绘制一个圆。假设圆的半径为r,则正方形的边长为2r。 3. 生成大量的随机点,并判断每个点是否落在圆内。可以使用random模块中的random()函数生成0到1之间的随机数,然后将其映射到正方形内的坐标。 4. 统计落在圆内的点的数量。可以使用一个计数器变量来记录落在圆内的点的数量。 5. 根据蒙特卡罗方法的原理,圆的面积与正方形的面积之比等于落在圆内的点的数量与总点数之比。由于正方形的面积为(2r)^2=4r^2,而圆的面积为πr^2,所以可以得到以下公式: π/4 = 落在圆内的点的数量 / 总点数 6. 根据上述公式,可以计算出π的近似: π ≈ 4 * (落在圆内的点的数量 / 总点数) 下面是一个使用Python实现蒙特卡罗方法计算圆周率的示例代码: ```python import random def estimate_pi(num_points): points_inside_circle = 0 total_points = num_points for _ in range(num_points): x = random.uniform(-1, 1) y = random.uniform(-1, 1) distance = x**2 + y**2 if distance <= 1: points_inside_circle += 1 pi_approximation = 4 * (points_inside_circle / total_points) return pi_approximation # 调用函数进行估计 approx_pi = estimate_pi(1000000) print("Approximation of pi:", approx_pi) ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值