圆周率的几种计算方法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
圆周率π的几种计算方法
姓名李至佳
学号 06205013
专业基础数学
摘要:本文简要的介绍了圆周率的起源及其计算方法,正是圆周率这个数的特殊性,致使从古到今许多数学家为之奉献毕生的经历来研究π的精确值。因此,用什么样的方法计算π使其值更加精确,这是一个很值得研究的问题。
关键词:圆周率,计算方法,正多边形,连分数
一、很早以前就有了π
从人类祖先的祖先诞生在这个地球上算起,经历了几千万年的时间。我们看见的太阳几乎总是圆的,而月亮由于地球的遮挡,有圆有缺。
椭圆、抛物线,双曲线等都是很晚才发现的曲线。地球诞生之前,太阳就是圆形的。月亮大概是和地球同时诞生的.
在使用工具和火不久,人类对太阳和月亮,或者对动物和鱼类的眼睛是圆的,也就是说对圆这种形状一定感到很奇妙。远古,数刚诞生时,肯定只在1和许多个之间有区别。而且,很早以前,就只考虑1和2这两个数。以后因为1个人有2只脚和2只手,2个人就有4只脚和4只手,1头家畜有4只脚,2头家畜有8只脚,等等。不久,就知道了比例的概念。
到了这个阶段人们自然关顾圆周的长度与圆的直径之间一定的比例常数。尽管圆有大有小,但对一个圆来说,其周长l与直径d之间的比例常数就是圆周率π
二、π的几种计算方法
有一个关于圆周率的歌谣,盛行于古代:“山巅一寺一壶酒,尔乐苦煞吾,把酒吃,酒杀尔,杀不死,乐而乐。” 圆周率是圆的周长与直径之比,表示的是一个常数,符号是希腊字母。人们为了计算圆周率,公元前便开始对它进行计算。魏晋时期刘徽曾于公元263年用割圆术的方法求到3.14,这被称为“徽率”。
在公元460年,祖冲之应用了刘徽的割圆术(也就是下面提到的正多边形的方法),算得圆周率为3.1415926。祖冲之所求的π值,保持了1000多年的世界纪录。
1596年,荷兰数学家鲁道夫经过长期的努力和探索,把π值推算到15位小数,打破了祖冲之长达1000多年的纪录,后来他本人又把这个数推进到35位。
18世纪初,圆周率达到72位。19世纪时,圆周率又求到140位、200位、500位。1873年,威廉欣克用了几十年时间,将π值算到707位。
到了1946年,世界上第一台电子计算机(ENIAC)问世美国,有人在计算机上用了70个小时,算出圆周率达到2035位。1955年达到10 017位,1962年达到10万位。1973年达到100万位,1981年日本数学家把它推算到200万位。1990年美国数学家继续新的计算,将值π推到新的顶点4.8亿位。
经过长时间艰苦的计算,π值只是个近似值,这是一个永不循环的数学计算,也是数学史上的马拉松。
下面介绍几种计算π的方法:
(一) 公元前利用正多边形计算π
公元前1650年,埃及人著的兰德纸草书中提出π=(4/3) 3=3.1604。但是对π的第一次科学的尝试应归功于阿基米德。 阿基米德计算π值是采用内接和外切正多边形的方法。数学上一般把它称为计算机的古典方法。
在公元前3世纪,古希腊的数学非常发达,为了使得数学计算简便,人们选一个以长度为直径的圆。这样圆的周长在任何内接正多边形的周长和任何外切正多边形的周长之间。这样就容易得到π的上下界,因为计算内接和外切正多边形的财长比较简单。阿基米德也掌握了这一原理,他从内接和外切严六边形开始,按照这个方法逐次进行下去,就得出12、24、38、96边的内拉和外切正多边形的财长,他利用这一方法最后得到π值在223/71,22/7之间,取值为3.14。这一方法和数值发表在他的论文集《圆的量度》中。
我国古代第一个把求圆周率近似值的方法提高到理论高度上来认识的是刘微。他独立地创造了“割圆术”,并系统而严密地用内接正多边形来求得圆周率的近似值,他从内接正六边形算起,计算到圆内接正192边形的面积,从而得出3.141024<<π3.142704这一值,后来他沿着这一思路继续前进,一直算到圆内接正3072边形时,得到了π=3927/1250,π的值为3.14159。这是当时得到的最精确的取值。 南北朝时期,我国的大数学家祖冲之采用刘徽的割圆术,一直算到圆内接正24576边形,从而推得: 3.1415926<π<
3.1415927 这一成果记载在他的著作《缀术》中。可惜的是,这本书已经失传。为了应用方便,祖冲之对圆周率还给出了两个分数值355/113和22/7,前者称之为“密率”,后者称之为“给率”。其中“密率”355/133是一个很有趣的数字,分母分子恰好是三个最小奇数的重复,既整齐美观、又便于记忆。355/113=3+4 2/(7 2+8 2)也是很巧妙的组合。它与的π实际值相对误差只有0.00000009 。
(二)连分数计算π
用连分数计算π的人不多,要多次展开。首创连分数的是一个叫盖托蒂的数学家。布朗开罗(1620—1684)得到的表达式为
254928122 (92)
4112π+++++=+
这个式子源于下式 1
19422524928122...
π=++++++ 在一定范围内计算上式,先采用繁分数形式。
119422524928122
π=+++++ 再计算
8148185492170982682,2222858585
2585536212526612268268268
92687734226612661+⨯++
==+==⨯++==⨯+= 266177342661103951773477347734
++== 再由
1039547734
π
= 可得
4773430936 2.97601039510395
π⨯=
==⋅⋅⋅ 因为在展开式中取的项数有限,所以π值没有超过3。
由上可见,计算量很大,是古人对计算感兴趣吗?对现在的年轻人来讲,这是枯燥无味的,古人也许因为娱乐或兴趣而高兴这么干下去。
(三)一些计算圆周率的经典的常用公式
1、1593年,韦达给出
2
2π= 这一不寻常的公式是π的最早分析表达式。甚至在今天,这个公式的优美也会令我们赞叹不已。它表明仅仅借助数字2,通过一系列的加、乘、除和开平方就可算出 π值。
2、沃利斯1650年给出:
22446688213345577π
⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅=⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅
3、Machin 公式
1
11645239
arctg arctg π=- 357211(1)35721
n n x x x x arctgx x n --=-+-+⋅⋅⋅+-- 这个公式由英国天文学教授John Machin 于1706年发现。他利用这个公式计算到了100位的圆周率。Machin 公式每计算一项可以得到14位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。
还有很多类似于Machin 公式的反正切公式。在所有这些公式中,