圆周率的几种计算方法.docx

《圆周率的几种计算方法.docx》由会员分享，可在线阅读，更多相关《圆周率的几种计算方法.docx（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、圆周率的几种计算方法摘要：本文简要的介绍了圆周率的起源及其计算方法，正是圆周率这个数的 特殊性，致使从古到今许多数学家为之奉献毕生的经历来研究的精确值。 因此，用什么样的方法计算乃使其值更加精确，这是一个很值得研究的问题。关键词：圆周率，计算方法，正多边形，连分数一、很早以前就有了从人类祖先的祖先诞生在这个地球上算起，经历了几千万年的时间。我 们看见的太阳几乎总是圆的，而月亮由于地球的遮挡，有圆有缺。椭圆、抛物线，双曲线等都是很晚才发现的曲线。地球诞生之前，太阳就 是圆形的。月亮大概是和地球同时诞生的.在使用工具和火不久，人类对太阳和月亮，或者对动物和鱼类的眼睛是圆 的,也就是说对圆这种形状一

2、定感到很奇妙。远古，数刚诞生时，肯定只在1 和许多个之间有区别。而且，很早以前，就只考虑1和2这两个数。以后因 为1个人有2只脚和2只手，2个人就有4只脚和4只手，1头家畜有4只脚, 2头家畜有8只脚，等等。不久，就知道了比例的概念。到了这个阶段人们自然关顾圆周的长度与圆的直径之间一定的比例常 数。尽管圆有大有小，但对一个圆来说，其周长/与直径d之间的比例常数就 是圆周率万二、乃的几种计算方法有一个关于圆周率的歌谣，盛行于古代：“山巅一寺一壶酒，尔乐苦煞 吾，把酒吃，酒杀尔，杀不死，乐而乐。”圆周率是圆的周长与直径之比， 表示的是一个常数，符号是希腊字母。人们为了计算圆周率，公元前便开始 对它

3、进行计算。魏晋时期刘徽曾于公元263年用割圆术的方法求到3. 14,这 被称为“徽率”。在公元460年，祖冲之应用了刘徽的割圆术(也就是下面提到的正多边 形的方法)，算得圆周率为3. 1415926c祖冲之所求的乃值，保持了 1000多 年的世界纪录。1596年，荷兰数学家鲁道夫经过长期的努力和探索，把乃值推算到15 位小数，打破了祖冲之长达1000多年的纪录，后来他本人又把这个数推进到 35位。18世纪初，圆周率达到72位。19世纪时，圆周率又求到140位、200 位、500位。1873年，威廉欣克用了几十年时间，将Ji值算到707位。到了 1946年，世界上第一台电子计算机(ENIAC)问

4、世美国，有人在计算 机上用了 70个小时，算出圆周率达到2035位。1955年达到10 017位，1962 年达到10万位。1973年达到100万位，1981年日本数学家把它推算到200 万位。1990年美国数学家继续新的计算，将值乃推到新的顶点4. 8亿位。经过长时间艰苦的计算，值只是个近似值，这是一个永不循环的数学 计算，也是数学史上的马拉松。下面介绍几种计算乃的方法：(一)公元前利用正多边形计算不公元前1650年，埃及人著的兰德纸草书中提出乃=(4/3) 3=3. 1604o但是对乃的第一次科学的尝试应归功于阿基米德。阿基米德计算万值是采用 内接和外切正多边形的方法。数学上一般把它称为计

5、算机的古典方法。在公元前3世纪，古希腊的数学非常发达，为了使得数学计算简便，人 们选一个以长度为直径的圆。这样圆的周长在任何内接正多边形的周长和任 何外切正多边形的周长之间。这样就容易得到的上下界，因为计算内接和 外切正多边形的财长比较简单。阿基米德也掌握了这一原理，他从内接和外 切严六边形开始，按照这个方法逐次进行下去，就得出12、24、38、96边的 内拉和外切正多边形的财长，他利用这一方法最后得到乃值在223/71,22/7 之间，取值为3.14。这一方法和数值发表在他的论文集圆的量度中。我国古代第一个把求圆周率近似值的方法提高到理论高度上来认识的是 刘微。他独立地创造了 “割圆术”，并

6、系统而严密地用内接正多边形来求得 圆周率的近似值，他从内接正六边形算起，计算到圆内接正192边形的面积, 从而得出3. 141024VV3. 142704这一值，后来他沿着这一思路继续前进， 一直算到圆内接正3072边形时,得到了 =3927/1250,%的值为3.内接9。 这是当时得到的最精确的取值。南北朝时期，我国的大数学家祖冲之采用刘 徽的割圆术，一直算到圆内接正24576边形，从而推得：3. 1415926 3. 1415927这一成果记载在他的著作缀术中。可惜的是，这本书已经失 传。为了应用方便，祖冲之对圆周率还给出了两个分数值355/113和22/7, 前者称之为“密率”，后者称之

7、为“给率”。其中“密率” 355/133是一个 很有趣的数字，分母分子恰好是三个最小奇数的重复，既整齐美观、又便于 记忆。355/113=3+4 2/(7 2+8 2)也是很巧妙的组合。它与的不实际值相对误 差只有 0. 00000009 o(二)连分数计算乃用连分数计算不的人不多，要多次展开。首创连分数的是一个叫盖托蒂 的数学家。布朗开罗(16201684)得到的表达式为712+27 25乙+?492+邛2+.这个式子源于下式1V2 + 2 +图711= 1 +42 +在一定范围内计算上式，先采用繁分数形式。7111尸十9一2 +西一2 +2+上2 +迎2再计算c 814 + 8185c 4

8、9x2170 + 98268222858585C 25x85 536 + 2125 2661268268268。9x268 773426612661 2661 7734 + 2661 10395 773477347734再由103957 - 7734 可得=2.97604x7734 3093671 =1039510395因为在展开式中取的项数有限，所以值没有超过3。由上可见，计算量很大，是古人对计算感兴趣吗？对现在的年轻人来讲, 这是枯燥无味的，古人也许因为娱乐或兴趣而高兴这么干下去。（三）一些计算圆周率的经典的常用公式1、1593年，韦达给出2 V2 42 + 2，2 +也 +行 ZZ 71

9、222这一不寻常的公式是万的最早分析表达式。甚至在今天，这个公式的优美 也会令我们赞叹不已。它表明仅仅借助数字2,通过一系列的加、乘、除和 开平方就可算出乃值。2、沃利斯1650年给出:71 2244668,8工一 133455773、Machin 公式7i - 1 barctg 4arctg12392一1357XXX arctsx = x1357这个公式由英国天文学教授John Machin于1706年发现。他利用这个公 式计算到了 100位的圆周率。Machin公式每计算一项可以得到14位的十进 制精度。因为它的计算过程中被乘数和被除数都不大于长整数，所以可以很 容易地在计算机上编程实现。

10、还有很多类似于Machin公式的反正切公式。在所有这些公式中，Machin公式似乎是最快的了。虽然如此，如果要计算更多的位数，比如几千 万位，Machin公式就力不从心了。4、Ramanujan 公式下面介绍的算法，在PC机上计算大约一天时间，就可以得到圆周率的过 亿位的精度。这些算法用程序实现起来比较复杂。因为计算过程中涉及两个 大数的乘除运算，要用FFT(Fast Fourier Transfonn)算法。FFT可以将两个 大数的乘除运算时间由0(n2)缩短为0 (nlog (n)o兀=-2碎n=0980144,?(n!)4994,?(M!(1103 + 263907?)1914年，印度数

11、学家Srinivasa Ramanujan在他的论文里发表了一系列共14 条圆周率的计算公式，这是其中之一。这个公式每计算一项可以得到8位的 十进制精度。1985年Gosper用这个公式计算到了圆周率的17, 500, 000位。5、乌依塔6、牛顿 )兀-7、欧拉计算的方法还有很多本文仅列具几种常见并且典型的计算方法。最后 不得不提一句的是，为什么对于的位数的竞争会持续不断呢？这也应该是 相当一部分人的问题吧，可能是因为分数诞生的缘故吧！参考文献：1崛场芳数.万的奥秘M.北京：科学出版社，1989.2 . . . 4K 2005-1-17.3陈仁政著.说不尽的冗M.北京：高教出版社，1992.