π

圓周率

圓周率(Pi)是圓的周長與直徑的比值,一般用希臘字母π表示,是一個在數學及物理學中普遍存在的數學常數。π也等於圓形之面積與半徑平方之比。是精確計算圓周長、圓面積、球體積等幾何形狀的關鍵值。在分析學里,π可以嚴格地定義為滿足sin x = 0的最小正實數x。

圓周率用希臘字母 π(讀作pài)表示,是一個常數(約等於3.141592654),是代表圓周長和直徑的比值。它是一個無理數,即無限不循環小數。在日常生活中,通常都用3.14代表圓周率去進行近似計算。而用十位小數3.141592654便足以應付一般計算。即使是工程師或物理學家要進行較精密的計算,充其量也只需取值至小數點後幾百個位。

1965年,英國數學家約翰·沃利斯(John Wallis)出版了一本數學專著,其中他推導出一個公式,發現圓周率等於無窮個分數相乘的積。2015年,羅切斯特大學的科學家們在氫原子能級的量子力學計算中發現了圓周率相同的公式。

2021年8月18日,圓周率π計算到小數點后62.8萬億位,創下該常數迄今最精確值記錄。

歷史發展


實驗時期

一塊古巴比倫石匾(約產於公元前1900年至1600年)清楚地記載了圓周率 = = 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位小數大關,他利用了如下公式:
其中arctan x可由泰勒級數算出。類似方法稱為“梅欽類公式”。
斯洛維尼亞數學家Jurij Vega於1789年得出π的小數點后首140位,其中只有137位是正確的。這個世界紀錄維持了五十年。他利用了梅欽於1706年提出的數式。
到1948年英國的弗格森(D. F. Ferguson)和美國的倫奇共同發表了π的808位小數值,成為人工計算圓周率值的最高紀錄。

計算機時代

電子計算機的出現使π值計算有了突飛猛進的發展。1949年,美國製造的世上首部電腦-ENIAC(Electronic
圓周率
圓周率
Numerical Integrator And Computer)在阿伯丁試驗場啟用了。次年,里特韋斯納、馮紐曼和梅卓普利斯利用這部電腦,計算出π的2037個小數位。這部電腦只用了70小時就完成了這項工作,扣除插入打孔卡所花的時間,等於平均兩分鐘算出一位數。五年後,IBM NORC(海軍兵器研究計算機)只用了13分鐘,就算出π的3089個小數位。科技不斷進步,電腦的運算速度也越來越快,在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年祖沖之中國73.1415926<π<3.1415927
499年阿耶波多印度3π=3.1416
598年婆羅摩笈多印度1π=3.162277...
800年花拉子米烏茲別克3π=3.1416
12世紀婆什迦羅第二印度4π=3.14156
1220年斐波那契義大利3π=3.141818
1400年Madhava10π=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π=127位小數(只有112位正確)
1723年建部賢弘π=41位小數
1730年Kamataπ=25位小數
1734年萊昂哈德·歐拉引入希臘字母π並肯定其普及性
1739年松永良弼π=50位小數
1761年約翰·海因里希·蘭伯特證明π是無理數
1775年歐拉指出π可能是超越數
1794年Jurij Vegaπ=140位小數(只有136位正確)
1794年阿德里安-馬里·勒讓德-
1841年Rutherfordπ=208位小數(只有152位正確)
1844年Zacharias Dase及Strassnitzkyπ=200位小數
1847年Thomas Clausenπ=248位小數
1853年Lehmannπ=261位小數
1853年William Rutherfordπ=440位小數
1855年Richterπ=500位小數
1874年William Shanksπ=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“藍色基因/p”超級電腦π的60,000,000,000,000位二進位小數 
註:上表 正確位數是指小數點后的位數。

算準記錄

小數點后位數首次算準者首次算準時間
1巴比倫人前20世紀
2-3阿基米德前3世紀(距離上次1700年)
4-5劉徽263年(距離上次563年以上)
6-7祖沖之480年(距離上次217年)
8-10Madhava1400年(距離上次920年)
11-16Jamshid Masud Al Kashi1424年(距離上次24年)
17-20魯道夫·范·科伊倫1596年(距離上次172年)
21-321615年(距離上次19年)
33-35威理博·司乃耳,范·科伊倫的學生1621年(距離上次6年)
36-71Abraham Sharp1699年(距離上次78年)
72-100John Machin1706年(距離上次7年)
101-112De Lagny1719年(距離上次13年)
113-136Jurij Vega1794年(距離上次75年)
137-152Rutherford1841年(距離上次47年)
153-200Zacharias Dase及Strassnitzky1844年(距離上次3年)
201-248Thomas Clausen1847年(距離上次3年)
249-261Lehmann1853年(距離上次6年)
262-440William Rutherford1853年(距離上次0年)
441-500Richter1855年(距離上次2年)
501-527William Shanks1874年(距離上次19年)
528-620D. F. Ferguson1946年(距離上次72年)
621-7101947年(距離上次1年)
711-8081947年(距離上次0年)
備註:這裡只列出人工計算的最高記錄,808位

記號


是第十六個希臘字母的小寫。這個符號,亦是希臘語 περιφρεια (表示周邊,地域,圓周等意思)的首字母。1706年英國數學家威廉·瓊斯(William Jones ,1675-1749)最先使用“π”來表示圓周率。1736年,瑞士大數學家歐拉也開始用 π表示圓周率。從此,π便成了圓周率的代名詞。
要注意不可把 和其大寫Π混用,後者是指連乘的意思。

公式


π[圓周率]
π[圓周率]
圓周率(π )一般定義為一個圓形的周長(C )與直徑( d)之比: ,或直接定義為單位圓的周長的一半。由相似圖形的性質可知,對於任何圓形,的值都是一樣,這樣就定義出常數π。
注意:將 定義為單位圓的周長的一半是有意義的,這是因為從現代數學的角度來看,直徑為d、半徑為r的圓的周長C由以下積分給出:
上式中令,由定積分的換元法可得:
其中 是單位圓周的周長(C的表達式中取r=1即得)。若定義,則,與我們熟知的周長公式相符。
而半徑為r的圓的面積S由以下積分給出:
令,由定積分的換元法可得:
其中 是單位圓的面積(S的表達式中取r=1即得)。利用分部積分法
於是,
π[圓周率]
π[圓周率]
因此,我們得到關係式:
這樣一來也得到了我們熟知的圓面積公式
第二個做法是,以圓形半徑為邊長作一正方形,然後把圓形面積和此正方形面積的比例定為,即圓形之面積與半徑平方之比。
定義圓周率不一定要用到幾何概念,比如,我們可以定義π 為滿足
的最小正實數。
這裡的正弦函數定義為冪級數

特性


把圓周率的數值算得這麼精確,實際意義並不大。現代科技領域使用的圓周率值,有十幾位已經足夠了。如果以39位精度的圓周率值,來計算宇宙的大小,誤差還不到一個原子的體積。以前的人計算圓周率,是要探究圓周率是否循環小數。自從1761年蘭伯特證明了圓周率是無理數,1882年林德曼證明了圓周率是超越數后,圓周率的神秘面紗就被揭開了。
π在許多數學領域都有非常重要的作用。

幾何

平面圖形周長面積
圓環
扇形
註:① C 為周長, S 為面積, l 為弧長; d 為直徑, r為半徑(內圓半徑), R為外圓半徑, 為圓心度數。
②周長、弧長用長度單位,面積用面積單位。
立體圖形表面積體積
圓柱
圓錐
註:① C 為底面周長, 為底面積, 為側面積, S 為表面積, V為體積;d 為底面直徑, r 為底面半徑, h為高。
②底面周長用長度單位,表面積(含底面積和側面積)用面積單位,體積用體積單位或容積單位。

代數

π是個無理數,即不可表達成兩個整數之比,是由瑞士科學家約翰·海因里希·蘭伯特於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日為圓周率近似日(英國式日期記作,看成圓周率的近似分數)
有數學家認為應把"真正的圓周率"定義為2π,並將其記為τ(發音:tau)。
2019年3月14日,谷歌宣布日裔前谷歌工程師愛瑪(EmmaHarukaIwao)在谷歌雲平台的幫助下,計算到圓周率小數點后31.4萬億位,準確的說是31415926535897位,比2016年創下的紀錄又增加數萬億位。據了解,愛瑪的團隊使用了一個名為ycruncher的程序,能將π計算到小數點后數萬億位。該程序由谷歌雲平台計算引擎上運行的25個虛擬機驅動。而2016年紀錄的創造者皮特(PeterTrueb)是用一台電腦計算出來的。這項計算需要170TB的數據,與整個美國國會圖書館印刷藏品數據量大致相同,愛瑪經過大約4個月的計算才打破了此前的世界紀錄。
2021年8月18日,瑞士研究人員使用一台超級計算機,歷時108天,將著名數學常數圓周率π計算到小數點后62.8萬億位,創下該常數迄今最精確值記錄。