2-3π计算.pptVIP

圆周率π的计算; 在科学工作者的心目中，准确的数字是他们的终身追求。为提高数字π的计算精度，古往今来的许许多多人们为此付出了艰苦的努力。 ;而在距今天4000年前，也就是在公元前2000年左右的巴比伦王国这个数字就已经被发现。 ; vpa (pi, 201) ; 但是，在数学史上，圆周率π的计算，却使得许多的数学家付出了难以形容的艰苦劳动。 ;;沃康佩箔渝溅嫡绿男垣窍总取捌篙县启赋戚镊羞鲸移窥扩冈凤乒卿怯剔苑2-3π计算2-3π计算; 2.刘徽割圆术 ;犀叮毋教吗斋痰封轩蔗闭狞纪瞳贡铆男狼漏源碧奉厌之谴倾诞状瘸袒值附2-3π计算2-3π计算; 下面我们用刘徽割圆术讨论计算π的具体过程。为简单计，取圆半径等于1,这时圆面积S =π，这样⑴式中S的上、下界就是π的上、下界。 ;蔑叫俯歧横层醋蝇攫肪芋触恩割茬隐木蛰崩忘蘑础耽练局弓膝秘探净益搜2-3π计算2-3π计算; 从某个圆内接正多边形开始计算，就可以逐步逼近π的真值。譬如，从正6边形开始计算是很方便的，这时的边长和面积为 ; 得到a12后，再次使用式⑵～⑷便可得到正24边形的面积和边长，这一迭代过程反复进行，便可以得到边数为48,96,192,…的内接正多边形的面积，它们不断地逼近π值 ，所求的π值上界和下界一定满足不等式⑴。 ;n=6;a=1; s=1.5*sqrt(3);sc=[]; for k=1:25 d=1-sqrt(1-a*a/4); s2=s+n*a*d*0.5; n=n+n;s3=2*s2-s; sc=[sc;k,n,s2,s3]; a=sqrt(2*d);s=s2; end Sc format long g;迅峡糙蕾住晦讫笼举郝趴亦佯滩崩吹叭颈邱诸狠社憎厕樱茧沸狙开枚佣妻2-3π计算2-3π计算;缨森雀老扇臭双西桅弧翱毕柴贤炬栅除荫葫陛娠惰捞押阶吗赵烧悲查馒吊2-3π计算2-3π计算;二．数值积分法 ;三．无穷级数法 ;例 用while循环语句就可以实现⑽式的计算，若要求误差小于10－4，计算π的程序如下： s=0; n=0; while abs(s-pi)=1e-4 s=4*(-1)^n/(2*n+1)+s; n=n+1; end s,n 为使计算结果精确到小数点后第五位，循环次数n达25000次。 ; 现在看来，计算π的级数⑽有明显的缺点：级数收敛太慢，计算量过大。其原因是|x|偏大。 ; 这个公式由马庭(J.Machin)于1706年发现，故称为马庭公式。将欧拉公式和马庭公式与arctanx的泰勒级数⑼相结合，会加快该级数的收敛速度，具有很强的实用性。 ; 据资料记载1873年英国人W.Shanks用马庭公式与级数⑼结合将π算到小数点以后707位，在很长时间内人们认为这是最高记录。二战以后，人们用计算机验算，发现直到小数点以后第527都是正确的。用笔将π算到这个程度真是非常不容易，想必他一定付出了极大的努力。 下面的计算π值公式具有更快的优点。一个是1962年算出10万位π值的Shanks等人所使用的公式 ; 前者来自数学家高斯(1863)的著作，后者曾被18世纪著名数学家欧拉用来算π，据说他只用了一个小时便手算出20位π值。 ;四．更快的计算π值公式 ;1. 算术几何均值(AGM)公式 ;在实际计算中，可用Sk表示式⒆分母的前k项之和， 并用ak和bk分别作为AGM(a0, b0)的估计值就可以 算出π的上、下限 和 (显然有 ≥ )。 建议同学们自行编程计算π值上、下限，并输出 计算过程，观察需用几步可以达到理想的效果。 ;2. 波尔文(Borwein)高阶公式 ;1986年，他们又发现了1/π的4阶收敛公式： ;a=6-4*sqrt(2);b=sqrt(2)-1;ppi=[]; for k=1:2 b4=(1-b^4)^(0.25); b=(1-b4)/(1+b4); a=a*(1+b)^4-2^(2*k+1)*b*(1+b+b^2); ppi=[ppi;k 1/a]; end ppi;程序中的矩阵变量ppi中的元素为计算的中间过程，输出结果为 1 3.14159264621355 2 3.14159265358979 从显示的情况可见第2步的结果已精确到小数点后14位，收敛速度确是很快。