用蒙特卡洛法求圆周率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),也称统计模拟方法,一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。
用蒙特卡洛法求圆周率pi相关推荐
- 用随机投点法计算pi值matlab,用蒙特卡洛法求圆周率PI的近似值,其算法如下:使用随机函数Rnd,随机生成两个[0,1)之间的随机数,由此构成一个投点坐标(x,y),然后根据x...
(2) 实现上述功能的VB程序如下,请在划线处填入合适代码. Private Sub Command1_Click() Dim m1 As String, m2 As String, sum As S ...
- 蒙特卡洛法求圆周率 c语言,c++蒙特卡洛法求圆周率
在linux下用 g++ pi.cpp -o pi -lm 编译. //蒙特卡洛法求圆周率 #include #include #include #include using namespace st ...
- (C++)按照公式求圆周率pi的近似值的两种方法
方法一 #include<cstdio> #include<cmath> //求圆周率pi的近似值int main(){double PI=0,mu=1;int i=1;whi ...
- 求圆周率pi的近似值
题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求输出的结果总宽度占10位,其中小数部分为8位. 程序中使用浮点型数据时,请定义为双精度doubl ...
- 25593 Problem G 例题5-7 求圆周率pi的近似值
问题 G: 例题5-7 求圆周率pi的近似值 时间限制: 1 Sec 内存限制: 12 MB 题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求 ...
- 1067 例题5-7 求圆周率pi的近似值
题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求输出的结果总宽度占10位,其中小数部分为8位. 程序中使用浮点型数据时,请定义为双精度doubl ...
- 蒙特卡洛求圆周率c语言代码,蒙特卡洛法求圆周率
利用蒙特卡洛算法求圆周率是一个概率的方法,关于这方面的内容很多,而且也很容易理解,更多具体分析过程可以参考如下文章: http://blog.csdn.net/zzhflyqh/article/det ...
- 求圆周率pi的怪异程序 祖冲之圆周率算法绝技之迷
据说华罗庚曾讲过一个故事,说:有个教书先生喜欢喝酒,一天,他叫学生背圆周率,自己却提壶酒到山上庙里找老和尚喝酒去了.有个聪明的学生把圆周率编了个打油诗"山巅一寺一壶酒,尔乐苦煞吾,把酒吃:酒 ...
- 【python】蒙特卡洛法计算圆周率pi
描述 蒙特卡洛(Monte Carlo)方法是由数学家冯·诺伊曼提出的,诞生于上世纪40年代美国的"曼哈顿计划".蒙特卡洛是一个地名,位于赌城摩纳哥,象征概率.蒙特卡洛方法的原理是 ...
最新文章
- java property异常_Java常见的异常和解决的办法
- IcmpBackDoor
- PTA 1067 Sort with Swap(0, i) (25 分)(思维)
- PyQt5 笔记1 -- 环境配置
- Hadoop使用常见问题以及解决方法
- 支付宝支付-提现到个人支付宝
- Spring boot 默认静态资源路径与手动配置访问路径的方法
- 奇异递归模板模式(Curiously Recurring Template Pattern,CRTP)
- ASP.NET Core Web API
- 因为贫贱,所以分手,我错了吗?
- 【白皮书下载】2020中国社交零售白皮书-BCGx腾讯.pdf
- 终于来了!Pyston v2.0 发布,速度比 Python 快 20%!
- 抄底王牌tdx 副图(无未來函数)
- Tricks(三十八)—— 在不计算欧式距离的前提下判断点到两点的距离哪个更近
- 使用XenApp/XenDesktop发布App-V应用——安装App-V服务器
- 机械--NX2007(UG)有限元分析教程2--装配体
- 每年的节假日数据处理
- Java爬取喜马拉雅非付费音频【优化】
- wav转mp3的最简单方法
- 年度成绩大赏|2021,StreamNative 扬帆起航