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

蒲丰（Buffon）掷针实验 　　（3）由几何概率知道针和直线相交的概率为p=2L/πd，取m/n为p的近似值，则 　　特别取针的长度L为d/2时，π=n/m。 蒲丰（Buffon）掷针实验 　　（3）由几何概率知道针和直线相交的概率为p=2L/πd，取m/n为p的近似值，则 　　特别取针的长度L为d/2时，π=n/m。 拉马努金（Ramanujan)公式 　　目前，计算pi的一个极其有效的公式为 　　这个级数收敛得非常快，级数每增加一项，可提高大约8位小数的精度。 拉马努金（Ramanujan)公式 　　1985年，数学家比尔.高斯帕依使用这个公式在计算机上算出了pi的1750万位小数。这个神奇的公式归功于印度年轻的传奇数学家拉马努金(Ramanujan,1887-1929). 拉马努金（Ramanujan)公式 　　另一个经过改进的计算公式为： 　　级数每增加一项，可提高14位小数的精度。 迭代公式 　　迭代公式1： 1989年，BorWein发现了下列收敛于1/pi的迭代公式： 主 页 上一页 下一页 He Renbin * 数学实验之八 ——　的近似计算 　 在本次试验中，我们将追溯关于圆周率　的计算历程。通过对割圆术、韦达公式、级数加速法、迭代法等计算方法的介绍和计算体验，感受数学思想和数学方法的发展过程，提高对极限和级数收敛性及收敛速度的综合认识，同时使我们看到数学家对科学真理的永无止境的追求。 实验目的 主 页 上一页 下一页 主要内容 四、利用级数计算 二、韦达(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 0 6.0000 1.0000 3.0000 3.4019 12.0000 0.5176 3.1058 3.2117 24.0000 0.2611 3.1326 3.1594 48.0000 0.1308 3.1394 3.1461 96.0000 0.0654 3.1410 3.1427 割圆术的意义 　　刘徽创立的割圆术，其意义不仅在于计算出了Pi的近似值，而且还在于提供了一种研究数学的方法。这种方法相当于今天的“求积分”，后者经16世纪英国的牛顿和德国的莱布尼茨