Loading AI tools
數學常數 来自维基百科,自由的百科全书
圓周率是數學常數,為圓的周長和其直徑的比,近似值約3.14159265,常用符號表示。
各式各樣的數 |
基本 |
延伸 |
其他 |
是無理數,不能用分數表示出來(即它的小數部分是無限不循環小數),但近似等有理數。學界認為π的數字序列在統計上是隨機分布,但迄今未能證明。此外,π還是超越數,它不是任何有理係數多項式的根,化圓為方的問題不可能用尺規作圖解決。
幾個文明古國很早就須計算出π的精確值以便於生產的計算。西元5世紀,中國劉宋數學家祖沖之用幾何方法將圓周率計算到小數點後7位。大約同時,印度數學家也將圓周率計算到小數點後5位。史上首條π的精確無窮級數公式(即π的萊布尼茨公式)直到約1000年後才由印度數學家發現。[5][6]微積分出現,π的位數很快計到數百位,足以滿足任何科學工程的計算需求。在20和21世紀,計算機技術快速發展,π的計算精度急速提高。截至2024年3月,π的十進制精度已達105萬億位。[7]幾乎所有科學研究對π的精度要求都不超過幾百位,當前計算π的值主要都為打破記錄、測試超級計算機的計算能力和高精度乘法算法。[2]:17[8]
π的定義涉及圓,在三角學和幾何學的許多公式,特別是廣泛應用在圓形、球形或橢球形相關公式中。[9]在近代數學分析裡,π改由實數系統譜性質中的特徵值或週期來定義,其他數學領域如數論、統計以及幾乎所有物理學領域均有出現,π的廣泛用途使它成為科學界內外最廣為人知的數學常數。幾本專門介紹π的書籍經已出版,圓周率日(3月14日)和π值計算突破記錄也往往會成為報紙的新聞頭條。[10]此外,背誦π值的世界記錄已達10萬位。[11]
數學家用小寫希臘字母表示圓周和其直徑之比,有時也將其拼寫為「Pi」,來自希臘語「περίμετρος」(周長)的首字母。[12]英語π的發音與英文單詞「Pie」(/paɪ/,西式餡餅)相同。[13]π的小寫字母(或其無襯線體)在數學要和表示連乘積的大寫Π相區分開。
關於選擇符號π的原因,請參見引入π符號一節。
無論圓的大小如何,比值為恆值。如果圓的直徑變為原先的二倍,周長也變為二倍,比值不變。π目前的定義暗地用了歐幾里得幾何的一些定理,雖然圓的定義可擴展到任意曲面(即非歐幾里得幾何),但這些圓不符合定律。[2]
這裡,圓的周長指其圓周的弧長,弧長這概念可以不依賴幾何學,而是用微積分學的極限來定義。[14]例如,若想計算笛卡兒坐標系中單位圓上半部分的弧長,需要用到積分:[15]
上述積分是由卡爾·魏爾斯特拉斯於1841年對π的積分定義。[16]
π這些依賴周長、且暗地依賴積分的定義如今在文獻中並不常見。雷默特(Remmert (1991))解釋說現代教微積分時,大學一般將微分學課程安排在積分學課程之前,所以不依賴於後者的π的定義就很有必要了。其中一種定義由理查·巴爾策提出,[17]由愛德蒙·蘭道推廣,[18]其表述如下:π是兩倍於能使餘弦函數等於零的最小正數。[2][15][19]餘弦函數可以由獨立於幾何之外的冪級數[20]定義,或者使用微分方程式的解來定義。[19]
在相似的啟發下,π可以用關於複變數的複指數函數來定義。複指數類似餘弦函數,可用多種方式定義。令函數值為一的複數集合是如下所示的(虛)等差數列:
基於同樣想法但更抽象的定義運用了精巧的拓撲學和代數學概念,用以下定理描述:[22]存在唯一的從加法模數整數組成的實數群R/Z到絕對值為1的複數組成的乘法群的連續同態(拓撲學概念,指在拓撲空間之間的一種態射)。數字π定義為此同態派生的模的一半。[23]
周長固定,圓會圍成最大面積,π同樣表述為等周不等式中出現的常數(乘四分之一)。此外,在很多其他緊密相關的方程式中,π作為某些幾何或者物理過程的特徵值出現;詳見下文。
π是無理數,無法表示成兩整數之比的形式(形如的分數常用來近似表達π,但是沒有任何普通分數(指整數的比)可以取到π的精確值)。[2]:5由於是無理數,故可表示為無限不循環小數。有多種方法能證明π是無理數,這些證明也都要用到微積分學和反證法。可以用有理數來近似的程度還無法準確得知(稱為無理性度量),不過估計其無理性度量比e或ln(2)的要大,但是小於萊歐維爾數的無理性度量[24]。
統計隨機性檢定,包括正規數檢定,可驗證的位數沒有明顯的固定模式。的小數中任意固定長度的序列(如3位數000,001……999)出現機率都相同[25]。不過有關π是正規數的猜想既無證明,亦無否證[2]:22-23[25]。
電腦出現後可生成大量π的不同位數,並統計分析之。金田康正詳細統計分析了π的十進制數字,並驗證了其分布正規:例如,假說檢定0到9十個數的出現頻率,找不到有特定重複規律的證據[2]:22, 28–30。根據無限猴子定理,任何任意長度、由隨機內容組成的子序列看起來都有可能像不隨機生成。因此,就算π的小數序列通過了隨機性統計測試,其中也可能有幾位的數字看起來似有規律可循而非隨機數,例如π的十進制寫法在小數第762位後開始出現了連續六個9[2]:3。
不僅是無理數,還是超越數,即不是任何有理係數多項式的根。(比方說,試圖解有限項方程式來求的值)[26][註 1]
的超越性衍生出一些重要的結果:不能經有限次四則運算和開平方運算有理數來獲得,因此不是規矩數。換言之,尺規作圖作不出長度為的線段,也就不可能用尺規方法做出與已知圓面積相等的正方形。後者即為有名的化圓為方問題,該問題早在古典時代即已提出,曾困擾人數千年之久[27][28]。直至今天,依然有民間數學愛好者聲稱他們解決了這問題[29]。
像所有無理數一樣無法表示成分數,但等全部無理數都能表示成一系列叫連分數的連續分數形式:
在這連分數的任意一點截斷化簡,都能得到π的近似值;前四位近似值是3、、、。這些數在歷史上是π最廣為人知且廣為使用的幾個近似值。用以上方式得出的的近似值要比任何有相同或更小的整數分母的其他整數分數近似值更接近π。[30]π是超越數,據定義來說它不是代數數,又因此不可能是二次無理數;是故π不能表示為循環連分數。儘管的簡單連分數沒有表現出任何其他明顯規律,[31]數學家發現了數條廣義連分數能表示π,例如:[32]
圓周率近似值包括:
其他進位制的近似值
任何複數(以為例)都可以表示為一組實數對:極坐標系用實數表示半徑,代表複數平面上複數離原點的距離;實數則表示夾角,即這條半徑(複數平面上複數與原點的連線)與正實數軸經順時針轉動的夾角。這樣一來,就可寫成[33]
複分析中,歐拉公式將三角函數與復指數函數糅合在一起[34]:
歐拉公式確立了的複指數與複數平面上以原點為圓心的單位圓上的點之間的關係,而且當時,歐拉公式就能改寫為歐拉恆等式的形式:
歐拉等式亦可用於求出方程式的個不同複數根(這些根叫做次單位根」[36]),可以根據以下公式求得:
常出現在有關幾何的問題中。然而,不少和幾何無關的問題也可看到的身影。
在許多用處中都會以特徵值形式出現。例如理想的振動弦問題可以建模為函數在單位區間的圖形,固定邊界值為。弦振動的模態會是微分方程式的,此處λ是相關的特徵值。受史特姆-萊歐維爾理論限制,只能是一些特定的數值。而即為一個特徵值,因為函數滿足邊界條件及微分方程式[37]。
是上述方程式的最小特徵值,也和弦振動的基本模式有關。一種讓弦振動的方式是提供弦能量,能量會滿足維廷格函數不等式[38],其中提到若函數使得,且和都是平方可積函數,則以下的不等式成立:
此例中等號成立的條件恰好是為倍數的時候。因此似乎是維爾丁格不等式的最佳常數,也是最小的特徵值(根據雷利商數的計算方式)
在更高維度的分析也有類似的角色,出現在其他類似問題的特徵值中。就如以上所述,的一項特點是等周定理中的最佳常數:周長為的平面若爾當曲線,所圍面積滿足以下的不等式
及,故等號成立的條件是曲線為圓形[39]。
圓周率π也和龐加萊不等式的最佳常數有關[40],是一維及二維的狄氏能量特徵向量最佳值中最小,會出現在許多經典的物理現象中,例如經典的位勢論[41][42][43]。其一維的情形即為維廷格不等式。
圓周率π也是傅立葉轉換的重要常數,傅立葉轉換屬於積分轉換,將實數線上有複數值、可積分的函數,轉換為以下形式:
傅立葉轉換有幾種不同的寫法,但不論怎麼寫,傅立葉轉換及反傅立葉轉換中,一定會有某處出現。不過上述的定義是最經典的,因為其描述了L2空間中唯一的么正算符,也是空間到空間的代數同態[44]。
不確定性原理也用到。不確定性原理提出了可以將函數在空間及在頻域中局部化程度的下限,用傅立葉轉換的方式表示:
物理的結果,有關量子力學中同時觀測位置及動量的不確定性,見下文。傅立葉分析中出現π是史東-凡紐曼定理的結果,證實了海森伯群的薛定諤表示是唯一[45]。
高斯積分是對高斯函數在整條實數軸上的積分,即函數下方與X軸圍成的面積,其結果為,
此積分的計算可以先計算對整條實數軸的積分的平方,通過轉換笛卡兒坐標系為極坐標系從而求得
其他計算方法可參閱高斯積分。高斯函數更一般的形式為,求一般形式的高斯積分均可通過換元積分法轉化為求的積分。
另外,當高斯函數為以下形式時,它則是平均數為和標準差為的常態分布的機率密度函數[46]:
這函數是機率密度函數,函數下方與X軸圍成的面積必須為1,令和即可轉換得出。機率論與統計學領域經常使用常態分布來作為複雜現象的簡單模型:例如科學家通常假設大多數試驗觀測值的隨機誤差都是服從常態分布[47]。
機率論與統計學中的中央極限定理解釋了常態分布以及的核心作用,這定理本質上是聯繫著的譜特徵與海森堡不確定性原理相關的特徵值,並且在不確定性原理中有
這裡的與分別為位置與動量的標準差,是約化普朗克常數,而不等式的等號若且唯若粒子的波函數為高斯函數使成立[48]。
同樣地,作為唯一獨特的常數使得高斯函數等於其自身的傅立葉轉換,此時的高斯函數形式為[49]。根據豪(Howe)的說法,建立傅立葉分析基本定理的「全部工作(whole business)」簡化為高斯積分。
圓周率在遠古時期(西元前一千紀)已估算至前兩位(3.1)。有些埃及學家聲稱,遠至古王國時期時期的古埃及人已經用作為圓周率的約數[50][註 2],但這說法受到質疑。[52][53][54][55]
最早有記載的對圓周率估值在古埃及和巴比倫出現,兩估值都與圓周率的正確數值相差不到百分之一。巴比倫曾出土一塊西元前1900至1600年的泥板,泥板上的幾何學陳述暗示人們當時把圓周率視同(等於3.125)。[2]:167埃及的萊因德數學紙草書(鑑定撰寫年份為西元前1650年,但抄自一份西元前1850年的文本)載有用作計算圓面積的公式,該公式中圓周率等於(≈3.1605)。[2]:167
西元前4世紀的《百道梵書》的天文學運算把(≈3.139,精確到99.91%)用作圓周率估值[56]。西元前150年前其他印度文獻把圓周率視為(≈3.1622)[2]:169。
第一條有紀錄、嚴謹計算π數值的演算法是用正多邊形的幾何算法,在西元前250年由希臘數學家阿基米德發明。[2]:170這算法用了有一千年之久,因而有時π亦稱阿基米德常數。[2]:175、205阿基米德的算法是在計算圓的外切正六邊形及內接正六邊形的邊長,以此計算的上限及下限,之後再將六邊形變成十二邊形,繼續計算邊長,一直計到正96邊形為止。他根據多邊形的邊長證明(也就是)[57]。阿基米德得到的上限也造成常見誤解,認為就等於[2]:171。在西元前150年,希臘羅馬的科學家克勞狄烏斯·托勒密在《天文學大成》一書中提到π的數值是3.1416,可能來自阿基米德,也可能來自阿波羅尼奧斯。[2]:176[58]數學家在1630年利用多邊形的方式計算π到第39位小數,一直到1699年,其他數學家才利用無窮級數的方式打破其紀錄,計算到第71位小數[59]。
中國歷史上,的數值有3[60]、3.1547(西元前一世紀)、(西元前100年,數值約3.1623)及(第三世紀,數值約3.1556)[2]:176–177。大約在西元265年,曹魏數學家劉徽創立割圓術,用正3072邊形計算出π的數值為3.1416。[61][2]:177他後來又發明了較快的算法,利用邊數差兩倍的正多邊形,其面積的差值會形成等比數列,其公比為的原理,配合96邊形算出π的值為3.14。[61]祖沖之在西元480年利用割圓術計算12288邊形邊長,得到π的值在3.1415926和3.1415927之間。他同時提出了π的約率和密率。在之後的八百年內,這都是π最準確的估計值。[2]:178為紀念祖沖之對圓周率發展的貢獻,日本數學家三上義夫將這推算值命名為「祖沖之圓周率」,簡稱「祖率」。[62]
印度天文學家阿耶波多在西元499年的著作《阿里亞哈塔曆書》中使用了3.1416的數值。[2]:179斐波那契在大約1220年用獨立於阿基米德多邊形法,計算出3.1418[2]:180。義大利作家但丁·阿利吉耶里用的數值則是。[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。
16及17世紀時,開始改用無窮級數的方式去計π。無窮級數是一組無窮數列的和[2]:185–191。無窮級數讓數學家可以計算出比阿基米德以及其他用幾何方式計算的數學家更準確的結果。[2]:185–191雖然詹姆斯·格雷果里及戈特弗里德·萊布尼茨等歐洲數學家利用無窮數列計算π而使得該方法為大家所知,但這種方法最早是由印度科學家在大約1400到1500年之間發現。[2]:185-186[67]第一個記載用無窮級數計算π的人是約西元1500年左右時,印度天文學家尼拉卡莎·薩默亞士在他的著作《系統匯編》中用梵語詩所記錄。[68]當時沒有這數列對應的證明,而證明出現在另一本較晚的印度作品《基本原理》,年代約在西元1530年。尼拉卡莎將該數列歸功於更早期的印度數學家桑加馬格拉馬的馬德哈瓦(1350–1425)。[68]相關的無窮級數有許多,包括有關、及的,現在稱為馬德哈瓦數列或π的萊布尼茨公式[68]。瑪達瓦在1400年用無窮級數計算π到第11位小數,但在1430年一位波斯數學家卡西利用多邊形算法否定了他算的結果[69]。
歐洲發現的第一條無窮項圓周率公式是無窮乘積(和一般用來計算π的無窮級數不同),由法國科學家弗朗索瓦·韋達在1593年發現[2]:187[71]:
約翰·沃利斯在1655年發現了沃利斯乘積,是歐洲發現的第二條無窮項圓周率公式[2]:187:
微積分學由英國科學家艾薩克·牛頓及德國數學家戈特弗里德·萊布尼茨在1660年代發明,許多計π的無窮級數出現。牛頓自己就用反正弦()數列在1655年或1666年將π近似到第15位小數,後來寫到「我很羞愧告訴你我為了計算它用了多少數字,我當時沒有做其他事。」[70]
蘇格蘭數學家詹姆斯·格雷果里在1671年發現了馬德哈瓦公式,萊布尼茨也在1674年發現:[2]:188–189[72]
這公式即為格雷果里-萊布尼茨公式,在時數值為。[72]1699年時英國數學家亞伯拉罕·夏普用格雷果里-萊布尼茨公式,在時計算,計算到π的第71位小數,打破由多邊形算法得到的第39位小數的記錄。[2]:189格雷果里-萊布尼茨公式在時非常簡單,但收斂到最終值的速度非常慢,現在不會再用此公式來計π。[2]:156
約翰·梅欽在1706年用格雷果里-萊布尼茨級數產生了可以快速收斂的公式:[2]:192–193
梅欽用這公式計到π第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
隨著一項一項的值加入總和中,只要項次夠多,總和最後會慢慢接近π。不過此數列的收斂速度很慢,要到50萬項之後,才會精確到π的第五位小數[74]。
尼拉卡莎在15世紀發展了π的另一條無窮級數,收斂速度比格雷果里-萊布尼茨公式快很多:[75]
以下比較兩條級數的收斂速率:
π的無窮級數 | 第1項 | 前2項 | 前3項 | 前4項 | 前5項 | 收斂到 |
---|---|---|---|---|---|---|
4.0000 | 2.6666… | 3.4666… | 2.8952… | 3.3396… | 3.1415… | |
3.0000 | 3.1666… | 3.1333… | 3.1452… | 3.1396… |
計算前五項後,格雷果里-萊布尼茨級數的和跟π的誤差為0.2,而尼拉卡莎級數和的誤差為0.002。尼拉卡莎級數收斂快很多,也甚為適合用來計π的值。收斂更快的級數有梅欽類公式及楚德諾夫斯基算法,後者每計一項就可以得到14位正確的小數位[73]。
並非所有和π有關的研究都旨在提高計算它的準確度。1735年,歐拉解決了巴塞爾問題,建立了所有平方數倒數和與π的關係。之後歐拉發現了歐拉乘積公式,得到了π、質數的重要關聯,對日後黎曼ζ函數的研究影響深遠。[76]
1761年,瑞士數學家約翰·海因里希·朗伯用正切函數的無窮連分數表達式證明了π是無理數。[2]:5[77]1794年,法國數學家阿德里安-馬里·勒壤得證明了也是無理數。1882年,德國數學家費迪南德·馮·林德曼證明了對任何非零代數數,都是超越數,該結論後來由魏爾斯特拉斯推廣為林德曼-魏爾斯特拉斯定理。據此定理和歐拉公式,π只能是超越數,進而證實了勒壤得和歐拉提出的π超越性猜想。[2]:196[78]哈代在其著作《數論導引》中則稱此證明在提出後,經過希爾伯特、施瓦次和其他一些人化簡過。[79]
在用π專指「圓周率」之前,希臘字母即已用於幾何概念中[2]:166。威廉·奧特雷德在1647年起在《數學之鑰》(Clavis Mathematicae)就已經用及(對應p和d的希臘字母)來表示圓的周長及直徑的比例。
威廉·瓊斯在他1706年出版的《新數學導論》(A New Introduction to the Mathematics)提到了,是目前已知最早專門用希臘字母表示圓周和其直徑比例的人[80]。這希臘字母第一次出現是在書中討論一塊半徑1的圓時提到「其圓周長一半()」。瓊斯選用可能因它是希臘文「周邊」一詞「περιφέρεια」的首字母[81]。不過瓊斯提到,他那些有關的算式出自「真正聰明的約翰·梅欽先生」,人們推測在瓊斯之前,約翰·梅欽就已開始用表示圓周率[2]:166。
瓊斯在1706年開始使用此希臘字母,但直到萊昂哈德·歐拉在其1736年出版的《力學》中開始使用之後,其他數學家才紛紛開始用指代圓周率。在此之前,數字家可能用像c或p之類的字母代表圓周率[2]:166。歐拉與歐洲其他數學家間時常互相寫信來往,的用法迅速傳播開來[2]:166。1748年歐拉在他的《無窮小分析引論》再一次提到了,寫道:「簡潔起見,我們將此數字寫為,等於半徑為1的圓周長的一半。」這表示方式之後也推展到整片西方世界[2]:166。
二十世紀中期計算機技術發展、革新再次引發了計算π位數的熱潮。美國數學家約翰·倫奇及李維·史密斯在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下式即為一例,其中用到了模方程式:
這無窮級數收斂速度遠快於絕大多數反正切數列,包括梅欽公式。[2]:104第一位使用拉馬努金公式計算π並取得進展的是比爾·高斯珀,他在1985年算得了小數點後一千七百萬位。[2]:104, 206拉馬努金公式開創了現代數值近似算法的先河,此後波爾文兄弟和楚德諾夫斯基兄弟進一步發展了這類算法。[2]:110–111後者於1987年提出了楚德諾夫斯基公式,如下所示:
此公式每計算一項就能得到π的約14位數值[88],因而用於突破圓周率的數位的計算。利用這公式,楚德諾夫斯基兄弟於1989年算得π小數點後10億(109)位,法布里斯·貝拉於2009年算得2.7千億(2.7×1012)位,亞歷山大·易和近藤滋在2011年算得一萬億(1013)位。[2]:110–111, 206[89][90]類似的公式還有拉馬努金-佐藤級數。
2006年,加拿大數學家西蒙·普勞夫利用PSLQ整數關係算法[91]按照以下模版生成了幾條計算π的新公式:
統計模擬法是以機率統計理論為指導的一類非常重要的計數方法,經大量重複試驗計算事件發生頻率,按照大數法則(即當試驗次數充分大時,頻率充分接近機率)可以求得的近似值[93]。 布芬(Buffon)投針問題就是其中一項實例:長度的針隨機往畫滿間距的平行線的平面上拋擲次, 如果針與平行直線相交次,充分大就可根據以下公式算出的近似值[94]:
用統計模擬法計的另一例子是隨機往內切四分之一圓的正方形內拋擲大量點,落在四分之一圓內的點的數量與拋擲點的總量的比值會近似於。[2]:39–40[95]
此外還可用隨機漫步試驗,並用統計模擬法計算值,如拋擲一枚均勻的硬幣次,並記錄正面朝上的次數,所得結果中,正面朝上的次數服從二項分布且
因為硬幣均勻,所以N次試驗中每次試驗結果相互獨立。由此可定義一系列獨立的隨機變數,當拋擲結果為正面時否則為-1,且且取何值有相同機率(即,正面朝上和背面朝上的機率相同)。對隨機變數求和可得
設k為「硬幣正面朝上的次數」減去「硬幣反面朝上的次數」,即可得到。轉換式子,得,因此
可證明[96],
並且當N變大時,的值會漸近於,因此當N充分大時可根據以下公式算出的近似值:[97]
和其他計算值的方法相比,蒙地卡羅方法收斂速度很慢,而且無論實驗多少次,都無從得知的估值已經精確到第幾位。因此,當追求速度或精度時,蒙地卡羅方法不適合用來估計。[2]:43[98]
1995年引入的兩條算法開闢了研究的新途徑。因為每計算出一位數字,該數就會像流過閥門的水一樣不會再出現在後續的計算過程中,這種新進算法叫閥門算法。[2]:77–84[99]這就與無窮級數及迭代算法形成對比——無窮級數和迭代算法自始至終的每一步計算都會涉及到之前所有步驟計算出的中間值。[2]:77–84
1995年,美國數學家斯坦·瓦格納和斯坦利·拉比諾維茨(Stanley Rabinowitz)發明了一種簡單的閥門算法[99][2]:77[100],其運算速度類似arctan演算法,但速度比迭代算法慢[2]:77。
貝利-波爾溫-普勞夫公式(BBP)是另一條閥門算法,屬於一種位數萃取演算法。1995年,西蒙·普勞夫等人發現[2]:117, 126–128[101]
這公式和其他公式不同,可以計算的任何十六進小數位,而不用計算前面全部小數位[2]:117, 126–128。十六進數位可計算得到特定二進數位;想要得到八進制數位的話,計算一、兩位十六進小數即可。目前也已發現一些這種演算法的變體,不過還沒有發現針對十進制、可以快速生成特定小數位的位數萃取演算法[102]。位數萃取演算法的一項重要用途是用來確認聲稱是計算到小數位數的新紀錄:若有聲稱是新紀錄的計算結果出現,先將十進制的數值轉換到十六進制,再用貝利-波爾溫-普勞夫公式去確認最後一些位數(用亂數決定),若這些位數都對,就能有一定把握認為此計算結果是對的[90]。
1998年到2000年間,分布式計算計畫PiHex用貝拉公式(貝利-波爾溫-普勞夫公式的一種變體)計算第1015位,結果是0[2]:20[103]。2010年9月,有雅虎員工用公司的Apache Hadoop應用程式在上千台電腦計算π在2×1015位開始往後256位,其第2×1015位剛好也是0[104]。
利用伽瑪函數計算
伽瑪函數,,可以被用作計算圓周率。
與圓密切相關,出現在許多幾何學和三角學的公式中(特別是與圓、橢圓和球體相關的那些)。 此外,也出現在其他學科的重要公式中,比如統計學、物理學,傅立葉分析和數論的公式。
出現在基於圓的幾何圖形(如橢圓、球、圓錐與環面)的面積、體積公式中。下面是一些用到π的常見公式:[9]
上述公式是n維球的體積與其邊界((n−1)維球的球面)的表面積的特殊情況,具體將在後文給出解釋。
描述由圓生成的圖形的周長、面積或體積的定積分常涉及π。例如,表示半徑為1的半圓的面積的積分為[105]
的積分表示上半圓(此處的平方根由畢氏定理得出),從-1到1的積分可用來計算計算半圓與x 軸間的面積。
三角函數要用到角,而數學家常用弧度作角度單位。π在弧度制起重要作用,數學家將周角,即360度定義為2π度。[106]由這條定義可得,180度=π弧度,1度=弧度。[106]因此,常用的三角函數的週期為的倍數;例如,正弦和餘弦週期為π,[107]任何角度和任何整數都有
常數出現在將平面微分幾何及其拓撲學聯繫起來的高斯-博內定理中。具體來說,如果緊曲面Σ的高斯曲率為,那麼有
其中是該曲面的歐拉示性數,是整數。[108]例如,曲率為1(也就是說其曲率半徑也為1,對於球面而言此時的曲率半徑與半徑重合)的球面的表面積。球面的歐拉特徵數可以通過其同源組計算,其結果為2。於是,便得出
即為半徑為1的球面的表面積公式。
向量分析是與向量場的性質有關的微積分的分支,並有許多物理用途,例如用在電磁學中。位於三維笛卡兒坐標系原點的點源的牛頓位勢為[110]
表示位於距原點的單位質量(或電荷)的勢能,而是維度常數。在這裡由表示的場可以是(牛頓)引力場或(庫侖)電場,是位勢的負梯度:
特殊情況有庫侖定律和牛頓萬有引力定律。高斯定律表明,通過包含原點的任何平滑、簡單、封閉、可定向曲面的場的向外通量等於:
標準形式會將的這因子吸收到常數中,但這種說法表明了它必須出現在「某處」。此外,是單位球面的表面積,但並沒有假設是球面。然而,作為散度定理的結果,由於遠離原點的區域是真空(無源的),只有中的表面的同調類與計算積分有關,因此可以由相同同調類中的任何方便的表面代替,特別是球形,因為球面坐標可以用於計算積分。
高斯定律的結果之一是位勢的負拉普拉斯算子等於狄拉克δ函數的倍:
通過摺積就能得到物質(或電荷)的更一般分布,給出卜瓦松方程式
其中是分布函數。
常數在與愛因斯坦場方程式中的四維勢起類似的作用,愛因斯坦方程式是形成廣義相對論基礎的一條基本公式,並且把引力的基本相互作用描述為物質和能量引起的時空彎曲的結果:[111]
是里奇曲率張量,是純量曲率,是度量張量,是宇宙學常數,是萬有引力常數,是真空中的光速,而是應力-能量張量。愛因斯坦方程式的左邊是度量張量的拉普拉斯算子的非線性模擬,並化簡(reduce)至在弱域的極限,而右邊是分布函數的模擬乘以。
在複分析中,沿複數平面若爾當曲線的圍道積分是研究解析函數的重要手段之一。簡化版的柯西積分公式表明,對任何若爾當曲線內任一點,以下圍道積分給出:[112]
該命題是柯西積分定理的直接推論,後者表明上述圍道積分在圍道的同倫轉換下保持不變,因而沿任一曲線的積分和沿以為圓心的圓周積分的結果相同。更為一般地,該公式對不通過點的任意可求長曲線都成立,但等式右邊要乘以曲線關於該點的卷繞數。
一般形式的柯西積分公式建立了全純函數在若爾當曲線上的值與曲線內任意點處值的關係:[113][114]
柯西積分定理是留數定理的一項特例。根據留數定理,在區域內除去有限個點解析的亞純函數在邊界上的圍道積分與函數在這些點的留數之和滿足:
階乘函數的值等於所有小於等於的正整數之積,它的定義域只包含非負整數。Γ函數則是階乘的推廣。它在複數平面的右半平面定義為:
再利用解析延拓可以將它的定義域擴展到除去非正整數的整塊複數域。當自變數取正整數時,函數給出階乘;當自變數取半整數時,計算結果含有。例如,。[115]
根據魏爾施特拉斯分解定理,函數可分解為如下的無窮乘積:[116]
是歐拉-馬斯刻若尼常數。利用該分解公式和函數在的值,亦可以證明沃利斯乘積式。函數和黎曼ζ函數、函數行列式的恆等式存在關聯,其中扮演著重要的角色。
函數常用於計算維歐氏空間中n 維球的體積和n 維球面的表面積。對維歐氏空間中半徑為維球,其體積和表面積滿足:[117]
兩者還滿足如下的關係式:
當很大,用函數可得到階乘的近似公式,稱斯特靈公式[118],等價於:
斯特靈近似的幾何應用之一是埃爾哈特體積猜想。將維歐幾里得空間的單體記作,則表示該單體的所有面擴大。於是
這是僅含一點晶格點之凸體體積的(最佳)上界[119]。
黎曼ζ函數 在數學的許多領域均有應用。當自變數 ,可寫作
找到這無窮級數的解析解是數學界著名的「巴塞爾問題」。1735年,歐拉解決了這問題,他得到該無窮級數等於[76]。歐拉的結論可推導出數論中一項結果,即兩隨機整數互質(無公因數)的機率為 [2]:41–43[120]。整數可由質數整除的機率為(例如,連續7個正整數只有一個可以7整除),任取兩隨機整數都能以質數整除的機率為,至少有一數不能整除的機率則為。又,一隨機整數能否以兩不同質數整除是相互獨立事件,兩隨機整數互質的機率可以表示成關於所有質數的無窮乘積[121]
這結論可結合隨機數生成器,用統計模擬法計的近似值。[2]:43
巴塞爾問題的結論意味著幾何導出量的數值與質數的分布有著深刻的關聯。巴塞爾問題是谷山-志村定理的一種特殊情況,是安德烈·韋伊對玉河數的猜想的一項特例,即猜想一個這種形式的算術量關於所有質數的無窮乘積能夠等於一個幾何量——某局部對稱空間體積的倒易。巴塞爾問題中,這空間是雙曲3-流形SL2(R)/SL2(Z)。[122]
函數同樣滿足黎曼方程式的公式,其中用到了和伽瑪公式:
除此之外, 函數導數也滿足
最終的結果是可以從諧振子泛函行列式中求得。這泛函行列式可以無窮乘積展開式計算,而且這種方法等價於沃利斯乘積公式。[123]這種方法可用於量子力學,尤其是玻爾模型中的變分。[124]
週期函數的傅立葉級數很自然出現了。週期函數即實數的小數部分所構成群上的函數。傅立葉分解指出,上的複值函數可表示為無窮多個的酉特徵的線性疊加之和。也就是說,到圓群(模為1的複數組成的乘法群)的映射是連續群同態。的特徵都有的形式,是一條定理。
有唯一的特徵值,直到複共軛,那是一群同態。在圓群用哈爾測度,常數是這特徵值的拉東-尼科迪姆導數值的一半。其他的特徵值的導數值為的正整數倍。[23]因此,常數是獨特的數字,以至於配備了其哈爾測度的群,有對於整數倍的點陣的龐特里亞金對偶性[126]。這是卜瓦松和公式的一維版本。
常數與模形式和Θ函數密切相關——比如,橢圓曲線中的j變量就很大程度涉及楚德諾夫斯基算法(一種快速計算π的方法)。
模形式是以在上半平面的全純函數的在模群(或其子群,是的一格)下的變換特性歸納。Θ函數便是一例:
它是一種名為雅可比形式的模形式,[127]有時以諾姆表達。
常數是特殊常數,它會使雅可比函數形成自守式,即該函數會以特定方式變換。有若干恆等式在所有自守式下成立。,例如:
它使得必然在離散海森伯群下以表示(representation)變換。一般模形式和其他函數也包含,這也是根據史東–馮紐曼定理。[127]
柯西分布的香農熵等於,也含。
柯西分布在位勢論中扮演著重要的角色因為它是最簡單的福斯坦堡測度和與在半平面上做布朗運動相關聯的經典卜瓦松核[128]。共軛諧波函數以及希爾伯特轉換與卜瓦松核的漸近線有關。希爾伯特轉換是由奇異積分的柯西主值給出的積分轉換
常數是唯一的(正)歸一化因子因此定義了一個在實數軸上的平方可積分實值函數的希爾伯特空間上的線性複結構[129]。 和傅立葉轉換一樣,希爾伯特轉換就其在希爾伯特空間的轉換特性而言可以完全特徵化。直到歸一化,它是唯一的與正膨脹對易且與實數軸的所有反射反對易有界線性算子[130]。常數是唯一能使這轉換么正的歸一化因子。
大衛·波(David Boll)在1991年發現在曼德博集合碎形也有π出現[131]。他檢查在曼德博集合在位置的特性。若考慮坐標在「頸部」的點,而趨近零,在發散之前迭代的次數和相乘,會趨近。若是在右側尖點處附近的點也會有類似的特性:在發散之前迭代的次數和的平方根相乘,也會趨近[131][132]。
與圓以及球坐標系關係密切,即使不是物理常數,也常出現在描述宇宙的基本原則方程式中。比方說,經典力學領域的簡單公式給出長L的單擺小幅擺動的近似週期,為地球引力加速度常數。[133]
海森堡不確定性原理是量子力學的基本公式,表明測量粒子時,其位置不確定度()與動量不確定度()不可能同時達到任意小(為普朗克常數):[134]。
近似三這特性,和電子偶素的半衰期相對較長有密切的聯繫。其半衰期的倒數和精細結構常數的關係為[135],為電子質量。
許多結構工程的公式也有,例如歐拉推導的挫曲公式說明了長度為、截面二次軸矩為I的細長形物體,在不挫曲的條件下可以承受的最大軸向負載[136]:
流體動力學的斯托克斯定律中也有。斯托克斯定律是半徑約為的小球體在黏度的流體中以速度運動時會受到的阻力滿足[137]:
在理想狀態下,河的曲折程度(河道本身的長度與源頭到入海口的比值)隨著時間的推移逐漸趨向於。河流外邊際的快速水流彎曲會使河流內邊際加倍侵蝕,河道變得更彎曲,整條河彎折更厲害。然而,這股彎折勁兒最終會導致河流折回一開始彎折的地方,導致「短路」,並形成河跡湖。這兩種相反因素使河道長度與源頭到入海口的比值的平均值為π。[138][139]
π文字學(或譯作圓周率的語言學)是指記住的大量位值[2]:44–45,並將其世界紀錄載於金氏世界紀錄大全的做法。維爾·美拿(Rajveer Meena)於2015年3月21日在印度於9小時27分鐘內背誦了7萬位的π,創下金氏世界紀錄大全認證的世界紀錄。[140]2006年,日本退休工程師原口證自稱已經背誦了十萬位小數,但他未獲金氏世界紀錄大全認證。[141]
常用於記憶π的一項技巧是背誦以單詞長度代表數值的故事或詩歌:第一單詞有三字母,第二單詞有一字母,第三單詞有四字母,第四單詞有一字母,第五單詞有五字母,如此類推。早期例子是英國科學家詹姆士·金斯設計的詩歌:「How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics.」[2]:44–45這類詩歌有時在英文中稱為「piem」。除了英文,用於記憶π的詩歌亦有不同語言的版本[2]:44-45。但是,記憶的人一般並不以詩歌記憶來創下紀錄,而是用如記憶數字規律或軌跡法的方法。[142]
有好幾位作家仿照上述記憶技巧,用的數值創作了新型的約束寫作方式,當中單詞長度須符合的數值。《The Cadaeic Cadenza》以上述技巧包含了前3835位的值[143],一本標準長度的書《Not a Wake》有一萬單詞,其中各單詞亦代表了的一位。[144]
也許因為的公式很簡短而且四處可見, 比其他數學常數在流行文化中更常見[註 3]。
在2008年由英國公開大學及英國廣播公司聯合製作的記錄片《數學的故事》於2008年十月由英國廣播公司第四台播放。影片講述了英國數學家馬庫斯·杜·索托伊在到訪印度研究當地三角學的貢獻時,展示出歷史上π最精確的計算公式的資訊圖形。[147]
巴黎的科學博物館發現宮有間圓形房間叫「房」,牆上刻有的707位數,數字貼在圓頂狀的天花板上,由大型的木製字符組成。數值是1853年由英國數學家威廉·尚克思計算出來,但是該結果於第528位後開始出現謬誤,在1946年發現,1949年修正。[148][2]:50
卡爾·薩根的小說《接觸未來》中則暗示說,宇宙的創造者在π的數字中暗藏了一則資訊。[149]π的數字也用在凱特·布希所出的專輯Aerial中的《Pi》的歌詞裡。[150]
美國人在3月14日慶祝圓周率日,此節日在學生中很流行。[151]一些自稱「數學極客」的人常常用與其數位來創作一些數學或技術圈內人士才能領會到的笑話。麻省理工學院則有幾句包含「3.14159」的大學歡呼口號。[152]2015年的圓周率日格外重要,按照美式寫法,當天的日期時間3/14/15 9:26:53
較其他圓周率日包含更多位數的。[153]
在北電網絡於2011年舉行的技術專利拍賣會上,谷歌用了一些包含在內的數學或科學常數來競價。[154]
在1958年,阿爾伯特·伊格爾提議將換成τ(tau)以便簡化公式。在此定義為的兩倍[155]。然而,沒有任何其他作者曾這樣使用過。有些人使用不同的值,。[156]這些人稱不論是作為弧度制下圓周長的1轉還是作為弧長與半徑的比值(而不是與直徑的比值)都比自然,也能因此簡化許多公式。[157][158]有媒體報道稱,因為的值大小約為6.28,現已有人在6月28日慶祝「節」,並吃「兩個派」;[159]然而,主流數學界還並未使用。[160]
1897年,有業餘美國數學家試圖藉印第安納州議會來通過後世所謂印第安納圓周率法案的法案。這法案試圖以法律命令強制規定數學常數之值而臭名遠播。該法案描述化圓為方的方法,並間接提到了的錯誤值,例如3.2。該法案通過了印第安納州眾議院的表決,但參議院否決之。[2]:211–212[161][162]