怎样才能精确计算圆周率?

怎样才能精确计算圆周率?
怎样才能精确计算圆周率?

怎样才能精确计算圆周率?
圆周率的计算方法
古人计算圆周率,一般是用割圆法.即用圆的内接或外切正多边形来逼近圆的周长.Archimedes用正96边形得到圆周率小数点后3位的精度；刘徽用正3072边形得到5位精度；Ludolph Van Ceulen用正262边形得到了35位精度.这种基于几何的算法计算量大,速度慢,吃力不讨好.随着数学的发展,数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式.下面挑选一些经典的常用公式加以介绍.除了这些经典公式外,还有很多其它公式和由这些经典公式衍生出来的公式,就不一一列举了.
1、 Machin公式
[这个公式由英国天文学教授John Machin于1706年发现.他利用这个公式计算到了100位的圆周率.Machin公式每计算一项可以得到1.4位的十进制精度.因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现.
Machin.c 源程序
还有很多类似于Machin公式的反正切公式.在所有这些公式中,Machin公式似乎是最快的了.虽然如此,如果要计算更多的位数,比如几千万位,Machin公式就力不从心了.下面介绍的算法,在PC机上计算大约一天时间,就可以得到圆周率的过亿位的精度.这些算法用程序实现起来比较复杂.因为计算过程中涉及两个大数的乘除运算,要用FFT(Fast Fourier Transform)算法.FFT可以将两个大数的乘除运算时间由O(n2)缩短为O(nlog(n)).
2、 Ramanujan公式
1914年,印度数学家Srinivasa Ramanujan在他的论文里发表了一系列共14条圆周率的计算公式,这是其中之一.这个公式每计算一项可以得到8位的十进制精度.1985年Gosper用这个公式计算到了圆周率的17,500,000位.
1989年,David & Gregory Chudnovsky兄弟将Ramanujan公式改良成为：
这个公式被称为Chudnovsky公式,每计算一项可以得到15位的十进制精度.1994年Chudnovsky兄弟利用这个公式计算到了4,044,000,000位.Chudnovsky公式的另一个更方便于计算机编程的形式是：
3、AGM(Arithmetic-Geometric Mean)算法
Gauss-Legendre公式：
这个公式每迭代一次将得到双倍的十进制精度,比如要计算100万位,迭代20次就够了.1999年9月Takahashi和Kanada用这个算法计算到了圆周率的206,158,430,000位,创出新的世界纪录.
4、Borwein四次迭代式：
这个公式由Jonathan Borwein和Peter Borwein于1985年发表,它四次收敛于圆周率.
这个公式简称BBP公式,由David Bailey,Peter Borwein和Simon Plouffe于1995年共同发表.它打破了传统的圆周率的算法,可以计算圆周率的任意第n位,而不用计算前面的n-1位.这为圆周率的分布式计算提供了可行性.1997年,Fabrice Bellard找到了一个比BBP快40％的公式：

圆周长除以直径
通常圆的内接正多边形的周长除以直径
所以原内接正多边形的变数越多，求得的值越接近于π
古代的π值都是这样计算的

1、新世界纪录
圆周率的最新计算纪录由日本人金田康正的队伍所创造。他们于2002年算出π值1,241,100,000,000 位小数，这一结果打破了他们于1999年9月18日创造的206,000,000,000位小数的世界纪录。
2、个人计算圆周率的世界纪录
在一个现场解说验证活动中，一名59岁日本老人Akira Haraguchi将圆周率π算到了小数点后的83431位，...

全部展开

1、新世界纪录
圆周率的最新计算纪录由日本人金田康正的队伍所创造。他们于2002年算出π值1,241,100,000,000 位小数，这一结果打破了他们于1999年9月18日创造的206,000,000,000位小数的世界纪录。
2、个人计算圆周率的世界纪录
在一个现场解说验证活动中，一名59岁日本老人Akira Haraguchi将圆周率π算到了小数点后的83431位，这名孜孜不倦的59岁老人向观众讲解了长达13个小时，最终获得认同。这一纪录已经被收入了Guinness世界大全中。据报道，此前的纪录是由一名日本学生于1995年计算出的，当时的精度是小数点后的42000位。

收起

圆周率三种公式:http://elephant.linux.net.cn/articles/pi.php?lang=zh
圆周率的计算方法
古人计算圆周率，一般是用割圆法。即用圆的内接或外切正多边形来逼近圆的周长。Archimedes用正96边形得到圆周率小数点后3位的精度；刘徽用正3072边形得到5位精度；Ludolph Van Ceulen用正262边形得到了35位精度。这种...

全部展开

圆周率三种公式:http://elephant.linux.net.cn/articles/pi.php?lang=zh
圆周率的计算方法
古人计算圆周率，一般是用割圆法。即用圆的内接或外切正多边形来逼近圆的周长。Archimedes用正96边形得到圆周率小数点后3位的精度；刘徽用正3072边形得到5位精度；Ludolph Van Ceulen用正262边形得到了35位精度。这种基于几何的算法计算量大，速度慢，吃力不讨好。随着数学的发展，数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式。下面挑选一些经典的常用公式加以介绍。除了这些经典公式外，还有很多其他公式和由这些经典公式衍生出来的公式，就不一一列举了。
http://www.jason314.com/palgorithm.htm
这个公式由英国天文学教授John Machin于1706年发现。他利用这个公式计算到了100位的圆周率。Machin公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数，所以可以很容易地在计算机上编程实现。
Machin.c 源程序
还有很多类似于Machin公式的反正切公式。在所有这些公式中，Machin公式似乎是最快的了。虽然如此，如果要计算更多的位数，比如几千万位，Machin公式就力不从心了。下面介绍的算法，在PC机上计算大约一天时间，就可以得到圆周率的过亿位的精度。这些算法用程序实现起来比较复杂。因为计算过程中涉及两个大数的乘除运算，要用FFT(Fast Fourier Transform)算法。FFT可以将两个大数的乘除运算时间由O(n2)缩短为O(nlog(n))。
*********************************
http://www.sdxcx.com/a/a278/text/278zs_10.htm
**********************************************
在我国，首先是由数学家刘徽得出较精确的圆周率。公元263年前后，刘徽提出著名的割圆术，得出 π ＝3.14，通常称为“徽率”，他指出这是不足近似值。虽然他提出割圆术的时间比阿基米德晚一些，但其方法确有着较阿基米德方法更美妙之处。割圆术仅用内接正多边形就确定出了圆周率的上、下界，比阿基米德用内接同时又用外切正多边形简捷得多。另外，有人认为在割圆术中刘徽提供了一种绝妙的精加工办法，以致于他将割到192边形的几个粗糙的近似值通过简单的加权平均，竟然获得具有4位有效数字的圆周率 π ＝3927/1250 ＝3.1416。而这一结果，正如刘徽本人指出的，如果通过割圆计算得出这个结果，需要割到3072边形。这种精加工方法的效果是奇妙的。这一神奇的精加工技术是割圆术中最为精彩的部分，令人遗憾的是，由于人们对它缺乏理解而被长期埋没了。
恐怕大家更加熟悉的是祖冲之所做出的贡献吧。对此，《隋书·律历志》有如下记载：“宋末，南徐州从事祖冲之更开密法。以圆径一亿为丈，圆周盈数三丈一尺四寸一分五厘九毫二秒七忽，朒数三丈一尺四寸一分五厘九毫二秒六忽，正数在盈朒二限之间。密率：圆径一百一十三，圆周三百五十五。约率，圆径七，周二十二。”
这一记录指出，祖冲之关于圆周率的两大贡献。其一是求得圆周率
3.1415926 ＜ π ＜ 3.1415927
其二是，得到 π 的两个近似分数即：约率为22／7；密率为355／113。
他算出的 π 的8位可靠数字，不但在当时是最精密的圆周率，而且保持世界记录九百多年。以致于有数学史家提议将这一结果命名为“祖率”。
这一结果是如何获得的呢？追根溯源，正是基于对刘徽割圆术的继承与发展，祖冲之才能得到这一非凡的成果。因而当我们称颂祖冲之的功绩时，不要忘记他的成就的取得是因为他站在数学伟人刘徽的肩膀上的缘故。后人曾推算若要单纯地通过计算圆内接多边形边长的话，得到这一结果，需要算到圆内接正12288边形，才能得到这样精确度的值。祖冲之是否还使用了其它的巧妙办法来简化计算呢？这已经不得而知，因为记载其研究成果的著作《缀术》早已失传了。这在中国数学发展史上是一件极令人痛惜的事。
中国发行的祖冲之纪念邮票
祖冲之的这一研究成果享有世界声誉：巴黎“发现宫”科学博物馆的墙壁上著文介绍了祖冲之求得的圆周率，莫斯科大学礼堂的走廊上镶嵌有祖冲之的大理石塑像，月球上有以祖冲之命名的环形山……
对于祖冲之的关于圆周率的第二点贡献，即他选用两个简单的分数尤其是用密率来近似地表示 π 这一点，通常人们不会太注意。然而，实际上，后者在数学上有更重要的意义。
密率与 π 的近似程度很好，但形式上却很简单，并且很优美，只用到了数字1、3、5。数学史家梁宗巨教授验证出：分母小于16604的一切分数中，没有比密率更接近 π 的分数。在国外，祖冲之死后一千多年，西方人才获得这一结果。
可见，密率的提出是一件很不简单的事情。人们自然要追究他是采用什么办法得到这一结果的呢？他是用什么办法把圆周率从小数表示的近似值化为近似分数的呢？这一问题历来为数学史家所关注。由于文献的失传，祖冲之的求法已不为人知。后人对此进行了各种猜测。

收起