首页 - 常识问答 - 正文

π的计算方法有哪些?

日期:2023-4-80 次浏览

1

π的计算方法有哪些?

1、马青公式 π=16arctan1/5-4arctan1/239 这个公式由英国天文学教授约翰·马青于1706年发现.他利用这个公式计算到了100位的圆周率.马青公式每计算一项可以得到1.4位的十进制精度.因为它的计算过程中被乘数和被除数都不大于长整数。

所以可以很容易地在计算机上编程实现.还有很多类似于马青公式的反正切公式.在所有这些公式中,马青公式似乎是最快的了.虽然如此,如果要计算更多的位数,比如几千万位,马青公式就力不从心了.

2、拉马努金公式 1914年,印度天才数学家拉马努金在他的论文里发表了一系列共14条圆周率的计算公式.这个公式每计算一项可以得到8位的十进制精度.1985年Gosper用这个公式计算到了圆周率的17,500,000位.

扩展资料:

【圆周率的历史】

π=Pài(π=Pi)

古希腊欧几里德《几何原本》(约公元前3世纪初)中提到圆周率是常数,中国古算书《周髀算经》( 约公元前2世纪)中有“径一而周三”的记载,也认为圆周率是常数.历史上曾采用过圆周率的多种近似值,早期大都是通过实验而得到的结果,如古埃π及纸草书(约公元前1700)中取pi=(4/3)^4≈3.1604 。

第一个用科学方法寻求圆周率数值的人是阿基米德,他在《圆的度量》(公元前3世纪)中用圆内接和外切正多边形的周长确定圆周长的上下界,从正六边形开始,逐次加倍计算到正96边形,得到(3+(10/71))

本回答被网友采纳

2

π的计算方法

求pi可以直接用arctan的taylor展开.arctan(x)=x-x^3/3+x^5/5-x^7/7....因为arctan(1)=pi/4,所以pi=4*(1-1/3+1/5-1/7....);这样就可以直接求出pi了,为何要产生随机角度呢?另外几个求pi的方法1:阿基米德的逼近方法(套在两个多边行里的圆,2*pi介于两多边行周长).2:莱布尼茨法,泰勒展开特例(arctan法)3:牛顿的积分法.

3

π的计算方法

求pi可以直接用arctan的taylor展开.arctan(x)=x-x^3/3+x^5/5-x^7/7....因为arctan(1)=pi/4,所以pi=4*(1-1/3+1/5-1/7....);这样就可以直接求出pi了,为何要产生随机角度呢?另外几个求pi的方法1:阿基米德的逼近方法(套在两个多边行里的圆,2*pi介于两多边行周长).2:莱布尼茨法,泰勒展开特例(arctan法)3:牛顿的积分法.

4

5

π有几种算法能得到?

3.1415926535897932384626 古人计算圆周率,一般是用割圆法。即用圆的内接或外切正多边形来逼近圆的周长。Archimedes用正96边形得到圆周率小数点后3位的精度;刘徽用正3072边形得到5位精度;Ludolph Van Ceulen用正262边形得到了35位精度。这种基于几何的算法计算量大,速度慢,吃力不讨好。随着数学的发展,数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式。下面挑选一些经典的常用公式加以介绍。除了这些经典公式外,还有很多其他公式和由这些经典公式衍生出来的公式,就不一一列举了。 Machin公式 这个公式由英国天文学教授John Machin于1706年发现。他利用这个公式计算到了100位的圆周率。Machin公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。 Machin.c 源程序 还有很多类似于Machin公式的反正切公式。在所有这些公式中,Machin公式似乎是最快的了。虽然如此,如果要计算更多的位数,比如几千万位,Machin公式就力不从心了。下面介绍的算法,在PC机上计算大约一天时间,就可以得到圆周率的过亿位的精度。这些算法用程序实现起来比较复杂。因为计算过程中涉及两个大数的乘除运算,要用FFT(Fast Fourier Transform)算法。FFT可以将两个大数的乘除运算时间由O(n2)缩短为O(nlog(n))。 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%的公式: 3.1415926<3.1415927

6

-

以上是关于π的计算公式的问答