计算圆周率引发的思考——Python

通过公式计算圆周率：

pi = 0 #圆周率未知
N = 100
for k in range(N):
    pi += 1/pow(16,k)*(
              4/(8*k+1) - 2/(8*k+4) -
              1/(8*k+5) - 1/(8*k+6))
print("圆周率值是: {}".format(pi))

>>>圆周率值是: 3.141592653589793

如果不知道公式要求一个圆面积呢？？？

蒙特卡罗方法也称统计模拟方法

• 通过随机撒点，当撒点够多，通过撒点总数与某区域撒点数的比值可以算出区域面积

from random import random
from time import perf_counter
Points = 1000*1000           # 设定撒点数量为 1千万
hits = 0.0                  # 设定区域内的落点为0
start = perf_counter()      # 程序计时测量时间
for i in range(1, Points+1):
    x, y = random(), random()        # 点坐标 0-1的小数值
    dist = pow(x ** 2 + y ** 2, 0.5) # 点到圆心的距离公式
    if dist <= 1.0:         # 判断点在圆心内外
        hits = hits + 1     # 在区域内 数量➕1
pi = 4 * (hits/Points)       # (hits/DARTS)计算出来是4/pi 的面积 所以*4
print("圆周率值是: {}".format(pi))
print("运行时间是: {:.5f}s".format(perf_counter() - start))

》圆周率值是: 3.142064
》运行时间是: 1.27685s

这是一种有数学思维转向计算思维的模式

• 数学思维：找到公式，利用公式求解
• 计算思维：抽象一种过程，用计算机自动化求解

上述代码 

• - 使用time库的计时方法获得程序运行时间
• - 改变撒点数量，理解程序运行时间的分布

 

思考：如果是其他不规则图形是否也可以用撒点计算出面积呢？