数学实验之8.pi的近似计算.pptVIP

利用级数计算Pi 1、莱布尼茨级数（1674年发现） 利用级数计算Pi 　　1844年，数学家达什在没有计算机的情况下利用此式算出了Pi的前200位小数。使用误差估计式 计算一下要精确到Pi的200位小数需要取级数的多少项？ 利用级数计算Pi 2、欧拉的两个级数（1748年发现） 这两个级数收敛也非常缓慢，计算时实用价值不大。 利用级数计算Pi 3、基于arctan x的级数 对泰勒级数 即为莱布尼茨级数，直接使用时收敛速度极慢，必须考虑加速算法。 利用级数计算Pi 　　观察级数可知，x的值越接近于0，级数收敛越快。由此可以考虑令 利用级数计算Pi 因此，β＝4α－pi/4非常接近0。 利用级数计算Pi 加速效果非常明显！ 蒙特卡罗（Monte Carlo）法 　　单位圆的面积等于Pi，使用蒙特卡罗法，即用随机投点的方法来求出这个面积Pi的近似值。具体方法如下： 　　在平面直角坐标系中，以O(0,0)，A(1,0)，C(1,1)，B(0,1)为四个顶点作一个正方形，其面积S＝1。以原点O为圆心的单位圆在这个正方形内的部分是圆心角为直角的扇形，面积为S1＝Pi/4。 蒙特卡罗（Monte Carlo）法 　　在这个正方形内随机地投入n个点，设其中有m个点落在单位扇形内。则 ? 随机投点如何来实现？ 蒲丰（Buffon）掷针实验 　　另一种用蒙特卡罗法来计算Pi的方法是 1777年法国数学家蒲丰（Buffon）提出的随机掷针实验。其步骤如下： 　　（1）取一张白纸，在上面画出许多间距为d的等距平行线。 　　（2）取一根长度为l(ld)的均匀直针，随机地向画有平行线的纸上掷去，一共掷n次。观察针和直线相交的次数m。 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 教学课件 数学实验之八 ——　的近似计算 　 在本次试验中，我们将追溯关于圆周率　的计算历程。通过对割圆术、韦达公式、级数加速法、迭代法等计算方法的介绍和计算体验，感受数学思想和数学方法的发展过程，提高对极限和级数收敛性及收敛速度的综合认识，同时使我们看到数学家对科学真理的永无止境的追求。 实验目的 主 页 上一页 下一页 主要内容 四、利用级数计算 二、韦达(VieTa)公式 三、数值积分方法 一、割圆术 六、拉马努金(Ramanujan)公式 五、蒙特卡罗(Monte Carlo)法 实验指导 　　 π是使人们最经常使用的数学常数。人们对π的研究已经持续了2500多年。在今天，这种探索还在继续…… 　　世界上数学家们一致公认：“历史上一个国家计算圆周率的准确度，可以作为衡量这个国家当时数学水平的一个标志。” 实验指导 π值——算法美的追求 π作为圆周率的符号，是由著名数学家欧勒于公元1737年首先使用的。古代的希伯来人，在描述所罗门庙宇中的“熔池”时曾经这样写道：“池为圆形，对径为十腕尺，池高为五腕尺，其周长为三十腕尺。”可见，古希伯来人认为圆周率等于3。不过，那时的建筑师们，似乎没有人不明白，圆周长与直径的比要比3大一些。 　　公元前3世纪古希腊大数学家阿基米德求出了223/71π 22/7。 “割圆术”中学问多 我国2000多年前的《周髀算经》称“周三径一”，这是π的第一个近似值，叫做“古率”。 　　据说，汉代大科学家、文学家张衡，有“圆周率一十之面”的推算。清代李潢考证这句话意思为π≈sqrt(10)。 　　魏晋间刘徽由圆内接正六边形依次倍增到正192边形，计算周长与值径之比，得 　　　　3.141024 π3.142704 实际应用时取3.14，或分数值157/50。 “割圆术”中学问多 他的割圆术已含有无限逼近的极限思想，这是比求π值更可宝贵的。从方法上说，他得到了重要的“刘徽不等式”。 　　设圆内接正n边形的边长为an，圆内接正n边形的面积为Sn。根据勾股定理，边长有如下递推公式： “割之弥细，失之弥少，割之又 割，则与圆合体而无所失矣。” 面积与边长有如下关系： 圆面积S与多边形的面积Sn之间有如下关系： 刘徽不等式 借助于计算机来完成刘徽的工作： a(1)=sqrt(3);b(1)=3*sqrt(3)/2; for i=2:6 a(i)=sqrt(2-sqrt(4-a(i-1)^2)); b(i)=3*2^(i-2)*a(i); 　c(i)=2*b(i)-b(i-1); end n=[3,6,12,24,48,96]; size(b) result=[n;a;b] 刘徽不等式 result‘＝ ans = 3.0000 1.7321 2.5981