Python实现 圆周率的计算（嵩天老师）

示例6：圆周率的计算（公式+蒙特卡洛方法）

圆周率的计算公式

代码1：（公式计算）

#CalPiV1.py
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))

ps：

\表示换行，与写在一行效果无异

简单的循环，程序不再逐行解释

代码2：（蒙特卡洛方法）

#CaLP1V2.py
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("圆周率值是：{}".format(pi))
print("运行时间是：{:.5f}s".format(perf_counter()-start))

 逐行分析：

#CaLP1V2.py
from random import random
从random库中引入random函数

from time import perf_counter
从time库中引入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)
由于x，y都为正，所以仅计算了4分之一的圆，计算圆面积需要乘4

print("圆周率值是：{}".format(pi))
print("运行时间是：{:.5f}s".format(perf_counter()-start))

举一反三

1.数学思维：找到公式，利用公式求解
2.计算思维：抽象一种过程，用计算机自动化求解
-使用time库的计时方法获得程序运行时间
-改变撒点数量，理解程序运行时间得分布
-初步掌握简单的的程序性能分析方法
一、不求解圆周率，而是某个特定图形得面积
二、在工程计算中寻找蒙特卡洛方法的应用场景

Python random库的使用_华泽的花的博客-CSDN博客

python（time库的使用）_华泽的花的博客-CSDN博客

在这里附上random库和time库的使用方法，如在逐行解释程序时有不理解的地方，欢迎来学习一下