π值计算公式.doc

π值计算公式 1. Machin公式 ?????? 这个公式由英国天文学教授John Machin于1706年发现。他利用这个公式计算到了100位的PI值。Machin公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数，所以可以很容易地在计算机上编程实现。??? 还有很多类似于Machin公式的反正切公式。 2．迭代公式 （1）二阶公式 （2）Borwein四阶公式 Ramanujan公式 ??? ??? 1914年，印度数学家Srinivasa Ramanujan在他的论文里发表了一系列共14条PI的计算公式，这是其中之一。这个公式每计算一项可以得到8位的十进制精度。1985年Gosper用这个公式计算到了PI的17,500,000位。??? 1989年，David Gregory Chudnovsky兄弟将Ramanujan公式改良成为：??? ???这个公式被称为Chudnovsky公式，每计算一项可以得到15位的十进制精度。1994年Chudnovsky兄弟利用这个公式计算到了4,044,000,000位。Chudnovsky公式的另一个更方便于计算机编程的形式是： ??? AGM(Arithmetic-Geometric Mean)算法 ??? 初值：??? ???重复计算：??? ??? 最后计算：??? ??? 这个公式每迭代一次将得到双倍的十进制精度，比如要计算100万位，迭代20次就够了。1999年9月Takahashi和Kanada用这个算法计算到了PI的206,158,430,000位，创出新的世界纪录。 ??? Borwein四次迭代式： ??? 初值：??? ??? 重复计算：??? ??? ??? 最后计算：??? ??? 这个公式由Jonathan Borwein和Peter Borwein于1985年发表，它四次收敛于PI。 Bailey-Borwein-Plouffe算法 ??? ??? 这个公式简称BBP公式，由David Bailey, Peter Borwein和Simon Plouffe于1995年共同发表。它打破了传统的PI的算法，可以计算PI的任意第n位，而不用计算前面的n-1位。这为PI的分布式计算提供了可行性。1997年，Fabrice Bellard找到了一个比BBP快40％的公式： ??? π值的计算公式(或方法)可找到上百种甚至更多一些，仅此一项就可写成一本书。但不是每个公式都可用于高精度计算 。