- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PI值计算方法1
计算圆周率的一些公式???? -|waruqi 发表于 2005-12-8 9:24:00
Machin公式
??????? 这个公式由英国天文学教授John Machin于1706年发现。他利用这个公式计算到了100位的圆周率。Machin公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。
Machin公式的反正切公式:
pi/4=arctg(1/2)+arctg(1/5)+ arctg(1/8) 1844.达塞利
= arctg(1/2)+ arctg(1/3)
=2 arctg(1/3)+ arctg(1/7)
=12 arctg(1/18)+8 arctg(1/57)-5 arctg(1/239)
在所有这些公式中,Machin公式似乎是最快的了。虽然如此,如果要计算更多的位数,比如几千万位,Machin公式就力不从心了。下面介绍的算法,在PC机上计算大约一天时间,就可以得到圆周率的过亿位的精度。这些算法用程序实现起来比较复杂。因为计算过程中涉及两个大数的乘除运算,要用FFT(Fast Fourier Transform)算法。FFT可以将两个大数的乘除运算时间由O(n2)缩短为O(nlog(n))。
(FFT算法不在此文讲诉)
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公式的另一个更方便于计算机编程的形式是:
???
AGM(Arithmetic-Geometric Mean)
??? Gauss-Legendre公式:
??? 初值:??????? 重复计算:??????? 最后计算:??????? 这个公式每迭代一次将得到双倍的十进制精度,比如要计算100万位,迭代20次就够了。1999年9月Takahashi和Kanada用这个算法计算到了圆周率的206,158,430,000位,创出新的世界纪录。
??? Borwein四次迭代式:
??? 初值:??????? 重复计算:??????????? 最后计算:??????? 这个公式由Jonathan Borwein和Peter Borwein于1985年发表,它四次收敛于圆周率。
Bailey-Borwein-Plouffe算法
???
??? 这个公式简称BBP公式,由David Bailey, Peter Borwein和Simon Plouffe于1995年共同发表。它打破了传统的圆周率的算法,可以计算圆周率的任意第n位,而不用计算前面的n-1位。这为圆周率的分布式计算提供了可行性。1997年,Fabrice Bellard找到了一个比BBP快40%的公式:???
(此上文为转载并改编)
计算pi的另一些公式:
1、? 1620-1684)
4/pi=(3*3*5*5*7*7*9*9*…)/(2*4*4*6*6*8*8*10*.,..)
并由数学家瓦利斯于1655年变换为连分数:
4/pi=1+ 1^2
——————————————
2+ 3^2
——————————
2+ 5^2
——————
2+………
?
2、pi=3+ 1
————————————————
7+ 1
————————————
15+ 1
————————
1+ 1
——————
292+……..
=[3,7,15,1,292,1,1,1,2,1,3,1,14,2,1,1,2,2,2,2,1,84,2,……..]
(南巴特.1770)
由此得近似数:
3/1,22/7,33/106,355/113,1039932/33102,10
文档评论(0)