圆周率

[yuán zhōu lǜ]

播报讨论上传视频

数学词汇

展开7个同名词条

查看我的收藏

0有用+1

本词条由中国科学院大学本科部、中国科学院自然科学史研究所参与编辑并审核，经科普中国·科学百科认证。

圆周率，是指圆周长与直径的比值，即圆周率=圆周长÷直径，也等于圆面积与半径平方之比，即圆周率=圆面积÷半径2，一般用希腊字母π表示。

圆周率是一个在数学及物理学中普遍存在的数学常数，是精确计算圆周长、圆面积、球体积等几何形状的关键值。在分析学里，π可以严格地定义为满足sinx=0的最小正数x。

圆周率是一个无理数，即无限不循环小数。在日常生活中，通常都用3.14代表圆周率去进行近似计算。而用九位小数3.141592654便足以应付一般计算。即使是工程师或物理学家要进行较精密的计算，充其量也只需取值至小数点后几百位^[1]。

1665年，英国数学家约翰·沃利斯（John Wallis）出版了一本数学专著^[24]，其中他推导出一个公式，发现圆周率等于无穷个分数相乘的积。2015年，罗切斯特大学的科学家们在氢原子能级的量子力学计算中发现了圆周率相同的公式^[2]。

2021年8月17日，瑞士研究人员使用一台超级计算机历时108天将圆周率计算到小数点后62.8万亿位，创下该常数迄今最精确值纪录^[22]。

中文名: 圆周率
外文名: Pi；Ratio of circumference to itsdiameter
符号表示: π

近似值: 22/7（约率）、355/113（密率）
属性: 无理数
计算位数: 105万亿位（2024年3月14日）^[22]^[32]

历史发展

播报

编辑

实验时期

一块古巴比伦石匾（约产于公元前1900年至公元前1600年）清楚地记载了圆周率=25/8=3.125^[4]。同一时期的古埃及文物，莱因德数学纸草书（Rhind Mathematical Papyrus）也表明圆周率等于分数16/9的平方，约等于3.1605^[4]。埃及人似乎在更早的时候就知道圆周率了。英国作家John Taylor（1781—1864）在其名著《金字塔》（《The Great Pyramid: Why was it built, and who built it?》）中指出，造于公元前2500年左右的胡夫金字塔和圆周率有关。例如，金字塔的周长和高度之比等于圆周率的两倍，正好等于圆的周长和半径之比。公元前800至600年成文的古印度宗教巨著《百道梵书》（Satapatha Brahmana^[29]）显示了圆周率等于分数339/108，约等于3.139^[5]。

几何法时期

古希腊作为古代几何王国对圆周率的贡献尤为突出。古希腊大数学家阿基米德（公元前287年—公元前212年）开创了人类历史上通过理论计算圆周率近似值的先河。阿基米德从单位圆出发，先用内接正六边形求出圆周率的下界为3，再用外接正六边形并借助勾股定理求出圆周率的上界小于4。接着，他对内接正六边形和外接正六边形的边数分别加倍，将它们分别变成内接正12边形和外接正12边形，再借助勾股定理改进圆周率的下界和上界。他逐步对内接正多边形和外接正多边形的边数加倍，直到内接正96边形和外接正96边形为止。最后，他求出圆周率的下界和上界分别为223/71和22/7，并取它们的平均值3.141851为圆周率的近似值。阿基米德用到了迭代算法和两侧数值逼近的概念，称得上是“计算数学”的鼻祖。

中国古算书《周髀算经》（约公元前2世纪）的中有“径一而周三”的记载，意即取

^[6]。汉朝时，张衡得出

，即

（约为3.162）。这个值不太准确，但它简单易理解^[7]。

公元263年，中国数学家刘徽用“割圆术”计算圆周率，他先从圆内接正六边形，逐次分割一直算到圆内接正192边形。他说：“割之弥细，所失弥少，割之又割，以至于不可割，则与圆周合体而无所失矣。”这包含了求极限的思想。刘徽给出π=3.141024的圆周率近似值，刘徽在得圆周率=3.14之后，将这个数值和晋武库中汉王莽时代制造的铜制体积度量衡标准嘉量斛的直径和容积检验，发现3.14这个数值还是偏小。于是继续割圆到1536边形，求出3072边形的面积，得到令自己满意的圆周率

。

公元480年左右，南北朝时期的数学家祖冲之进一步得出精确到小数点后7位的结果，给出不足近似值3.1415926和过剩近似值3.1415927，还得到两个近似分数值，密率

和约率

。密率是个很好的分数近似值，要取到

才能得出比

略准确的近似^[8]。（参见丢番图逼近）

在之后的800年里祖冲之计算出的π值都是最准确的。其中的密率在西方直到1573年才由德国人奥托（Valentinus Otho）得到，1625年发表于荷兰工程师安托尼斯（Metius）的著作中，欧洲称之为Metius' number。

约在公元530年，印度数学大师阿耶波多算出圆周率约为

。婆罗摩笈多采用另一套方法，推论出圆周率等于10的算术平方根。

阿拉伯数学家卡西在15世纪初求得圆周率17位精确小数值，打破祖冲之保持近千年的纪录。德国数学家鲁道夫·范·科伊伦（Ludolph van Ceulen）于1596年将π值算到20位小数值，后投入毕生精力，于1610年算到小数后35位数，该数值被用他的名字称为鲁道夫数。

分析法时期

这一时期人们开始利用无穷级数或无穷连乘积求π，摆脱可割圆术的繁复计算。无穷乘积式、无穷连分数、无穷级数等各种π值表达式纷纷出现，使得π值计算精度迅速增加。

第一个快速算法由英国数学家梅钦（John Machin）提出，1706年梅钦计算π值突破100位小数大关，他利用了如下公式：^[9]

其中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位^[25]。

2024年3月15日，据美国趣味科学网站报道，在国际圆周率日，总部位于美国加州的计算机存储公司Solidigm发布声明称，该公司已将圆周率Pi（π）计算到小数点后约105万亿位，打破此前100万亿位的世界纪录^[28]。

计算历史

日期	计算者	国籍	正确位数	详细纪录
前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年	弗朗索瓦·韦达	法国^[10]	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	^[11]
2011年	IBM“蓝色基因” 超级电脑			π2的前60,000,000,000,000位二进制小数^[12]

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

算准记录

小数点后位数	首次算准者	首次算准时间
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）最先使用“π”来表示圆周率^[13]。1736年，瑞士大数学家欧拉也开始用

表示圆周率。从此，

便成了圆周率的代名词^[14]。

要注意不可把

和其大写Π混用，后者是指连乘的意思。

公式

播报

编辑

圆周率

圆周率（

）一般定义为一个圆形的周长（

）与直径（

）之比：

，或直接定义为单位圆的周长的一半。由相似图形的性质可知，对于任何圆形，

的值都是一样，这样就定义出常数

。

注意：将

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

即

上式中令

，由定积分的换元法可得：

其中

是单位圆周的周长（C的表达式中取r=1即得）。若定义

，则

，与周长公式相符。

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

令

，由定积分的换元法可得：

其中

是单位圆的面积（S的表达式中取r=1即得）。利用分部积分法，

于是，

因此，得到关系式：

这样一来也得到了圆面积公式

第二个做法是，以圆形半径为边长作一正方形，然后把圆形面积和此正方形面积的比例定为

，即圆形之面积与半径平方之比。

定义圆周率不一定要用到几何概念，比如，定义

为满足

的最小正实数

。

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

圆周率的计算方法

播报

编辑

传统计算圆周率的方法依赖于几何构造、简单的数列展开或解析方法，例如割圆术、泰勒级数法等。现代计算圆周率的方法主要依赖于复杂的迭代算法、快速收敛的级数与强大的计算机的计算能力，例如查德诺夫斯基算法 (Chudnovsky Algorithm)、高斯-勒让德算法等。

割圆术

割圆术的基本思想是通过在圆内外分别嵌入正多边形，并逐步增加多边形的边数，使其逐渐接近圆的周长。由于多边形的周长可以精确计算，通过这种方法可以获得圆周长的上界和下界，从而逐步逼近圆周率。如图所示，在一个圆中嵌入一个内接的正多边形（如正六边形）。正六边形的边长等于圆的半径，并可通过几何关系递推得到边数更多的多边形的边长。设正n边形的边长为

，推导边的关系如下：

BC即内接正12边形的边长。归纳得到

与

之间的递推式^[31]：

所以，由

可计算出所有正内接正多边形的边长

,... 进而计算相应的周长，得到一列递增的数列，圆的周长就是数列的极限值^[31]。用圆周长除以圆直径就可以得到π的值。

泰勒级数法

的泰勒展开式为^[30]：

带入x=1，上述式子可化简为：

已知

，因此有：

即：

这个式子也被称为莱布尼茨级数(Gregory-Leibniz series)。这个公式非常简单，但收敛速度很慢，需要计算非常多的项才能得到准确的结果。

马青(Machin)公式

用

的泰勒级数来计算π收敛速度太慢，实际应用中不实用。要使收敛速度较快，|x| 应当远小于1^[30]。例如

就收敛得较快。利用

，可以通过倍角公式得到：

由于

，所以

。

这虽然能够得到一个近似值，但为了进一步修正，需要计算他们的误差

^[30]。

由

，

可得：

因此：

于是得到：

最终推导得出：

这个公式被称为马青（Machin）公式^[30]。

高斯-勒让德算法

这是一个快速收敛的迭代算法，可以非常高效地计算 π。它使用以下的迭代过程：

初始条件：

迭代公式：

经过几轮迭代后，π 的近似值为：

这个算法的收敛速度非常快，每次迭代大约可以提供倍增的精度。

查德诺夫斯基算法 (Chudnovsky Algorithm)

这是目前用于计算 π 最常用的算法之一，它是由两位数学家格雷戈里·查德诺夫斯基和大卫·查德诺夫斯基兄弟提出的。它基于快速收敛的级数，并且特别适合于大规模的高精度计算。公式如下：

这种方法可以非常快速地逼近 π，被用于大规模 π 的计算。

应用

播报

编辑

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

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

几何

平面图形	周长	面积
圆
圆环
扇形
注：① 为周长，为面积，为弧长；为直径，为半径（内圆半径），为外圆半径，为圆心角度数。 ②周长、弧长用长度单位，面积用面积单位。

立体图形	表面积	体积
圆柱



圆锥
圆锥
注：① 为底面周长，为底面积，为侧面积，为表面积，为体积；为底面直径，为底面半径，为高。 ②底面周长用长度单位，表面积（含底面积和侧面积）用面积单位，体积用体积单位或容积单位。

代数

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

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

数学分析

莱布尼茨(Leibniz)公式：

华里士(Wallis)公式：

当

时：

可以化为：

高斯积分：

斯特林公式

欧拉公式：

π的连分数表示：

数论

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

。

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

。

一个任意整数平均可用

个方法写成两个完全数之和。

概率论

在平面上画有一组间距为d的平行线，求一根长度为l（l＜d）的针任意投掷在这个平面上，求此针与平行线中和任一条相交的概率。这就是布丰投针问题。1777年，布丰自己本人证明了这个概率为

^[27]^[30]。

统计学

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

物理学

海森堡不确定性原理：

相对论的场方程：

国际圆周率日

播报

编辑

2011年，国际数学协会正式宣布，将每年的3月14日设为国际数学节，来源则是中国古代数学家祖冲之的圆周率。^[16]

国际圆周率日可以追溯至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位开始就算错了^[17]。

在谷歌公司2005年的一次公开募股中，共集资四十多亿美元，A股发行数量是14,159,265股，这当然是由π小数点后的位数得来。（顺便一提，谷歌公司2004年的首次公开募股，集资额为$2,718,281,828，与数学常数e有关^[18]）

排版软件TeX从第三版之后的版本号为逐次增加一位小数，使之越来越接近π的值：3.1，3.14，……当前的最新版本号是3.1415926^[19-20]。

每年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）^[21]。

2019年3月14日，谷歌宣布日裔前谷歌工程师爱玛（Emma HarukaIwao）在谷歌云平台的帮助下，计算到圆周率小数点后31.4万亿位，准确的说是31415926535897位，比2016年创下的纪录又增加数万亿位。据了解，爱玛的团队使用了一个名为ycruncher的程序，能将π计算到小数点后数万亿位。该程序由谷歌云平台计算引擎上运行的25个虚拟机驱动。而2016年纪录的创造者皮特（Peter Trueb）是用一台电脑计算出来的。这项计算需要170TB的数据，与整个美国国会图书馆印刷藏品数据量大致相同，爱玛经过大约4个月的计算才打破了此前的世界纪录^[3]。

2020年，一个名为北阿拉巴马慈善计算的非营利组织的创始人蒂莫西·穆利肯使用个人电脑，将数值计算到小数点后50万亿位，耗时303天^[23]。

2021年8月17日，美国趣味科学网站报道，瑞士研究人员使用一台超级计算机，历时108天，将著名数学常数圆周率π计算到小数点后62.8万亿位，创下该常数迄今最精确值记录^[22]。

2022年6月，谷歌云（Google Cloud）在新闻稿中表示，谷歌的云服务已经打破了其在2019年创造的纪录，计算出了 100 万亿位圆周率数字，上一次，谷歌云把圆周率精确到了 31.4 万亿位^[26]。

2024年的国际圆周率日（3月14日），总部位于美国加州的计算机存储公司Solidigm发布声明称，该公司已将圆周率Pi（π）计算到小数点后约105万亿位，打破此前100万亿位的世界纪录^[32]。

新手上路

成长任务编辑入门编辑规则本人编辑

我有疑问

内容质疑在线客服官方贴吧意见反馈

投诉建议

举报不良信息未通过词条申诉投诉侵权信息封禁查询与解封

圆周率

目录