蒙特·卡罗算法Python实例解析

一 背景

此算法诞生的背景是:

  1. 曼哈顿计划,有极大的计算需求。
  2. 计算机刚开始发展,最适合做计算。

蒙特卡洛算法理论基础是概率论,实际就是暴力计算逼近理想结果。正是在以上两个背景下,它刚好得到了极大的应用和发展。

二 概念

蒙特·卡罗算法,也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。与它对应的是确定性算法。蒙特·卡罗方法在金融工程学,宏观经济学,计算物理学(如粒子输运计算、量子热力学计算、空气动力学计算)等领域应用广泛。

三 实例

上面叙述的是定义,我来描述一个例子:蒙特卡罗法计算圆周率

1.问题描述

  • 一个圆半径R,它有一个外切正方形边长2R。
    可以知道:
  • 圆面积=Pi R^2
  • 正方形面积 2R x 2R=4R^2
  • 从这个正方形内随机抽取一个点,对这个点的要求是在正方形内任意一点的概率平均分布。
    那么这个点在圆以内的概率大概就是pi*R2/4R2=pi/4
  • 生成若干个这样的点,利用平面上两点间距离公式计算这个点到圆心的距离来判断是否在圆内。
    当我们使用足够多的点来进行统计时,我们得到的概率值十分接近pi/4
    这样就可以得到pi值

2.代码

import random 
import math
def main():
    print '请输入迭代的次数:'
    n=int(raw_input())   #n是随机的次数  
    total=0   #total是所有落入圆内的随机点
    for i in xrange(n):
        x=random.random()
        y=random.random()
        if math.sqrt(x**2+y**2)<1.0:   #判断是否落入圆内
            total+=1
    mypi=4.0*total/n   #得到Pi值
    print '迭代次数是',n,'Pi的值是:',mypi
    print '数学pi:',math.pi
    print '误差是:',abs(math.pi-mypi)/math.pi   #计算误差
    
main()

3.结果

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 一. 增强学习简介 1.1 什么是增强学习? 机器学习的算法可以分为三类:监督学习,非监督学习和增强学习。 增强学...
    阿阿阿阿毛阅读 31,507评论 0 25
  • //Spark入门(Python版) - 文章 - 伯乐在线http://blog.jobbole.com/862...
    葡萄喃喃呓语阅读 3,517评论 0 1
  • 近日,为促进员工交流,加强团队建设,调整员工身心健康,佳联集团郭董事长携全体员工及员共开启了一段特别的身心灵旅行—...
    宜日宜酵阅读 3,283评论 0 0
  • 所以,我奔跑不像无定向的,我斗拳,不像打空气的。我是攻克己身,叫身服我,恐怕我传福音给别人,自己反被弃绝了。 理解...
    小王瑞阅读 11,857评论 0 0
  • 感赏今晚自己做自己的事,没有去聚焦女儿做作业是否认真,有没有磨蹭。女儿数学作业不会的不让我教,说要反映她的真实水平...
    勿忘我瑶阅读 916评论 2 3