圆周率π数学建模论文2.docVIP

圆周率的近似计算方法问题 摘要: 本文探讨了分别用割圆术、蒙特卡罗两种方法近似计算圆周率的方法，并且借助数学软件matlab完成了相应的程序编写和结果分析，对圆周率的近似计算做了简单的探讨。 割圆术就是等分圆周，依次连接等分点组成圆的内接正多边形。显然，圆的内接正多边形的边数越多，其面积与圆面积之差的绝对值越小。当增加的边数趋于无穷时，则该正多边形的面积就与圆的面积没有差别了。通过迭代逐次逼近的办法可以近似求得π的值。 蒙特卡罗的方法就是利用单位圆与边长为1的正方形面积之比来计算π的近似值。具体思想如下：图二为单位圆的1/4为一个扇形，它是边长为1的正方形的一部分。考虑扇形面积在正方形面积中所占的比例k，得出其结果为π/4,然后乘以4就可以得到π的值。 关键词： 割圆术 迭代逼近法 蒙特卡罗 Matlab 1.引言 在科学工作者的心目中，准确的数字是他们的终身追求。为提高数字π的计算精度，古往今来的许许多多人们为此付出了艰苦的努力。3.14，对这样一个古老的数字，人们已经是再熟悉不过了。随着科学技术的发展，尤其是计算机的广泛应用，对获得高精度的π已非难事，利用专业的数学软件，很轻松就可以获得小数点后几百位的数值。尤其是在近几年，利用计算机可以将π精确到小数点后万亿次，计算π的精确值已经成为测试或者检验超级计算机各项性能一种方法，引发新的概念、方法和思想，产生新的问题提供了广阔的前景。下文就针对π的求解进行一系列的探讨和分析。 模型假设： 1.近似逼近得到的π在保留有限位小数后可以认为是精确值。 2.投点时落在扇形每一个位置的机会是均等的。 符号说明： 圆的面积 S 内接正n边形的面积 内接正2n边形的面积 圆半径 R 弦长AB 弦心距0F CF的长度 落在扇形内的点数 m 总投点数 n 2.问题的分析 π的求解方法，自古以来就有前人开始研究，到现在方法已经很多很多，对于本问题的求解，选用了刘辉经典的割圆术和蒙特卡罗方法中的随机投点方法，借助计算机，快速便捷的对π进行近似计算。 割圆术就是等分圆周，依次连接等分点组成圆的内接正多边形。显然，圆的内接正多边形的边数越多，其面积与圆面积之差的绝对值越小。当增加的边数趋于无穷时，则该正多边形的面积就与圆的面积没有差别了。根据这个原理，画出它的几何关系图形，列出对应的几何关系，通过计算，可以算的π的上、下界。随着边数的增加，通过迭代的思想，逐步的逼近π的精确值。 蒙特卡罗的方法就是利用单位圆与边长为1的正方形面积之比来计算π的近似值。具体思想如下：图二为单位圆的1/4为一个扇形，它是边长为1的正方形的一部分。考虑扇形面积在正方形面积中所占的比例k，得出其结果为π/4,然后乘以4就可以得到π的值。在实现的过程中，利用数学软件，产生均匀分布的随机数，这些数落在圆内的点数与总共点数的比值在乘以4就是近似求得的π值。但是这种方法随机性很强，会产生随机误差，进而导致算出的结果不一定准确，为了减少误差，可以采用多次求去平均值的方法，这样得到的结果更加接近真是值。 3.模型的建立与求解 方法1 . 割圆术 1.模型的建立 割圆术就是等分圆周，依次连接等分点组成圆的内接正多边形。显然，圆的内接正多边形的边数越多，其面积与圆面积之差的绝对值越小。当增加的边数趋于无穷时，则该正多边形的面积就与圆的面积没有差别了。 图1 记圆的面积为S，其内接正n边形的面积为，图1中AB为正n边形的一边，当分点加倍后，面积为。将图中两个四边形AOBC和ABED的面积分别记为和，有 从而： （1） 上式就是刘徽用于计算π的圆面积不等式。 下面我们用刘徽割圆术讨论计算π的具体过程。为简单计，取圆半径等于1,这时圆面积S =π，这样（1）式中S的上、下界就是π的上、下界。 记 （2） 从图1可见，多边形由n边变为2n之后，其面积在原来的基础上增加了n/2个矩形ABED的面积，有 （3） 新的正2n边形的边长为： （4） 上述的式（1）～（4）就是计算π的数学模型。 2.模型的求解 从某个圆内接正多边形开始计算，逐次迭代，就可以逐渐逼近π的值。这里从正6边形开始计算，可以得到，，依次使用计算公式就可以得到正12边形的边长和面积，正24边形的边长和面积，等等。得到的这些面积不断的逼近π的值，所求的π的上、下界均满足不等式（1）。根据题意，用matlab编写程序，如下： format long g n=6;a=1; s=1.5*sqrt(3);sc=[]; for k=1