圓周率 π 的歷史溯源

圣經(jīng)里有一段不太為人熟知的經(jīng)文，是這么寫(xiě)的：
“他又鑄了一個(gè)銅海，樣式是圓的，從這邊到那邊共十肘，高五肘，圍腰三十肘?！保ā读型跫o(jì)上》7:23）

這段內(nèi)容在《歷代志下》4:2 中也能找到，記載的是公元前 950 年左右建造所羅門(mén)圣殿時(shí)的計(jì)數(shù)規(guī)格。這里提到的圓周率 π 值是 3，顯然不算精確，即便放在當(dāng)時(shí)也不算——要知道更早之前，埃及人和巴比倫人就已經(jīng)算出了更接近真實(shí)值的結(jié)果：埃及人的 π 約為 ，巴比倫人的則是 。不過(guò)這里也要替所羅門(mén)的工匠們說(shuō)句公道話，他們描述的“銅海”本就是某種巨大的銅器，這種大物件既不可能、也沒(méi)必要追求極高的幾何精度。

早期 π 值：從測(cè)量到理論

“圓的周長(zhǎng)與直徑的比值是恒定的一個(gè)數(shù)”這個(gè)事實(shí)，人們?cè)缇椭懒?，具體起源已無(wú)從考證。最早的 π 值（包括圣經(jīng)里的 3），幾乎都是通過(guò)測(cè)量得到的。比如公元前 1650 年左右的埃及《萊因德數(shù)學(xué)紙草書(shū)》（Rhind Papyrus）中，就有充分證據(jù)表明當(dāng)時(shí)人們用 作為 π 的近似值。

而第一個(gè)通過(guò)理論計(jì)算得到 π 值的，大概是敘拉古的阿基米德（Archimedes of Syracuse，公元前 287-前 212 年）。他算出的結(jié)果是：

在簡(jiǎn)單說(shuō)說(shuō)他的證明思路前，有個(gè)點(diǎn)值得注意：這里用到的不等式，其實(shí)暗含了相當(dāng)高的數(shù)學(xué)技巧。阿基米德很清楚——直到現(xiàn)在還有很多人不清楚——π 并不等于 ，他也從沒(méi)宣稱自己找到了 π 的精確值。如果我們?nèi)∷麅蓚€(gè)邊界值的平均值，得到的結(jié)果是 3.1418，誤差僅約 0.0002。

阿基米德的計(jì)算思路

假設(shè)有一個(gè)半徑為 1 的圓，在圓內(nèi)作一個(gè)邊數(shù)為 的正內(nèi)接多邊形，其半周長(zhǎng)記為 ；再作一個(gè)同樣邊數(shù)的正外切多邊形，半周長(zhǎng)記為 。（下面是 時(shí)的示意圖，即邊數(shù)為 12 的多邊形）

這種做法會(huì)定義出兩個(gè)序列：一個(gè)是遞增序列 ，另一個(gè)是遞減序列 ，而這兩個(gè)序列的極限都是 π。

用三角函數(shù)表示的話，這兩個(gè)半周長(zhǎng)可以寫(xiě)成：
，
其中 。同樣，下一個(gè)邊數(shù)（ ）的多邊形半周長(zhǎng)滿足：
，

通過(guò)基礎(chǔ)三角函數(shù)運(yùn)算，能推導(dǎo)出兩個(gè)關(guān)鍵關(guān)系式：

阿基米德從初始值 、 開(kāi)始，先用公式(1)算出 ，再用公式(2)算出 ，接著又用(1)算 ，用(2)算 ，一直算到 和 。最終他得出結(jié)論： 。

這里要明白一點(diǎn)：我們剛才用三角函數(shù)描述，其實(shí)不符合歷史——阿基米德當(dāng)時(shí)沒(méi)有代數(shù)和三角函數(shù)符號(hào)，只能純靠幾何方法推導(dǎo)(1)和(2)；并且他連我們現(xiàn)在用的十進(jìn)制記數(shù)法都沒(méi)有，單靠公式(1)(2)計(jì)算 和 ，本身就不是件容易事。

咱們不該覺(jué)得 “他怎么咋算到 96 邊形就停了”，反而該驚嘆“居然能硬生生算到 96 邊形”—— 要知道在沒(méi)有任何現(xiàn)代數(shù)學(xué)工具的情況下，這已經(jīng)是超了不起的成就了。

歷代數(shù)學(xué)家的 π 值探索

理論上，只要不斷增加多邊形邊數(shù)，就能算出更精確的 π 值。后來(lái)確實(shí)有不少人這么做，比如：

• 托勒密（約公元 150 年）：3.1416

• 祖沖之（430-501 年）： （約 3.1415929）

• 花拉子米（約公元 800 年）：3.1416

• 卡西（約 1430 年）：精確到 14 位小數(shù)

• 韋達(dá)（1540-1603 年）：9 位小數(shù)

• 魯曼（1561-1615 年）：17 位小數(shù)

• 范·塞倫（約 1600 年）：35 位小數(shù)

除了祖沖之，這些改進(jìn)都沒(méi)有理論上的突破，靠的只是更強(qiáng)的計(jì)算耐力。有意思的是，就像其他科學(xué)領(lǐng)域一樣，在公元 400 到 1400 年這一千年里，π 值計(jì)算的領(lǐng)先地位是從歐洲轉(zhuǎn)移到了東方。

花拉子米生活在巴格達(dá)，順便說(shuō)一句，“算法（algorithm）”這個(gè)詞就源自他的名字，而他某本著作標(biāo)題里的“代數(shù)學(xué)（al jabr）”，則成了“代數(shù)（algebra）”一詞的來(lái)源?？ㄎ魃畹牡胤礁繓|，在撒馬爾罕。

文藝復(fù)興后的 π 值公式革命

歐洲文藝復(fù)興最終帶來(lái)了全新的數(shù)學(xué)世界，其中一個(gè)早期成果就是出現(xiàn)了計(jì)算 π 的數(shù)學(xué)公式。最早的公式之一是沃利斯（1616-1703 年）提出的：

而最知名的公式之一，是這個(gè)無(wú)窮級(jí)數(shù)：

這個(gè)級(jí)數(shù)有時(shí)會(huì)歸功于萊布尼茨（1646-1716 年），但實(shí)際上最早是詹姆斯·格雷戈里（1638-1675 年）發(fā)現(xiàn)的。

這些公式既驚艷又出人意料——右邊都是純算術(shù)形式，而 π 最初是從幾何里來(lái)的。它們向我們展示了無(wú)窮過(guò)程能帶來(lái)多么驚人的結(jié)果，也為現(xiàn)代數(shù)學(xué)的豐富內(nèi)涵指明了方向。

不過(guò)從計(jì)算 π 值的角度看，這兩個(gè)公式其實(shí)沒(méi)什么用。比如格雷戈里級(jí)數(shù)，要想得到 4 位精確小數(shù)，誤差得小于 0.00005，這意味著需要約 10000 項(xiàng)才能做到。但格雷戈里還證明了一個(gè)更通用的結(jié)果：

只要令 ，就能得到前面那個(gè) π 的級(jí)數(shù)。但如果利用 ，就能得到一個(gè)收斂快得多的級(jí)數(shù)：

這個(gè)級(jí)數(shù)的第 10 項(xiàng)是 ，小于 0.00005，也就是說(shuō)，只需要 9 項(xiàng)就能得到至少 4 位精確小數(shù)。

還有個(gè)更好的辦法，是利用這個(gè)公式：

然后把 和 分別代入公式(3)，計(jì)算兩個(gè)級(jí)數(shù)的和就行。

顯然，如果能找到這樣的公式：

其中 和 都是大數(shù)，收斂速度會(huì)快得多。1706 年，馬青（Machin）就找到了這樣一個(gè)公式：

其實(shí)只要會(huì)證明公式(4)，證明(5)也沒(méi)什么額外難度，就是計(jì)算過(guò)程麻煩點(diǎn)。當(dāng)然，能想出這個(gè)公式本身，就是另一回事了。

計(jì)算 π 的“執(zhí)念”與錯(cuò)誤

有了這樣的公式，計(jì)算 π 的唯一難點(diǎn)就只剩過(guò)程的枯燥了。不用說(shuō)，還真有人愿意花大量時(shí)間和精力在這種單調(diào)又無(wú)用的追求上。其中一位是英國(guó)的尚克斯（Shanks），他用梅欽公式算出了 π 的 707 位小數(shù)，并在 1873 年發(fā)表了自己多年的計(jì)算成果。尚克斯能留名，原因還挺特別，我們馬上就會(huì)說(shuō)到。

先看看 π 值計(jì)算精度的提升歷程：

• 1699 年：夏普用格雷戈里的成果算出 71 位精確小數(shù)

• 1701 年：梅欽用改進(jìn)方法算出 100 位，后來(lái)的人也沿用他的思路：

• 1719 年：德拉尼算出 112 位精確小數(shù)

• 1789 年：維加算出 126 位，1794 年又算到 136 位

• 1841 年：盧瑟福算出 152 位，1853 年算到 440 位

• 1873 年：尚克斯算出 707 位，但其中只有前 527 位是對(duì)的

尚克斯知道 π 是無(wú)理數(shù)（irrational number），因?yàn)樘m伯特（Lambert）早在 1761 年就證明了這一點(diǎn)。就在尚克斯完成計(jì)算后不久，林德曼（Lindemann）證明了 π 是超越數(shù)（transcendental number）——也就是說(shuō)，π 不是任何整系數(shù)多項(xiàng)式方程的解。其實(shí)林德曼的這個(gè)結(jié)果，也證明了“化圓為方”是不可能的：既然 π 是超越數(shù)，就不可能用尺規(guī)作圖畫(huà)出一個(gè)和給定圓面積相等的正方形。

尚克斯的計(jì)算結(jié)果公布后沒(méi)多久，德摩根（De Morgan）就發(fā)現(xiàn)了一個(gè)奇怪的統(tǒng)計(jì)現(xiàn)象：在這 707 位小數(shù)里，數(shù)字 7 出現(xiàn)的次數(shù)異常少。他在 1872 年的《悖論集錦》里提到了這一點(diǎn)，但這個(gè)疑問(wèn)一直懸到 1945 年才解開(kāi)——弗格森（Ferguson）發(fā)現(xiàn)尚克斯在第 528 位小數(shù)處算錯(cuò)了，后面的所有數(shù)字自然也都錯(cuò)了。1949 年，人們用計(jì)算機(jī)算出了 π 的 2000 位小數(shù)。在這次及之后所有計(jì)算機(jī)算出的結(jié)果里，數(shù)字 7 的出現(xiàn)次數(shù)都和預(yù)期差不多，而且到目前為止，π 的小數(shù)序列通過(guò)了所有隨機(jī)性統(tǒng)計(jì)測(cè)試。

π 符號(hào)的由來(lái)

我們得說(shuō)說(shuō) π 這個(gè)符號(hào)是怎么來(lái)的。1647 年，奧特雷德（Oughtred）用“ ”表示圓的直徑與周長(zhǎng)的比值；1697 年，大衛(wèi)·格雷戈里（David Gregory）用“ ”表示圓的周長(zhǎng)與半徑的比值。而第一個(gè)用 π 表示現(xiàn)在這個(gè)含義（周長(zhǎng)與直徑的比值）的，是威爾士數(shù)學(xué)家威廉·瓊斯（William Jones）——1706 年，他在著作中寫(xiě)道“3.14159 等等= ”。1737 年，歐拉（Euler）采用了這個(gè)符號(hào)，之后 π 很快成為標(biāo)準(zhǔn)表示法。

關(guān)于 π 的趣味歷史

最后，我們?cè)僬f(shuō)說(shuō) π 值計(jì)算過(guò)程中兩個(gè)有趣的統(tǒng)計(jì)故事，首先是布豐投針實(shí)驗(yàn)（Buffon's needle experiment）。如果有一組等距的平行直線（間距為 1），把一根長(zhǎng)度 的針隨機(jī)扔到這組直線上，針與直線相交的概率是 。有不少人試著用這種方法算 π，其中最出名的是拉澤里尼（Lazzerini）在 1901 年的實(shí)驗(yàn)——他扔了 34080 次針，算出

巧合的是，這個(gè)值正是祖沖之算出的結(jié)果。這個(gè)結(jié)果準(zhǔn)得有點(diǎn)可疑，而“34080 次”這個(gè)奇怪的次數(shù)也暴露了問(wèn)題?？系?tīng)枺↘endall）和莫蘭（Moran）指出，只要在“最佳時(shí)機(jī)”停止實(shí)驗(yàn)，就能得到一個(gè)理想值；但如果事先定好扔針次數(shù)，這種方法算 π 其實(shí)非常不準(zhǔn)。他們還調(diào)侃說(shuō)，不如找一塊圓形木頭，用卷尺量它的周長(zhǎng)和直徑，結(jié)果會(huì)更準(zhǔn)。

說(shuō)到“可疑實(shí)驗(yàn)”，格里奇曼（Gridgeman）在一篇嘲諷拉澤里尼等人的論文里，還搞了個(gè)有趣的惡作劇：他選了一根長(zhǎng)度 的針，扔了兩次，有一次和直線相交，然后根據(jù)公式計(jì)算：

算出 ，這個(gè)結(jié)果看起來(lái)還挺可信——當(dāng)然，他根本沒(méi)當(dāng)真。

π 引發(fā)的爭(zhēng)議與鬧劇

1934 年，居然有人以 π 的定義為借口，對(duì)著名數(shù)學(xué)家埃德蒙·蘭道（Edmund Landau）發(fā)起了種族攻擊，這事聽(tīng)起來(lái)簡(jiǎn)直難以置信。那一年，蘭道在哥廷根出版的教材里，用了現(xiàn)在相當(dāng)常見(jiàn)的方法定義 π： 是 1 到 2 之間使 的 值。這個(gè)定義引發(fā)了一場(chǎng)學(xué)術(shù)爭(zhēng)議，最終導(dǎo)致蘭道被免去哥廷根大學(xué)的教授職位。

比貝爾巴赫（Bieberbach）雖然是一位知名數(shù)論學(xué)家，但本人的種族主義觀點(diǎn)讓其聲名狼藉，他試圖這樣解釋免職蘭道的原因：
“哥廷根學(xué)生群體勇敢地反對(duì)偉大數(shù)學(xué)家埃德蒙·蘭道，歸根結(jié)底是因?yàn)檫@個(gè)人在研究和教學(xué)中‘非德國(guó)式’的風(fēng)格，讓德國(guó)人無(wú)法忍受。一個(gè)民族既然意識(shí)到另一個(gè)種族正試圖強(qiáng)加外來(lái)思想，就必須拒絕接受異文化的教師?！?/p>

G·H·哈代（G·H·Hardy）立刻在一篇公開(kāi)評(píng)論中反駁比貝爾巴赫，針對(duì)這種“非德國(guó)式 π 定義”的說(shuō)法：
“我們中有很多人，無(wú)論是英國(guó)人還是德國(guó)人，在戰(zhàn)爭(zhēng)期間都說(shuō)過(guò)一些言不由衷、現(xiàn)在回想起來(lái)會(huì)后悔的話。為了自己的地位焦慮，害怕跟不上愚蠢的潮流，不顧一切想不被落下，這些理由或許情有可原，哪怕不算英勇。但比貝爾巴赫教授的聲望，讓這些理由無(wú)法解釋他的言論——我只能得出一個(gè)不太善意的結(jié)論：他是真的相信這些話?！?/p>

不光德國(guó)因?yàn)?π 出了問(wèn)題，美國(guó)也曾因?yàn)?π 的取值引發(fā)過(guò)激烈的政治爭(zhēng)論。1897 年，印第安納州眾議院全票通過(guò)了一項(xiàng)法案，宣稱要確立一個(gè)“新的數(shù)學(xué)真理”：
“印第安納州議會(huì)頒布法令：經(jīng)發(fā)現(xiàn)，圓的面積與‘等于圓周長(zhǎng)四分之一的線段所構(gòu)成的正方形’的面積之比，等于矩形（這里實(shí)指正方形）的面積與其一邊構(gòu)成的正方形的面積之比?！保?897 年《印第安納州眾議院第 246 號(hào)法案》第一節(jié)）

還好印第安納州參議院更理智，把這項(xiàng)法案無(wú)限期擱置了！

關(guān)于 π 的開(kāi)放問(wèn)題

1. 數(shù)字 0-9 在 π 的小數(shù)展開(kāi)中，是否每個(gè)都出現(xiàn)無(wú)窮多次？

2. 布勞威爾的問(wèn)題：π 的小數(shù)展開(kāi)中，是否存在某一處連續(xù)出現(xiàn) 1000 個(gè) 0？

3. π 在十進(jìn)制下是“簡(jiǎn)單正規(guī)數(shù)”嗎？也就是說(shuō)，長(zhǎng)期來(lái)看，每個(gè)數(shù)字出現(xiàn)的頻率是否相等？

4. π 在十進(jìn)制下是“正規(guī)數(shù)”嗎？也就是說(shuō)，長(zhǎng)期來(lái)看，任意長(zhǎng)度的數(shù)字塊出現(xiàn)的頻率是否相等？

5. π 是“絕對(duì)正規(guī)數(shù)”嗎？也就是說(shuō)，長(zhǎng)期來(lái)看，在任意進(jìn)制下，任意長(zhǎng)度的數(shù)字塊出現(xiàn)的頻率是否相等？這個(gè)概念是博雷爾（Borel）在 1909 年提出的。

6. 另一個(gè)關(guān)于正規(guī)性的問(wèn)題：我們知道 π 不是有理數(shù)，所以它的小數(shù)不會(huì)從某一處開(kāi)始循環(huán)。但如果 π 是正規(guī)數(shù)，那么“314159265358979…”這前 100 萬(wàn)位數(shù)字，總會(huì)在某個(gè)位置再次出現(xiàn)；即便 π 不是正規(guī)數(shù)，這種情況也可能發(fā)生。事實(shí)果真如此嗎？如果是，會(huì)從哪一位開(kāi)始？（注：截至 2000 萬(wàn)位小數(shù)，出現(xiàn)過(guò)的最長(zhǎng)匹配是“31415926”，出現(xiàn)過(guò)兩次）

最后給大家一個(gè)記憶 π 小數(shù)展開(kāi)的口訣——每個(gè)單詞的字母數(shù)，對(duì)應(yīng) π 的一位小數(shù)：
“How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics. All of thy geometry, Herr Planck, is fairly hard...”（中文大意：聽(tīng)了量子力學(xué)的沉重講座后，我多想來(lái)杯酒，當(dāng)然是含酒精的。普朗克先生，你的幾何課實(shí)在太難了……）

對(duì)應(yīng)的 π 值是：

本文譯自 MacTutor 網(wǎng)站，作者 J J O'Connor and E F Robertson。原文除另有標(biāo)注外遵循 CC BY-SA 4.0 國(guó)際許可協(xié)議。翻譯：【遇見(jiàn)數(shù)學(xué)】，譯文繼承原協(xié)議：可自由復(fù)制、改編，但需標(biāo)注原文來(lái)源、作者及本譯文；改編后分享需采用同類許可。

往期推薦閱讀

1、

2、

3、

4、

5、

微信公眾號(hào)：考研競(jìng)賽數(shù)學(xué)(ID: xwmath)大學(xué)數(shù)學(xué)公共基礎(chǔ)課程分享交流平臺(tái)！支持咱號(hào)請(qǐng)點(diǎn)贊分享！