无理数π的计算.docVIP

无理数π的计算 谈起数字，人们常常为文人以数字入诗词的才华叫绝。初唐的骆宾王被称作算博士。骆诗有秦塞重关一百二，汉家离宫三十六，写得颇有气势。苏东坡邦使臣的掌故众人皆知，苏东坡以四诗风雅颂回应三光日月星终让辽使称臣，传下一则数字对的佳话π的计算精度，古往今来的许许多多人们为此付出了艰苦的努力。 π，或者其近似值3.14，从小学的算术课本中，我们就认识了这个数字。而在距今天4000年前，也就是在公元前2000年左右的巴比伦王国这个数字就已经被发现。对这样一个古老的数字，人们已经是再熟悉不过了。现在，在计算机的帮助下，截止到1995年，人们已经将这个古老的常数精确地计算到了六十亿位以上。譬如，利用Matlab的命令，我们可以十分容易得到π的小数点后200位的值如下： vpa (pi, 201) ans = 3. 1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 5058223172 5359408128 4811174502 8410270193 8521105559 6446229489 5493038196 但是，在数学史上，圆周率π的计算，却使得许多的数学家付出了难以形容的艰苦劳动。 在计算机技术发展到今天的情况下，我们来介绍一点π的计算方法，并自己动手试着计算π的值，看看会得到什么样的结论。 一．古典方法 1.阿基米德的方法 通过使正多边形外切或内接于圆，将其边数逐步增多来计算圆周长度的方法，很早就用于圆周率的计算了。显然，若取圆的直径D = 1，则圆周长度就是圆周率π！ 希腊数学家阿基米德就用这种方法从计算圆内接正6边形的周长的开始，最后算到圆内接正96边形，发现圆周率大于。进而他又继续计算圆外切正96边形的周长，发现圆周率小于，就是说， 圆周率 ；用近似小数表示就是 3.140845… π 3.142857… 以后，许多学者就用过这种方法计算圆周率。阿基米德的这种方法只能正确地计算到小数点后第二位，久而久之圆周率就成了3.14。 通过计算圆内接或圆外切正多边形，历史上有些学者计算了圆周率π的近似值，有关资料提供了他们的计算结果。 计算人 多边形的边数 正确值 阿基米德 皮沙诺 比埃塔 罗马奴斯 鲁道夫 6×24 6×24 6×216 5×224 262 3.14 3.141 小数点后10位 小数点后15位 小数点后35位 2.刘徽割圆术 刘徽是我国古代著名的数学家，他的最大功绩是为《九章算术》作注。他的割圆术就是等分圆周，依次连接等分点组成圆的内接正多边形。显然，圆的内接正多边形的边数越多，其面积与圆面积之差的绝对值越小。当增加的边数趋于无穷时，则该正多边形的面积就与圆的面积没有差别了。 刘徽割圆术的具体做法如图1所示。 图1 刘徽割圆术示意图 记圆的面积为S，其内接正n边形的面积为Sn，图1中AB为正n边形的一边，当分点加倍后，面积为S2n。将图中两个四边形AOBC和ABED的面积分别记为SAOBC和SABED ，有 S n×SAOBC = S2n S n×(SAOBC + ) = S2n + (S2n － Sn) 从而 S2n S 2S2n － Sn ⑴ 上式就是刘徽用于计算π的圆面积不等式。 下面我们用刘徽割圆术讨论计算π的具体过程。为简单计，取圆半径等于1,这时圆面积S =π，这样⑴式中S的上、下界就是π的上、下界。记an = |AB|，有 hn = |OF| = dn = |EF| = 1－hn = 1－ ⑵ 从图1可见，多边形由n边变为2n之后，其面积在原来的基础上增加了n/2个矩形ABED的面积，有 S2n = Sn + ⑶ 新的正2n边形的边长为 a2n = = ⑷ 上述的式⑴～⑷就是我们计算π的数学模型。从某个圆内接正多边形开始计算，就可以逐步逼近π的真值。譬如，从正6边形开始计算是很方便的，这时的边长和面积为