3 π的计算

来源:互联网 发布:mac无法充电 灯不亮了 编辑:程序博客网 时间:2024/06/10 09:52

π的计算

计算思想:采用蒙特卡洛方法,利用随机实验求解问题。构建一个单位正方形和1/4圆,向正方形内随机抛撒点,计算抛撒的点落在1/4圆内的概率。随着抛撒点的增多,得到的值越精确。

输入:抛撒点的数量
处理:对于抛撒的每个点,计算点到圆心的距离,通过距离判断点在圆内还是圆外。统计落在圆内的点的数量。
输出:π值。

需要用到的库是math库和random库。同时引入计时程序衡量运行时间。

math库常用的函数如下图所示:

这里写图片描述

random库常用的函数如下图所示:

这里写图片描述

π的计算程序如下:

# pi.pyfrom random import randomfrom math import sqrtfrom time import clockDARTS = 1200#抛撒点的数量hits = 0clock()#计时for i in range(1,DARTS):    x, y = random(), random()#通过random函数产生随机坐标    dist = sqrt(x**2 + y**2)#计算坐标到原点的距离    if dist <= 1.0:#如果距离小于或等于1        hits = hits + 1#落入1/4圆的数量就加1pi = 4 * (hits/DARTS)#计算抛撒点落入1/4圆的概率print("Pi的值是 %s" % pi)print("程序运行时间是 %-5.5ss" % clock())