圓周率 π 的歷史溯源
圣經里有一段不太為人熟知的經文,是這么寫的:
“他又鑄了一個銅海,樣式是圓的,從這邊到那邊共十肘,高五肘,圍腰三十肘。”(《列王紀上》7:23)
這段內容在《歷代志下》4:2 中也能找到,記載的是公元前 950 年左右建造所羅門圣殿時的計數規格。這里提到的圓周率 π 值是 3,顯然不算精確,即便放在當時也不算——要知道更早之前,埃及人和巴比倫人就已經算出了更接近真實值的結果:埃及人的 π 約為 ,巴比倫人的則是 。不過這里也要替所羅門的工匠們說句公道話,他們描述的“銅海”本就是某種巨大的銅器,這種大物件既不可能、也沒必要追求極高的幾何精度。
早期 π 值:從測量到理論![]()
“圓的周長與直徑的比值是恒定的一個數”這個事實,人們早就知道了,具體起源已無從考證。最早的 π 值(包括圣經里的 3),幾乎都是通過測量得到的。比如公元前 1650 年左右的埃及《萊因德數學紙草書》(Rhind Papyrus)中,就有充分證據表明當時人們用 作為 π 的近似值。
而第一個通過理論計算得到 π 值的,大概是敘拉古的阿基米德(Archimedes of Syracuse,公元前 287-前 212 年)。他算出的結果是:
在簡單說說他的證明思路前,有個點值得注意:這里用到的不等式,其實暗含了相當高的數學技巧。阿基米德很清楚——直到現在還有很多人不清楚——π 并不等于 ,他也從沒宣稱自己找到了 π 的精確值。如果我們取他兩個邊界值的平均值,得到的結果是 3.1418,誤差僅約 0.0002。
阿基米德的計算思路![]()
假設有一個半徑為 1 的圓,在圓內作一個邊數為 的正內接多邊形,其半周長記為 ;再作一個同樣邊數的正外切多邊形,半周長記為 。(下面是 時的示意圖,即邊數為 12 的多邊形)
![]()
這種做法會定義出兩個序列:一個是遞增序列 ,另一個是遞減序列 ,而這兩個序列的極限都是 π。
用三角函數表示的話,這兩個半周長可以寫成:
,
其中 。同樣,下一個邊數( )的多邊形半周長滿足:
,
通過基礎三角函數運算,能推導出兩個關鍵關系式:
阿基米德從初始值 、 開始,先用公式(1)算出 ,再用公式(2)算出 ,接著又用(1)算 ,用(2)算 ,一直算到 和 。最終他得出結論: 。
這里要明白一點:我們剛才用三角函數描述,其實不符合歷史——阿基米德當時沒有代數和三角函數符號,只能純靠幾何方法推導(1)和(2);并且他連我們現在用的十進制記數法都沒有,單靠公式(1)(2)計算 和 ,本身就不是件容易事。
咱們不該覺得 “他怎么咋算到 96 邊形就停了”,反而該驚嘆“居然能硬生生算到 96 邊形”—— 要知道在沒有任何現代數學工具的情況下,這已經是超了不起的成就了。
歷代數學家的 π 值探索
理論上,只要不斷增加多邊形邊數,就能算出更精確的 π 值。后來確實有不少人這么做,比如:
托勒密(約公元 150 年):3.1416
祖沖之(430-501 年): (約 3.1415929)
花拉子米(約公元 800 年):3.1416
卡西(約 1430 年):精確到 14 位小數
韋達(1540-1603 年):9 位小數
魯曼(1561-1615 年):17 位小數
范·塞倫(約 1600 年):35 位小數
除了祖沖之,這些改進都沒有理論上的突破,靠的只是更強的計算耐力。有意思的是,就像其他科學領域一樣,在公元 400 到 1400 年這一千年里,π 值計算的領先地位是從歐洲轉移到了東方。
花拉子米生活在巴格達,順便說一句,“算法(algorithm)”這個詞就源自他的名字,而他某本著作標題里的“代數學(al jabr)”,則成了“代數(algebra)”一詞的來源。卡西生活的地方更靠東,在撒馬爾罕。
![]()
文藝復興后的 π 值公式革命
歐洲文藝復興最終帶來了全新的數學世界,其中一個早期成果就是出現了計算 π 的數學公式。最早的公式之一是沃利斯(1616-1703 年)提出的:
而最知名的公式之一,是這個無窮級數:
這個級數有時會歸功于萊布尼茨(1646-1716 年),但實際上最早是詹姆斯·格雷戈里(1638-1675 年)發現的。
這些公式既驚艷又出人意料——右邊都是純算術形式,而 π 最初是從幾何里來的。它們向我們展示了無窮過程能帶來多么驚人的結果,也為現代數學的豐富內涵指明了方向。
不過從計算 π 值的角度看,這兩個公式其實沒什么用。比如格雷戈里級數,要想得到 4 位精確小數,誤差得小于 0.00005,這意味著需要約 10000 項才能做到。但格雷戈里還證明了一個更通用的結果:
只要令 ,就能得到前面那個 π 的級數。但如果利用 ,就能得到一個收斂快得多的級數:
這個級數的第 10 項是 ,小于 0.00005,也就是說,只需要 9 項就能得到至少 4 位精確小數。
還有個更好的辦法,是利用這個公式:
然后把 和 分別代入公式(3),計算兩個級數的和就行。
顯然,如果能找到這樣的公式:
其中 和 都是大數,收斂速度會快得多。1706 年,馬青(Machin)就找到了這樣一個公式:
其實只要會證明公式(4),證明(5)也沒什么額外難度,就是計算過程麻煩點。當然,能想出這個公式本身,就是另一回事了。
計算 π 的“執念”與錯誤
有了這樣的公式,計算 π 的唯一難點就只剩過程的枯燥了。不用說,還真有人愿意花大量時間和精力在這種單調又無用的追求上。其中一位是英國的尚克斯(Shanks),他用梅欽公式算出了 π 的 707 位小數,并在 1873 年發表了自己多年的計算成果。尚克斯能留名,原因還挺特別,我們馬上就會說到。
先看看 π 值計算精度的提升歷程:
1699 年:夏普用格雷戈里的成果算出 71 位精確小數
1701 年:梅欽用改進方法算出 100 位,后來的人也沿用他的思路:
1719 年:德拉尼算出 112 位精確小數
1789 年:維加算出 126 位,1794 年又算到 136 位
1841 年:盧瑟福算出 152 位,1853 年算到 440 位
1873 年:尚克斯算出 707 位,但其中只有前 527 位是對的
尚克斯知道 π 是無理數(irrational number),因為蘭伯特(Lambert)早在 1761 年就證明了這一點。就在尚克斯完成計算后不久,林德曼(Lindemann)證明了 π 是超越數(transcendental number)——也就是說,π 不是任何整系數多項式方程的解。其實林德曼的這個結果,也證明了“化圓為方”是不可能的:既然 π 是超越數,就不可能用尺規作圖畫出一個和給定圓面積相等的正方形。
尚克斯的計算結果公布后沒多久,德摩根(De Morgan)就發現了一個奇怪的統計現象:在這 707 位小數里,數字 7 出現的次數異常少。他在 1872 年的《悖論集錦》里提到了這一點,但這個疑問一直懸到 1945 年才解開——弗格森(Ferguson)發現尚克斯在第 528 位小數處算錯了,后面的所有數字自然也都錯了。1949 年,人們用計算機算出了 π 的 2000 位小數。在這次及之后所有計算機算出的結果里,數字 7 的出現次數都和預期差不多,而且到目前為止,π 的小數序列通過了所有隨機性統計測試。
π 符號的由來
我們得說說 π 這個符號是怎么來的。1647 年,奧特雷德(Oughtred)用“ ”表示圓的直徑與周長的比值;1697 年,大衛·格雷戈里(David Gregory)用“ ”表示圓的周長與半徑的比值。而第一個用 π 表示現在這個含義(周長與直徑的比值)的,是威爾士數學家威廉·瓊斯(William Jones)——1706 年,他在著作中寫道“3.14159 等等= ”。1737 年,歐拉(Euler)采用了這個符號,之后 π 很快成為標準表示法。
關于 π 的趣味歷史![]()
最后,我們再說說 π 值計算過程中兩個有趣的統計故事,首先是布豐投針實驗(Buffon's needle experiment)。如果有一組等距的平行直線(間距為 1),把一根長度 的針隨機扔到這組直線上,針與直線相交的概率是 。有不少人試著用這種方法算 π,其中最出名的是拉澤里尼(Lazzerini)在 1901 年的實驗——他扔了 34080 次針,算出
巧合的是,這個值正是祖沖之算出的結果。這個結果準得有點可疑,而“34080 次”這個奇怪的次數也暴露了問題。肯德爾(Kendall)和莫蘭(Moran)指出,只要在“最佳時機”停止實驗,就能得到一個理想值;但如果事先定好扔針次數,這種方法算 π 其實非常不準。他們還調侃說,不如找一塊圓形木頭,用卷尺量它的周長和直徑,結果會更準。
說到“可疑實驗”,格里奇曼(Gridgeman)在一篇嘲諷拉澤里尼等人的論文里,還搞了個有趣的惡作劇:他選了一根長度 的針,扔了兩次,有一次和直線相交,然后根據公式計算:
算出 ,這個結果看起來還挺可信——當然,他根本沒當真。
π 引發的爭議與鬧劇
1934 年,居然有人以 π 的定義為借口,對著名數學家埃德蒙·蘭道(Edmund Landau)發起了種族攻擊,這事聽起來簡直難以置信。那一年,蘭道在哥廷根出版的教材里,用了現在相當常見的方法定義 π: 是 1 到 2 之間使 的 值。這個定義引發了一場學術爭議,最終導致蘭道被免去哥廷根大學的教授職位。
比貝爾巴赫(Bieberbach)雖然是一位知名數論學家,但本人的種族主義觀點讓其聲名狼藉,他試圖這樣解釋免職蘭道的原因:
“哥廷根學生群體勇敢地反對偉大數學家埃德蒙·蘭道,歸根結底是因為這個人在研究和教學中‘非德國式’的風格,讓德國人無法忍受。一個民族既然意識到另一個種族正試圖強加外來思想,就必須拒絕接受異文化的教師。”
G·H·哈代(G·H·Hardy)立刻在一篇公開評論中反駁比貝爾巴赫,針對這種“非德國式 π 定義”的說法:
“我們中有很多人,無論是英國人還是德國人,在戰爭期間都說過一些言不由衷、現在回想起來會后悔的話。為了自己的地位焦慮,害怕跟不上愚蠢的潮流,不顧一切想不被落下,這些理由或許情有可原,哪怕不算英勇。但比貝爾巴赫教授的聲望,讓這些理由無法解釋他的言論——我只能得出一個不太善意的結論:他是真的相信這些話。”
不光德國因為 π 出了問題,美國也曾因為 π 的取值引發過激烈的政治爭論。1897 年,印第安納州眾議院全票通過了一項法案,宣稱要確立一個“新的數學真理”:
“印第安納州議會頒布法令:經發現,圓的面積與‘等于圓周長四分之一的線段所構成的正方形’的面積之比,等于矩形(這里實指正方形)的面積與其一邊構成的正方形的面積之比。”(1897 年《印第安納州眾議院第 246 號法案》第一節)
還好印第安納州參議院更理智,把這項法案無限期擱置了!
關于 π 的開放問題
數字 0-9 在 π 的小數展開中,是否每個都出現無窮多次?
布勞威爾的問題:π 的小數展開中,是否存在某一處連續出現 1000 個 0?
π 在十進制下是“簡單正規數”嗎?也就是說,長期來看,每個數字出現的頻率是否相等?
π 在十進制下是“正規數”嗎?也就是說,長期來看,任意長度的數字塊出現的頻率是否相等?
π 是“絕對正規數”嗎?也就是說,長期來看,在任意進制下,任意長度的數字塊出現的頻率是否相等?這個概念是博雷爾(Borel)在 1909 年提出的。
另一個關于正規性的問題:我們知道 π 不是有理數,所以它的小數不會從某一處開始循環。但如果 π 是正規數,那么“314159265358979…”這前 100 萬位數字,總會在某個位置再次出現;即便 π 不是正規數,這種情況也可能發生。事實果真如此嗎?如果是,會從哪一位開始?(注:截至 2000 萬位小數,出現過的最長匹配是“31415926”,出現過兩次)
最后給大家一個記憶 π 小數展開的口訣——每個單詞的字母數,對應 π 的一位小數:
“How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics. All of thy geometry, Herr Planck, is fairly hard...”(中文大意:聽了量子力學的沉重講座后,我多想來杯酒,當然是含酒精的。普朗克先生,你的幾何課實在太難了……)
對應的 π 值是:
本文譯自 MacTutor 網站,作者 J J O'Connor and E F Robertson。原文除另有標注外遵循 CC BY-SA 4.0 國際許可協議。翻譯:【遇見數學】,譯文繼承原協議:可自由復制、改編,但需標注原文來源、作者及本譯文;改編后分享需采用同類許可。
往期推薦閱讀
1、
2、
3、
4、
5、
微信公眾號:考研競賽數學(ID: xwmath)大學數學公共基礎課程分享交流平臺!支持咱號請點贊分享!
![]()
特別聲明:以上內容(如有圖片或視頻亦包括在內)為自媒體平臺“網易號”用戶上傳并發布,本平臺僅提供信息存儲服務。
Notice: The content above (including the pictures and videos if any) is uploaded and posted by a user of NetEase Hao, which is a social media platform and only provides information storage services.