700字范文,内容丰富有趣,生活中的好帮手!
700字范文 > 用蒙特卡洛法求圆周率pi

用蒙特卡洛法求圆周率pi

时间:2020-04-07 08:42:23

相关推荐

用蒙特卡洛法求圆周率pi

思路

正方形的面积为:

圆的面积为:

圆的面积比上正方形的面积为:π / 4

使用蒙特卡洛法在正方形内随机撒点,落在圆内的点 / 落在正方形内的点(全部的点),就约等于圆的面积 / 正方形的面积 = π / 4

代码实现以第一象限的1/4圆为例。

代码

import randomdef count_pi(n):#这里用1/4圆写代码逻辑i = 0count = 0# n 为传入的总点数量while i < n:# 随机产生x,y坐标# random.random()用于生成一个0到1的随机符点数: 0 <= n < 1.0x = random.random()y = random.random()# 如果x平方 + y平方 < 1,说明在圆内if (pow(x, 2) + pow(y, 2)) < 1:count += 1i += 1# π的值为:4 * (落在圆内的点/总的点)return 4 * (count / n)pi = count_pi(100000)print(pi)

运行结果

运行100次,取平均值。

蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。