Python——随机法（蒙特卡罗方法）计算圆周率

基本概念

蒙特卡罗方法：蒙特卡罗方法于20世纪40年代美国在第二次世界大战中研制原子弹的“曼哈顿计划”计划的成员S.M.乌拉姆和J.冯·诺伊曼首先提出。数学家冯·诺伊曼用驰名世界的赌城—摩纳哥的Monte Carlo—来命名这种方法，为它蒙上了一层神秘色彩。在这之前，蒙特卡罗方法就已经存在。1777年，法国数学家布丰（Georges Louis Leclere de Buffon，1707—1788）提出用投针实验的方法求圆周率π。这被认为是蒙特卡罗方法的起源。

算法思想

当所求解问题是某种随机事件出现的概率，或者是某个随机变量的期望值时，通过某种"实验"的方法，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征，并将其作为问题的解。

问题分析

通过多次撒点，模拟概率，求面积，而是否在圆内，可以通过到圆心的距离求解，利用计算机的运算速度，可以很快求出圆周率。

撒点次数越多，圆周率就越准确。

源代码

#CalPI.py
from random import random
from math import sqrt
from time import process_time
DARTS=1000
hits=0.0
process_time()
for i in range(1,DARTS+1):x,y=random(),random()dist=sqrt(x**2+y**2)if(dist<=1.0):hits=hits+1
pi=4*(hits/DARTS)
print("PI值是{}.".format(pi))
print("运行时间是：{:.5f}s".format(process_time()))

10e6次

10e7次

10e8次

10e9次

10e10次

实在跑不了，Python太慢了

参考文章

https://baike.so.com/doc/5989374-6202341.html

https://blog.csdn.net/q1694222672/article/details/81985391

Python——随机法（蒙特卡罗方法）计算圆周率相关推荐

Python用蒙特卡罗方法计算圆周率近似值
实验目的: 1.理解蒙特Ÿ卡罗方法原理. 2.理解for循环本质与工作原理. 3.了解random模块中常用函数. 实验内容: 蒙特Ÿ卡罗方法是一种通过概率来得到问题近似解的方法,在很多领域都有重要的 ...
Python+numpy实现蒙特卡罗方法估计圆周率近似值
问题描述:使用蒙特卡罗方法估计圆周率近似值,具体描述详见以前发的文章蒙特.卡罗方法求解圆周率近似值原理与Python实现技术要点:Python扩展库numpy中的模块random可以批量生成特定范围 ...
【Python】蒙特卡罗方法计算圆周率及给定随机数种子
蒙特卡罗方法利用随机点分布,让随机点足够的多,用指定区域随机点数与总区域随机点数做比值来求取指定区域面积求解圆周率,在一个正方形中撒点,根据在1/4圆内点的数量占总撒点数的比例计算圆周率值. 代码 ...
Python项目实践：蒙特卡罗方法计算圆周率
一.数学思维 # CalPi.py n = 100 Pi = 0 for k in range(n):Pi += 1 / pow(16, k) * (4 / (8 * k + 1) - 2 / (8 ...
python：蒙特卡罗方法计算圆周率
求解圆周率可以采用蒙特卡罗方法,在一个正方形中撒点,根据在1/4圆内点的数量占总撒点数的比例计算圆周率值.‪‬‪‬‪‬‪‬‪‬‮‬‫‬‪‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‭‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‭ ...
使用蒙特卡罗方法计算圆周率
介绍这个东西网上已经有很多相关的内容了, 我只是对这个过程感兴趣, 所以自己来做了一个. 网上, 大多数纯粹计算的, 对于我这种视觉型的人来说, 还是要看一下比较好. 这里我加入了 matplotl ...
使用蒙特卡罗方法计算圆周率π Python
[问题描述] 假设有一个单位圆,其面积就是π, 单位圆的外接正方形的边长为2,故正方形的面积是4 在正方形内随机产生m个点,假设落在圆内的点的数量为n,则n/m近似等于圆与正方形的面积比也就是n/m ...
种子数是123，如何用蒙特卡罗方法计算圆周率
#pi.py from random import random,seed DARTS=eval(input()) seed(123) hits=0.0 for i in range (DARTS): ...
蒙特卡洛方法计算圆周率
import random import mathdef calpai():n = 10000000 #随机实验次数r = 1.0 #圆的半径a, b = (0.0, 0.0) #中心点x_neg, ...
蒙特卡罗（Monte Carlo）方法计算圆周率π
一.蒙特卡洛(Monte Carlo)方法简介蒙特卡洛是一个地名,位于赌城摩纳哥,象征概率.蒙特卡洛(Monte Carlo)方法是由大名鼎鼎的数学家冯·诺伊曼提出的,诞生于上世纪40年代美国的&q ...