圆周率
数学词汇
圆周率,是指周长直径比值,即圆周率=圆周长÷直径,也等于圆面积半径平方之比,即圆周率=圆面积÷半径2,一般用希腊字母π表示。
历史发展
实验时期
一块古巴比伦石匾(约产于公元前1900年至公元前1600年)清楚地记载了圆周率=25/8=3.125。同一时期的古埃及文物,莱因德数学纸草书(Rhind Mathematical Papyrus)也表明圆周率等于分数16/9的平方,约等于3.1605。埃及人似乎在更早的时候就知道圆周率了。英国作家John Taylor(1781—1864)在其名著《金字塔》(《The Great Pyramid: Why was it built, and who built it?》)中指出,造于公元前2500年左右的胡夫金字塔和圆周率有关。例如,金字塔的周长和高度之比等于圆周率的两倍,正好等于圆的周长和半径之比。公元前800至600年成文的古印度宗教巨著《百道梵书》(Satapatha Brahmana)显示了圆周率等于分数339/108,约等于3.139。
几何法时期
古希腊作为古代几何王国对圆周率的贡献尤为突出。古希腊大数学家阿基米德(公元前287年—公元前212年)开创了人类历史上通过理论计算圆周率近似值的先河。阿基米德从单位圆出发,先用内接正六边形求出圆周率的下界为3,再用外接正六边形并借助勾股定理求出圆周率的上界小于4。接着,他对内接正六边形和外接正六边形的边数分别加倍,将它们分别变成内接正12边形和外接正12边形,再借助勾股定理改进圆周率的下界和上界。他逐步对内接正多边形和外接正多边形的边数加倍,直到内接正96边形和外接正96边形为止。最后,他求出圆周率的下界和上界分别为223/71和22/7,并取它们的平均值3.141851为圆周率的近似值。阿基米德用到了迭代算法和两侧数值逼近的概念,称得上是“计算数学”的鼻祖。
中国古算书《周髀算经》(约公元前2世纪)的中有“径一而周三”的记载,意即取。汉朝时,张衡得出,即(约为3.162)。这个值不太准确,但它简单易理解。
公元263年,中国数学家刘徽用“割圆术”计算圆周率,他先从圆内接正六边形,逐次分割一直算到圆内接正192边形。他说:“割之弥细,所失弥少,割之又割,以至于不可割,则与圆周合体而无所失矣。”这包含了求极限的思想。刘徽给出π=3.141024的圆周率近似值,刘徽在得圆周率=3.14之后,将这个数值和晋武库中汉王莽时代制造的铜制体积度量衡标准嘉量斛的直径和容积检验,发现3.14这个数值还是偏小。于是继续割圆到1536边形,求出3072边形的面积,得到令自己满意的圆周率。
公元480年左右,南北朝时期的数学家祖冲之进一步得出精确到小数点后7位的结果,给出不足近似值3.1415926和过剩近似值3.1415927,还得到两个近似分数值,密率和约率。密率是个很好的分数近似值,要取到才能得出比略准确的近似。(参见丢番图逼近
在之后的800年里祖冲之计算出的π值都是最准确的。其中的密率在西方直到1573年才由德国人奥托(Valentinus Otho)得到,1625年发表于荷兰工程师安托尼斯(Metius)的著作中,欧洲称之为Metius' number。
约在公元530年,印度数学大师阿耶波多算出圆周率约为。婆罗摩笈多采用另一套方法,推论出圆周率等于10的算术平方根
阿拉伯数学家卡西在15世纪初求得圆周率17位精确小数值,打破祖冲之保持近千年的纪录。德国数学家鲁道夫·范·科伊伦(Ludolph van Ceulen)于1596年将π值算到20位小数值,后投入毕生精力,于1610年算到小数后35位数,该数值被用他的名字称为鲁道夫数。
分析法时期
这一时期人们开始利用无穷级数或无穷连乘积求π,摆脱可割圆术的繁复计算。无穷乘积式、无穷连分数、无穷级数等各种π值表达式纷纷出现,使得π值计算精度迅速增加。
第一个快速算法由英国数学家梅钦(John Machin)提出,1706年梅钦计算π值突破100位小数大关,他利用了如下公式:
其中arctanx可由泰勒级数算出。类似方法称为“梅钦类公式”。
斯洛文尼亚数学家Jurij Vega于1789年得出π的小数点后首140位,其中只有137位是正确的。这个世界纪录维持了50年。他利用了梅钦于1706年提出的数式。
到1948年英国的弗格森(D. F. Ferguson)和美国的伦奇共同发表了π的808位小数值,成为人工计算圆周率值的最高纪录。
计算机时期
电子计算机的出现使π值计算有了突飞猛进的发展。1949年,美国制造的世上首部电脑——ENIAC(ElectronicNumerical Integrator And Computer)在阿伯丁试验场启用了。次年,里特韦斯纳、冯纽曼和梅卓普利斯利用这部电脑,计算出π的2037个小数位。这部电脑只用了70小时就完成了这项工作,扣除插入打孔卡所花的时间,等于平均两分钟算出一位数。五年后,IBM NORC(海军兵器研究计算机)只用了13分钟,就算出π的3089个小数位。科技不断进步,电脑的运算速度也越来越快,在20世纪60年代至70年代,随着美、英、法的电脑科学家不断地进行电脑上的竞争,π的值也越来越精确。在1973年,Jean Guilloud和Martin Bouyer以电脑CDC 7600发现了π的第一百万个小数位。
在1976年,新的突破出现了。萨拉明(Eugene Salamin)发表了一条新的公式,那是一条二次收敛算则,也就是说每经过一次计算,有效数字就会倍增。高斯以前也发现了一条类似的公式,但十分复杂,在那没有电脑的时代是不可行的。这算法被称为布伦特-萨拉明(或萨拉明-布伦特)演算法,亦称高斯-勒让德演算法。
1989年美国哥伦比亚大学研究人员用克雷-2型(Cray-2)和IBM-3090/VF型巨型电子计算机计算出π值小数点后4.8亿位数,后又继续算到小数点后10.1亿位数。2010年1月7日——法国工程师法布里斯·贝拉将圆周率算到小数点后27000亿位。2010年8月30日——日本计算机奇才近藤茂利用家用计算机和云计算相结合,计算出圆周率到小数点后5万亿位。
2011年10月16日,日本长野县饭田市公司职员近藤茂利用家中电脑将圆周率计算到小数点后10万亿位,刷新了2010年8月由他自己创下的5万亿位吉尼斯世界纪录。当时56岁的近藤茂使用的是自己组装的计算机,从10月起开始计算,花费约一年时间刷新了纪录。
2022年3月14日是国际圆周率日。经吉尼斯世界纪录认证,目前π的最准确值,超过小数点后62,831,853,071,796位。
2024年3月15日,据美国趣味科学网站报道,在国际圆周率日,总部位于美国加州的计算机存储公司Solidigm发布声明称,该公司已将圆周率Pi(π)计算到小数点后约105万亿位,打破此前100万亿位的世界纪录。
计算历史
注:上表正确位数是指小数点后的位数。
算准记录
记号
是第十六个希腊字母的小写。这个符号,亦是希腊语περιφρεια(表示周边、地域、圆周等意思)的首字母。1706年英国数学家威廉·琼斯(William Jones,1675—1749)最先使用“π”来表示圆周率。1736年,瑞士大数学家欧拉也开始用表示圆周率。从此,便成了圆周率的代名词。
要注意不可把和其大写Π混用,后者是指连乘的意思。
公式
圆周率()一般定义为一个圆形的周长()与直径()之比:,或直接定义为单位圆的周长的一半。由相似图形的性质可知,对于任何圆形,的值都是一样,这样就定义出常数。
注意:将定义为单位圆的周长的一半是有意义的,这是因为从现代数学的角度来看,直径为d、半径为r的圆的周长C由以下积分给出:
上式中令,由定积分的换元法可得:
其中是单位圆周的周长(C的表达式中取r=1即得)。若定义,则,与周长公式相符。
而半径为r的圆的面积S由以下积分给出:
令,由定积分的换元法可得:
其中是单位圆的面积(S的表达式中取r=1即得)。利用分部积分法,
于是,
因此,得到关系式:
这样一来也得到了圆面积公式
第二个做法是,以圆形半径为边长作一正方形,然后把圆形面积和此正方形面积的比例定为,即圆形之面积与半径平方之比。
定义圆周率不一定要用到几何概念,比如,定义为满足
的最小正实数。
这里的正弦函数定义为幂级数
圆周率的计算方法
传统计算圆周率的方法依赖于几何构造、简单的数列展开或解析方法,例如割圆术、泰勒级数法等。现代计算圆周率的方法主要依赖于复杂的迭代算法、快速收敛的级数与强大的计算机的计算能力,例如查德诺夫斯基算法 (Chudnovsky Algorithm)、高斯-勒让德算法等。
割圆术
割圆术的基本思想是通过在圆内外分别嵌入正多边形,并逐步增加多边形的边数,使其逐渐接近圆的周长。由于多边形的周长可以精确计算,通过这种方法可以获得圆周长的上界和下界,从而逐步逼近圆周率。如图所示,在一个圆中嵌入一个内接的正多边形(如正六边形)。正六边形的边长等于圆的半径,并可通过几何关系递推得到边数更多的多边形的边长。推导边的关系如下:
BC即内接正12边形的边长。归纳得到与之间的递推式:
所以,由 可计算出所有正内接正多边形的边长 ,,,,... 进而计算相应的周长,得到一列递增的数列,圆的周长就是数列的极限值。用圆周长除以圆直径就可以得到π的值。
泰勒级数法
的泰勒展开式为:
带入x=1,上述式子可化简为:
已知,因此有:
即:
这个式子也被称为莱布尼茨级数(Gregory-Leibniz series)。这个公式非常简单,但收敛速度很慢,需要计算非常多的项才能得到准确的结果。
马青(Maqin)公式
用的泰勒级数来计算π收敛速度太慢,实际应用中不实用。要使收敛速度较快,|x| 应当远小于1。例如就收敛得较快。利用,可以通过倍角公式得到:
由于,所以。
这虽然能够得到一个近似值,但为了进一步修正,需要计算他们的误差 。
由,可得:
因此:
于是得到:
最终推导得出:
这个公式被称为马青(Maqin) 公式。
高斯-勒让德算法
这是一个快速收敛的迭代算法,可以非常高效地计算 π。它使用以下的迭代过程:
初始条件:
迭代公式:
经过几轮迭代后,π 的近似值为:
这个算法的收敛速度非常快,每次迭代大约可以提供倍增的精度。
查德诺夫斯基算法 (Chudnovsky Algorithm)
这是目前用于计算 π 最常用的算法之一,它是由两位数学家格雷戈里·查德诺夫斯基和大卫·查德诺夫斯基兄弟提出的。它基于快速收敛的级数,并且特别适合于大规模的高精度计算。公式如下:
这种方法可以非常快速地逼近 π,被用于大规模 π 的计算。
应用
把圆周率的数值算得这么精确,实际意义并不大。现代科技领域使用的圆周率值,有十几位已经足够了。如果以39位精度的圆周率值,来计算可观测宇宙(observable universe)的大小,误差还不到一个原子的体积。以前的人计算圆周率,是要探究圆周率是否循环小数。自从1761年兰伯特证明了圆周率是无理数,1882年林德曼证明了圆周率是超越数后,圆周率的神秘面纱就被揭开了。
π在许多数学领域都有非常重要的作用。
几何
代数
π是个无理数,即不可表达成两个整数之比,是由德国科学家约翰·海因里希·兰伯特于1761年证明的。1882年,林德曼(Ferdinand von Lindemann)更证明了π是超越数,即π不可能是任何整系数多项式的根。
圆周率的超越性否定了化圆为方这古老尺规作图问题的可能性,因所有尺规作图只能得出代数数,而超越数不是代数数。
数学分析
莱布尼茨(Leibniz)公式:
华里士(Wallis)公式:
当时:
可以化为:
高斯积分:
斯特林公式
欧拉公式:
π的连分数表示:
数论
两个任意自然数是互质概率是。
任取一个任意整数,该整数没有重复质因子的概率为。
一个任意整数平均可用个方法写成两个完全数之和。
概率论
在平面上画有一组间距为d的平行线,求一根长度为l(l<d)的针任意投掷在这个平面上,求此针与平行线中和任一条相交的概率。这就是布丰投针问题。1777年,布丰自己本人证明了这个概率为。
统计学
物理学
相对论的场方程:
国际圆周率日
2011年,国际数学协会正式宣布,将每年的3月14日设为国际数学节,来源则是中国古代数学家祖冲之的圆周率。
国际圆周率日可以追溯至1988年3月14日,旧金山科学博物馆的物理学家Larry Shaw,他组织博物馆的员工和参与者围绕博物馆纪念碑做3又1/7圈(22/7,π的近似值之一)的圆周运动,并一起吃水果派。之后,旧金山科学博物馆继承了这个传统,在每年的这一天都举办庆祝活动。
2009年,美国众议院正式通过一项无约束力决议,将每年的3月14日设定为“圆周率日”。决议认为,“鉴于数学和自然科学是教育当中有趣而不可或缺的一部分,而学习有关π的知识是一教孩子几何、吸引他们学习自然科学和数学的迷人方式……π约等于3.14,因此3月14日是纪念圆周率日最合适的日子。”
趣闻事件
历史上最马拉松式的人手π值计算,其一是德国的鲁道夫·范·科伊伦(Ludolph van Ceulen),他几乎耗尽了一生的时间,于1609年得到了圆周率的35位精度值,以至于圆周率在德国被称为Ludolphine number;其二是英国的威廉·山克斯(William Shanks),他耗费了15年的光阴,在1874年算出了圆周率的小数点后707位,并将其刻在了墓碑上作为一生的荣誉。可惜,后人发现,他从第528位开始就算错了。
谷歌公司2005年的一次公开募股中,共集资四十多亿美元,A股发行数量是14,159,265股,这当然是由π小数点后的位数得来。(顺便一提,谷歌公司2004年的首次公开募股,集资额为$2,718,281,828,与数学常数e有关)
排版软件TeX从第三版之后的版本号为逐次增加一位小数,使之越来越接近π的值:3.1,3.14,……当前的最新版本号是3.1415926。
每年3月14日为圆周率日。“终极圆周率日”则是1592年3月14日6时54分(因为其英式记法为“3/14/15926.54”)和3141年5月9日2时6分5秒(从前往后,3.14159265)。
7月22日为圆周率近似日(英国式日期记作22/7,看成圆周率的近似分数)。
有数学家认为应把“真正的圆周率”定义为2π,并将其记为τ(英文:tau)。
2019年3月14日,谷歌宣布日裔前谷歌工程师爱玛(Emma HarukaIwao)在谷歌云平台的帮助下,计算到圆周率小数点后31.4万亿位,准确的说是31415926535897位,比2016年创下的纪录又增加数万亿位。据了解,爱玛的团队使用了一个名为ycruncher的程序,能将π计算到小数点后数万亿位。该程序由谷歌云平台计算引擎上运行的25个虚拟机驱动。而2016年纪录的创造者皮特(Peter Trueb)是用一台电脑计算出来的。这项计算需要170TB的数据,与整个美国国会图书馆印刷藏品数据量大致相同,爱玛经过大约4个月的计算才打破了此前的世界纪录。
2020年,一个名为北阿拉巴马慈善计算的非营利组织的创始人蒂莫西·穆利肯使用个人电脑,将数值计算到小数点后50万亿位,耗时303天。
2021年8月17日,美国趣味科学网站报道,瑞士研究人员使用一台超级计算机,历时108天,将著名数学常数圆周率π计算到小数点后62.8万亿位,创下该常数迄今最精确值记录。
2022年6月,谷歌云(Google Cloud)在新闻稿中表示,谷歌的云服务已经打破了其在2019年创造的纪录,计算出了 100 万亿位圆周率数字,上一次,谷歌云把圆周率精确到了 31.4 万亿位。
2024年的国际圆周率日(3月14日),总部位于美国加州的计算机存储公司Solidigm发布声明称,该公司已将圆周率Pi(π)计算到小数点后约105万亿位,打破此前100万亿位的世界纪录。
最新修订时间:2024-12-21 16:37
目录
概述
历史发展
参考资料