凭直观推测或实物度量,来计算 π 值的实验方法所得到的结果是相当粗略的。
真正使圆周率计算建立在科学的基础上,首先应归功于阿基米德。他是科学地研究这一常数的第一个人,是他首先提出了一种能够借助数学过程而不是通过测量的、能够把 π 的值精确到任意精度的方法。由此,开创了圆周率计算的第二阶段。
圆周长大于内接正四边形而小于外切正四边形,因此 2√2 < π < 4 。
当然,这是一个差劲透顶的例子。据说阿基米德用到了正96边形才算出他的值域。
阿基米德求圆周率的更精确近似值的方法,体现在他的一篇论文《圆的测定》之中。在这一书中,阿基米德第一次创用上、下界来确定 π 的近似值,他用几何方法证明了"圆周长与圆直径之比小于 3+(1/7) 而大于 3 + (10/71) ",他还提供了误差的估计。重要的是,这种方法从理论上而言,能够求得圆周率的更准确的值。到公元150年左右,希腊天文学家托勒密得出 π =3.1416,取得了自阿基米德以来的巨大进步。
割圆术。不断地利用勾股定理,来计算正N边形的边长。
在我国,首先是由数学家刘徽得出较精确的圆周率。公元263年前后,刘徽提出著名的割圆术,得出 π =3.14,通常称为"徽率",他指出这是不足近似值。虽然他提出割圆术的时间比阿基米德晚一些,但其方法确有着较阿基米德方法更美妙之处。割圆术仅用内接正多边形就确定出了圆周率的上、下界,比阿基米德用内接同时又用外切正多边形简捷得多。另外,有人认为在割圆术中刘徽提供了一种绝妙的精加工办法,以致于他将割到192边形的几个粗糙的近似值通过简单的加权平均,竟然获得具有4位有效数字的圆周率 π =3927/1250 =3.1416。而这一结果,正如刘徽本人指出的,如果通过割圆计算得出这个结果,需要割到3072边形。这种精加工方法的效果是奇妙的。这一神奇的精加工技术是割圆术中最为精彩的部分,令人遗憾的是,由于人们对它缺乏理解而被长期埋没了。
恐怕大家更加熟悉的是祖冲之所做出的贡献吧。对此,《隋书·律历志》有如下记载:"宋末,南徐州从事祖冲之更开密法。以圆径一亿为丈,圆周盈数三丈一尺四寸一分五厘九毫二秒七忽,朒数三丈一尺四寸一分五厘九毫二秒六忽,正数在盈朒二限之间。密率:圆径一百一十三,圆周三百五十五。约率,圆径七,周二十二。"
这一记录指出,祖冲之关于圆周率的两大贡献。其一是求得圆周率
3.1415926 < π < 3.1415927
其二是,得到 π 的两个近似分数即:约率为22/7;密率为355/113。
他算出的 π 的8位可靠数字,不但在当时是最精密的圆周率,而且保持世界记录九百多年。以致于有数学史家提议将这一结果命名为"祖率"。
这一结果是如何获得的呢?追根溯源,正是基于对刘徽割圆术的继承与发展,祖冲之才能得到这一非凡的成果。因而当我们称颂祖冲之的功绩时,不要忘记他的成就的取得是因为他站在数学伟人刘徽的肩膀上的缘故。后人曾推算若要单纯地通过计算圆内接多边形边长的话,得到这一结果,需要算到圆内接正12288边形,才能得到这样精确度的值。祖冲之是否还使用了其它的巧妙办法来简化计算呢?这已经不得而知,因为记载其研究成果的著作《缀术》早已失传了。这在中国数学发展史上是一件极令人痛惜的事。
中国发行的祖冲之纪念邮票
祖冲之的这一研究成果享有世界声誉:巴黎"发现宫"科学博物馆的墙壁上著文介绍了祖冲之求得的圆周率,莫斯科大学礼堂的走廊上镶嵌有祖冲之的大理石塑像,月球上有以祖冲之命名的环形山……
对于祖冲之的关于圆周率的第二点贡献,即他选用两个简单的分数尤其是用密率来近似地表示 π 这一点,通常人们不会太注意。然而,实际上,后者在数学上有更重要的意义。
密率与 π 的近似程度很好,但形式上却很简单,并且很优美,只用到了数字1、3、5。数学史家梁宗巨教授验证出:分母小于16604的一切分数中,没有比密率更接近 π 的分数。在国外,祖冲之死后一千多年,西方人才获得这一结果。
可见,密率的提出是一件很不简单的事情。人们自然要追究他是采用什么办法得到这一结果的呢?他是用什么办法把圆周率从小数表示的近似值化为近似分数的呢?这一问题历来为数学史家所关注。由于文献的失传,祖冲之的求法已不为人知。后人对此进行了各种猜测。
让我们先看看国外历史上的工作,希望能够提供出一些信息。
1573年,德国人奥托得出这一结果。他是用阿基米德成果22/7与托勒密的结果377/120用类似于加成法"合成"的:(377-22) / (120-7) = 355/113。
1585年,荷兰人安托尼兹用阿基米德的方法先求得:333/106 < π < 377/120,用两者作为 π 的母近似值,分子、分母各取平均,通过加成法获得结果:3 ((15+17)/(106+120) = 355/113。
两个虽都得出了祖冲之密率,但使用方法都为偶合,无理由可言。
在日本,十七世纪关孝和重要著作《括要算法》卷四中求圆周率时创立零约术,其实质就是用加成法来求近似分数的方法。他以3、4作为母近似值,连续加成六次得到祖冲之约率,加成一百十二次得到密率。其学生对这种按部就班的笨办法作了改进,提出从相邻的不足、过剩近似值就近加成的办法,(实际上就是我们前面已经提到的加成法)这样从3、4出发,六次加成到约率,第七次出现25/8,就近与其紧邻的22/7加成,得47/15,依次类推,只要加成23次就得到密率。
钱宗琮先生在《中国算学史》(1931年)中提出祖冲之采用了我们前面提到的由何承天首创的"调日法"或称加权加成法。他设想了祖冲之求密率的过程:以徽率157/50,约率22/7为母近似值,并计算加成权数x=9,于是 (157 + 22×,9) / (50+7×9) = 355/113,一举得到密率。钱先生说:"冲之在承天后,用其术以造密率,亦意中事耳。"
另一种推测是:使用连分数法。
由于求二自然数的最大公约数的更相减损术远在《九章算术》成书时代已流行,所以借助这一工具求近似分数应该是比较自然的。于是有人提出祖冲之可能是在求得盈 二数之后,再使用这个工具,将3.14159265表示成连分数,得到其渐近分数:3,22/7,333/106,355/113,102573/32650…
最后,取精确度很高但分子分母都较小的355/113作为圆周率的近似值。至于上面圆周率渐近分数的具体求法,这里略掉了。你不妨利用我们前面介绍的方法自己求求看。英国李约瑟博士持这一观点。他在《中国科学技术史》卷三第19章几何编中论祖冲之的密率说:"密率的分数是一个连分数渐近数,因此是一个非凡的成就。"
我国再回过头来看一下国外所取得的成果。
1150年,印度数学家婆什迦罗第二计算出 π= 3927/1250 = 3.1416。1424年,中亚细亚地区的天文学家、数学家卡西著《圆周论》,计算了3×228=805,306,368边内接与外切正多边形的周长,求出 π 值,他的结果是:
π=3.14159265358979325
有十七位准确数字。这是国外第一次打破祖冲之的记录。
16世纪的法国数学家韦达利用阿基米德的方法计算 π 近似值,用 6×216正边形,推算出精确到9位小数的 π 值。他所采用的仍然是阿基米德的方法,但韦达却拥有比阿基米德更先进的工具:十进位置制。17世纪初,德国人鲁道夫用了几乎一生的时间钻研这个问题。他也将新的十进制与早的阿基米德方法结合起来,但他不是从正六边形开始并将其边数翻番的,他是从正方形开始的,一直推导出了有262条边的正多边形,约4,610,000,000,000,000,000边形!这样,算出小数35位。为了记念他的这一非凡成果,在德国圆周率 π 被称为"鲁道夫数"。但是,用几何方法求其值,计算量很大,这样算下去,穷数学家一生也改进不了多少。到鲁道夫可以说已经登峰造极,古典方法已引导数学家们走得很远,再向前推进,必须在方法上有所突破。
17世纪出现了数学分析,这锐利的工具使得许多初等数学束手无策的问题迎刃而解。 π 的计算历史也随之进入了一个新的阶段。.