python弹性碰撞次数圆周率_期末作业 - 作业部落 Cmd Markdown 编辑阅读器

期末作业

物基一班　　陈天懿

学号：2013301020146

计算圆周率 蒙特卡罗方法 设计物理模型计算圆周率

摘要：

圆周率 是圆的周长与直径的比值，数学家已严格证明 是一个无理数。自从发现以来,它便成为了数学与物理学等自然科学中的重要常数，在生产生活中占有极为重要的地位。因此，如何精确计算圆周率就一直是人们研究的重要课题。如今，在高度发展的电子计算机的帮助下，圆周率 的计算速度和精度都达到了前所未有的高度。

本文利用python，分别从几何逼近，概率计算，幂级数展开直接计算和设计物理模型计算这4个角度探讨并实现了 的理论计算和数值模拟。前几个方法较为基础，本文重点在于最后一个方法，它同时也是物理理论模型和数值模拟拟合的一次成功尝试

背景：

是第十六个希腊字母的小写。 这个符号，亦是希腊语 περιφρεια (表示周边，地域，圆周等意思)的首字母。1706年英国数学家威廉·琼斯(William Jones ，1675－1749)最先使用“”来表示圆周率 。1736年，瑞士大数学家欧拉也开始用 表示圆周率。从此， 便成了圆周率的代名词。此后的几百年间，数学家们对 值更为精确的计算从未停止过。

正文：

计算圆周率的四种方法：几何逼近——割圆术：

早在公元263年，中国数学家刘徽用“割圆术”计算圆周率，原理就是：当圆内接多边形的边数无限增加时，多边形的周长和面积可以无限接近圆。他先从圆内接正六边形，逐次分割一直算到圆内接正192边形。他说“割之弥细，所失弥少，割之又割，以至于不可割，则与圆周合体而无所失矣。”这已包含了求极限的思想。刘徽给出 的圆周率近似值。

公元480年左右，南北朝时期的数学家祖冲之进一步得出精确到小数点后7位的结果，给出不足近似值3.1415926和过剩近似值3.1415927，还得到两个近似分数值，密率 和约率 。在之后的800年里祖冲之计算出的π值都是最准确的。其中的密率在西方直到1573年才由德国人奥托(Valentinus Otho)得到，1625年发表于荷兰工程师安托尼斯(Metius)的著作中，欧洲称之为Metius' number。

已割圆术计算圆周率的具体思路如下：

如图：这是一个圆和其内接正多边形的一部分，易知：

当 趋向于无穷大时，圆周长

So

计算结果如下：

…………

虽然这种方法的效率并不低，但是有一点极为严重：采用这种方法，我们必须事先知道 的值以计算 本身，这显然存在因果上的矛盾，因此这只能作为验证 值的手段，而不能作为计算 的方法。

概率法:

蒙特 · 卡罗方法 (Monte Carlo method)，也称统计模拟方法，是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。它的基本思想是：当所求解问题是某种随机事件出现的概率，或者是某个随机变量的期望值时，通过某种“实验”的方法，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征，并将其作为问题的解。

计算圆周率时，我们可以利用几何概型的概率原理，通过