多语言展示
当前在线:152今日阅读:75今日分享:44

圆周率π的来历

介绍π的历史发展,以及π的计算方法。
方法/步骤
1

什么是π?圆周率,一般以π来表示,是一个在数学及物理学普遍存在的数学常数。它定义为圆形之周长与直径之比值。它圆周率π也等于圆形之面积与半径平方之比值。是精确计算圆周长、圆面积、球体积等几何形状的关键值。 在分析学上,π可以严格地定义为满足sin(x) = 0的最小正实数x。2011年6月部分学者认为圆周率定义不合理,要求改为6.28。π是第十六个希腊字母,本来它是和圆周率没有关系的,但大数学家欧拉从一七三六年开始,在书信和论文中都用π来表示圆周率。因为他是大数学家,所以人们也有样学样地用π来表示圆周率了。但π除了表示圆周率外,也可以用来表示其他事物,在统计学中也能看到它的出现。π=Pai(π=Pi)古希腊欧几里德《几何原本》(约公元前3世纪初)中提到圆周率是常数,中国古算书《周髀算经》( 约公元前2世纪)中有“径一而周三”的记载,也认为圆周率是常数。历史上曾采用过圆周率的多种近似值,早期大都是通过实验而得到的结果,如古埃及纸草书(约公元前1700)中取pi=(4/3)^4≒3.1604 。第一个用科学方法寻求圆周率数值的人是阿基米德,他在《圆的度量》(公元前3世纪)中用圆内接和外切正多边形的周长确定圆周长的上下界,从正六边形开始,逐次加倍计算到正96边形,得到(3+(10/71))<π<(3+(1/7)) ,开创了圆周率计算的几何方法(亦称古典方法,或阿基米德方法),得出精确到小数点后两位的π值。中国数学家刘徽在注释《九章算术》(263年)时只用圆内接正多边形就求得π的近似值,也得出精确到两位小数的π值,他的方法被后人称为割圆术。他用割圆术一直算到圆内接正192边形,得出π≈根号10(约为3.14)。

2

发展历史。古希腊作为古代几何王国对圆周率的贡献尤为突出。古希腊大数学家阿基米德(公元前287–212 年) 开创了人类历史上通过理论π计算圆周率近似值的先河。阿基米德从单位圆出发,先用内接正六边形求出圆周率的下界为3,再用外接正六边形并借助勾股定理求出圆周率的上界小于4。接着,他对内接正六边形和外接正六边形的边数分别加倍,将它们分别变成内接正12边形和外接正12边形,再借助勾股定理改进圆周率的下界和上界。他逐步对内接正多边形和外接正多边形的边数加倍,直到内接正96边形和外接正96边形为止。最后,他求出圆周率的下界和上界分别为223/71 和22/7, 并取它们的平均值3.141851 为圆周率的近似值。阿基米德用到了迭代算法和两侧数值逼近的概念,称得上是“计算数学”的鼻祖。南北朝时代著名数学家祖冲之进一步得出精确到小数点后7位的π值(约5世纪下半叶),得出圆周率π应该介于3.1315926和3.1415927之间,还得到两个近似分数值,密率355/113和约率22/7(分子/分母)。他的辉煌成就比欧洲至少早了近千年。其中的密率在西方直到1573才由德国人奥托得到,1625年发表于荷兰工程师安托尼斯的著作中,欧洲不知道是祖冲之先知道密率的,将密率错误的称之为安托尼斯率。阿拉伯数学家卡西在15世纪初求得圆周率17位精确小数值,打破祖冲之保持近千年的纪录。德国数学家柯伦于1596年将π值算到20位小数值,后投入毕生精力,于1610年算到小数后35位数,该数值被用他的名字称为鲁道夫数。无穷乘积式、无穷连分数、无穷级数等各种π值表达式纷纷出现,π值计算精度也迅速增加。1706年英国数学家梅钦计算π值突破100位小数大关。1873 年另一位英国数学家尚可斯将π值计算到小数点后707位,可惜他的结果从528位起是错的。到1948年英国的弗格森和美国的伦奇共同发表了π的808位小数值,成为人工计算圆周率值的最高纪录。相关教学电子计算机的出现使π值计算有了突飞猛进的发展。1949年美国马里兰州阿伯丁的军队弹道研究实验室首次用计算机(ENIAC)计算π值,一下子就算到2037位小数,突破了千位数。1989年美国哥伦比亚大学研究人员用克雷-2型和IBM-VF型巨型电子计算机计算出π值小数点后4.8亿位数,后又继续算到小数点后10.1亿位数,创下最新的纪录。2010年1月7日——法国一工程师将圆周率算到小数点后27000亿位。2010年8月30日——日本计算机奇才近藤茂利用家用计算机和云计算相结合,计算出圆周率到小数点后5万亿位。2011年10月16日,日本职员近藤茂利用家中电脑将圆周率计算到小数点后10万亿位,刷新了2010年8月由他自己创下的5万亿位吉尼斯世界纪录。今年56岁近藤茂使用的是自己组装的计算机,从去年10月起开始计算,花费约一年时间刷新了纪录。而如今计算机高速发展,人们虽然已经知道π是一个无理数,而且已经计算得越来越精准,而人们不管是工程测量、数学解题过程中,大部分都取前两位数,就是π≈3.14,也产生了圆周率日(3月14日)。

3

各地区的发展。在历史上,有不少数学家都对圆周率做出过研究,当中著名的有阿基米德(Archimedes ofSyracuse)、托勒密(Claudius Ptolemy)、张衡、祖冲之等。他们在自己的国家用各自的方法,辛辛苦苦地去计算圆周率的值。下面,就是世上各个地方对圆周率的研究成果。折叠亚洲中国,最初在《周髀算经》中就有“径一周三”的记载,取π值为3。魏晋时,刘徽曾用使正多边形的边数逐渐增加去逼近圆周的方法(即“割圆术”),求得π的近似值3.1416。汉朝时,张衡得出π的平方除以16等于5/8,即π等于10的开方(约为3.162)。虽然这个值不太准确,但它简单易理解,所以也在亚洲风行了一阵。 王蕃(229-267)发现了另一个圆周率值,这就是3.156,但没有人知道他是如何求出来的。公元5世纪,祖冲之和他的儿子以正24576边形,求出圆周率约为355/113,和真正的值相比,误差小于八亿分之一。这个纪录在一千年后才给打破。印度,约在公元530年,数学大师阿耶波多利用384边形的周长,算出圆周率约为√9.8684。婆罗门笈多采用另一套方法,推论出圆周率等于10的算术平方根。折叠欧洲斐波那契算出圆周率约为3.1418。韦达用阿基米德的方法,算出3.1415926535<π<3.1415926537他还是第一个以无限乘积叙述圆周率的人。(阿基米德,前287-212,古希腊数学家,从单位圆出发,先用内接六边形求出圆周率的下界是3,再用外接六边形结合勾股定理求出圆周率的上限为4,接着对内接和外界正多边形的边数加倍,分别变成了12边形,直到内接和外接96边形为止。最后他求出上界和下界分别为22╱7和223╱71,并取他们的平均值3.141851为近似值,用到了迭代算法和两数逼近的概念,称得算是计算的鼻祖。鲁道夫万科伦以边数多过32的多边形算出有35个小数位的圆周率。华理斯在1655年求出一道公式π/2=2×2×4×4×6×6×8×8...../3×3×5×5×7×7×9×9......欧拉发现的e的iπ次方加1等于0,成为证明π是超越数的重要依据。之后,不断有人给出反正切公式或无穷级数来计算π,在这里就不多说了。

4

计算历史古今中外,许多人致力于圆周率的研究与计算。为了计算出圆周率的越来越好的近似值,一代代的数学家为这个神秘的数贡献了无数的时间与心血。十九世纪前,圆周率的计算进展相当缓慢,十九世纪后,计算圆周率的世界纪录频频创新。整个十九世纪,可以说是圆周率的手工计算量最大的世纪。进入二十世纪,随着计算机的发明,圆周率的计算有了突飞猛进。借助于超级计算机,人们已经得到了圆周率的2,061亿位精度。历史上最马拉松式的计算,其一是德国的Ludolph Van Ceulen,他几乎耗尽了一生的时间,计算到圆的内接正262边形,于1609年得到了圆周率的35位精度值,以至于圆周率在德国被称为Ludolph数;其二是英国的威廉·山克斯,他耗费了15年的光阴,在1874年算出了圆周率的小数点后707位,并将其刻在了墓碑上作为一生的荣誉。可惜,后人发现,他从第528位开始就算错了。把圆周率的数值算得这么精确,实际意义并不大。现代科技领域使用的圆周率值,有十几位已经足够了。如果用鲁道夫算出的35位精度的圆周率值,来计算一个能把太阳系包起来的一个圆的周长,误差还不到质子直径的百万分之一。以前的人计算圆周率,是要探究圆周率是否循环小数。自从1761年兰伯特证明了圆周率是无理数,1882年林德曼证明了圆周率是超越数后,圆周率的神秘面纱就被揭开了。现在的人计算圆周率,多数是为了验证计算机的计算能力,还有,就是为了兴趣。计算方法古人计算圆周率,一般是用割圆法。即用圆的内接或外切正多边形来逼近圆的周长。阿基米德用正96边形得到圆周率小数点后3位的精度;刘徽用正3072边形得到5位精度;鲁道夫用正262边形得到了35位精度。这种基于几何的算法计算量大,速度慢,吃力不讨好。随着数学的发展,数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式。下面挑选一些经典的常用公式加以介绍。除了这些经典公式外,还有很多其它公式和由这些经典公式衍生出来的公式,就不一一列举了。马青公式π=16arctan1/5-4arctan1/239这个公式由英国天文学教授约翰·马青于1706年发现。他利用这个公式计算到了100位的圆周率。马青公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。还有很多类似于马青公式的反正切公式。在所有这些公式中,马青公式似乎是最快的了。虽然如此,如果要计算更多的位数,比如几千万位,马青公式就力不从心了。拉马努金公式1914年,印度天才数学家拉马努金在他的论文里发表了一系列共14条圆周率的计算公式。这个公式每计算一项可以得到8位的十进制精度。1985年Gosper用这个公式计算到了圆周率的17,500,000位。1989年,大卫·丘德诺夫斯基和格雷高里·丘德诺夫斯基兄弟将拉马努金公式改良,这个公式被称为丘德诺夫斯基公式,每计算一项可以得到15位的十进制精度。1994年丘德诺夫斯基兄弟利用这个公式计算到了4,044,000,000位。丘德诺夫斯基公式的另一个更方便于计算机编程的形式是:AGM(Arithmetic-Geometric Mean)算法。高斯-勒让德公式这个公式每迭代一次将得到双倍的十进制精度,比如要计算100万位,迭代20次就够了。1999年9月,日本的高桥大介和金田康正用这个算法计算到了圆周率的206,158,430,000位,创出新的世界纪录。波尔文四次迭代式这个公式由乔纳森·波尔文和彼得·波尔文于1985年发表的。bailey-borwein-plouffe算法这个公式简称BBP公式,由David Bailey, Peter Borwein和Simon Plouffe于1995年共同发表。它打破了传统的圆周率的算法,可以计算圆周率的任意第n位,而不用计算前面的n-1位。这为圆周率的分布式计算提供了可行性。丘德诺夫斯基公式这是由丘德诺夫斯基兄弟发现的,十分适合计算机编程,是目前计算机使用较快的一个公式。以下是这个公式的一个简化版本:莱布尼茨公式π/4=1-1/3+1/5-1/7+1/9-1/11+……最新纪录圆周率的最新计算纪录由日本筑波大学所创造。他们于2009年算出π值2,576,980,370,000 位小数,这一结果打破了由日本人金田康正的队伍于2002年创造的00位小数的世界纪录。法国软件工程师法布里斯-贝拉德日前宣称,他已经计算到了小数点后27,000亿位,从而成功打破了由日本科学家2009年利用超级计算机算出来的小数点后25779亿位的吉尼斯世界纪录。个人背诵圆周率的世界纪录11月20日,在位于陕西杨凌的西北农林科技大学,生命科学学院研究生吕超结束背诵圆周率之后,戴上了象征成功的花环。当日,吕超同学不间断、无差错背诵圆周率至小数点后6,7890位,此前,背诵圆周率的吉尼斯世界纪录为背诵小数点后42,195位。整个过程用时24小时04分。

5

PC机计算 PiFast目前PC机上流行的最快的圆周率计算程序是PiFast。它除了计算圆周率,还可以计算e和sqrt(2)。PiFast可以利用磁盘缓存,突破物理内存的限制进行超高精度的计算,最高计算位数可达240亿位,并提供基于Fabrice Bellard公式的验算功能。PC机上的最高计算记录最高记录:12,884,901,372位时间:2000年10月10日记录创造者:Shigeru Kondo所用程序:PiFast ver3.3机器配置:Pentium III 1G,1792M RAM,WindowsNT4.0,40GBx2(IDE,FastTrak66)计算时间:1884375秒(21天19时26分15秒)验算时间:29小时C++计算程序演示#include#include#include#include#define N 30015//SOURCE-CODE from Haoso.com//ReWeite & Debug by Codester//Dev C++ 5.9.2  using namespace std;void mult (int *a,int b,int *s){for(int i=N,c=0;i>=0;i--){int y=(*(a+i))*b+c;c=y/10;*(s+i)=y%10;}}void divi (int *a,int b,int *s){for(int i=0,c=0;i<=N;i++){int y=(*(a+i))+c*10;c=y%b;*(s+i)=y/b;}}void incr(int *a,int *b,int *s){for(int i=N,c=0;i>=0;i--){int y=(*(a+i))+(*(b+i))+c;c=y/10;*(s+i)=y%10;}}bool eqs(int *a,int *b){int i=0;while(((*(a+i))==(*(b+i)))&&(i<=N)) i++;return i>N;}int main(int argc, char *argv[]){system('title 圆周率计算');cout<<'■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■\n\n'<using namespace std;long long a=1000000, b, c=2800000, d, e, f[2801000], g;int main(){for( ;b-c; ) f[b++] =a/5;for( ; d=0, g=c*2; c-=14,printf('%.4d',e+d/a),e=d%a)for(b=c; d+=f*a,f =d%--g,d/=g--,--b; d*=b ) ;return 0;}注:在自己机器上运行CPU使用率一直在百分之六十运算结果在3万位左右

6

π的数值近似值π≈3.1415926535 8979323846倍数值    当1π=3.14时    2π=6.28    3π=9.42    4π=12.56    5π=15.7    6π=18.84    7π=21.98    8π=25.12    9π=28.26    10π=31.4

推荐信息