蒙特卡罗方法模拟计算圆周率_基于R语言的实现方法_彭滔.pdfVIP

ISSN1009-3044 E-mail:xsjl@ ComputerKnowledgeandTechnology电脑知识与技术 第10卷第17期 (2014年6月) ComputerKnowledgeandTechnology 电脑知识与技术 Vol.10,No.17,June2014 Tel:+86-55165690964 蒙特卡罗方法模拟计算圆周率——基于R语言的实现方法 彭滔，魏雷阳 (海南三亚理工职业学院，海南三亚572022) 摘要：蒙特卡罗方法是一种新型计算方法，它需要真实的随机数，在统计学方面有强有力的应用，随着高性能计算机变得 越来越便宜，此方面变得愈发普遍，该文通过设计一个随机试验，建立pi值与试验次数的联系方程，使用R语言来模拟计 算PI值。 关键词：蒙特卡罗方法；随机模拟；R语言 中图分类号：TP311 文献标识码：A 文章编号：1009-3044(2014)17-4038-02 ANewWaytoComputeSimulatelyPIwithMonteCarloMethodinRLanguage PENGTao,WEILei-yang (SanyaPolythenicCollege,Sanya572022,China) Abstract: MonteCarlomethodsisanewwayofcomputationalalgorithms,basicallyMonteCarlosimulationmethodsrequire trulyrandomnumbers,MonteCarloisapowerfulnumericaltechniqueusefulforsolvingproblemsespeciallyinstatistics.When high-speedcomputersbecomescheap,themethodhasgainedinimportanceandpopularity,whenacirclewithitsbounding squareisdrawn,letrandomlygenerateapointinthesqure,theprobabilitythatthepointfallsinthecircleiscirclearea/squarear⁃ ea=PI/4,wecangetpifromtheequation.WefinishthestochasticsimulationwithRlanguageinthearticle. Keywords:MonteCarlomethods;stochasticsimulation;Rlanguange 1蒙特卡罗方法简介 蒙特卡罗方法是一种随机模拟方法，以概率和统计理论方法为基础，使用随机数（或更常见的伪随机数）来解决计算问题的方 法。其早期的雏形存在于17-18世纪，1777年法国数学家普丰（GeorgesLouisLecleredeBuffon，1707—1788）提出用投针实验的方 法求圆周率π，统计史上称为普丰投针问题。 在平面上画有一组间距为a的平行线，将一根长度为l（la）的针任意掷在这个平面上，求此针与平行线中任一条相交的概率。 普丰证明了，这个概率是： p=2l/（πa) π为圆周率。 几个世纪以来，先后有无数人重复了这个试验，比较有名气的试验见下表。 表1 试验者 时间 投掷次数 相交次数 圆周率估计值 Wolf 1850年 5000 2532 3.1596 Smith 1855年