pi的计算实验报告.docx

班级：02 姓名：张海洋 学号:9 圆周率〔选做〕 要求：先查资料看看古人是怎样计算 π 的，再对 π 的各种计算方法进展争论和争论〔收敛速度等〕，并给出不同算法算出 π 的小数点后第 10000 位的数字是什么，你觉得该数字应当是多少 第一局部： 圆周率简介 圆周率是指平面上圆的周长与直径之比 〔ratio of the circumference of a circle to the diameter〕 。用符号 π〔读音：pài〕表示。中国古代有圆率、圆率、周等名称。它是一个常数〔约等于〕。它是一个无理数，即无限不循环小数。在日常生活中，通常都用代表圆周率去进展近似计算。而用十位小数便足以应付一般计算。即使是工程师或物理学家要进展较周密的计算，充其量也只需取值至小数点后几百个位。计算圆周率的方法 “历史上一个国家所算得的圆周率的准确程度，可以作为衡量这个国家产时数学进展水平的指标。” 历史上最马拉松式的计算，其一是德国的 Ludolph Van Ceulen，他几乎耗尽了一生的时间，计算到圆的内接正 262 边形，于 1609 年得到了圆周率的 35 位精度值，以至于圆周率在德国被称为 Ludolph 数；其二是英国的 William Shanks，他消耗了 15 年的光阴，在1874 年算出了圆周率的小数点后 707 位。惋惜，后人觉察，他从第 528 位开头就算错了。把圆周率的数值算得这么准确，实际意义并不大。现代科技领域使用的圆周率值，有十几位已经足够了。假设用 Ludolph Van Ceulen 算出的 35 位精度的圆周率值，来计算一个能把太阳系包起来的一个圆的周长，误差还不到质子直径的百万分之一。以前的人计算圆周率，是要探究圆周率是否循环小数。自从 1761 年 Lambert 证明白圆周率是无理数， 1882 年 Lindemann 证明白圆周率是超越数后，圆周率的奇特面纱就被揭开了。 在中国，公元 263 年前后，刘徽提出著名的 “割圆术”求出了比较准确的圆周率。他觉察：当圆内接正多边形的边数不断增加后，多边形的周长会越来越靠近圆周长，而多边形的面积也会越来越靠近圆面积。于是，刘徽利用正多边形面积和圆面积之间的关系，从正六边形开头，逐步把边数加倍：正十二边形、正二十四边形，正四十八边形……， 始终到正三○七二边形，算出圆周率等于三点一四一六，将圆周率的精度提高到小数点后第四位。在刘徽争论的根底上，祖冲之进一步地进展，经过既漫长又烦琐的计算，始终算到圆内接正 24576 边形，而得到一个结论： ＜ π ＜ 同时得到π 的两个近似分 数：约率为 22／7；密率为 355／113。他算出的 π 的 8 位牢靠数字，不但在当时是最周密的圆周率，而且保持世界记录九百多年。以致于有数学史家提议将这一结果命名为“祖率”。 现在的人计算圆周率, 多数是为了验证计算机的计算力量，还 有，就是为了兴趣。 其次局部： 古人计算圆周率方法 古人计算圆周率，一般是用割圆法。即用圆的内接或外切正多边形来靠近圆的周长。Archimedes 用正 96 边形得到圆周率小数点后 3 位的精度；刘徽用正 3072 边形得到 5 位精度；Ludolph Van Ceulen 用正 262 边形得到了 35 位精度。这种基于几何的算法计算量大，速度慢，吃力不讨好。随着数学的进展，数学家们在进展数学争论时有意无意地觉察了很多计算圆周率的公式。下面选择一些经典的常用公式加以介绍。除了这些经典公式外，还有很多其它公式和由这些经典公式衍生出来的公式，就不一一列举了 1. Machin 公式 ? ? 16 arctg 1 5  4arg tg  1 239 x 3 5 7x xarctg ( x ) ? x ? ? ? ? L ? ( ? 1) x 3 5 7x x x 2 n ? 1 3 5 7 2 n ? 1 这个公式由英国天文学教授 John Machin 于 1706 年觉察。他利用这个公式计算到了100 位的圆周率。Machin 公式每计算一项可以得到位的十进制精度。由于它的计算过程中被乘数和被除数都不大于长整数，所以可以很简洁地在计算机上编程实现。 还有很多类似于 Machin 公式的反正切公式。在全部这些公式中，Machin 公式似乎是最快的了。虽然如此，假设要计算更多的位数，比方几千万位，Machin 公式就力不从心了。下面介绍的算法，在 PC 机上计算大约一天时间，就可以得到圆周率的过亿位的精度。这些算法用程序实现起来比较简单。由于计算过程中涉及两个大数的乘除运算，要用FFT(Fast Fourier Transform)算法。FFT 可以将两个大数的乘除运算时间由 O(n2)缩短为 O(nlog(n))。 2、 Ram