計算圓周率一直是令人著迷的議題,從古埃及至今,無數專家學者乃至業餘數學家前仆後繼地投入,1777 年,Georges-Louis Leclerc,Comte de Buffon (1707—1788) 提出嶄新的途徑,將圓周率這等幾何問題出發的計算,巧妙地以機率統計原理來表示,自此,開創使用隨機數值處理典型數學表示的先河,我們就來看看傳奇性的 Buffon 投針試驗。
法國數學家、科學家、《自然史》作者,也是風格家的 Comte de Buffon 在 1777 年某日,邀請賓客齊聚大廳,共襄盛舉一次試驗活動。古稀之年的 Buffon 鋪好一張白紙,其上預先畫好了一條條等距的平行線,接著取出一大把質量均等、長度為平行線間距一半的小針,待賓客就座後,Buffon 發言道:
「煩請各位將這些小針一根一根扔往白紙上,並且告知扔下的針是否與紙上平行線相交」
客隨主意,雖摸不著頭緒,但也一個個加入了試驗的行列。一把小針扔完了,把它撿起來又扔,而 Buffon 則在一旁不停地記數著,忙碌了將近一個鐘頭。最後,Buffon 高聲宣佈:
「各位賓客,依據我的紀錄,剛才的投針結果,共投針 2212 次,其中與平行線相交有 704 次。而總數 2212 與相交數 704 的比值為 3.142。」
說到這裡,Buffon 故作停頓,神秘張望賓客,接著說:
「這就是圓周率π的近似值!」
Buffon 利用平凡不過的除法,計算出圓周率的近似值,並宣稱投針的數目越多,圓周率的近似值將會越精準,這就是數學史上著名的 Buffon 投針問題,記載於其著作《機率算術試驗》(1777 年),此外,Buffon 也試著將機率應用於審判場合,比方說,若能對每個審判員規定某個足以理解真相或說出真相機會的數值,即可算出法庭作出正確判決的機會,換言之,就是「審判的概率」(Probabilite des jugements)。
圓周率π在這種看似雜亂的場合出現,實在出乎意料。一個直觀的理解途徑可透過物理上的對應,取一根鐵絲,將其彎成一個圓圈,適度剪裁使其直徑恰等於平行線間距離 d。於是乎,對於這個圓圈來說,無論如何扔下,都將和平行線有兩個交點。也就是說,若圓圈扔下的次數為 n,那麼,相交的交點總數必為 2n。接著,我們展開物理的形變,將圓圈拉開、拉直,這樣就成為長度為 πd 的鐵絲,再將這條鐵絲扔下,與平行線相交的情況就複雜許多,由於 1 < π <>πd : (1/2)d ≈ 2n : m
這也是 Buffon 投針試驗中所作的參數配置,約分後可得漂亮的式子:
π ≈ n / m
在古典數學中,求圓周率之值是幾何問題,而 Buffon 卻以此拍案叫絕的方式,以機率方法打通兩個看似風馬牛不相及的領域,成為幾何概率的典型例子。
近代科學的發展下,原本壁壘分明的個別人文、科學、哲學思想領域走向空前的大融通,匯流而成當代種種巨大變革,一如 Buffon 首次打破機率與幾何學的藩籬。數學領域的變遷也受到這等啟蒙,1904 年,R·Chartres 甚至提出另一種表示法:若寫下任意兩個整數,測這兩者互質的機率為 6 / π^2。
經過幾百年的演繹與探討,Buffon 投針試驗逐漸演化為一種數值方法的前身:「蒙地卡羅方法」(Monte Carlo method),也就是透過利用亂數取樣 (random sampling) 模擬來解決數學問題。第二次世界大戰期間,Monte Carlo 方法被系統性地應用於科學研究中,誕生了 MANIAC (Mathematical Analyzer, Numerical Integrator and Computer),而 Stanislaw Ulam、John von Neumann、Nicholas Metropolis、Enrico Fermi 等人發展法一種基於樣本統計的方法,來解決關於在原子彈設計中,中子隨機擴散問題和 Schrodinger 等式的特徵值估計問題。該方法的原理最初是 Stanislaw Ulam 闡述的,後來由 John von Neumann 深入研究,於 1949 年發表一篇名為 "The Monte Carlo method" 的論文而聞名,當然,到了進入電腦時代,這個方法才得以由原本手動產生亂數來解決問題,變成實際性的數值方法。
Monte Carlo 方法是由 Nicholas Metropolis 所命名,取自其亂數機率有如賭博一般,而恰似北非最西側的摩洛哥首都 Monte Carlo,也就是知名賭城,種種奇豔動人的賭場生活寫照。所有具有隨機效應的過程,均可能以 Monte Carlo 方法大量模擬單一事件,並藉由統計上平均值,獲得某設定條件下實際最可能測量值,更廣泛來說,自然界裏的布朗運動、電波的噪音、基因的突變、交通即時路況等等,無處不含有隨機的變化,均有可適用的場合。
參考資料:
Wikipedia:
布豐投針問題
Read more...