圓周率（簡(jiǎn)稱(chēng)π）是數(shù)學(xué)中最神秘又熟悉的數(shù)字之一。無(wú)論是在建筑、工程，還是航天、物理，π都發(fā)揮著至關(guān)重要的作用。簡(jiǎn)單的說(shuō)，圓周率就是一個(gè)圓的周長(zhǎng)與其直徑之比?？雌饋?lái)要得到這個(gè)數(shù)字很簡(jiǎn)單，只要量一下一個(gè)圓的周長(zhǎng)，用這個(gè)周長(zhǎng)除以直徑，圓周率就得到了。
但困難的是，這個(gè)圓用什么尺子才能夠?qū)⒅荛L(zhǎng)量的十分精確呢？從古到今，幾千年過(guò)去了，這個(gè)傷腦筋的難題也沒(méi)有完全得到解決。而得到的結(jié)果就是，圓周率是一個(gè)無(wú)理數(shù)，意味著它的小數(shù)部分無(wú)限多，永無(wú)止境，且沒(méi)有任何循環(huán)規(guī)律。
那么這個(gè)小數(shù)位數(shù)到底有多少呢？沒(méi)人知道，因?yàn)槠袢藗冎恢朗菬o(wú)限的，無(wú)限就是永無(wú)止境的。幾千年前，古人們就試圖用手工的方法來(lái)測(cè)量出精確的圓周長(zhǎng)，他們用尺規(guī)和幾何方法手工計(jì)算π，得到了一個(gè)大致的數(shù)據(jù)，最終精確到了小數(shù)點(diǎn)后九位數(shù)，即3.141615926。這是咱們老祖宗祖沖之得到的，可以自豪一下。
今天，隨著科學(xué)的進(jìn)步，超級(jí)計(jì)算機(jī)將π推算到100萬(wàn)億位，但也只能說(shuō)明它還是個(gè)近似值，依然沒(méi)有窮盡。100萬(wàn)億位是什么概念？就是一個(gè)人如果從出生就開(kāi)始算數(shù)，不吃不喝24小時(shí)不停的算，平均每秒鐘算兩個(gè)數(shù)的話，要算158萬(wàn)多年。也就是說(shuō)，一個(gè)人從原始時(shí)代猿人老祖宗開(kāi)始，世世代代接著算，算到現(xiàn)在也還沒(méi)算完。
那么，圓周率真的有這么多位嗎，人類(lèi)開(kāi)始時(shí)如何得到的呢？我們一起來(lái)了解。
古人的智慧：用割圓術(shù)計(jì)算π
無(wú)論是西方還是東方，早在遠(yuǎn)古時(shí)代就有那么一群智者對(duì)圓周率有了興趣，并開(kāi)始研究。東西合璧殊途同歸，得到了圓周率的近似值。具體說(shuō)來(lái)，有如下一些典型代表：
阿基米德的“多邊形逼近法”
最早系統(tǒng)計(jì)算π的數(shù)學(xué)家之一是古希臘的阿基米德（約公元前287—前212年）。他的想法很簡(jiǎn)單：
- 畫(huà)一個(gè)圓，然后在圓里面畫(huà)一個(gè)正六邊形，再在外面畫(huà)一個(gè)正六邊形。
- 計(jì)算這兩個(gè)六邊形的周長(zhǎng)，就可以得到一個(gè)π的上下限。
- 然后，把六邊形的邊數(shù)加倍，變成十二邊形，再算一次周長(zhǎng)。
- 繼續(xù)加倍，變成二十四邊形、四十八邊形……，直到逼近圓的真實(shí)周長(zhǎng)。
阿基米德用這種方法，把π的值估算在3.1408 到 3.1429 之間，這在沒(méi)有計(jì)算工具的時(shí)代已經(jīng)是非常驚人的成就！
祖沖之：計(jì)算到全球領(lǐng)先1000年的精度
祖沖之是南北朝時(shí)期（公元5世紀(jì)）中國(guó)古代數(shù)學(xué)家，他改進(jìn)了割圓術(shù)，將一個(gè)正元分割出24576個(gè)邊，通過(guò)對(duì)這24576邊形的測(cè)量計(jì)算，得到了比阿基米德的精度更高的π近似值：3.1415926到3.1415927，這個(gè)精度在世界范圍內(nèi)領(lǐng)先了1000年！
但問(wèn)題是，割圓術(shù)計(jì)算π非常慢。如果想要更精確的π，就得畫(huà)邊數(shù)更多的多邊形，計(jì)算量成倍增長(zhǎng)，最終會(huì)變得難以承受。所以，祖沖之之后的1000年間，π的精度再也難以提升。到了近代，數(shù)學(xué)家們開(kāi)始尋找更高效的計(jì)算方法，讓?duì)械挠?jì)算速度大幅提升。
數(shù)學(xué)公式的力量：不用畫(huà)圖，也能精確計(jì)算π，且比割圓術(shù)快很多
古人計(jì)算π是用形狀逼近圓，而現(xiàn)代數(shù)學(xué)家則用公式直接計(jì)算π，這種方法比割圓術(shù)快得多，且更精準(zhǔn)。比較有代表性的公式有：
18世紀(jì)：馬青公式（Machin's Formula）
1706年，英國(guó)數(shù)學(xué)家**約翰·馬青（John Machin）**提出了一個(gè)快速計(jì)算π的公式：
π=16tan??1(1/5)?4tan??1(1/239)\pi = 16 \tan^{-1} (1/5) - 4 \tan^{-1} (1/239)π=16tan?1(1/5)?4tan?1(1/239)
這個(gè)公式可以讓數(shù)學(xué)家用級(jí)數(shù)展開(kāi)的方式快速計(jì)算π，而不需要使用古代的割圓術(shù)方法。從18世紀(jì)到19世紀(jì)，數(shù)學(xué)家們不斷改進(jìn)這類(lèi)公式，比如：
- 約翰遜（John W. Wrench）等人手工計(jì)算π到808位（1946年），當(dāng)時(shí)是世界紀(jì)錄。
- 手工計(jì)算一次 π 的100位數(shù)值，大約需要數(shù)個(gè)月的時(shí)間，計(jì)算錯(cuò)誤率較高。
高斯-勒讓德算法（Gauss-Legendre Algorithm）——倍增精度
到了1970年代，數(shù)學(xué)家發(fā)現(xiàn)了一種指數(shù)級(jí)提高精度的方法：
- 先選兩個(gè)數(shù)，一個(gè)代表圓的半徑，另一個(gè)代表多邊形的周長(zhǎng)。
- 不斷用數(shù)學(xué)公式調(diào)整這兩個(gè)數(shù)，使它們?cè)絹?lái)越接近真實(shí)的π值。
- 這個(gè)方法每計(jì)算一次，π的精度就能翻倍，計(jì)算速度比割圓術(shù)快得多。
但在現(xiàn)代計(jì)算機(jī)問(wèn)世之前，即便有了更好的公式，對(duì)π計(jì)算的精度和速度比割圓術(shù)快了許多，傳統(tǒng)的手工計(jì)算依然是很復(fù)雜緩慢的。一直到1948年，英國(guó)的弗格森和美國(guó)的倫奇共同發(fā)表了π的808位小數(shù)值，成為人工計(jì)算圓周率值的最高紀(jì)錄 ?，F(xiàn)代計(jì)算機(jī)問(wèn)世，對(duì)π的計(jì)算才有了質(zhì)的飛躍。
計(jì)算機(jī)計(jì)算π：效率突飛猛進(jìn)
1949年：第一次使用計(jì)算機(jī)計(jì)算π
1949年，美國(guó)人約翰·馮·諾依曼（John von Neumann）和尼古拉斯·梅特羅波利斯（Nicholas Metropolis）首次使用電子計(jì)算機(jī)計(jì)算π。他們使用的計(jì)算機(jī)是ENIAC（電子數(shù)字積分計(jì)算機(jī)），計(jì)算到了2037位，打破了人類(lèi)歷史上所有的手工計(jì)算紀(jì)錄。
計(jì)算效率提升對(duì)比：
- 手工計(jì)算808位 π 需要幾個(gè)月，但ENIAC 用70小時(shí)就算出了2037位，速度提升了數(shù)百倍。
- ENIAC每秒可計(jì)算5000次加法，遠(yuǎn)超人工手算。
這次計(jì)算標(biāo)志著計(jì)算機(jī)在數(shù)學(xué)計(jì)算中的首次突破，也讓?duì)械挠?jì)算邁入了現(xiàn)代計(jì)算時(shí)代。
1960年代 - 1980年代：計(jì)算機(jī) π 計(jì)算突破百萬(wàn)位
隨著計(jì)算機(jī)技術(shù)的發(fā)展，數(shù)學(xué)家開(kāi)始使用更加高效的算法，比如高斯-勒讓德算法（Gauss-Legendre Algorithm）。這使得π的計(jì)算速度指數(shù)級(jí)增長(zhǎng)：
計(jì)算效率提升對(duì)比：
- 從1949年 ENIAC 計(jì)算2037位需要70小時(shí)，到1982年 CRAY-1 計(jì)算800萬(wàn)位僅需5小時(shí)，速度提升了10萬(wàn)倍以上。
1990年代至今：超級(jí)計(jì)算機(jī)計(jì)算π
1987年，數(shù)學(xué)家楚德諾夫斯基兄弟提出了一種超快的計(jì)算方法：
1π=12∑k=0∞(?1)k(6k)!(545140134k+13591409)(3k)!(k!)3(640320)3k+3/2\frac{1}{\pi} = 12 \sum_{k=0}^{\infty} \frac{(-1)^k (6k)! (545140134k + 13591409)}{(3k)! (k!)^3 (640320)^{3k+3/2}}π1=12k=0∑∞(3k)!(k!)3(640320)3k+3/2(?1)k(6k)!(545140134k+13591409)
這個(gè)公式雖然看起來(lái)很復(fù)雜，但它的優(yōu)勢(shì)是——計(jì)算速度極快！
- 每次計(jì)算，π的精度可以增加100萬(wàn)位！
- 這是目前計(jì)算機(jī)計(jì)算π時(shí)最常用的方法。
進(jìn)入21世紀(jì)后，數(shù)學(xué)家開(kāi)始使用更加高效的算法（比如楚德諾夫斯基算法）和超級(jí)計(jì)算機(jī)，π的計(jì)算速度和精度進(jìn)一步提升：
計(jì)算效率提升對(duì)比：
手工計(jì)算1000位 π 需要數(shù)月，計(jì)算機(jī)可以在秒級(jí)完成。隨著數(shù)學(xué)公式的優(yōu)化和計(jì)算機(jī)的不斷升級(jí)，計(jì)算速度越來(lái)越快：
- 1949年 ENIAC 計(jì)算2037位 π 需要70小時(shí)，2022年超級(jí)計(jì)算機(jī)計(jì)算100萬(wàn)億位僅需157天。
- 計(jì)算位數(shù)從2000位增長(zhǎng)到100萬(wàn)億位，增長(zhǎng)了50億倍，計(jì)算時(shí)間卻只增長(zhǎng)了1000倍左右，說(shuō)明計(jì)算機(jī)的計(jì)算能力大幅提升。
現(xiàn)代的π值，計(jì)算機(jī)已經(jīng)完全取代了手工計(jì)算。截至2023年，瑞士的研究團(tuán)隊(duì)已經(jīng)用超級(jí)計(jì)算機(jī)將π計(jì)算到100萬(wàn)億（101?）位，打破了歷史記錄。今天，我們計(jì)算π的主要瓶頸已經(jīng)不再是數(shù)學(xué)公式，而是計(jì)算機(jī)硬件的速度和存儲(chǔ)容量。未來(lái)，隨著量子計(jì)算的發(fā)展，π的計(jì)算可能會(huì)變得更加高效。
計(jì)算這么多位的π值到底有什么用呢？
其實(shí)，日常生活中，我們用的π通常不會(huì)超過(guò)3.1416，就連NASA計(jì)算航天器軌道時(shí)，也只用到15位（3.14159265358979）。但是，計(jì)算超高精度的π仍然有很多意義：
- 測(cè)試超級(jí)計(jì)算機(jī)的性能：計(jì)算π需要大量計(jì)算資源，能測(cè)試計(jì)算機(jī)的處理能力和穩(wěn)定性。
- 數(shù)學(xué)研究：數(shù)學(xué)家想知道，π的無(wú)窮小數(shù)部分是否真的完全隨機(jī)，或者是否有隱藏的數(shù)學(xué)規(guī)律。
- 科學(xué)工程：在某些精密科學(xué)研究（比如量子計(jì)算、黑洞模擬）中，需要超高精度的π值。
盡管我們可能永遠(yuǎn)不需要全部的100萬(wàn)億位π，但計(jì)算它的過(guò)程，推動(dòng)了數(shù)學(xué)、計(jì)算機(jī)科學(xué)和工程學(xué)的發(fā)展。
π的精確值計(jì)算仍在繼續(xù)，未來(lái)沒(méi)有止境
因?yàn)棣惺莻€(gè)無(wú)理數(shù)，因此是小數(shù)點(diǎn)是不循環(huán)無(wú)限多的，這就決定了π值的精確值求索永無(wú)止境。從阿基米德的割圓術(shù)到超級(jí)計(jì)算機(jī)的楚德諾夫斯基算法，人類(lèi)對(duì)π的追求已經(jīng)持續(xù)了2000多年。今天，π的計(jì)算精度已經(jīng)遠(yuǎn)超實(shí)際需求，但科學(xué)家們?nèi)匀辉诓粩嗵魬?zhàn)極限，不僅僅是為了計(jì)算π，而是為了推動(dòng)數(shù)學(xué)和計(jì)算機(jī)科學(xué)的發(fā)展。
由此，π不僅僅是一個(gè)數(shù)學(xué)常數(shù)，它是人類(lèi)智慧和科技進(jìn)步的象征。未來(lái)，人類(lèi)將會(huì)用更先進(jìn)的方法，把π計(jì)算得更遠(yuǎn)，或許能揭開(kāi)它更深層的秘密！
時(shí)空通訊原創(chuàng)文章，請(qǐng)尊重作者版權(quán)，感謝閱讀。
熱門(mén)跟貼