圓周率 圓周率 識別 種類 無理數 超越數 符號 π {\displaystyle \pi } 位數 數列編號 A000796 性質 定義 π = C d . {\displaystyle \pi ={\frac {C}{d}}.} ,其中 C {\displaystyle C} 為圓周長、 d {\displaystyle d} 為直徑 π = ∫ − 1 1 d x 1 − x 2 . {\displaystyle \pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.} 連分數 π = 3 + 1 7 + 1 15 + 1 1 + 1 292 + 1 1 + 1 1 + 1 1 + ⋱ {\displaystyle \pi =3+\textstyle {\frac {1}{7+\textstyle {\frac {1}{15+\textstyle {\frac {1}{1+\textstyle {\frac {1}{292+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\ddots }}}}}}}}}}}}}}} 以此為根 的多項式或函數 e i x + 1 = 0 {\displaystyle e^{ix}+1=0} 表示方式 值 π ≈ {\displaystyle \pi \approx } 3.14159265無窮級數 π = ∑ k = 0 ∞ 4 ( − 1 ) k 2 k + 1 {\displaystyle \pi =\sum _{k=0}^{\infty }{\frac {4(-1)^{k}}{2k+1}}} 二进制 11.00100100 0011 1111 0110 … [ 1] 十进制 3.14159265 3589 7932 3846 … 十六进制 3.243F6A88 85A3 08D3 1319 … [ 2] :242 六十进制 3;8,29,44,0,47,25,53,7,24,57,36…[ 3] [ 4]
圓周率 是数学常数 ,为圆 的周长 和其直径 的比 ,近似值 约3.14159265,常用符号 π {\displaystyle \pi } 表示。
π {\displaystyle \pi } 是无理数 ,不能用分数 表示出来(即它的小数 部分是无限不循环小数 ),但近似 22 7 {\textstyle {\frac {22}{7}}} 等有理数。學界認為π的数字序列在统计上 是随机分布 ,但迄今未能证明。此外,π还是超越数 ,它不是任何有理 系数 多项式 的根 ,化圆为方 的问题不可能用尺规作图 解决。
几个文明古国 很早就須计算出π的精确值以便于生产的计算。西元5世纪,中國劉宋 数学家祖冲之 用几何方法将圆周率计算到小数点后7位。大约同时,印度 数学家也将圆周率计算到小数点后5位。史上首條π的精确无穷级数 公式(即π的莱布尼茨公式 )直到约1000年后才由印度数学家发现。[ 5] [ 6] 微積分 出現,π的位數很快計到數百位,足以滿足任何科學工程的計算需求。在20和21世纪,计算机技术 快速发展,π的计算精度急速提高。截至2024年3月,π的十进制 精度已達105万亿位。[ 7] 几乎所有科学研究对π的精度要求都不超过几百位,当前计算π的值主要都为打破记录、测试超级计算机 的计算能力和高精度乘法算法 。[ 2] :17 [ 8]
π的定义涉及圆,在三角学 和几何学 的许多公式,特别是广泛应用在圆形、球形或椭球形相關公式中。[ 9] 在近代數學分析 裡,π改由實數 系統譜性質中的特征值 或週期 來定義,其他數學領域 如數論 、統計 以及幾乎所有 物理學 領域均有出現,π的广泛用途使它成为科学界内外最广为人知的数学常数 。几本专门介绍π的书籍经已出版,圆周率日(3月14日)和π值计算突破记录也往往会成为报纸的新闻头条。[ 10] 此外,背诵 π值的世界记录已达10萬位。[ 11]
直徑為一的圓的周長是π(3.14159265...) 数学家用小写希腊字母 π {\displaystyle \pi } 表示圆周和其直径之比,有时也将其拼写为“Pi ”,来自希腊语“περίμετρος ”(周长)的首字母。[ 12] 英语π的发音与英文单词“Pie”( ,西式馅饼 )相同。[ 13] π的小写字母(或其无衬线体 )在数学要和表示连乘积 的大写Π相区分开。
关于选择符号π的原因,请参见引入π符号 一节。
圆周长略大于其直径的三倍;精确的比例称为 π {\displaystyle \pi } π常用定义为圆 的周长 C {\displaystyle C} 与直径 d {\displaystyle d} 的比值:[ 2] :8
π = C d {\displaystyle \pi ={\frac {C}{d}}} 。 无论圆的大小如何,比值 C d {\displaystyle {\frac {C}{d}}} 为恒值。如果圆的直径变为原先的二倍,周长也变为二倍,比值 C d {\displaystyle {\frac {C}{d}}} 不变。π目前的定义暗地用了欧几里得几何 的一些定理,虽然圆的定义可扩展到任意曲面(即非欧几里得几何 ),但这些圆不符合定律 π = C d {\displaystyle \pi ={\frac {C}{d}}} 。[ 2]
这里,圆的周长指其圆周的弧长 ,弧长这概念可以不依赖几何学,而是用微积分 学的极限 来定义。[ 14] 例如,若想计算笛卡儿坐标系 中单位圆 x 2 + y 2 = 1 {\displaystyle x^{2}+y^{2}=1} 上半部分的弧长,需要用到积分 :[ 15]
π = ∫ − 1 1 d x 1 − x 2 . {\displaystyle \pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.} 上述积分是由卡尔·魏尔斯特拉斯 于1841年对π的积分定义。[ 16]
π这些依赖周长、且暗地依赖积分 的定义如今在文献中并不常见。雷默特(Remmert (1991 ))解释说现代教微积分時,大学一般将微分学 课程安排在积分学课程之前,所以不依赖于后者的π的定义就很有必要了。其中一种定义由理查·巴爾策 提出,[ 17] 由愛德蒙·蘭道 推广,[ 18] 其表述如下:π是两倍于能使余弦函数 等于零的最小正数。[ 2] [ 15] [ 19] 余弦函数可以由独立于几何之外的幂级数 [ 20] 定义,或者使用微分方程 的解来定义。[ 19]
在相似的启发下,π可以用关于复变量 z {\displaystyle z} 的复指数 函数 exp ( z ) {\displaystyle \exp(z)} 来定义。复指数类似余弦函数,可用多种方式定义。令函数 exp ( z ) {\displaystyle \exp(z)} 值为一的复数集合是如下所示的(虚)等差數列:
{ … , − 2 π i , 0 , 2 π i , 4 π i , … } = { 2 π k i | k ∈ Z } {\displaystyle \{\dots ,-2\pi i,0,2\pi i,4\pi i,\dots \}=\{2\pi ki|k\in \mathbb {Z} \}} , 并且其中包括独特的正实数 π {\displaystyle \pi } 。[ 15] [ 21]
基于同样想法但更抽象的定义运用了精巧的拓扑学 和代数学 概念,用以下定理描述:[ 22] 存在唯一的从加法模数整数组成的实数群R /Z 到绝对值为1的复数组成的乘法群的连续同态 (拓扑学概念,指在拓扑空间 之间的一种态射 )。数字π定义为此同态 派生的模的一半。[ 23]
周长固定,圆会围成最大面积,π同樣表述为等周不等式 中出现的常数(乘四分之一)。此外,在很多其他紧密相关的方程中,π作为某些几何或者物理过程的特征值 出现;详见下文 。
π是无理数 ,无法表示成两整数之比 的形式(形如 22 7 {\textstyle {\frac {22}{7}}} 的分数常用来近似表达π,但是没有任何普通分数 (指整数的比)可以取到π的精确值)。[ 2] :5 由于 π {\displaystyle \pi } 是无理数,故可表示为无限不循环小数。有多种方法能证明π是无理数 ,这些证明也都要用到微积分学 和反证法 。 π {\displaystyle \pi } 可以用有理数 来近似的程度還無法準確得知(稱為無理性度量 ),不過估計其無理性度量比e 或ln(2) 的要大,但是小於刘维尔数 的無理性度量[ 24] 。
統計隨機性 检验,包括正规数 检验,可验证 π {\displaystyle \pi } 的位數沒有明顯的固定模式。 π {\displaystyle \pi } 的小数中任意固定长度的序列(如3位數000,001……999)出現機率都相同[ 25] 。不過有關π 是正规数 的猜想既無證明,亦無证伪 [ 2] :22-23 [ 25] 。
電腦出現後可生成大量π 的不同位数,并統計分析之。金田康正 詳細統計分析了π 的十進制數字,并验证了其分布正规:例如,假設檢定 0到9十個數的出現頻率,找不到有特定重复规律的證據[ 2] :22, 28–30 。根據無限猴子定理 ,任何任意長度、由隨機內容組成的子序列看起來都有可能像不隨機生成。因此,就算π 的小数序列通過了隨機性統計測試,其中也可能有幾位的數字看起來似有规律可循而非隨機数,例如π 的十進制写法在小數第762位后开始出现了連續六個9 [ 2] :3 。
由于π 是超越數 ,不能利用尺规作图 化圓為方 。 π {\displaystyle \pi } 不仅是无理数,还是超越数 ,即 π {\displaystyle \pi } 不是任何有理系数多项式 的根 。(比方说,试图解有限项方程 x 5 120 − x 3 6 + x = 0 {\textstyle {\frac {x^{5}}{120}}-{\frac {x^{3}}{6}}+x=0} 来求 π {\displaystyle \pi } 的值)[ 26] [ 註 1]
π {\displaystyle \pi } 的超越性衍生出一些重要的结果: π {\displaystyle \pi } 不能经有限次四则运算和开平方运算有理数来获得,因此不是规矩数 。换言之,尺规作图 作不出长度为 π {\displaystyle \pi } 的线段,也就不可能用尺规方法做出与已知圆面积相等的正方形。后者即为有名的化圓為方 问题,该问题早在古典时代 即已提出,曾困扰人数千年之久[ 27] [ 28] 。直至今天,依然有民间数学爱好者 声称他们解决了这问题[ 29] 。
π {\displaystyle \pi } 像所有无理数一样无法表示成分数 ,但 π {\displaystyle \pi } 等全部无理数都能表示成一系列叫连分数 的连续分数形式:
π = 3 + 1 7 + 1 15 + 1 1 + 1 292 + 1 1 + 1 1 + 1 1 + ⋱ {\displaystyle \pi =3+\textstyle {\frac {1}{7+\textstyle {\frac {1}{15+\textstyle {\frac {1}{1+\textstyle {\frac {1}{292+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\ddots }}}}}}}}}}}}}}} 在这连分数的任意一点截断化简,都能得到π的近似值;前四位近似值是3、 22 7 {\displaystyle {\frac {22}{7}}} 、 333 106 {\displaystyle {\frac {333}{106}}} 、 355 113 {\displaystyle {\frac {355}{113}}} 。这些数在历史上是π最广为人知且广為使用的几个近似值。用以上方式得出的 π {\displaystyle \pi } 的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近π。[ 30] π是超越数,据定义来说它不是代數數 ,又因此不可能是二次無理數 ;是故π不能表示为循环连分数 。尽管 π {\displaystyle \pi } 的简单连分数没有表现出任何其他明显规律,[ 31] 数学家發現了数條广义连分数能表示π,例如:[ 32]
π = 4 1 + 1 2 2 + 3 2 2 + 5 2 2 + 7 2 2 + 9 2 2 + ⋱ = 3 + 1 2 6 + 3 2 6 + 5 2 6 + 7 2 6 + 9 2 6 + ⋱ = 4 1 + 1 2 3 + 2 2 5 + 3 2 7 + 4 2 9 + ⋱ {\displaystyle \pi =\textstyle {\cfrac {4}{1+\textstyle {\frac {1^{2}}{2+\textstyle {\frac {3^{2}}{2+\textstyle {\frac {5^{2}}{2+\textstyle {\frac {7^{2}}{2+\textstyle {\frac {9^{2}}{2+\ddots }}}}}}}}}}}}=3+\textstyle {\frac {1^{2}}{6+\textstyle {\frac {3^{2}}{6+\textstyle {\frac {5^{2}}{6+\textstyle {\frac {7^{2}}{6+\textstyle {\frac {9^{2}}{6+\ddots }}}}}}}}}}=\textstyle {\cfrac {4}{1+\textstyle {\frac {1^{2}}{3+\textstyle {\frac {2^{2}}{5+\textstyle {\frac {3^{2}}{7+\textstyle {\frac {4^{2}}{9+\ddots }}}}}}}}}}} 圆周率近似值 包括:
整数 :3 分数 (依准确度顺序排列):13 / 4 、16 / 5 、19 / 6 、22 / 7 、179 / 57 、267 / 85 、333 / 106 、355 / 113 、52163 / 16604 、53228 / 16943 、55358 / 17621 、57843 / 18412 、60328 / 19203 、103993 / 33102 、245850922 / 78256779 [ 30] (选自 A063674 及 A063673 。) 小數 (整数后首80位):3.1415926535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899... [ 2] :240 (另见 A000796 ) 其他进位制的近似值
二进制 (整数后首48位):11.00100100 0011 1111 0110 1010 1000 1000 1000 0101 1010 0011… 十六进制 (整数后首20位):3.243F6A88 85A3 08D3 1319… [ 2] :242 六十进制 (整数后首20位):3;8,29,44,0,47,25,53,7,24,57,36,17,43,4,29,7,10,3,41,17…[ 3] [ 4] 欧拉公式 给出了e 的复指数与复平面 上以原點 为圆心的单位圆 上的点 之间的关系。 任何复数 (以 z {\displaystyle z} 为例)都可以表示为一组实数 对:极坐标系 用实数 r {\displaystyle r} 表示半径 ,代表复平面 上复数 z {\displaystyle z} 离原點的距离;实数 φ {\displaystyle \varphi } 则表示夹角,即这条半径(复平面上复数 z {\displaystyle z} 与原点的连线)与正实轴经顺时针转动的夹角。这样一来, z {\displaystyle z} 就可写成[ 33]
z = r ⋅ ( cos φ + i sin φ ) {\displaystyle z=r\cdot (\cos \varphi +i\sin \varphi )} ,这里 i {\displaystyle i} 代表虛數單位 ,即 i 2 {\displaystyle i^{2}} =-1。 复分析 中,欧拉公式 将三角函数 与复指数函数 糅合在一起[ 34] :
e i φ = cos φ + i sin φ {\displaystyle e^{i\varphi }=\cos \varphi +i\sin \varphi } ,这里数学常数e 是自然對數 的底数。 欧拉公式确立了 e {\displaystyle e} 的复指数与复平面上以原点为圆心的单位圆 上的点之间的关系,而且当 φ = π {\displaystyle \varphi =\pi } 时,欧拉公式就能改写为歐拉恆等式 的形式:
e i π + 1 = 0 {\displaystyle {{{e}^{{i}\,{\pi }}}+{1}}=0} 。此等式亦稱“最奇妙的数学公式”(英語:the most remarkable formula in mathematics ),全因它将五个最基本的数学常数简洁联系起来[ 34] [ 35] 。 欧拉等式亦可用于求出方程 z n = 1 {\displaystyle z^{n}=1} 的 n {\displaystyle n} 个不同复数根(这些根叫做 n {\displaystyle n} 次单位根 ”[ 36] ),可以根据以下公式求得:
e 2 π i k / n ( k = 0 , 1 , 2 , … , n − 1 ) {\displaystyle e^{2\pi ik/n}\qquad (k=0,1,2,\dots ,n-1)} 。 震盪弦的泛音 是二次微分的本徵函數 ,會形成泛音列 。對應的本徵值會形成由π 整數倍組成的等差数列 π {\displaystyle \pi } 常出现 在有关几何的问题中。然而,不少和几何无关的问题也可看到 π {\displaystyle \pi } 的身影。
π {\displaystyle \pi } 在許多用處中都會以特征值 形式出現。例如理想的振動弦 問題可以建模為函數 f {\displaystyle f} 在單位區間 [ 0 , 1 ] {\displaystyle [0,1]} 的圖形,固定邊界值 为 f ( 0 ) = f ( 1 ) = 0 {\displaystyle f(0)=f(1)=0} 。弦振動的模態會是微分方程 的 f n ( x ) + λ 2 f ( x ) = 0 {\displaystyle f^{n}(x)+\lambda ^{2}f(x)=0} ,此處λ 是相關的特徵值。受施图姆-刘维尔理论 限制, λ {\displaystyle \lambda } 只能是一些特定的數值。而 λ = π {\displaystyle \lambda =\pi } 即為一個特征值,因為函數 f ( x ) = sin ( π x ) {\displaystyle f(x)=\sin(\pi x)} 滿足邊界條件及微分方程 λ = π {\displaystyle \lambda =\pi } [ 37] 。
依照第一代开尔文男爵威廉·汤姆森 所述的一篇傳說,古迦太基 城的外形是等周長問題的一項解(Thompson 1894 )。這些包圍著海的區域由迦太基女王狄多 所圍,城不靠海的邊界須用指定大小的牛皮圍住,後來是將牛皮剪成小段 π {\displaystyle \pi } 是上述方程的最小特征值,也和弦振動的基本模式 有關。一種讓弦振動的方式是提供弦能量,能量會滿足維廷格函數不等式 [ 38] ,其中提到若函數 f : [ 0 , 1 ] → C {\displaystyle f:[0,1]\rightarrow \mathbb {C} } 使得 f ( 0 ) = f ( 1 ) = 0 {\displaystyle f(0)=f(1)=0} ,且 f {\displaystyle f} 和 f ′ {\displaystyle f'} 都是平方可積函數 ,則以下的不等式成立:
π 2 ∫ 0 1 | f ( x ) | 2 d x ≤ ∫ 0 1 | f ′ ( x ) | 2 d x , {\displaystyle \pi ^{2}\int _{0}^{1}|f(x)|^{2}\,dx\leq \int _{0}^{1}|f'(x)|^{2}\,dx,} 此例中等號成立的條件恰好是 f {\displaystyle f} 為 sin ( π x ) {\displaystyle \sin(\pi x)} 倍數的時候。因此 π {\displaystyle \pi } 似乎是維爾丁格不等式的最佳常數,也是最小的特征值(根據雷利商數 的計算方式)
π {\displaystyle \pi } 在更高維度的分析也有類似的角色,出現在其他類似問題的特徵值中。就如以上所述 , π {\displaystyle \pi } 的一項特點是等周定理 中的最佳常數:周長為 P {\displaystyle P} 的平面若尔当曲线 ,所圍面積 A {\displaystyle A} 滿足以下的不等式
4 π A ≤ P 2 {\displaystyle 4\pi A\leq P^{2}} , A = π r 2 {\displaystyle A=\pi r^{2}} 及 P = 2 π r {\displaystyle P=2\pi r} ,故等號成立的條件是曲线為圓形[ 39] 。
圓周率π 也和庞加莱不等式 的最佳常數有關[ 40] , π {\displaystyle \pi } 是一維及二維的狄氏能量 特征向量 最佳值中最小,會出現在許多經典的物理現象中,例如經典的位势论 [ 41] [ 42] [ 43] 。其一維的情形即為維廷格不等式。
圓周率π 也是傅里叶变换 的重要常數,傅里叶变换屬於积分变换 ,將實數線上有複數值、可積分的函數,轉換為以下形式:
f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x . {\displaystyle {\hat {f}}(\xi )=\int _{-\infty }^{\infty }f(x)e^{-2\pi ix\xi }\,dx.} 傅里叶变换有幾種不同的寫法,但不論怎麼寫,傅里叶变换及反傅里叶变换中,一定會有某處出現 π {\displaystyle \pi } 。不過上述的定義是最經典的,因為其描述了L 2 空間 中唯一的幺正算符 ,也是 L 1 {\displaystyle L^{1}} 空間到 L ∞ {\displaystyle L^{\infty }} 空間的代數同態 [ 44] 。
不确定性原理 也用到 π {\displaystyle \pi } 。不确定性原理提出了可以將函數在空間及在頻域中局部化程度的下限,用傅立葉轉換的方式表示:
∫ − ∞ ∞ x 2 | f ( x ) | 2 d x ∫ − ∞ ∞ ξ 2 | f ^ ( ξ ) | 2 d ξ ≥ ( 1 4 π ∫ − ∞ ∞ | f ( x ) | 2 d x ) 2 {\displaystyle \int _{-\infty }^{\infty }x^{2}|f(x)|^{2}\,dx\ \int _{-\infty }^{\infty }\xi ^{2}|{\hat {f}}(\xi )|^{2}\,d\xi \geq \left({\frac {1}{4\pi }}\int _{-\infty }^{\infty }|f(x)|^{2}\,dx\right)^{2}} 。 物理的結果,有關量子力学 中同時觀測位置及動量的不確定性,見下文 。傅立葉分析 中出現π 是史東-凡紐曼定理 的結果,證實了海森伯群 的薛定諤表示 是唯一[ 45] 。
高斯函数 f ( x ) = e − x 2 {\textstyle f(x)=e^{-x^{2}}} 的图像,函数下方与X轴围成的阴影部分面积为 π {\textstyle {\sqrt {\pi }}} 。 高斯积分 是对高斯函数 e − x 2 {\displaystyle e^{-x^{2}}} 在整条实轴上的积分,即函数下方与X轴围成的面积,其结果为 π {\displaystyle {\sqrt {\pi }}} ,
∫ − ∞ ∞ e − x 2 d x = π {\displaystyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}} 此积分的计算可以先计算 f ( x ) = e − x 2 {\displaystyle f(x)=e^{-x^{2}}} 对整条实轴的积分的平方,通过转换笛卡尔坐标系 为极坐标系 从而求得
( ∫ − ∞ ∞ e − x 2 d x ) 2 = ∬ R 2 e − ( x 2 + y 2 ) d x d y = ∫ 0 2 π ∫ 0 ∞ e − r 2 r d r d θ = π {\displaystyle \left(\int _{-\infty }^{\infty }e^{-x^{2}}\,dx\right)^{2}=\iint _{\mathbf {R} ^{2}}e^{-(x^{2}+y^{2})}\,dxdy=\int _{0}^{2\pi }\int _{0}^{\infty }e^{-r^{2}}r\,dr\,d\theta =\pi } 其他计算方法可参阅高斯积分 。高斯函数更一般的形式为 f ( x ) = a exp − ( x − b ) 2 2 c 2 {\textstyle f(x)=a\exp {\frac {-(x-b)^{2}}{2c^{2}}}} ,求一般形式的高斯积分均可通过换元积分法 转化为求 f ( x ) = e − x 2 {\textstyle f(x)=e^{-x^{2}}} 的积分。
另外,当高斯函数为以下形式时,它则是平均数 为 μ {\displaystyle \mu } 和標準差 为 σ {\displaystyle \sigma } 的正态分布 的機率密度函數 [ 46] :
f ( x ) = 1 σ 2 π exp − ( x − μ ) 2 2 σ 2 {\displaystyle f(x)={1 \over \sigma {\sqrt {2\pi }}}\,\exp {\frac {-(x-\mu )^{2}}{2\sigma ^{2}}}} 这函数是概率密度函数,函数下方与X轴围成的面积必须为1,令 μ = 0 {\displaystyle \mu =0} 和 σ = 1 {\displaystyle \sigma =1} 即可变换得出 ∫ − ∞ ∞ e − x 2 d x = π {\textstyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}} 。概率论 与统计学 领域经常使用正态分布来作为复杂现象的简单模型:例如科学家通常假设大多数试验观测值的随机误差都是服从正态分布[ 47] 。
由一维布朗运动 的反正弦定律,可以通过试验正信号相对于负信号领先权过零点的分布反过来推算π 概率论与统计学中的中心极限定理 解释了正态分布以及 π {\displaystyle \pi } 的核心作用,这定理本质上是联系着 π {\displaystyle \pi } 的谱特征 与海森堡 不确定性原理 相关的特征值,并且在不确定性原理中有
σ x σ p ≥ ℏ / 2 {\displaystyle \sigma _{x}\sigma _{p}\geq \hbar /2} , 这里的 σ x {\displaystyle \sigma _{x}} 與 σ p {\displaystyle \sigma _{p}} 分別為位置與動量的標準差 , ℏ {\displaystyle \hbar } 是約化普朗克常数 ,而不等式的等号当且仅当粒子的波函数为高斯函数使成立[ 48] 。
同样地, π {\displaystyle \pi } 作为唯一独特的常数使得高斯函数等于其自身的傅里叶变换 ,此时的高斯函数形式为 f ( x ) = e − π x 2 {\displaystyle f(x)=e^{-\pi x^{2}}} [ 49] 。根据豪(Howe )的说法,建立傅里叶分析基本定理的“全部工作(whole business)”简化为高斯积分。
圓周率在远古时期(西元前一千纪 )已估算至前两位(3.1)。有些埃及學家 聲稱,遠至古王國時期 時期的古埃及人已經用 22 7 {\textstyle {\frac {22}{7}}} 作為圓周率的約數[ 50] [ 註 2] ,但這說法受到質疑。[ 52] [ 53] [ 54] [ 55]
最早有記載的对圓周率估值在古埃及 和巴比伦 出现,兩估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一塊西元前1900至1600年的泥板 ,泥板上的幾何學陳述暗示人们当时把圓周率視同 25 8 {\textstyle {\frac {25}{8}}} (等於3.125)。[ 2] :167 埃及的莱因德数学纸草书 (鉴定撰寫年份為西元前1650年,但抄自一份西元前1850年的文本)載有用作計算圓面積的公式,该公式中圓周率等于 ( 16 9 ) 2 {\textstyle ({\frac {16}{9}})^{2}} (≈3.1605)。[ 2] :167
西元前4世紀的《百道梵書 》的天文學運算把 339 108 {\textstyle {\frac {339}{108}}} (≈3.139,精确到99.91%)用作圓周率估值[ 56] 。西元前150年前其他印度文獻把圓周率視為 10 {\textstyle {\sqrt {10}}} (≈3.1622)[ 2] :169 。
π 可以透過計算圓的外切多邊形及內接多邊形周長來估算 第一條有紀錄、嚴謹計算π 數值的演算法是用正多邊形的幾何算法,在西元前250年由希臘數學家阿基米德 發明。[ 2] :170 這算法用了有一千年之久,因而有時π 亦稱阿基米德常數。[ 2] :175、205 阿基米德的算法是在計算圓的外切正六邊形及內接正六邊形的邊長,以此計算 π {\displaystyle \pi } 的上限及下限,之後再將六邊形變成十二邊形,繼續計算邊長,一直計到正96邊形為止。他根據多邊形的邊長證明 223 71 < π < 22 7 {\textstyle {\frac {223}{71}}<\pi <{\frac {22}{7}}} (也就是 3.1408 < π < 3.1429 {\textstyle 3.1408<\pi <3.1429} )[ 57] 。阿基米德得到的上限 22 7 {\textstyle {\frac {22}{7}}} 也造成常見誤解,認為 π {\displaystyle \pi } 就等於 22 7 {\textstyle {\frac {22}{7}}} [ 2] :171 。在西元前150年,希臘羅馬的科學家克劳狄乌斯·托勒密 在《天文学大成 》一書中提到π 的數值是3.1416,可能來自阿基米德,也可能來自阿波罗尼奥斯 。[ 2] :176 [ 58] 數學家在1630年利用多邊形的方式計算π 到第39位小數,一直到1699年,其他數學家才利用無窮級數的方式打破其紀錄,計算到第71位小數[ 59] 。
阿基米德發展了用多邊形近似π 的計算方式 中国历史 上, π {\displaystyle \pi } 的數值有3[ 60] 、3.1547(西元前一世紀)、 10 {\displaystyle {\sqrt {10}}} (西元前100年,數值約3.1623)及 142 45 {\textstyle {\frac {142}{45}}} (第三世紀,數值約3.1556)[ 2] :176–177 。大約在西元265年,曹魏 數學家刘徽 創立割圆术 ,用正3072邊形計算出π 的數值為3.1416。[ 61] [ 2] :177 他後來又發明了較快的算法,利用邊數差兩倍的正多邊形,其面積的差值會形成等比數列,其公比為 1 4 {\textstyle {\frac {1}{4}}} 的原理,配合96邊形算出π的值為3.14。[ 61] 祖冲之 在西元480年利用割圆术計算12288邊形邊長,得到π的值在3.1415926和3.1415927之间。他同时提出了π的约率 22 7 {\textstyle {\frac {22}{7}}} 和密率 355 113 {\textstyle {\frac {355}{113}}} 。在之後的八百年內,這都是π 最準確的估計值。[ 2] :178 為紀念祖沖之 對圓周率發展的貢獻,日本 數學家三上義夫 將這推算值命名為“祖沖之圓周率”,簡稱“祖率”。[ 62]
印度天文學家阿耶波多 在西元499年的著作《阿里亞哈塔曆書 》中使用了3.1416的數值。[ 2] :179 斐波那契 在大約1220年用獨立於阿基米德多邊形法,計算出3.1418[ 2] :180 。義大利作家但丁·阿利吉耶里 用的數值則是 3 + 2 10 ≈ 3.14142 {\textstyle 3+{\frac {\sqrt {2}}{10}}\approx 3.14142} 。[ 2] :180
波斯天文學家卡西 在1424年利用3×228 邊的多邊形,計算到六十進制的第9位小數,相當十進制的第16位小數。[ 63] [ 64] 這一突破成為當時的紀錄,延續了約180年。[ 65] 法國數學家弗朗索瓦·韦达 在1579年用3×217 邊形計算到第9位小數[ 65] ,佛蘭芒數學家阿德里安·范·羅門 在1593年計算到第15位小數[ 65] 。荷蘭數學家鲁道夫·范·科伊伦 在1596年計算到第20位小數,他之後又計算到第35位小數(因此在二十世紀初之前,圓周率在德國會稱為鲁道夫數)。[ 2] :182–183 荷蘭科學家威理博·司乃耳 在1621年計算到第34位小數[ 2] :183 ,而奧地利天文學家克里斯托夫·格林伯格 在1630年用1040 邊形計算到第38位小數[ 66] ,至今這仍是利用多邊形算法可以達到最準確的結果[ 2] :183 。
比較幾條曾用來計π 的無窮級數的收斂情形。Sn 是只取前n項的近似值。每張圖都是對應前一張圖的陰影部份,然後放大橫軸10倍。(點擊察看細節) 16及17世紀時,開始改用無窮级数 的方式去計π。無窮级数是一組無窮數列 的和[ 2] :185–191 。無窮级数讓數學家可以計算出比阿基米德 以及其他用幾何方式計算的數學家更準確的結果。[ 2] :185–191 雖然詹姆斯·格雷果里 及戈特弗里德·莱布尼茨 等歐洲數學家利用無窮數列計算π 而使得该方法为大家所知,但这种方法最早是由印度 科學家在大約1400到1500年之間發現。[ 2] :185-186 [ 67] 第一個记载用無窮级数計算π 的人是约西元1500年左右时,印度天文學家尼拉卡莎·薩默亞士 在他的著作《系統匯編 》中用梵語詩所記錄。[ 68] 當時沒有這數列對應的證明,而證明出現在另一本較晚的印度作品《基本原理 》,年代約在西元1530年。尼拉卡莎將該數列歸功於更早期的印度數學家桑加馬格拉馬的馬德哈瓦 (1350–1425)。[ 68] 相關的無窮级数有許多,包括有關 sin {\displaystyle \sin } 、 tan {\displaystyle \tan } 及 cos {\displaystyle \cos } 的,現在稱為馬德哈瓦數列 或π的莱布尼茨公式 [ 68] 。瑪達瓦在1400年用無窮级数計算π 到第11位小數,但在1430年一位波斯數學家卡西 利用多邊形算法否定了他算的結果[ 69] 。
艾萨克·牛顿 利用無窮级数計算π 到第15位,後來寫道:「我很羞愧的告訴你我為了這個計算用了多少個數字。」[ 70] 歐洲發現的第一條無窮項圓周率公式 是無窮乘積 (和一般用來計算π 的無窮級數不同),由法國科學家弗朗索瓦·韦达 在1593年發現[ 2] :187 [ 71] :
2 π = 2 2 ⋅ 2 + 2 2 ⋅ 2 + 2 + 2 2 ⋯ {\displaystyle {\frac {2}{\pi }}={\frac {\sqrt {2}}{2}}\cdot {\frac {\sqrt {2+{\sqrt {2}}}}{2}}\cdot {\frac {\sqrt {2+{\sqrt {2+{\sqrt {2}}}}}}{2}}\cdots } 約翰·沃利斯 在1655年發現了沃利斯乘积 ,是歐洲發現的第二條無窮項圓周率公式[ 2] :187 :
π 2 = 2 1 ⋅ 2 3 ⋅ 4 3 ⋅ 4 5 ⋅ 6 5 ⋅ 6 7 ⋅ 8 7 ⋅ 8 9 ⋯ {\displaystyle {\frac {\pi }{2}}={\frac {2}{1}}\cdot {\frac {2}{3}}\cdot {\frac {4}{3}}\cdot {\frac {4}{5}}\cdot {\frac {6}{5}}\cdot {\frac {6}{7}}\cdot {\frac {8}{7}}\cdot {\frac {8}{9}}\cdots } 微积分学 由英國科學家艾萨克·牛顿 及德國數學家戈特弗里德·莱布尼茨 在1660年代發明,許多計π 的無窮級數出現。牛頓自己就用反正弦 ( arcsin {\displaystyle \arcsin } )數列在1655年或1666年將π 近似到第15位小數,後來寫到「我很羞愧告訴你我為了計算它用了多少數字,我當時沒有做其他事。」[ 70]
蘇格蘭數學家詹姆斯·格雷果里 在1671年發現了馬德哈瓦公式,莱布尼茨也在1674年發現:[ 2] :188–189 [ 72]
arctan z = z − z 3 3 + z 5 5 − z 7 7 + ⋯ {\displaystyle \arctan z=z-{\frac {z^{3}}{3}}+{\frac {z^{5}}{5}}-{\frac {z^{7}}{7}}+\cdots } 這公式即為格雷果里-莱布尼茨公式,在 z = 1 {\displaystyle z=1} 時數值為 π 4 {\displaystyle {\frac {\pi }{4}}} 。[ 72] 1699年時英國數學家亚伯拉罕·夏普 用格雷果里-莱布尼茨公式,在 z = 1 3 {\displaystyle z={\frac {1}{\sqrt {3}}}} 時計算,計算到π的第71位小數,打破由多邊形算法得到的第39位小數的记录。[ 2] :189 格雷果里-莱布尼茨公式在 z = 1 {\displaystyle z=1} 時非常簡單,但收斂 到最終值的速度非常慢,現在不会再用此公式來計π。[ 2] :156
約翰·梅欽 在1706年用格雷果里-莱布尼茨級數產生了可以快速收斂的公式:[ 2] :192–193
π 4 = 4 arctan 1 5 − arctan 1 239 {\displaystyle {\frac {\pi }{4}}=4\,\arctan {\frac {1}{5}}-\arctan {\frac {1}{239}}} 梅欽用這公式計到π第100位小數[ 2] :72–74 後來其他數學家也發展了一些類似公式,現在稱為梅欽類公式 ,創下了許多計算π位數的紀錄。[ 2] :72–74 在進入電腦時代時,梅欽類公式仍然是耳熟能详可以計算π的公式,而且在约250年的时间里,很多有關π位數的紀錄都是梅欽類公式所得,比如在1946年時由達尼爾·弗格森(Daniel Ferguson )用這類公式計到第620位小數,是沒有計算設備輔助的最佳紀錄。[ 2] :192–196, 205
1844年,計算天才扎卡里亞斯·達斯 在德國數學家卡爾·弗里德里希·高斯 的要求下以梅欽類公式心算了π的200位小數,並創下紀錄。[ 2] :194-196 英國數學家威廉·謝克斯 花了15年的時間計算π 到小數707位,不過第528位小數出錯,後面的小數也都不正確。[ 2] :194–196
有些π 的無窮級數收斂 的比其他級數要快,數學家一般會選用收斂速度較快的級數,可以在較少的計算量下計算π,且達到需要的準確度[ 73] [ 2] :15–17, 70–72, 104, 156, 192–197, 201–202 。以下是π 的莱布尼茨公式 :[ 2] :69–72
π = 4 1 − 4 3 + 4 5 − 4 7 + 4 9 − 4 11 + 4 13 − ⋯ {\displaystyle \pi ={\frac {4}{1}}-{\frac {4}{3}}+{\frac {4}{5}}-{\frac {4}{7}}+{\frac {4}{9}}-{\frac {4}{11}}+{\frac {4}{13}}-\cdots } 隨著一項一項的值加入總和中,只要項次夠多,總和最後會慢慢接近π。不過此數列的收斂速度很慢,要到50萬項之後,才會精確到π的第五位小數[ 74] 。
尼拉卡莎在15世紀發展了π的另一條無窮級數,收斂速度比格雷果里-萊布尼茨公式快很多:[ 75]
π = 3 + 4 2 × 3 × 4 − 4 4 × 5 × 6 + 4 6 × 7 × 8 − 4 8 × 9 × 10 + ⋯ {\displaystyle \pi =3+{\frac {4}{2\times 3\times 4}}-{\frac {4}{4\times 5\times 6}}+{\frac {4}{6\times 7\times 8}}-{\frac {4}{8\times 9\times 10}}+\cdots } 以下比較兩條級數的收斂速率:
π的無窮級數 第1項 前2項 前3項 前4項 前5項 收斂到 π = 4 1 − 4 3 + 4 5 − 4 7 + 4 9 − 4 11 + 4 13 ⋯ . {\displaystyle \pi ={\frac {4}{1}}-{\frac {4}{3}}+{\frac {4}{5}}-{\frac {4}{7}}+{\frac {4}{9}}-{\frac {4}{11}}+{\frac {4}{13}}\cdots .} 4.0000 2.6666… 3.4666… 2.8952… 3.3396… 3.1415… π = 3 + 4 2 × 3 × 4 − 4 4 × 5 × 6 + 4 6 × 7 × 8 ⋯ . {\displaystyle \pi ={3}+{\frac {4}{2\times 3\times 4}}-{\frac {4}{4\times 5\times 6}}+{\frac {4}{6\times 7\times 8}}\cdots .} 3.0000 3.1666… 3.1333… 3.1452… 3.1396…
計算前五項後,格雷果里-萊布尼茨級數的和跟π的誤差為0.2,而尼拉卡莎級數和的誤差為0.002。尼拉卡莎級數收斂快很多,也甚為適合用來計π的值。收斂更快的級數有梅欽類公式 及楚德诺夫斯基算法 ,後者每計一項就可以得到14位正確的小數位[ 73] 。
并非所有和π有关的研究都旨在提高计算它的准确度。1735年,欧拉 解决了巴塞尔问题 ,建立了所有平方数倒数 和与π的关系。之后欧拉发现了欧拉乘积公式 ,得到了π、素数 的重要關聯,對日後黎曼ζ函數 的研究影響深遠。[ 76]
π 2 6 = 1 1 2 + 1 2 2 + 1 3 2 + 1 4 2 + ⋯ {\displaystyle {\frac {\pi ^{2}}{6}}={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots } 1761年,瑞士数学家约翰·海因里希·朗伯 用正切函数 的无穷连分数表达式证明 了π是無理數 。[ 2] :5 [ 77] 1794年,法国数学家阿德里安-马里·勒让德 证明了 π 2 {\displaystyle \pi ^{2}} 也是无理数。1882年,德国数学家费迪南德·冯·林德曼 证明了对任何非零代数数 α {\displaystyle \alpha } , e α {\displaystyle e^{\alpha }} 都是超越数,该结论后来由魏尔斯特拉斯 推广为林德曼-魏尔斯特拉斯定理 。据此定理和欧拉公式,π只能是超越數 ,進而证实了勒让德和欧拉提出的π超越性猜想。[ 2] :196 [ 78] 哈代 在其著作《数论导引》中则称此证明在提出後,經過希尔伯特 、施瓦兹 和其他一些人化简过。[ 79]
萊昂哈德·歐拉 在他1736年到1748年的作品中開始用希臘字母π 表示圓周率,數學界也開始廣為使用 在用π 专指“圆周率”之前,希腊字母即已用於幾何概念中[ 2] :166 。威廉·奥特雷德 在1647年起在《數學之鑰》(Clavis Mathematicae)就已經用 π {\displaystyle \pi } 及 δ {\displaystyle \delta } (對應p和d的希臘字母)來表示圓的周長及直徑的比例。
威廉·琼斯 在他1706年出版的《新數學導論》(A New Introduction to the Mathematics )提到了 π {\displaystyle \pi } ,是目前已知最早专门用希臘字母 π {\displaystyle \pi } 表示圓周和其直徑比例的人[ 80] 。這希臘字母第一次出现是在书中討論一塊半徑1的圓時提到「其圓周長一半( π {\displaystyle \pi } )」。琼斯選用 π {\displaystyle \pi } 可能因它是希臘文“周边”一词“περιφέρεια ”的首字母[ 81] 。不過琼斯提到,他那些有關 π {\displaystyle \pi } 的算式出自「真正聰明的約翰·梅欽先生」,人们推測在瓊斯之前,約翰·梅欽 就已开始用 π {\displaystyle \pi } 表示圓周率[ 2] :166 。
瓊斯在1706年開始使用此希臘字母,但直到萊昂哈德·歐拉 在其1736年出版的《力學 》中開始使用之后,其他数学家才纷纷开始用 π {\displaystyle \pi } 指代圆周率。在此之前,數字家可能用像c 或p 之類的字母代表圓周率[ 2] :166 。歐拉與歐洲其他數學家间时常互相写信来往, π {\displaystyle \pi } 的用法迅速傳播开来[ 2] :166 。1748年歐拉在他的《无穷小分析引论 》再一次提到了 π {\displaystyle \pi } ,写道:「簡潔起見,我們將此數字寫為 π {\displaystyle \pi } , π {\displaystyle \pi } 等於半徑為1的圓周長的一半。」这表示方式之後也推展到整片西方世界 [ 2] :166 。
约翰·冯·诺伊曼 所屬的團隊是用數位計算機ENIAC 來計π 的第一隊 高斯-勒让德算法 :
一開始設定
a 0 = 1 b 0 = 1 2 t 0 = 1 4 p 0 = 1 {\displaystyle \scriptstyle a_{0}=1\quad b_{0}={\frac {1}{\sqrt {2}}}\quad t_{0}={\frac {1}{4}}\quad p_{0}=1} 迭代計算: a n + 1 = a n + b n 2 b n + 1 = a n b n {\displaystyle \scriptstyle a_{n+1}={\frac {a_{n}+b_{n}}{2}}\quad \quad b_{n+1}={\sqrt {a_{n}b_{n}}}}
t n + 1 = t n − p n ( a n − a n + 1 ) 2 p n + 1 = 2 p n {\displaystyle \scriptstyle t_{n+1}=t_{n}-p_{n}(a_{n}-a_{n+1})^{2}\quad \quad p_{n+1}=2p_{n}} 則π 的估計值為
π ≈ ( a n + b n ) 2 4 t n {\displaystyle \scriptstyle \pi \approx {\frac {(a_{n}+b_{n})^{2}}{4t_{n}}}} 二十世紀中期计算机技术发展、革新再次引发了計算π 位數的熱潮。美國數學家约翰·伦奇 及李維·史密斯在1949年用桌上型計算機計算到1120位[ 2] :205 。同年,喬治·韋斯納(George Reitwiesner)及约翰·冯·诺伊曼 帶領的團隊利用反三角函数 (arctan)的無窮級數 ,用ENIAC 計算到了小數後2037位,花了70小時的電腦工作時間[ 82] 。這紀錄後來多次由其他透過arctan級數计算出的結果打破(1957年到7480位小數,1958年到第一萬位數,1961年到第十萬位小數),直到1973年,小数点后第一百萬位小數經已算出[ 2] :197 。
1980年代有两项發明加速計算了π。第一项是發现了新的迭代法去计π 的值,計算速度比無窮級數快很多;另一项是發现了可以快速計算大數字乘積的乘法演算法 [ 2] :15–17 。電腦大部分的工作時間都是在計乘法,這類演算法對現代計π 格外重要[ 2] :131 。這類演算法包括嘉良對馬(Karatsuba)算法 、譚曲(Toom-Cook)乘法 及以傅里叶变换 為基礎的乘法演算法(傅里叶乘法)[ 2] :132, 140 。
迭代演算法最早是在1975年至1976年间分别由美國物理學家尤金·薩拉明 及奧地利科學家理查·布蘭特 独立提出[ 2] :87 。這两條演算法没有依赖無窮級數來計算。迭代會重覆特定計算,将前一次的計算結果作为這一次的輸入值,使得計算結果漸漸的趨近理想值。此方式的原始版本其實是在160年前由卡爾·弗里德里希·高斯 提出,現在稱為算术-几何平均数算法 (AGM法)或高斯-勒让德算法 [ 2] :87 。薩拉明及布蘭特都曾修改之,这算法也稱為薩拉明-布蘭特演算法。
迭代演算法收斂速度比無窮級數快很多,在1980年代以後廣為使用。無窮級數隨著項次的增加,一般來說正確的位數也會增加幾位,但迭代演算法每計算多一次,正確位數會呈几何级数增长。例如薩拉明-布蘭特演算法每計算多一次,正確位數會是之前的二倍。1984年加拿大人喬納森·波温 及彼得·波温 提出迭代演算法,每計算多一次,正確位數會是之前的四倍,1987年時有另一條迭代演算法,每計算多一次,正確位數會是之前的五倍[ 83] 。日本數學家金田康正 使用的演算法在1955年及2002年間創下了若干項紀錄[ 84] 。不過迭代演算法的快速收斂也有其代價,需要的記憶體明顯比無窮級數多[ 84] 。
當數學家發現新的算法、電腦變得普及时,π 的已知小數位急剧增加。注意垂直坐标使用了对数坐标 。 一般而言,π值并不需要过于精确便能够满足大部分数学运算的需求。按照約·安(Jörg Arndt)及古里斯佗夫·希奴(Christoph Haenel )的计算,39位精確度已可将可觀測宇宙 圆周的精确度準確至一粒原子大小,足以運算絕大多數宇宙学 的计算需求[ 85] 。尽管如此,和π有關的成就往往成為世界各地的新聞頭條;部分人出于對破紀錄的冲动,依然奋力算出π小数点后上千甚至上百萬位[ 2] :17–19 [ 86] [ 87] 。此外也有測試超级计算机 、測試数值分析算法(包括高精度乘法算法 )等實際好處。純粹數學這领域也能计算π的位数评定其隨機度[ 2] :18 。
斯里尼瓦瑟·拉马努金 的肖像,他在印度独立工作时提出了许多计算π的新颖数列。 现代计算π的程序不仅局限于迭代算法。20世纪80与90年代,出现了可用来计算π的新无穷级数,其收敛速度可与迭代算法媲美,而又有着复杂度、内存密集度更低的优势。[ 84] 印度数学家斯里尼瓦瑟·拉马努金 是这方面的先驱,他在1914年发表了许多与π相关的公式,这些公式十分新颖,极为优雅而又颇具数学深度,收敛速度也非常快。[ 2] :103–104 下式即为一例,其中用到了模方程 :
1 π = 2 2 9801 ∑ k = 0 ∞ ( 4 k ) ! ( 1103 + 26390 k ) k ! 4 ( 396 4 k ) . {\displaystyle {\frac {1}{\pi }}={\frac {2{\sqrt {2}}}{9801}}\sum _{k=0}^{\infty }{\frac {(4k)!(1103+26390k)}{k!^{4}(396^{4k})}}.} 这无穷级数收敛速度远快于绝大多数反正切数列,包括梅钦公式。[ 2] :104 第一位使用拉马努金公式计算π并取得进展的是比尔·高斯珀 ,他在1985年算得了小数点后一千七百万位。[ 2] :104, 206 拉马努金公式开创了现代数值近似算法的先河,此后波尔文兄弟和楚德诺夫斯基兄弟 进一步发展了这类算法。[ 2] :110–111 后者于1987年提出了楚德诺夫斯基公式 ,如下所示:
1 π = 12 640320 3 / 2 ∑ k = 0 ∞ ( 6 k ) ! ( 13591409 + 545140134 k ) ( 3 k ) ! ( k ! ) 3 ( − 640320 ) 3 k . {\displaystyle {\frac {1}{\pi }}={\frac {12}{640320^{3/2}}}\sum _{k=0}^{\infty }{\frac {(6k)!(13591409+545140134k)}{(3k)!(k!)^{3}(-640320)^{3k}}}.} 此公式每计算一项就能得到π的约14位数值[ 88] ,因而用於突破圆周率的数位的计算。利用这公式,楚德诺夫斯基兄弟于1989年算得π小数点后10亿(109 )位,法布里斯·贝拉 于2009年算得2.7千亿(2.7×1012 )位,亚历山大·易和近藤滋在2011年算得一万亿(1013 )位。[ 2] :110–111, 206 [ 89] [ 90] 类似的公式还有拉马努金-佐藤级数 。
2006年,加拿大数学家西蒙·普勞夫 利用PSLQ整数关系算法 [ 91] 按照以下模版生成了几條计算π的新公式:
π k = ∑ n = 1 ∞ 1 n k ( a q n − 1 + b q 2 n − 1 + c q 4 n − 1 ) {\displaystyle \pi ^{k}=\sum _{n=1}^{\infty }{\frac {1}{n^{k}}}\left({\frac {a}{q^{n}-1}}+{\frac {b}{q^{2n}-1}}+{\frac {c}{q^{4n}-1}}\right)} , q {\displaystyle q} 为e π {\displaystyle \pi } , k {\displaystyle k} 是奇数 , a , b , c {\displaystyle a,b,c} 是普勞夫计算出的有理常数。[ 92]
蒙特卡洛方法 基于随机试验结果计算
π {\displaystyle \pi } 的近似值
統計模擬法 是以概率统计理论为指导的一类非常重要的计数方法,經大量重复试验计算事件发生频率,按照大数定律 (即当试验次数充分大时,频率充分接近概率)可以求得 π {\displaystyle \pi } 的近似值[ 93] 。 布芬(Buffon)投針問題 就是其中一項實例:长度 l {\displaystyle l} 的针随机往画满间距 t ( l ≤ t ) {\displaystyle t\left(l\leq t\right)} 的平行线的平面上抛掷 n {\displaystyle n} 次, 如果针与平行直线相交 m {\displaystyle m} 次, n {\displaystyle n} 充分大就可根据以下公式算出 π {\displaystyle \pi } 的近似值[ 94] :
π ≈ 2 n ℓ m t {\displaystyle \pi \approx {\frac {2n\ell }{mt}}} 用統計模擬法计 π {\displaystyle \pi } 的另一例子是随机往内切四分之一圆的正方形内抛掷大量点,落在四分之一圆内的点的数量与抛掷点的总量的比值会近似于 π 4 {\textstyle {\frac {\pi }{4}}} 。[ 2] :39–40 [ 95]
此外还可用随机游走 试验,并用統計模擬法计算 π {\displaystyle \pi } 值,如抛掷一枚均匀的硬币 N {\displaystyle N} 次,并记录正面朝上的次数,所得结果中,正面朝上的次数 n N {\displaystyle n_{N}} 服从二項分佈 且
Pr ( n N = m ) = ( N m ) ( 1 2 ) m ( 1 2 ) N − m {\displaystyle \Pr(n_{N}=m)={\binom {N}{m}}({\frac {1}{2}})^{m}({\frac {1}{2}})^{N-m}} 因为硬币均匀,所以N 次试验中每次试验结果相互独立。由此可定义一系列独立的随机变量 X k ( k = 1 , 2 , … ) {\displaystyle X_{k}\left(k=1,2,\ldots \right)} ,当抛掷结果为正面时 X k = 1 {\displaystyle X_{k}=1} 否则为-1,且 X k = ± 1 {\displaystyle X_{k}=\pm 1} 且取何值有相同概率(即,正面朝上和背面朝上的概率相同)。对随机变量 X k ( k = 1 , 2 , … , N ) {\displaystyle X_{k}\left(k=1,2,\ldots ,N\right)} 求和可得
W N = ∑ k = 1 N X k {\displaystyle W_{N}=\sum _{k=1}^{N}X_{k}} 设k 为“硬币正面朝上的次数”减去“硬币反面朝上的次数”,即可得到 m − ( N − m ) = k {\displaystyle m-\left(N-m\right)=k} 。变换式子,得 m = N + k 2 {\displaystyle m={\frac {N+k}{2}}} ,因此
Pr ( W N = k ) = ( N N + k 2 ) 1 2 N {\displaystyle \Pr(W_{N}=k)={\binom {N}{\frac {N+k}{2}}}{\frac {1}{2^{N}}}} ,其中 k = − N , − N + 2 , − N + 4 , … , N − 2 , N {\displaystyle k=-N,-N+2,-N+4,\ldots ,N-2,N} 。 可证明[ 96] ,
E ( W N ) = 0 {\displaystyle E(W_{N})=0} , E ( W N 2 ) = N {\displaystyle E(W_{N}^{2})=N} ,以及 E ( | W N | ) = ( N ⌈ N / 2 ⌉ ⌈ N / 2 ⌉ 2 N − 1 ) = { ( N − 1 ) ! ! ( N − 2 ) ! ! , 若 N 偶, N ! ! ( N − 1 ) ! ! , 若 N 奇。 {\textstyle E(|W_{N}|)={\binom {N}{\left\lceil {N/2}\right\rceil {\frac {\left\lceil {N/2}\right\rceil }{2^{N-1}}}}}={\begin{cases}{\frac {(N-1)!!}{(N-2)!!}},&{\text{若 }}N{\text{偶,}}\\{\frac {N!!}{(N-1)!!}},&{\mbox{若 }}N{\mbox{奇。}}\end{cases}}} 并且当N 变大时, E ( | W N | ) {\textstyle E\left(\left\vert W_{N}\right\vert \right)} 的值会渐近于 2 N π {\textstyle {\sqrt {\frac {2N}{\pi }}}} ,因此当N 充分大时可根据以下公式算出 π {\displaystyle \pi } 的近似值:[ 97]
π ≈ 2 N | W N | 2 {\displaystyle \pi \approx {\frac {2N}{|W_{N}|^{2}}}} 和其他计算 π {\displaystyle \pi } 值的方法相比,蒙特卡洛方法收敛速度很慢,而且无论实验多少次,都无从得知 π {\displaystyle \pi } 的估值已经精确到第几位。因此,当追求速度或精度时,蒙特卡洛方法不适合用来估计 π {\displaystyle \pi } 。[ 2] :43 [ 98]
1995年引入的兩條算法开辟了研究 π {\displaystyle \pi } 的新途径。因为每计算出一位数字,該數就会像流过阀门的水一样不会再出现在后续的计算过程中,这种新進算法叫阀门算法 。[ 2] :77–84 [ 99] 这就与无穷级数及迭代算法形成对比——无穷级数和迭代算法自始至终的每一步计算都会涉及到之前所有步骤计算出的中间值。[ 2] :77–84
1995年,美國數學家斯坦·瓦格纳 和斯坦利·拉比諾維茨(Stanley Rabinowitz )发明了一种簡單的阀门算法[ 99] [ 2] :77 [ 100] ,其運算速度類似arctan演算法,但速度比迭代算法慢[ 2] :77 。
贝利-波尔温-普劳夫公式 (BBP)是另一條阀门算法,屬於一种位數萃取演算法 。1995年,西蒙·普勞夫 等人發現[ 2] :117, 126–128 [ 101]
π = ∑ k = 0 ∞ 1 16 k ( 4 8 k + 1 − 2 8 k + 4 − 1 8 k + 5 − 1 8 k + 6 ) {\displaystyle \pi =\sum _{k=0}^{\infty }{\frac {1}{16^{k}}}\left({\frac {4}{8k+1}}-{\frac {2}{8k+4}}-{\frac {1}{8k+5}}-{\frac {1}{8k+6}}\right)} 這公式和其他公式不同,可以計算 π {\displaystyle \pi } 的任何十六进 小數位,而不用計算前面全部小數位[ 2] :117, 126–128 。十六进数位可计算得到特定二进数位;想要得到八进制 数位的话,计算一、两位十六进小數即可。目前也已發現一些這種演算法的變體,不過還沒有发现針對十進制、可以快速生成特定小數位的位數萃取演算法[ 102] 。位數萃取演算法的一項重要用途是用來確認聲稱是計算到 π {\displaystyle \pi } 小數位數的新紀錄:若有聲稱是新紀錄的計算結果出現,先將十進制的數值轉換到十六進制,再用贝利-波尔温-普劳夫公式去確認最後一些位數(用亂數決定),若這些位數都對,就能有一定把握认为此計算結果是对的[ 90] 。
1998年到2000年間,分布式计算 計畫PiHex 用貝拉公式 (贝利-波尔温-普劳夫公式的一種變體)計算 π {\displaystyle \pi } 第1015 位,結果是0[ 2] :20 [ 103] 。2010年9月,有雅虎 員工用公司的Apache Hadoop 應用程式在上千台電腦計算π在2×1015 位开始往后256位,其第2×1015 位剛好也是0[ 104] 。
利用伽瑪函數計算
伽瑪函數, Γ ( n ) = ( n − 1 ) ! {\displaystyle \Gamma (n)=(n-1)!} ,可以被用作計算圓周率。
π = Γ ( 1 2 ) 2 {\displaystyle \pi =\Gamma ({\frac {1}{2}})^{2}}
π {\displaystyle \pi } 与圆密切相关,出现在许多几何学和三角学的公式 中(特别是与圆、椭圆和球体相关的那些)。 此外, π {\displaystyle \pi } 也出现在其他学科的重要公式中,比如统计学、物理学,傅立叶分析和数论的公式。
圆的面积等于 π {\displaystyle \pi } 乘以阴影部分面积。 π {\displaystyle \pi } 出现在基于圆的几何图形(如椭圆 、球 、圆锥 与环面 )的面积、体积公式中。下面是一些用到π的常见公式:[ 9]
半径 r {\displaystyle r} 的圆周长 2 π r {\displaystyle 2\pi r} 。 半径 r {\displaystyle r} 的圆面积 π r 2 {\displaystyle \pi r^{2}} 。 半径 r {\displaystyle r} 的球体积 4 3 π r 3 {\textstyle {\frac {4}{3}}\pi r^{3}} 。 半径 r {\displaystyle r} 的球面面积 4 π r 2 {\displaystyle 4\pi r^{2}} 。 上述公式是n 维球 的体积与其边界((n −1)维球的球面 )的表面积的特殊情况,具体将在后文 给出解释。
描述由圆生成的图形的周长、面积或体积的定积分 常涉及π。例如,表示半径为1的半圆的面积的积分为[ 105]