生活

圆周率

圆的周长与直径的比值

中文名:圆周率 外文名:Ratio of circumference to diameter;Pi 别名: 近似值:22/7(约率)、355/113(密率) 符号表示:π 计算位数:62.8万亿位(2021年8月17日)
圆周率介绍
圆周率(Pi)是圆的周长与直径的比值,一般用希腊字母π表示,是一个在数学及物理学中普遍存在的数学常数。π也等于圆形之面积与半径平方之比,是精确计算圆周长、圆面积、球体积等几何形状的关键值。在分析学里,π可以严格地定义为满足sinx=0的最小正实数x。圆周率用希腊字母π(读作[paɪ])表示,是一个常数(约等于3.141592654),是代表圆周长和直径的比值。它是一个无理数,即无限不循环小数。在日常生活中,通常都用3.14代表圆周率去进行近似计算。而用十位小数3.141592654便足以应付一般计算。即使是工程师或物理学家要进行较精密的计算,充其量也只需取值至小数点后几百个位。1965年,英国数学家约翰·沃利斯(John Wallis)出版了一本数学专著,其中他推导出一个公式,发现圆周率等于无穷个分数相乘的积。2015年,罗切斯特大学的科学家们在氢原子能级的量子力学计算中发现了圆周率相同的公式。2021年8月17日,美国趣味科学网站报道,瑞士研究人员使用一台超级计算机,历时108天,将著名数学常数圆周率π计算到小数点后62.8万亿位,创下该常数迄今最精确值记录。[1]

发展历史

实验时期

一块古巴比伦石匾(约产于公元前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月起开始计算,花费约一年时间刷新了纪录。

计算历史

日期

计算者

正确位数

详细纪录

前20世纪

未知

1

π=3.125

前20世纪

未知

1

π=3.160493...

前12世纪

未知

-

π=3

前6世纪中

圣经列王记上7章23节

-

π=3

前3世纪

阿基米德

3

π=3.1418

公元前20年

维特鲁威

1

π=3.125

公元前50年-公元前23年

刘歆

1

π=3.1547

130年

张衡

1

π=3.162277…

150年

未知

3

π=3.141666…

250年

王蕃

1

π=3.155555…

263年

刘徽

5

π=3.14159

480年

祖冲之

7

3.1415926

499年

阿耶波多

3

π=3.1416

598年

婆罗摩笈多

1

π=3.162277…

800年

花拉子米

3

π=3.1416

12世纪

婆什迦罗第二

4

π=3.14156

1220年

斐波那契

3

π=3.141818

1400年

Madhava

10

π=3.14159265359

1424年

Jamshid Masud Al Kashi

16

1573年

Valentinus Otho

6

1593年

弗朗索瓦·韦达

9

1593年

Adriaan van Roomen

15

1596年

鲁道夫·范·科伊伦

20

1615年

32

1621年

威理博·司乃耳,范·科伊伦的学生

35

1665年

牛顿

16

1699年

Abraham Sharp

71

1700年

关孝和

10

1706年

John Machin

100

1706年

William Jones

引入希腊字母π

1719年

De Lagny

112

得出127位

前112位正确

1723年

建部贤弘

41

1730年

Kamata

25

1734年

莱昂哈德·欧拉

引入希腊字母π并肯定其普及性

1739年

松永良弼

50

1761年

约翰·海因里希·兰伯特

证明π是无理数

1775年

欧拉

指出π可能是超越数

1794年

Jurij Vega

136

得出140位小数

前136位正确

1794年

阿德里安-马里·勒让德

-

1841年

Rutherford

152

得出208位小数

前152位正确

1844年

Zacharias Dase及Strassnitzky

200

1847年

Thomas Clausen

248

1853年

Lehmann

261

1853年

William Rutherford

440

1855年

Richter

500

1874年

William Shanks

527

得出707位小数

前527位正确

1882年

Lindemann

证明π是超越数

1946年

D. F. Ferguson

620

1947年

710

1947年

808

1949年

J. W. Wrench爵士和L. R. Smith

2,037

首次使用计算机

1955年

J. W. Wrench爵士及L. R. Smith

3,089

1957年

G.E.Felton

7,480

1958年

Francois Genuys

10,000

1958年

G.E.Felton

10,020

1959年

Francois Genuys

16,167

1961年

IBM 7090

晶体管计算机

20,000

1961年

J. W. Wrench, Jr,及L. R. Smith

100,000

1966年

250,000

1967年

500,000

1974年

1,000,000

1981年

金田康正

2,000,000

1982年

4,000,000

1983年

8,000,000

1983年

16,000,000

1985年

Bill Gosper

17,000,000

1986年

David H. Bailey

29,000,000

1986年

金田康正

33,000,000

1986年

67,000,000

1987年

134,000,000

1988年

201,000,000

1989年

楚诺维斯基兄弟

480,000,000

1989年

535,000,000

1989年

金田康正

536,000,000

1989年

楚诺维斯基兄弟

1,011,000,000

1989年

金田康正

1,073,000,000

1992年

2,180,000,000

1994年

楚诺维斯基兄弟

4,044,000,000

1995年

金田康正和高桥大介

4,294,960,000

1995年

6,000,000,000

1996年

楚诺维斯基兄弟

8,000,000,000

1997年

金田康正和高桥大介

51,500,000,000

1999年

68,700,000,000

1999年

206,000,000,000

2002年

金田康正的队伍

1,241,100,000,000

2009年

高桥大介

2,576,980,370,000

2009年

法布里斯·贝拉

2,699,999,990,000

2010年

近藤茂

5,000,000,000,000

2011年

IBM“蓝色基因”

超级电脑

π的前60,000,000,000,000位二进制小数

计算历史

注:上表正确位数是指小数点后的位数。

算准记录

小数点后位数

首次算准者

首次算准时间

1

巴比伦人

前20世纪

2-3

阿基米德

前3世纪(距离上次1700年)

4-5

刘徽

263年(距离上次563年以上)

6-7

祖冲之

480年(距离上次217年)

8-10

Madhava

1400年(距离上次920年)

11-16

Jamshid Masud Al Kashi

1424年(距离上次24年)

17-20

鲁道夫·范·科伊伦

1596年(距离上次172年)

21-32

1615年(距离上次19年)

33-35

威理博·司乃耳, 范·科伊伦的学生

1621年(距离上次6年)

36-71

Abraham Sharp

1699年(距离上次78年)

72-100

John Machin

1706年(距离上次7年)

101-112

De Lagny

1719年(距离上次13年)

113-136

Jurij Vega

1794年(距离上次75年)

137-152

Rutherford

1841年(距离上次47年)

153-200

Zacharias Dase及Strassnitzky

1844年(距离上次3年)

201-248

Thomas Clausen

1847年(距离上次3年)

249-261

Lehmann

1853年(距离上次6年)

262-440

William Rutherford

1853年(距离上次0年)

441-500

Richter

1855年(距离上次2年)

501-527

William Shanks

1874年(距离上次19年)

528-620

D. F. Ferguson

1946年(距离上次72年)

621-710

1947年(距离上次1年)

711-808

1947年(距离上次0年)

备注:这里只列出人工计算的最高记录,808位

算准记录

记号

是第十六个希腊字母的小写。这个符号,亦是希腊语περιφρεια(表示周边、地域、圆周等意思)的首字母。1706年英国数学家威廉·琼斯(William Jones,1675—1749)最先使用“π”来表示圆周率。1736年,瑞士大数学家欧拉也开始用表示圆周率。从此,便成了圆周率的代名词。

要注意不可把和其大写Π混用,后者是指连乘的意思。

公式

圆周率()一般定义为一个圆形的周长()与直径()之比:,或直接定义为单位圆的周长的一半。由相似图形的性质可知,对于任何圆形,的值都是一样,这样就定义出常数

注意:将定义为单位圆的周长的一半是有意义的,这是因为从现代数学的角度来看,直径为d、半径为r的圆的周长C由以下积分给出:

上式中令,由定积分的换元法可得:

其中是单位圆周的周长(C的表达式中取r=1即得)。若定义,则,与我们熟知的周长公式相符。

而半径为r的圆的面积S由以下积分给出:

,由定积分的换元法可得:

其中是单位圆的面积(S的表达式中取r=1即得)。利用分部积分法,

于是,

因此,我们得到关系式:,这样一来也得到了我们熟知的圆面积公式

第二个做法是,以圆形半径为边长作一正方形,然後把圆形面积和此正方形面积的比例定为,即圆形之面积与半径平方之比。

定义圆周率不一定要用到几何概念,比如,我们可以定义为满足的最小正实数

这里的正弦函数定义为幂级数

特性

把圆周率的数值算得这么精确,实际意义并不大。现代科技领域使用的圆周率值,有十几位已经足够了。如果以39位精度的圆周率值,来计算可观测宇宙(observable universe)的大小,误差还不到一个原子的体积。以前的人计算圆周率,是要探究圆周率是否循环小数。自从1761年兰伯特证明了圆周率是无理数,1882年林德曼证明了圆周率是超越数后,圆周率的神秘面纱就被揭开了。

π在许多数学领域都有非常重要的作用。

代数

π是个无理数,即不可表达成两个整数之比,是由德国科学家约翰·海因里希·兰伯特于1761年证明的 。1882年,林德曼(Ferdinand von Lindemann)更证明了π是超越数,即π不可能是任何整系数多项式的根。

圆周率的超越性否定了化圆为方这古老尺规作图问题的可能性,因所有尺规作图只能得出代数数,而超越数不是代数数。

数学分析

Leibniz定理:

wallis公式:

高斯积分:

斯特林公式:

欧拉公式:

π的连分数表示:

数论

两个任意自然数是互质的概率是

任取一个任意整数,该整数没有重复质因子的概率为

一个任意整数平均可用个方法写成两个完全数之和。

概率论

设我们有一个以平行且等距木纹铺成的地板,随意抛一支长度比木纹之间距离小的针,求针和其中一条木纹相交的概率。这就是布丰投针问题。1777年,布丰自己解决了这个问题——这个概率值是1/π。

统计学

正态分布的概率密度函数:

物理学

海森堡不确定性原理:

相对论的场方程:

国际圆周率日

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万亿位,创下该常数迄今最精确值记录。

内容声明

1、本网站为开放性注册平台,以上所有展示信息均由会员自行提供,内容的真实性、准确性和合法性均由发布会员负责,本网站对此不承担任何法律责任。

2、网站信息如涉嫌违反相关法律规定或侵权,请发邮件至599385753@qq.com删除。

Copyright © 趣爱秀