蒙特卡罗方法 python 实现2-阿里云开发者社区

开发者社区> 人工智能> 正文
登录阅读全文

蒙特卡罗方法 python 实现2

简介:

如果不考虑作图,这里的两个例子可以改写成下面的样子:

求圆周率

import random

'''
蒙特卡罗模拟

投点法计算圆周率
'''


# 投点游戏
def play_game():
    # 圆
    r = 1.0         # 半径
    a, b = (0., 0.) # 圆心
    
    # 正方形区域边界
    x_min, x_max = a-r, a+r
    y_min, y_max = b-r, b+r
    
    # 在 正方形 区域内随机投点
    x = random.uniform(x_min, x_max) # 均匀分布
    y = random.uniform(y_min, y_max)
    
    # 计算点到圆心距离
    d = (x-a)**2 + (y-b)**2
    
    # 根据落在圆内与否,返回1,0(为方便计数)
    return [0, 1][d<r]


# 游戏次数
n = 100000

# 蒙特卡罗方法,模拟 n 次游戏
pi = 4.0 * sum((play_game() for _ in range(n))) / n  # 4 倍频率近似于圆周率
print('pi: ', pi)

求定积分

import random

'''
蒙特卡罗模拟

投点法计算函数 y=x^2在[0,1]内的定积分
'''

# 函数 y=x^2
def f(x):
    return x**2

# 投点游戏
def play_game():
    
    # 矩形区域边界
    x_min, x_max = 0, 1
    y_min, y_max = 0, 1
    
    # 在 矩形 区域内随机投点
    x = random.uniform(x_min, x_max) # 均匀分布
    y = random.uniform(y_min, y_max)
    
    # 根据点落在函数 y=x^2图像下方与否,返回1,0(为方便计数)
    return [0, 1][y<f(x)]


# 游戏次数
n = 100000

# 蒙特卡罗方法,模拟 n 次游戏
integral = sum((play_game() for _ in range(n))) / n  # 频率近似于定积分
print('integral: ', integral)
本文转自罗兵博客园博客,原文链接:http://www.cnblogs.com/hhh5460/p/6718971.html,如需转载请自行联系原作者

版权声明:本文内容由阿里云实名注册用户自发贡献,版权归原作者所有,阿里云开发者社区不拥有其著作权,亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容,填写侵权投诉表单进行举报,一经查实,本社区将立刻删除涉嫌侵权内容。

分享: