13 python 随机数random库的使用 蒙特卡罗方法计算π

13 python 随机数random库的使用

1.常用随机数函数

在这里插入图片描述

  • 设定随机数种子,可以保证程序再次执行的可重复性,不会因为随机数影响最终的结果,如果不适用种子,产生的随机数完全随机,程序的结果完全随机

2.扩展随机数函数

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

3.随机数应用实例:蒙特卡罗方法求π

1.数学原理:
在这里插入图片描述
2.利用无穷级数计算π值
在这里插入图片描述

pi=0
N=100 #这里面无穷项用很大的数代替
for k in range(N):
	pi+=1/16**k*(4/(8*k+1)-2/(8*k+4)-1/(8*k+5)-1/(8*k+6))
print(pi)

3.蒙特卡罗方法实现

# -*- coding: utf-8 -*-
from random import random
from time import perf_counter
DARTS=10000000
hits=0.0 #记录有多少个在圆内部的点
start=perf_counter()
for i in range(DARTS+1):
	x,y=random(),random()
	if (x**2+y**2)<=1:
		hits+=1
pi=4*hits/DARTS
print(pi)
print('运行时间为{:.5f}s'.format(perf_counter()-start))

猜你喜欢

转载自blog.csdn.net/bj_zhb/article/details/104891708