python实现蒙特卡罗方法教程
作者：laozhang 发布时间：2023-01-29 16:36:02
标签：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个随机数，可以看出两个输出结果相差并不大。
0
投稿猜你喜欢
- 这篇文章详细的介绍了Oracle数据库SQL语句性能调整的基本原则，具体内容请参考下文。一、问题的提出在应用系统开发初期，由于开发数据库数据
本文用到的文件的下载地址百度网盘链接: https://pan.baidu.com/s/1tmpdEfAZKff5TOMAitUXqQ提取码- 一 概述Python生态系统的一些核心基础数据分析库：NumPy：支持大量的维度数组与矩阵运算，此外也针对数组运算提供大量的数学函数库，包含
- Python中列表（list）的实现其实是一个数组，当要查找某一个元素的时候时间复杂度是O(n)，使用list.index()方法，但是随着
- 前言：其实很简单！举一个简单的例子！这是一个简单的显示时刻的温度折线图：from matplotlib import pyplot as p
- 1 数据导出 python manage.py dumpdata不指定 appname 时默认为导出所有的apppython manage.
- python的库一般都用pip安装。但是有时候也会出现在线安装失败的情况，如下图安装PIL模块时报错：这时候可以采取离线安装的方式；一、首先
- 很多时候关心的是优化SELECT 查询，因为它们是最常用的查询，而且确定怎样优化它们并不总是直截了当。相对来说，将数据装入数据库是直截了当的
- 本文主要介绍了OpenCV实现背景分离(证件照背景替换)，具有一定的参考价值，感兴趣的可以了解一下实现原理图像背景分离是常见的图像处理方法之
- # _*_ coding:utf-8 _*_# name gefile.pyimport osimport statimport socke
- pip是常用的python包管理工具，用python的同学，都离不开pip ～～第一种 在系统自带的python2.7的路径下pip是pyt
- 下面通过实例代码给大家分享Python切片操作去除字符串首尾的空格的方法，具体内容如下所示：#利用切片操作，实现一个trim()函数，去除字
- 我听说在ASP里，可以不用DSN来连接EXCEL，请问如何实现？确实ASP提供了替代DSN的其他方法，甚至可以在运行时建立连接看看下面的范例
- 视频观看视频Python 类和对象面向对象编程是编写表示现实世界中的事物和情景的类，并基于这些类来创建对象。编写类时，你定义一大类对象都有的
- 有了ORM，我们就可以把Web App需要的3个表用Model表示出来：import time, uuidfrom transwarp.db
- 本文实例讲述了Django2 连接MySQL及model测试。分享给大家供大家参考，具体如下：参考：https://www.jb51.net
- 一：简介由paramiko是用python语言写的一个模块，遵循SSH2协议，支持以加密和认证的方式，进行远程服务器的连接。由于使用的是py
- Python 3 利用 Dlib 19.7 实现摄像头人脸检测特征点标定0.引言利用python开发，借助Dlib库捕获摄像头中的人脸，进行
- 一、mysqlcheck简介mysqlcheck客户端可以检查和修复MyISAM表。它还可以优化和分析表。mysqlcheck的功能类似my
- 原来工作中曾经碰到过UL列表里一些异常的表现,加上昨天看到了http://bbs.blueidea.com/thread-2984871-1
