各种各样的数 基本 N ⊆ Z ⊆ Q ⊆ R ⊆ C {\displaystyle \mathbb {N} \subseteq \mathbb {Z} \subseteq \mathbb {Q} \subseteq \mathbb {R} \subseteq \mathbb {C} }
正數 R + {\displaystyle \mathbb {R} ^{+}} 自然数 N {\displaystyle \mathbb {N} } 正整數 Z + {\displaystyle \mathbb {Z} ^{+}} 小数 有限小数 无限小数 循环小数有理数 Q {\displaystyle \mathbb {Q} } 代數數 A {\displaystyle \mathbb {A} } 实数 R {\displaystyle \mathbb {R} } 複數 C {\displaystyle \mathbb {C} } 高斯整數 Z [ i ] {\displaystyle \mathbb {Z} [i]}
负数 R − {\displaystyle \mathbb {R} ^{-}} 整数 Z {\displaystyle \mathbb {Z} } 负整數 Z − {\displaystyle \mathbb {Z} ^{-}} 分數 單位分數 二进分数 規矩數 無理數 超越數 虚数 I {\displaystyle \mathbb {I} } 二次無理數 艾森斯坦整数 Z [ ω ] {\displaystyle \mathbb {Z} [\omega ]}
延伸 二元数 四元數 H {\displaystyle \mathbb {H} } 八元数 O {\displaystyle \mathbb {O} } 十六元數 S {\displaystyle \mathbb {S} } 超實數 ∗ R {\displaystyle ^{*}\mathbb {R} } 大實數 上超實數
雙曲複數 雙複數 複四元數共四元數 超复数 超數 超現實數
其他 質數 P {\displaystyle \mathbb {P} } 可計算數 基數 阿列夫數 同餘 整數數列 公稱值
規矩數 可定义数 序数 超限数p 進數 数学常数
圓周率 π = 3.14159265 {\displaystyle \pi =3.14159265} … 自然對數的底 e = 2.718281828 {\displaystyle e=2.718281828} … 虛數單位 i = − 1 {\displaystyle i={\sqrt {-{1}}}} 無限大 ∞ {\displaystyle \infty }
π {\displaystyle \pi } 是无理数,不能用分数表示出来(即它的小数部分是无限不循环小数),但近似22 7 {\textstyle {\frac {22}{7}}} 等有理数。學界認為π {\displaystyle \pi } 的数字序列在统计上是随机分布,但迄今未能证明。此外,π {\displaystyle \pi } 还是超越数,亦即它不是任何有理 系数多项式的根;化圆为方的问题不可能用尺规作图解决。
几个文明古国很早就須计算出π {\displaystyle \pi } 的精确值以便于生产的计算。西元5世纪,中國劉宋数学家祖冲之用几何方法将圆周率计算到小数点后7位。大约同时,印度 数学家也将圆周率计算到小数点后5位。史上首條π {\displaystyle \pi } 的精确无穷级数公式(即π的莱布尼茨公式)直到约1000年后才由印度数学家发现。微積分出現,π {\displaystyle \pi } 的位數很快計到數百位,足以滿足任何科學工程的計算需求。在20和21世纪,计算机技术 快速发展,π {\displaystyle \pi } 的计算精度急速提高。截至2024年3月,π {\displaystyle \pi } 的十进制精度已達105万亿位。几乎所有科学研究对π {\displaystyle \pi } 的精度要求都不超过几百位,当前计算π {\displaystyle \pi } 的值主要都为打破记录、测试超级计算机的计算能力和高精度乘法算法。:17
π {\displaystyle \pi } 的定义涉及圆,在三角学和几何学 的许多公式,特别是广泛应用在圆形、球形或椭球形相關公式中。在近代數學分析裡,π {\displaystyle \pi } 改由實數系統譜性質中的特征值或週期來定義,其他數學領域如數論、統計 以及幾乎所有物理學 領域均有出現,π {\displaystyle \pi } 的广泛用途使它成为科学界内外最广为人知的数学常数。几本专门介绍π {\displaystyle \pi } 的书籍经已出版,圆周率日(3月14日)和π {\displaystyle \pi } 值计算突破记录也往往会成为报纸的新闻头条。此外,背诵π {\displaystyle \pi } 值的世界记录已达10萬位。
基本概念
名称 数学家用小写希腊字母π {\displaystyle \pi } 表示圆周和其直径之比,有时也将其拼写为“Pi ”,来自希腊语“περίμετρος ”(周长)的首字母。英语π {\displaystyle \pi } 的发音与英文单词“Pie”( ,西式馅饼)相同。π的小写字母(或其无衬线体)在数学要和表示连乘积的大写Π相区分开。
关于选择符号π {\displaystyle \pi } 的原因,请参见引入π符号 一节。
定义 π {\displaystyle \pi } 常用定义为圆的周长C {\displaystyle C} 与直径 d {\displaystyle d} 的比值::8
π = C d {\displaystyle \pi ={\frac {C}{d}}} 。无论圆的大小如何,比值C d {\displaystyle {\frac {C}{d}}} 为恒值。如果圆的直径变为原先的二倍,周长也变为二倍,比值C d {\displaystyle {\frac {C}{d}}} 不变。π {\displaystyle \pi } 目前的定义暗地用了欧几里得几何 的一些定理,虽然圆的定义可扩展到任意曲面(即非欧几里得几何 ),但这些圆不符合定律π = C d {\displaystyle \pi ={\frac {C}{d}}} 。
这里,圆的周长指其圆周的弧长,弧长这概念可以不依赖几何学,而是用微积分学的极限来定义。例如,若想计算笛卡儿坐标系中单位圆x 2 + y 2 = 1 {\displaystyle x^{2}+y^{2}=1} 上半部分的弧长,需要用到积分:
π = ∫ − 1 1 d x 1 − x 2 . {\displaystyle \pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.} 上述积分是由卡尔·魏尔斯特拉斯于1841年对π {\displaystyle \pi } 的积分定义。
π {\displaystyle \pi } 这些依赖周长、且暗地依赖积分的定义如今在文献中并不常见。雷默特(Remmert (1991) )解释说现代教微积分時,大学一般将微分学课程安排在积分学课程之前,所以不依赖于后者的π的定义就很有必要了。其中一种定义由理查·巴爾策 提出,由愛德蒙·蘭道推广,其表述如下:π {\displaystyle \pi } 是两倍于能使余弦函数等于零的最小正数。余弦函数可以由独立于几何之外的幂级数定义,或者使用微分方程的解来定义。
在相似的启发下,π {\displaystyle \pi } 可以用关于复变量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 } 。
基于同样想法但更抽象的定义运用了精巧的拓扑学 和代数学概念,用以下定理描述:存在唯一的从加法模数整数组成的实数群R /Z 到绝对值为1的复数组成的乘法群的连续同态(拓扑学概念,指在拓扑空间之间的一种态射)。数字π定义为此同态派生的模的一半。
周长固定,圆会围成最大面积,π同樣表述为等周不等式中出现的常数(乘四分之一)。此外,在很多其他紧密相关的方程中,π作为某些几何或者物理过程的特征值出现;详见下文 。
无理及正规性 π {\displaystyle \pi } 是无理数,无法表示成两整数之比 的形式(形如22 7 {\textstyle {\frac {22}{7}}} 的分数常用来近似表达π {\displaystyle \pi } ,但是没有任何普通分数(指整数的比)可以取到π {\displaystyle \pi } 的精确值)。:5 由于π {\displaystyle \pi } 是无理数,故可表示为无限不循环小数。有多种方法能证明π是无理数,这些证明也都要用到微积分学和反证法。π {\displaystyle \pi } 可以用有理数 来近似的程度還無法準確得知(稱為無理性度量),不過估計其無理性度量比e {\displaystyle e} 或ln ( 2 ) {\displaystyle \ln(2)} 的要大,但是小於刘维尔数的無理性度量。
統計隨機性 检验,包括正规数检验,可验证π {\displaystyle \pi } 的位數沒有明顯的固定模式。π {\displaystyle \pi } 的小数中任意固定长度的序列(如3位數000,001……999)出現機率都相同。不過有關π {\displaystyle \pi } 是正规数的猜想既無證明,亦無证伪:22-23 。
電腦出現後可生成大量π {\displaystyle \pi } 的不同位数,并統計分析之。金田康正詳細統計分析了π {\displaystyle \pi } 的十進制數字,并验证了其分布正规:例如,假設檢定0到9十個數的出現頻率,找不到有特定重复规律的證據:22, 28–30 。根據無限猴子定理,任何任意長度、由隨機內容組成的子序列看起來都有可能像不隨機生成。因此,就算π {\displaystyle \pi } 的小数序列通過了隨機性統計測試,其中也可能有幾位的數字看起來似有规律可循而非隨機数,例如π {\displaystyle \pi } 的十進制写法在小數第762位后开始出现了連續六個9: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 } 的值)
π {\displaystyle \pi } 的超越性衍生出一些重要的结果:π {\displaystyle \pi } 不能经有限次四则运算和开平方运算有理数来获得,因此不是规矩数。换言之,尺规作图作不出长度为π {\displaystyle \pi } 的线段,也就不可能用尺规方法做出与已知圆面积相等的正方形。后者即为有名的化圓為方问题,该问题早在古典时代 即已提出,曾困扰人数千年之久。直至今天,依然有民间数学爱好者声称他们解决了这问题。
连分式 π {\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 } 最广为人知且广為使用的几个近似值。用以上方式得出的π {\displaystyle \pi } 的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近π {\displaystyle \pi } 。π {\displaystyle \pi } 是超越数,据定义来说它不是代數數,又因此不可能是二次無理數;是故π {\displaystyle \pi } 不能表示为循环连分数。尽管π {\displaystyle \pi } 的简单连分数没有表现出任何其他明显规律,数学家發現了数條广义连分数能表示π {\displaystyle \pi } ,例如:
π = 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 {\displaystyle {\frac {13}{4}}} 、16 5 {\displaystyle {\frac {16}{5}}} 、19 6 {\displaystyle {\frac {19}{6}}} 、22 7 {\displaystyle {\frac {22}{7}}} 、179 57 {\displaystyle {\frac {179}{57}}} 、267 85 {\displaystyle {\frac {267}{85}}} 、333 106 {\displaystyle {\frac {333}{106}}} 、355 113 {\displaystyle {\frac {355}{113}}} 、52163 16604 {\displaystyle {\frac {52163}{16604}}} 、53228 16943 {\displaystyle {\frac {53228}{16943}}} 、55358 17621 {\displaystyle {\frac {55358}{17621}}} 、57843 18412 {\displaystyle {\frac {57843}{18412}}} 、60328 19203 {\displaystyle {\frac {60328}{19203}}} 、103993 33102 {\displaystyle {\frac {103993}{33102}}} 、245850922 78256779 {\displaystyle {\frac {245850922}{78256779}}} (选自 A063674 及 A063673。) 小數 (整数后首80位):3.1415926535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899... :240 (另见 A000796)其他进位制的近似值
二进制(整数后首48位):11.00100100 0011 1111 0110 1010 1000 1000 1000 0101 1010 0011… 十六进制(整数后首20位):3.243F6A88 85A3 08D3 1319… :242 六十进制(整数后首20位):3;8,29,44,0,47,25,53,7,24,57,36,17,43,4,29,7,10,3,41,17…
复数与欧拉恒等式 任何复数(以z {\displaystyle z} 为例)都可以表示为一组实数对:极坐标系用实数r {\displaystyle r} 表示半径,代表复平面上复数z {\displaystyle z} 离原點的距离;实数φ {\displaystyle \varphi } 则表示夹角,即这条半径(复平面上复数z {\displaystyle z} 与原点的连线)与正实轴经顺时针转动的夹角。这样一来,z {\displaystyle z} 就可写成
z = r ⋅ ( cos φ + i sin φ ) {\displaystyle z=r\cdot (\cos \varphi +i\sin \varphi )} ,这里i {\displaystyle i} 代表虛數單位,即i 2 = − 1 {\displaystyle i^{2}=-1} 。复分析中,欧拉公式将三角函数与复指数函数糅合在一起:
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 ),全因它将五个最基本的数学常数简洁联系起来。欧拉等式亦可用于求出方程z n = 1 {\displaystyle z^{n}=1} 的n {\displaystyle n} 个不同复数根(这些根叫做n {\displaystyle n} 次单位根”),可以根据以下公式求得:
e 2 π i k n ( k = 0 , 1 , 2 , … , n − 1 ) {\displaystyle e^{\frac {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 } 。
π {\displaystyle \pi } 是上述方程的最小特征值,也和弦振動的基本模式 有關。一種讓弦振動的方式是提供弦能量,能量會滿足維廷格函數不等式,其中提到若函數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} ,故等號成立的條件是曲线為圓形。
圓周率π {\displaystyle \pi } 也和庞加莱不等式的最佳常數有關,π {\displaystyle \pi } 是一維及二維的狄氏能量 特征向量最佳值中最小,會出現在許多經典的物理現象中,例如經典的位势论。其一維的情形即為維廷格不等式。
圓周率π {\displaystyle \pi } 也是傅里叶变换的重要常數,傅里叶变换屬於积分变换,將實數線上有複數值、可積分的函數,轉換為以下形式:
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 }} 空間的代數同態。
不确定性原理也用到π {\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}} 。物理的結果,有關量子力学中同時觀測位置及動量的不確定性,見下文 。傅立葉分析中出現π {\displaystyle \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 } 的正态分布的機率密度函數:
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 }}} 。概率论与统计学领域经常使用正态分布来作为复杂现象的简单模型:例如科学家通常假设大多数试验观测值的随机误差都是服从正态分布。
概率论与统计学中的中心极限定理解释了正态分布以及π {\displaystyle \pi } 的核心作用,这定理本质上是联系着π {\displaystyle \pi } 的谱特征 与海森堡不确定性原理相关的特征值,并且在不确定性原理中有
σ x σ p ≥ ℏ 2 {\displaystyle \sigma _{x}\sigma _{p}\geq {\frac {\hbar }{2}}} ,这里的σ x {\displaystyle \sigma _{x}} 與σ p {\displaystyle \sigma _{p}} 分別為位置與動量的標準差,ℏ {\displaystyle \hbar } 是約化普朗克常数 ,而不等式的等号当且仅当粒子的波函数为高斯函数使成立。
同样地,π {\displaystyle \pi } 作为唯一独特的常数使得高斯函数等于其自身的傅里叶变换,此时的高斯函数形式为f ( x ) = e − π x 2 {\displaystyle f(x)=e^{-\pi x^{2}}} 。根据豪(Howe )的说法,建立傅里叶分析基本定理的“全部工作(whole business)”简化为高斯积分。
历史
远古时期 圓周率在远古时期(西元前一千纪 )已估算至前两位(3.1)。有些埃及學家聲稱,遠至古王國時期時期的古埃及人已經用22 7 {\textstyle {\frac {22}{7}}} 作為圓周率的約數,但這說法受到質疑。
最早有記載的对圓周率估值在古埃及 和巴比伦出现,兩估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一塊西元前1900至1600年的泥板,泥板上的幾何學陳述暗示人们当时把圓周率視同25 8 {\textstyle {\frac {25}{8}}} (等於3.125)。:167 埃及的莱因德数学纸草书(鉴定撰寫年份為西元前1650年,但抄自一份西元前1850年的文本)載有用作計算圓面積的公式,该公式中圓周率等于( 16 9 ) 2 {\textstyle ({\frac {16}{9}})^{2}} (≈3.1605)。:167
西元前4世紀的《百道梵書 》的天文學運算把339 108 {\textstyle {\frac {339}{108}}} (≈3.139,精确到99.91%)用作圓周率估值。西元前150年前其他印度文獻把圓周率視為10 {\textstyle {\sqrt {10}}} (≈3.1622):169 。
割圆时代 第一條有紀錄、嚴謹計算π {\displaystyle \pi } 數值的演算法是用正多邊形的幾何算法,在西元前250年由希臘數學家阿基米德發明。:170 這算法用了有一千年之久,因而有時π {\displaystyle \pi } 亦稱阿基米德常數。: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} )。阿基米德得到的上限22 7 {\textstyle {\frac {22}{7}}} 也造成常見誤解,認為π {\displaystyle \pi } 就等於22 7 {\textstyle {\frac {22}{7}}} :171 。在西元前150年,希臘羅馬的科學家克劳狄乌斯·托勒密在《天文学大成》一書中提到π {\displaystyle \pi } 的數值是3.1416,可能來自阿基米德,也可能來自阿波罗尼奥斯。:176 數學家在1630年利用多邊形的方式計算π {\displaystyle \pi } 到第39位小數,一直到1699年,其他數學家才利用無窮級數的方式打破其紀錄,計算到第71位小數。
中国历史 上,π {\displaystyle \pi } 的數值有3、3.1547(西元前一世紀)、10 {\displaystyle {\sqrt {10}}} (西元前100年,數值約3.1623)及142 45 {\textstyle {\frac {142}{45}}} (第三世紀,數值約3.1556):176–177 。大約在西元265年,曹魏 數學家刘徽創立割圆术,用正3072邊形計算出π 的數值為3.1416。:177 他後來又發明了較快的算法,利用邊數差兩倍的正多邊形,其面積的差值會形成等比數列,其公比為1 4 {\textstyle {\frac {1}{4}}} 的原理,配合96邊形算出π {\displaystyle \pi } 的值為3.14。祖冲之在西元480年利用割圆术計算12288邊形邊長,得到π {\displaystyle \pi } 的值在3.1415926和3.1415927之间。他同时提出了π的约率22 7 {\textstyle {\frac {22}{7}}} 和密率355 113 {\textstyle {\frac {355}{113}}} 。在之後的八百年內,這都是π 最準確的估計值。:178 為紀念祖沖之對圓周率發展的貢獻,日本 數學家三上義夫將這推算值命名為“祖沖之圓周率”,簡稱“祖率”。
印度天文學家阿耶波多在西元499年的著作《阿里亞哈塔曆書》中使用了3.1416的數值。:179 斐波那契在大約1220年用獨立於阿基米德多邊形法,計算出3.1418:180 。義大利作家但丁·阿利吉耶里用的數值則是3 + 2 10 ≈ 3.14142 {\textstyle 3+{\frac {\sqrt {2}}{10}}\approx 3.14142} 。:180
波斯天文學家卡西在1424年利用3×228 邊的多邊形,計算到六十進制的第9位小數,相當十進制的第16位小數。這一突破成為當時的紀錄,延續了約180年。法國數學家弗朗索瓦·韦达在1579年用3×217 邊形計算到第9位小數,佛蘭芒數學家阿德里安·范·羅門在1593年計算到第15位小數。荷蘭數學家鲁道夫·范·科伊伦在1596年計算到第20位小數,他之後又計算到第35位小數(因此在二十世紀初之前,圓周率在德國會稱為鲁道夫數)。:182–183 荷蘭科學家威理博·司乃耳在1621年計算到第34位小數:183 ,而奧地利天文學家克里斯托夫·格林伯格 在1630年用1040 邊形計算到第38位小數,至今這仍是利用多邊形算法可以達到最準確的結果:183 。
无穷级数 16及17世紀時,開始改用無窮级数的方式去計π {\displaystyle \pi } 。無窮级数是一組無窮數列的和:185–191 。無窮级数讓數學家可以計算出比阿基米德以及其他用幾何方式計算的數學家更準確的結果。:185–191 雖然詹姆斯·格雷果里及戈特弗里德·莱布尼茨等歐洲數學家利用無窮數列計算π 而使得该方法为大家所知,但这种方法最早是由印度 科學家在大約1400到1500年之間發現。:185-186 第一個记载用無窮级数計算π {\displaystyle \pi } 的人是约西元1500年左右时,印度天文學家尼拉卡莎·薩默亞士 在他的著作《系統匯編 》中用梵語詩所記錄。當時沒有這數列對應的證明,而證明出現在另一本較晚的印度作品《基本原理 》,年代約在西元1530年。尼拉卡莎將該數列歸功於更早期的印度數學家桑加馬格拉馬的馬德哈瓦 (1350–1425)。相關的無窮级数有許多,包括有關sin {\displaystyle \sin } 、tan {\displaystyle \tan } 及cos {\displaystyle \cos } 的,現在稱為馬德哈瓦數列 或π的莱布尼茨公式。瑪達瓦在1400年用無窮级数計算π {\displaystyle \pi } 到第11位小數,但在1430年一位波斯數學家卡西利用多邊形算法否定了他算的結果。
歐洲發現的第一條無窮項圓周率公式是無窮乘積(和一般用來計算π {\displaystyle \pi } 的無窮級數不同),由法國科學家弗朗索瓦·韦达在1593年發現:187 :
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年發現了沃利斯乘积,是歐洲發現的第二條無窮項圓周率公式: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年代發明,許多計π {\displaystyle \pi } 的無窮級數出現。牛頓自己就用反正弦(arcsin {\displaystyle \arcsin } )數列在1655年或1666年將π {\displaystyle \pi } 近似到第15位小數,後來寫到「我很羞愧告訴你我為了計算它用了多少數字,我當時沒有做其他事。」
蘇格蘭數學家詹姆斯·格雷果里在1671年發現了馬德哈瓦公式,莱布尼茨也在1674年發現::188–189
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}}} 。1699年時英國數學家亚伯拉罕·夏普用格雷果里-莱布尼茨公式,在z = 1 3 {\displaystyle z={\frac {1}{\sqrt {3}}}} 時計算,計算到π的第71位小數,打破由多邊形算法得到的第39位小數的记录。:189 格雷果里-莱布尼茨公式在z = 1 {\displaystyle z=1} 時非常簡單,但收斂到最終值的速度非常慢,現在不会再用此公式來計π。:156
約翰·梅欽在1706年用格雷果里-莱布尼茨級數產生了可以快速收斂的公式::192–193
π 4 = 4 arctan 1 5 − arctan 1 239 {\displaystyle {\frac {\pi }{4}}=4\,\arctan {\frac {1}{5}}-\arctan {\frac {1}{239}}} 梅欽用這公式計到π {\displaystyle \pi } 第100位小數:72–74 後來其他數學家也發展了一些類似公式,現在稱為梅欽類公式,創下了許多計算π {\displaystyle \pi } 位數的紀錄。:72–74 在進入電腦時代時,梅欽類公式仍然是耳熟能详可以計算π {\displaystyle \pi } 的公式,而且在约250年的时间里,很多有關π {\displaystyle \pi } 位數的紀錄都是梅欽類公式所得,比如在1946年時由達尼爾·弗格森(Daniel Ferguson )用這類公式計到第620位小數,是沒有計算設備輔助的最佳紀錄。:192–196, 205
1844年,計算天才扎卡里亞斯·達斯 在德國數學家卡爾·弗里德里希·高斯的要求下以梅欽類公式心算了π {\displaystyle \pi } 的200位小數,並創下紀錄。:194-196 英國數學家威廉·謝克斯 花了15年的時間計算π {\displaystyle \pi } 到小數707位,不過第528位小數出錯,後面的小數也都不正確。:194–196
收敛速度 有些π {\displaystyle \pi } 的無窮級數收斂的比其他級數要快,數學家一般會選用收斂速度較快的級數,可以在較少的計算量下計算π {\displaystyle \pi } ,且達到需要的準確度:15–17, 70–72, 104, 156, 192–197, 201–202 。以下是π {\displaystyle \pi } 的莱布尼茨公式::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萬項之後,才會精確到π {\displaystyle \pi } 的第五位小數。
尼拉卡莎在15世紀發展了π {\displaystyle \pi } 的另一條無窮級數,收斂速度比格雷果里-萊布尼茨公式快很多:
π = 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 } 以下比較兩條級數的收斂速率:
π {\displaystyle \pi } 的無窮級數 第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…
計算前五項後,格雷果里-萊布尼茨級數的和跟π {\displaystyle \pi } 的誤差為0.2,而尼拉卡莎級數和的誤差為0.002。尼拉卡莎級數收斂快很多,也甚為適合用來計π {\displaystyle \pi } 的值。收斂更快的級數有梅欽類公式及楚德诺夫斯基算法,後者每計一項就可以得到14位正確的小數位。
无理与超越性 并非所有和π {\displaystyle \pi } 有关的研究都旨在提高计算它的准确度。1735年,欧拉解决了巴塞尔问题,建立了所有平方数倒数和与π {\displaystyle \pi } 的关系。之后欧拉发现了欧拉乘积公式,得到了π {\displaystyle \pi } 、素数的重要關聯,對日後黎曼ζ函數的研究影響深遠。
π 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年,瑞士数学家约翰·海因里希·朗伯用正切函数的无穷连分数表达式证明了π {\displaystyle \pi } 是無理數。:5 1794年,法国数学家阿德里安-马里·勒让德证明了π 2 {\displaystyle \pi ^{2}} 也是无理数。1882年,德国数学家费迪南德·冯·林德曼证明了对任何非零代数数α {\displaystyle \alpha } ,e α {\displaystyle e^{\alpha }} 都是超越数,该结论后来由魏尔斯特拉斯推广为林德曼-魏尔斯特拉斯定理。据此定理和欧拉公式,π只能是超越數,進而证实了勒让德和欧拉提出的π超越性猜想。:196 哈代在其著作《数论导引》中则称此证明在提出後,經過希尔伯特、施瓦兹和其他一些人化简过。
引入π {\displaystyle \pi } 符号 在用π {\displaystyle \pi } 专指“圆周率”之前,希腊字母即已用於幾何概念中:166 。威廉·奥特雷德在1647年起在《數學之鑰》(Clavis Mathematicae)就已經用π {\displaystyle \pi } 及δ {\displaystyle \delta } (對應p和d的希臘字母)來表示圓的周長及直徑的比例。
威廉·琼斯在他1706年出版的《新數學導論》(A New Introduction to the Mathematics )提到了π {\displaystyle \pi } ,是目前已知最早专门用希臘字母π {\displaystyle \pi } 表示圓周和其直徑比例的人。這希臘字母第一次出现是在书中討論一塊半徑1的圓時提到「其圓周長一半(π {\displaystyle \pi } )」。琼斯選用π {\displaystyle \pi } 可能因它是希臘文“周边”一词“περιφέρεια ”的首字母。不過琼斯提到,他那些有關π {\displaystyle \pi } 的算式出自「真正聰明的約翰·梅欽先生」,人们推測在瓊斯之前,約翰·梅欽就已开始用π {\displaystyle \pi } 表示圓周率:166 。
瓊斯在1706年開始使用此希臘字母,但直到萊昂哈德·歐拉在其1736年出版的《力學 》中開始使用之后,其他数学家才纷纷开始用π {\displaystyle \pi } 指代圆周率。在此之前,數字家可能用像c 或p 之類的字母代表圓周率:166 。歐拉與歐洲其他數學家间时常互相写信来往,π {\displaystyle \pi } 的用法迅速傳播开来:166 。1748年歐拉在他的《无穷小分析引论》再一次提到了π {\displaystyle \pi } ,写道:「簡潔起見,我們將此數字寫為π {\displaystyle \pi } ,π {\displaystyle \pi } 等於半徑為1的圓周長的一半。」这表示方式之後也推展到整片西方世界 :166 。
现代数值近似
计算机时代与迭代算法 高斯-勒让德算法:
一開始設定
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位:205 。同年,喬治·韋斯納(George Reitwiesner)及约翰·冯·诺伊曼帶領的團隊利用反三角函数(arctan)的無窮級數,用ENIAC計算到了小數後2037位,花了70小時的電腦工作時間。這紀錄後來多次由其他透過arctan級數计算出的結果打破(1957年到7480位小數,1958年到第一萬位數,1961年到第十萬位小數),直到1973年,小数点后第一百萬位小數經已算出:197 。
1980年代有两项發明加速計算了π。第一项是發现了新的迭代法去计π 的值,計算速度比無窮級數快很多;另一项是發现了可以快速計算大數字乘積的乘法演算法:15–17 。電腦大部分的工作時間都是在計乘法,這類演算法對現代計π 格外重要:131 。這類演算法包括嘉良對馬(Karatsuba)算法、譚曲(Toom-Cook)乘法及以傅里叶变换為基礎的乘法演算法(傅里叶乘法):132, 140 。
迭代演算法最早是在1975年至1976年间分别由美國物理學家尤金·薩拉明 及奧地利科學家理查·布蘭特 独立提出:87 。這两條演算法没有依赖無窮級數來計算。迭代會重覆特定計算,将前一次的計算結果作为這一次的輸入值,使得計算結果漸漸的趨近理想值。此方式的原始版本其實是在160年前由卡爾·弗里德里希·高斯提出,現在稱為算术-几何平均数算法(AGM法)或高斯-勒让德算法:87 。薩拉明及布蘭特都曾修改之,这算法也稱為薩拉明-布蘭特演算法。
迭代演算法收斂速度比無窮級數快很多,在1980年代以後廣為使用。無窮級數隨著項次的增加,一般來說正確的位數也會增加幾位,但迭代演算法每計算多一次,正確位數會呈几何级数增长。例如薩拉明-布蘭特演算法每計算多一次,正確位數會是之前的二倍。1984年加拿大人喬納森·波温 及彼得·波温 提出迭代演算法,每計算多一次,正確位數會是之前的四倍,1987年時有另一條迭代演算法,每計算多一次,正確位數會是之前的五倍。日本數學家金田康正使用的演算法在1955年及2002年間創下了若干項紀錄。不過迭代演算法的快速收斂也有其代價,需要的記憶體明顯比無窮級數多。
计算π {\displaystyle \pi } 的意义 一般而言,π {\displaystyle \pi } 值并不需要过于精确便能够满足大部分数学运算的需求。按照約·安(Jörg Arndt)及古里斯佗夫·希奴(Christoph Haenel )的计算,39位精確度已可将可觀測宇宙 圆周的精确度準確至一粒原子大小,足以運算絕大多數宇宙学 的计算需求。尽管如此,和π有關的成就往往成為世界各地的新聞頭條;部分人出于對破紀錄的冲动,依然奋力算出π小数点后上千甚至上百萬位:17–19 。此外也有測試超级计算机、測試数值分析算法(包括高精度乘法算法 )等實際好處。純粹數學這领域也能计算π的位数评定其隨機度:18 。
快速收敛级数 现代计算π {\displaystyle \pi } 的程序不仅局限于迭代算法。20世纪80与90年代,出现了可用来计算π {\displaystyle \pi } 的新无穷级数,其收敛速度可与迭代算法媲美,而又有着复杂度、内存密集度更低的优势。印度数学家斯里尼瓦瑟·拉马努金是这方面的先驱,他在1914年发表了许多与π相关的公式,这些公式十分新颖,极为优雅而又颇具数学深度,收敛速度也非常快。: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})}}.} 这无穷级数收敛速度远快于绝大多数反正切数列,包括梅钦公式。:104 第一位使用拉马努金公式计算π并取得进展的是比尔·高斯珀 ,他在1985年算得了小数点后一千七百万位。:104, 206 拉马努金公式开创了现代数值近似算法的先河,此后波尔文兄弟和楚德诺夫斯基兄弟 进一步发展了这类算法。: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位数值,因而用於突破圆周率的数位的计算。利用这公式,楚德诺夫斯基兄弟于1989年算得π小数点后10亿(109 )位,法布里斯·贝拉于2009年算得2.7千亿(2.7×1012 )位,亚历山大·易和近藤滋在2011年算得一万亿(1013 )位。:110–111, 206 类似的公式还有拉马努金-佐藤级数 。
2006年,加拿大数学家西蒙·普勞夫利用PSLQ整数关系算法 按照以下模版生成了几條计算π的新公式:
π 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} 是普勞夫计算出的有理常数。
統計模擬法