思路

正方形的面积为:
圆的面积为:
圆的面积比上正方形的面积为:π / 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相关推荐

  1. 用随机投点法计算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 ...

  2. 蒙特卡洛法求圆周率 c语言,c++蒙特卡洛法求圆周率

    在linux下用 g++ pi.cpp -o pi -lm 编译. //蒙特卡洛法求圆周率 #include #include #include #include using namespace st ...

  3. (C++)按照公式求圆周率pi的近似值的两种方法

    方法一 #include<cstdio> #include<cmath> //求圆周率pi的近似值int main(){double PI=0,mu=1;int i=1;whi ...

  4. 求圆周率pi的近似值

    题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求输出的结果总宽度占10位,其中小数部分为8位. 程序中使用浮点型数据时,请定义为双精度doubl ...

  5. 25593 Problem G 例题5-7 求圆周率pi的近似值

    问题 G: 例题5-7 求圆周率pi的近似值 时间限制: 1 Sec  内存限制: 12 MB 题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求 ...

  6. 1067 例题5-7 求圆周率pi的近似值

    题目描述 用如下公式 求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加). 要求输出的结果总宽度占10位,其中小数部分为8位. 程序中使用浮点型数据时,请定义为双精度doubl ...

  7. 蒙特卡洛求圆周率c语言代码,蒙特卡洛法求圆周率

    利用蒙特卡洛算法求圆周率是一个概率的方法,关于这方面的内容很多,而且也很容易理解,更多具体分析过程可以参考如下文章: http://blog.csdn.net/zzhflyqh/article/det ...

  8. 求圆周率pi的怪异程序 祖冲之圆周率算法绝技之迷

    据说华罗庚曾讲过一个故事,说:有个教书先生喜欢喝酒,一天,他叫学生背圆周率,自己却提壶酒到山上庙里找老和尚喝酒去了.有个聪明的学生把圆周率编了个打油诗"山巅一寺一壶酒,尔乐苦煞吾,把酒吃:酒 ...

  9. 【python】蒙特卡洛法计算圆周率pi

    描述 蒙特卡洛(Monte Carlo)方法是由数学家冯·诺伊曼提出的,诞生于上世纪40年代美国的"曼哈顿计划".蒙特卡洛是一个地名,位于赌城摩纳哥,象征概率.蒙特卡洛方法的原理是 ...

最新文章

  1. java property异常_Java常见的异常和解决的办法
  2. IcmpBackDoor
  3. PTA 1067 Sort with Swap(0, i) (25 分)(思维)
  4. PyQt5 笔记1 -- 环境配置
  5. Hadoop使用常见问题以及解决方法
  6. 支付宝支付-提现到个人支付宝
  7. Spring boot 默认静态资源路径与手动配置访问路径的方法
  8. 奇异递归模板模式(Curiously Recurring Template Pattern,CRTP)
  9. ASP.NET Core Web API
  10. 因为贫贱,所以分手,我错了吗?
  11. 【白皮书下载】2020中国社交零售白皮书-BCGx腾讯.pdf
  12. 终于来了!Pyston v2.0 发布,速度比 Python 快 20%!
  13. 抄底王牌tdx 副图(无未來函数)
  14. Tricks(三十八)—— 在不计算欧式距离的前提下判断点到两点的距离哪个更近
  15. 使用XenApp/XenDesktop发布App-V应用——安装App-V服务器
  16. 机械--NX2007(UG)有限元分析教程2--装配体
  17. 每年的节假日数据处理
  18. Java爬取喜马拉雅非付费音频【优化】
  19. wav转mp3的最简单方法
  20. 年度成绩大赏|2021,StreamNative 扬帆起航

热门文章

  1. 【飞浆百度领航团零基础Python】学习笔记
  2. 后端: 互动直播之WebRTC服务开源技术选型
  3. 零基础自学网络安全/网络渗透攻防路线学习方法【建议收藏】
  4. React实现H5在线音乐播放器
  5. Oracle个人笔记
  6. Unity学习笔记:人物换装换武器
  7. 史上最纯净精简版win7系统封装过程,整个体积2.3G
  8. delphi random_random函数详解 ,随机函数 全
  9. Ribbon负载均衡策略、懒加载及饥饿加载
  10. Android单元测试系列(3)-Mock之PowerMock