眾所周知,圓周率是圓的周長與直徑的比值。在古代的數(shù)學(xué)史上,圓周率的研究和計算一定程度上反應(yīng)了當(dāng)時的數(shù)學(xué)水平。古希臘阿基米德,阿拉伯的卡西,古印度阿耶波多,古代中國祖沖之和劉徽等等數(shù)學(xué)家,都致力于圓周率的研究和計算,先后給出了圓周率的估值。劉徽等人使用的是割圓術(shù):使用內(nèi)接于圓的正多邊形逼近圓,多邊形的邊數(shù)越多,其周長與面積也越接近圓。思路很簡單,但其計算量是個不小的挑戰(zhàn)。也似乎在計算前,缺少了對計算的論證。不同于古代正多邊形的估算,現(xiàn)代借助計算機,已經(jīng)可以求得的圓周率達(dá)到了驚人的位數(shù),2019年3月14日,谷歌宣布圓周率已到小數(shù)點后31.4萬億位。那么現(xiàn)代計算的依據(jù)是什么呢?

圓周率是無理數(shù)
1737年,歐拉證明了e是無理數(shù);蘭伯特根據(jù)歐拉的工作證明:如果x是非零有理數(shù),那么和
都不是有理數(shù)。再根據(jù)此結(jié)果,由
,可得到
不是有理數(shù),因此圓周率是無理數(shù)。
圓周率是超越數(shù)
勒讓德首先猜測圓周率可能不是有理系數(shù)方程的根,勒讓德的猜測促進了無理數(shù)的分類。任何有理系數(shù)多項式代數(shù)方程的任何一個根叫做一個代數(shù)數(shù)。即方程的根叫做代數(shù)數(shù),其中
是有理數(shù)。而那些不是代數(shù)數(shù)的數(shù)叫做超越數(shù)。
1873年,法國數(shù)學(xué)家埃爾米特給出了自然常數(shù)的超越性的證明,1882年,德國數(shù)學(xué)家林德曼證明了
也是超越數(shù)。但是
是不是無理數(shù)、超越數(shù)還有待確定,歐拉常數(shù)
是不是無理數(shù)也同樣有待確定。
圓周率的計算公式
- 維加于1789年發(fā)表公式:
- 歐拉-維加公式:
- 歐拉公式:
- 克拉森公式:
- 盧瑟福公式(1841):
- 高斯公式:
- 達(dá)澤公式:
- 布賽爾公式:
- 艾斯克托公式(1896):
- 肖魯茲公式(1844):
- 斯特姆公式:
- 山克斯公式(1853):
- 赫頓:
- 馬庭(1706):
上述公式只是部分成果,而且除了正切反函數(shù)表示的公式外,也有一些使用正弦、余弦的反函數(shù)的圓周率公式。再使用正切反函數(shù)的級數(shù)展開式,通過計算機來計算圓周率。由于上述公式的收斂速度有快有慢,圓周率的計算會選擇收斂較快的公式。因此歐拉和馬庭公式使用的更多一些。
圓周率的級數(shù)表達(dá)
- 萊布尼茨(1673):
- 歐拉:
- 歐拉:
- 歐力斯(1656):
- 松永良弼:
- 牛頓:
- 牛頓:
- 烏衣塔:
- 夏普(1705):
- 馬庭(1706):
1673年,萊布尼茨找到圓周率的一個級數(shù)表示形式,結(jié)束了人工計算圓周率的工作,開啟了計算機計算圓周率的時代。當(dāng)然了,萊布尼茨的級數(shù)收斂速度是偏慢的。借助計算機,圓周率的計算位數(shù)也在不斷的刷新,很多時候都用來測試計算機的性能。
想要了解更多數(shù)學(xué)內(nèi)容,請關(guān)注公眾號:“究盡數(shù)學(xué)” 和 “究盡中學(xué)數(shù)學(xué)”