編者按:
柯普寧說,當(dāng)數(shù)學(xué)家導(dǎo)出方程式和公式,如同看到雕像、美麗的風(fēng)景,聽到優(yōu)美的曲調(diào)等等一樣而得到充分的快樂。南科大首屆數(shù)學(xué)文化節(jié)正在進(jìn)行中,為傳播數(shù)學(xué)文化,普及數(shù)學(xué)知識(shí),展示數(shù)學(xué)之美,打造學(xué)科與公益相融合的校園公益文化,基金會(huì)聯(lián)合數(shù)學(xué)系推出系列數(shù)學(xué)公益科普文章,以饗讀者。
圓周率的前世今生
撰文 | 付云皓
古代數(shù)學(xué)
——圓周率的認(rèn)識(shí)與估算
關(guān)于圓周率的定義有兩種,一種是周長(zhǎng)的定義,即
或者,
一種是面積的定義,即
從古希臘、古中國(guó)、古印度等古代文明留下的文獻(xiàn)可以推斷出,這些古代文明都認(rèn)為兩者是同樣的定義,即可以認(rèn)為它們都知道這樣一個(gè)公式:
在古代文明中,隨著對(duì)數(shù)學(xué)研究的深入,出現(xiàn)了要求得到圓周率的精確值或近似值的問題。如古希臘的三大幾何作圖問題中,圓化方問題是用尺規(guī)作圖作出一個(gè)與已知圓面積相等的正方形,實(shí)際上就是求的精確值的問題。
最早研究圓化方問題的是安納薩格拉斯(約前500-前428),公元前5世紀(jì)下半葉,開奧斯的希波克拉底解決了化圓為方(最初只有等腰直角三角形的圖形,通過勾股定理/畢達(dá)哥拉斯定理可以推廣到任意兩個(gè)月牙的情況)。
化圓為方
如圖,由勾股定理
故以AB為直徑的半圓面積等于分別以AC、CB為直徑的兩個(gè)半圓的面積之和。去掉公共部分后,即得到兩塊綠色的月牙型圖形面積之和等于三角形ABC的面積。
詭辯學(xué)派的代表人物安提豐(約公元前480-前411)首先提出了用圓內(nèi)接正多邊形的面積來逼近圓面積的方法,用正4,8,16,……邊形來逼近圓,而正多邊形化方是可以通過尺規(guī)作圖來實(shí)現(xiàn)的。雖然這實(shí)際上并沒有解決圓化方問題,但是他使用的這一“窮竭法”后來被阿基米德(公元前287-前212)等人發(fā)揚(yáng)光大,解決了求球體體積等一系列問題。安提豐也因此成為窮竭法的始祖。
考慮一個(gè)半徑為r的圓,我們需要找到一個(gè)與這個(gè)圓面積相同的矩形,即需要找到a使得
這里可以把a(bǔ)看成r和πr的比例中項(xiàng),當(dāng)然也可以直接想法求π的開平方。事實(shí)上,古希臘人早就掌握了求比例中項(xiàng)(或開平方)的幾何做法,具體方法是考慮直角三角形即可。
如下圖所示
用幾何法求π的開平方值
如果AD=πr,BD=r,先作AB中點(diǎn)O,然后以O(shè)為圓心,OA為半徑作圓,與過D且垂直于AB的直線交于點(diǎn)C,由射影定理即得
如果π是一個(gè)有理數(shù),那么能夠通過作圖做出πr,那么得到a也是沒有問題的了。即便π是無理數(shù),假如π形如
或者
也是可以作出的。不過,1882年德國(guó)數(shù)學(xué)家林德爾曼證明了π是超越數(shù)(即不是任何整系數(shù)代數(shù)方程的根),而另一方面,前人已經(jīng)證明尺規(guī)作圖能得到的所有長(zhǎng)度之間的比值都必須是代數(shù)數(shù)(即是某個(gè)整系數(shù)代數(shù)方程的根),這樣就對(duì)圓化方問題給出了否定的答案。
古印度的《繩法經(jīng)》中,遇到將圓轉(zhuǎn)化為等面積的正方形時(shí),使用的數(shù)據(jù)是正方形的邊長(zhǎng)為圓直徑的8/9或者
即認(rèn)為
或者
在中國(guó)古代,從《九章算術(shù)》開始就有了對(duì)圓周率的近似值。
術(shù)曰:半周半徑相乘得積步。
又術(shù)曰:周徑相乘,四而一。
又術(shù)曰:徑自相乘,三之,四而一。
又術(shù)曰:周自相乘,十二而一。
劉徽注《九章算術(shù)》寫到“方五斜七,圓三徑一”,即《九章算術(shù)》里的圓周率近似值為3。從西漢末年開始,新率陸續(xù)出現(xiàn),但仍然不精確,且沒有推算方法。
劉徽同樣以圓內(nèi)接正多邊形的面積來近似圓的面積,但與安提豐的思路不同,他思考的是求出圓面積的盡量精確的值。
若半徑為r的圓內(nèi)接正n邊形的邊長(zhǎng)是ln,那么
證明:
如下圖所示
劉徽的方法
設(shè)AB是以O(shè)為圓心的圓內(nèi)接正n邊形的一條邊,M是AB中點(diǎn),延長(zhǎng)OM交圓于C,則AC,CB都是圓內(nèi)接正2n邊形的一條邊。因此
故
設(shè)該圓的面積為S0,內(nèi)接正n邊形的面積為Sn,內(nèi)接正2n邊形的面積為S2n。由于
且
因此只要求出圓內(nèi)接正n邊形的邊長(zhǎng),就能輕松求出正2n邊形的面積。
正常情況下,要估計(jì)S0,除了求內(nèi)接正2n邊形的面積,還要求外切正2n邊形的面積,但是劉徽注意到
即
這是因?yàn)槿绻麅H看扇形OAB的話,S2n比Sn多的是三角形ACB的面積,而S0比S2n多的是兩個(gè)弓形的面積,兩個(gè)弓形的面積要小于三角形ACB的面積。這樣,劉徽就不需要求外切多邊形的面積了。
劉徽從正6邊形開始算到正192邊形,得出“徽率”3.14。
后來祖沖之沿用這個(gè)算法算到24576邊形(即12288邊形邊長(zhǎng)),得到
(朒nv4數(shù))3.1415926<π<3.1415927(盈數(shù))
按現(xiàn)在的角度來看,圓內(nèi)接正n邊形的邊長(zhǎng)即為
上面的等式恰好是一個(gè)三角恒等式。
在半角公式中,余弦的半角更加簡(jiǎn)單,即
用此迭代可以得到很小的角的余弦值,例如
然后計(jì)算正弦值也可以得到同樣的估計(jì)結(jié)果。
近代數(shù)學(xué)
——圓周率的逼近方式
從微積分的萌芽階段開始,由于三角函數(shù)的運(yùn)算更加靈活,數(shù)學(xué)家們得到了更多關(guān)于圓周率的逼近方式。
首先是Wallis公式。由于
故
因此
推論:
也就是說
這個(gè)結(jié)論在概率和數(shù)論中都有一些應(yīng)用。
Wallis公式的壞處是收斂得太慢了,可以看到不等式的左右兩邊的比例是2n:(2n+1),因此即便n取n=500,誤差也在千分之一左右。
由arctanx的展開可以得到如下的Leibniz級(jí)數(shù):
這個(gè)級(jí)數(shù)也有同樣的問題,收斂得太慢了。例如,寫到1/1001那一項(xiàng),誤差仍然是千分之一。
由上可知
故可以寫成級(jí)數(shù)的形式,并且收斂速度很快。因此如果能先把π寫成較小的x的反正切,收斂速度就會(huì)變快很多。馬青公式(Machin’s Formula)就是在此基礎(chǔ)上得到的一個(gè)收斂速度較快的公式。
設(shè)
那么
故
容易知道4α,β+(π/4)均在(0,π)中,因此有4α=β+(π/4),故
這就是馬青公式,再根據(jù)
可以得到較快收斂到圓周率的級(jí)數(shù)。
一些案例
——圓周率出現(xiàn)在其他領(lǐng)域
除了普通的級(jí)數(shù)之外,圓周率也在其它的領(lǐng)域里有出現(xiàn),這里舉兩個(gè)例子。
第一個(gè)例子是自然數(shù)平方和的倒數(shù)求和。歐拉求自然數(shù)平方和的倒數(shù)時(shí),使用的是下面的方法。
考慮
即
它的根為±π,±2π,±3π......
因此它應(yīng)與另一個(gè)無限乘積
的展開應(yīng)該相等。比較兩邊項(xiàng)x平方項(xiàng)的系數(shù)即得
故
當(dāng)然,上面的證明是不嚴(yán)謹(jǐn)?shù)?,使用Fourier級(jí)數(shù)可以給出一個(gè)較為嚴(yán)謹(jǐn)?shù)淖C明,這里略。
自然數(shù)平方和的倒數(shù)有如下的應(yīng)用。由于整數(shù)環(huán)上有唯一分解,也就是說任何一個(gè)正整數(shù)可以以唯一的方式分解成素?cái)?shù)的乘積(不計(jì)順序),因此如果假設(shè)p1,p2,...,pm,...是所有的素?cái)?shù),那么
另一方面,由級(jí)數(shù)求和公式易知
因此
故
現(xiàn)在,我們選擇一個(gè)非常大的正整數(shù)N,然后在1,2,...,N中隨機(jī)選擇兩個(gè)數(shù)a,b(允許相同),那么a,b互素的概率是多少呢?
注意a,b互素當(dāng)且僅當(dāng)它們不同時(shí)是任何一個(gè)素?cái)?shù)p的倍數(shù)。當(dāng)N相對(duì)于pm非常大時(shí),可以認(rèn)為a,b是pm倍數(shù)的概率均為1/pm,因此a,b不同時(shí)是pm倍數(shù)的概率為
又由于N當(dāng)相對(duì)于pm,p1非常大時(shí),1,2,...,N中隨機(jī)選擇一個(gè)數(shù),是否為pm的倍數(shù)與是否為p1的倍數(shù)是相互獨(dú)立事件。因此,當(dāng)N趨近于無窮時(shí),a,b互素的概率是
也就是在足夠大的范圍內(nèi)任取兩個(gè)自然數(shù),它們互素的概率趨近于6/(π^2)。
第二個(gè)例子是Buffon投針實(shí)驗(yàn)??紤]一組平行線,兩兩間距為d,如果將直徑為d的圓形鐵絲扔到平行線上,不論何種情況,圓都應(yīng)該與這組平行線恰有2個(gè)交點(diǎn)。
現(xiàn)在我們將鐵絲拉直再隨機(jī)拋擲到平行線上,直觀來看,鐵絲的“長(zhǎng)度”沒有變化,故交點(diǎn)的期望數(shù)不變。但是,現(xiàn)在鐵絲的長(zhǎng)度是πd,含有圓周率π,就可以利用它。由于現(xiàn)在鐵絲的長(zhǎng)度超過d,記錄交點(diǎn)的數(shù)目不容易。
我們可以將鐵絲變短一些,假設(shè)鐵絲的長(zhǎng)度是l<d,那么交點(diǎn)的期望值應(yīng)當(dāng)是
且交點(diǎn)至多有一個(gè),故鐵絲與平行線相交的概率為2l/πd。
由此看出,如果令l=d/2,那么鐵絲與平行線相交的概率恰為1/π。
當(dāng)然,上述的推理也是不嚴(yán)謹(jǐn)?shù)?,我們可以用積分的方法再推理一遍。
假設(shè)鐵絲與平行線的夾角為φ,那么鐵絲在垂直于平行線的方向的投影長(zhǎng)度為lsinφ。當(dāng)l<d時(shí),鐵絲與平行線相交的概率為lsinφ/d。
由于鐵絲是隨機(jī)拋擲到平行線上的,故可以認(rèn)為φ是均勻分布在【0,π/2】上,故鐵絲與平行線相交的概率為
作者簡(jiǎn)介
付云皓,南方科技大學(xué)數(shù)學(xué)系講師,主要研究方向?yàn)?span>圖論、教育數(shù)學(xué)。
親愛的小伙伴
你學(xué)會(huì)了嗎?
聯(lián)系客服