python实现蒙特卡罗方法教程
更新时间：2019年01月28日 08:54:07 投稿：laozhang
在本篇文章里小编给大家分享了关于python实现蒙特卡罗方法和知识点，有需要的朋友们学习下。
蒙特卡罗方法是一种统计模拟方法，由冯·诺依曼和乌拉姆提出，在大量的随机数下，根据概率估计结果，随机数据越多，获得的结果越精确。下面我们将用python实现蒙特卡罗方法。
1.首先我们做一个简单的圆周率的近似计算，在这个过程中我们要用到随机数，因此需要先使用import numpy as np导入numpy库。
2.代码实现：
import numpy as np total = 8000000 count = 0 for i in range(total): x = np.random.rand() y = np.random.rand() dis = (x**2+y**2)**0.5 if dis <= 1: count = count+1 PI = 4*count/total print(PI)
3.在上面的程序中我们用8000000个随机数进行投放，这样得到的结果会更精确一些，运行程序需要一定的时间，最终得到的结果如下
4.下面我们进行一项简单的应用，下图为我在画图工具中随便画的一个图，我们可以用蒙特卡罗方法来估算图中黑色部分的面积。
5.上面的图形是不规则的，我们只需知道在投放大量随机数的情况下，随机数在黑色部分出现的概率，再用总面积相乘即可估算黑色部分的面积。我们知道，黑色的rgb编码为(0,0,0)，所以需要统计rgb编码为(0,0,0)时随机数的投放概率即可。
6.代码实现：
from PIL import Image
import numpy as np
im = Image.open("C:/Users/21974/Desktop/handwrite2.PNG")
total = 9000000
count = 0
defin = 0
width = im.size[0]
height = im.size[1]
for i in range(total): #用蒙特卡罗方法获得估计值
x = np.random.randint(0, width-1)
y = np.random.randint(0, height-1)
k = im.getpixel((x, y))
if k[0]+k[1]+k[2] == 0:
count += 1
print(int(width*height*count/total))
for i in range(width): #用遍历获得准确值
for j in range(height):
k = im.getpixel((i, j))
if k[0] + k[1] + k[2] == 0:
defin += 1
print(defin)
上面的代码可分为两部分，第一个for后面是用蒙特卡罗方法获得的面积的估计值，第二个for后面是用遍历所有像素点的方法获得的面积的精确值，获得两个输出后进行对比。
我们在上面的程序中采用了9000000个随机数，可以看出两个输出结果相差并不大。
相关文章
@dataclass 是 Python 3.7 引入的一个装饰器,用于方便地定义符合数据类协议的类,本文主要介绍了python3中@dataclass的实现示例,感兴趣的可以了解一下2024-02-02
下面小编就为大家带来一篇python 实现上传图片并预览的3种方法(推荐)。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧2017-07-07
这篇文章主要为大家介绍了python数据可视化绘制世界人口地图的示例详解，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪2022-06-06
这篇文章主要介绍了Python3如何实现Win10桌面自动切换,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2020-08-08
这篇文章主要介绍了浅谈pymysql查询语句中带有in时传递参数的问题，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2020-06-06
python利用Appium实现自动控制移动设备并提取数据功能
这篇文章主要介绍了python利用Appium自动控制移动设备并提取数据,本文以控制抖音app滑动并获取抖音短视频发布者为例，通过实例代码给大家介绍的非常详细，需要的朋友可以参考下2021-09-09
python使用PIL和matplotlib获取图片像素点并合并解析
这篇文章主要介绍了python使用PIL和matplotlib获取图片像素点并合并解析,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2019-09-09
这篇文章主要为大家详细介绍了python实现文本界面网络聊天室，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-12-12
这篇文章主要介绍了使用Python更换外网IP的方法，，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧2018-07-07
这篇文章主要介绍了python实现在sqlite动态创建表的方法,涉及Python操作SQLite数据库创建数据表的技巧,具有一定参考借鉴价值,需要的朋友可以参考下2015-05-05
最新评论