python蒙特卡洛方法圆周率_蒙特卡罗方法 python 实现

蒙特卡罗(Monte Carlo)方法的精髓：用统计结果去计算频率，从而得到真实值的近似值。

一、求圆周率的近似值，采用投点法

import numpy as np

import matplotlib.pyplot as plt

from matplotlib.patches import Circle

# 投点次数

n = 10000

# 圆的信息

r = 1.0 # 半径

a, b = (0., 0.) # 圆心

# 正方形区域边界

x_min, x_max = a-r, a+r

y_min, y_max = b-r, b+r

# 在正方形区域内随机投点

x = np.random.uniform(x_min, x_max, n) # 均匀分布

y = np.random.uniform(y_min, y_max, n)

# 计算点到圆心的距离

d = np.sqrt((x-a)**2 + (y-b)**2)

# 统计落在圆内的点的数目

res = sum(np.where(d < r, 1, 0))

# 计算 pi 的近似值(Monte Carlo方法的精髓：用统计值去近似真实值)

pi = 4 * res / n

print('pi: ', pi)

# 画个图看看

fig = plt.figure()

axes = fig.add_subplot(111)

axes.plot(x, y,'ro',markersize = 1)

plt.axis('equal') # 防止图像变形

circle = Circle(xy=(a,b), radius=r, alpha=0.5)

axes.add_patch(circle)

plt.show()

效果图

二、求定积分(definite integral)的近似值，采用投点法

import numpy as np

import matplotlib.pyplot as plt

'''蒙特卡罗方法求函数 y=x^2 在[0,1]内的定积分(值)'''

def f(x):

return x**2

# 投点次数

n = 10000

# 矩形区域边界

x_min, x_max = 0.0, 1.0

y_min, y_max = 0.0, 1.0

# 在矩形区域内随机投点

x = np.random.uniform(x_min, x_max, n) # 均匀分布

y = np.random.uniform(y_min, y_max, n)

# 统计落在函数 y=x^2图像下方的点的数目

res = sum(np.where(y < f(x), 1, 0))

# 计算定积分的近似值(Monte Carlo方法的精髓：用统计值去近似真实值)

integral = res / n

print('integral: ', integral)

# 画个图看看

fig = plt.figure()

axes = fig.add_subplot(111)

axes.plot(x, y,'ro',markersize = 1)

plt.axis('equal') # 防止图像变形

axes.plot(np.linspace(x_min, x_max, 10), f(np.linspace(x_min, x_max, 10)), 'b-') # 函数图像

#plt.xlim(x_min, x_max)

plt.show()

效果图

python蒙特卡洛方法圆周率_蒙特卡罗方法 python 实现相关推荐

二维蒙特卡洛模拟居里温度_蒙特卡罗方法计算居里温度（上）
声明:旨在自学一些物理概念和计算方法后整理的笔记和心得,本文基于Kotze先生的Introduction to Monte Carlo methods for an Ising Model of a ...
matlab蒙特卡罗方法求体积_蒙特卡罗方法详细讲解与MATLAB实现.ppt
下面给出几个常用的α与的数值: ? 关于蒙特卡罗方法的误差需说明两点:第一,蒙特卡罗方法的误差为概率误差,这与其他数值计算方法是有区别的.第二,误差中的均方差σ是未知的,必须使用其估计值来代替,在计 ...
蒙特卡洛python求解派_Python实现蒙特卡罗方法
一.蒙特卡罗方法简介蒙特卡罗(Monte Carlo)方法:简单来说,蒙特卡洛的基本原理简单描述是先大量模拟,然后计算一个事件发生的次数,再通过这个发生次数除以总模拟次数,得到想要的结果,精髓就是: ...
python三重积分_蒙特卡罗方法。三重积分。Python。“+”的操作数父级不受支持...
我尝试用蒙特卡罗方法近似三重积分∫∫∫∫xyzdV,其中S=[0,1]×[0,1]×0,1]. 我有这个代码:from numpy import * import time from scipy.in ...
python 凯利公式_蒙特卡罗方法验证凯利公式
说明本文受知乎陈小米启发而写.有兴趣的朋友可以移步这里. 本文的代码完全是本人所撸. 问题描述假想一个游戏.赢的概率是60%,输的概率40%.入场费随意交.如果赢了获得2倍的入场费金额(1赔1), ...
values在python中的意思_相当于Python的values（）字典方法的Javascript
相当于Python的values()字典方法的Javascript 这个问题已经在这里有了答案: 如何获取Javascript对象的所有属性值(不知道键)? ...
matlab 蒙特卡洛树搜索算法_蒙特卡罗方法入门、蒙特卡洛树简介（转载）
本文通过五个例子,介绍蒙特卡罗方法(Monte Carlo Method). 一.概述蒙特卡罗方法是一种计算方法.原理是通过大量随机样本,去了解一个系统,进而得到所要计算的值. 它非常强大和灵活,又 ...
python怎么编辑文件_如何使用python中的方法对文件进行修改文件名
在使用python语言中的方法操作文件时,打开方法可以直接使用open,但是对文件重命名需要调用os模块中的方法,删除文件也是工具/原料 python 编辑器截图工具台式机方法/步骤 1 进入 ...
python模块下载过程_常用的python模块及安装方法
bsddb3:BerkeleyDB的连接组件 Cheetah-1.0:我比较喜欢这个版本的cheetah cherrypy:一个WEB framework ctypes:用来调用动态链接库 DBUti ...

python蒙特卡洛方法圆周率_蒙特卡罗方法 python 实现

python蒙特卡洛方法圆周率_蒙特卡罗方法 python 实现相关推荐

最新文章

热门文章