數(shù)學(xué)實(shí)驗之8.pi的近似計算

數(shù)學(xué)實(shí)驗之八——的近似計算MathematicalExperiments2021/6/27在本次試驗中，我們將追溯關(guān)于圓周率的計算歷程。通過對割圓術(shù)、韋達(dá)公式、級數(shù)加速法、迭代法等計算方法的介紹和計算體驗，感受數(shù)學(xué)思想和數(shù)學(xué)方法的發(fā)展過程，提高對極限和級數(shù)收斂性及收斂速度的綜合認(rèn)識，同時使我們看到數(shù)學(xué)家對科學(xué)真理的永無止境的追求。實(shí)驗?zāi)康闹黜撋弦豁撓乱豁?021/6/27主要內(nèi)容

四、利用級數(shù)計算二、韋達(dá)(VieTa)公式

三、數(shù)值積分方法

一、割圓術(shù)六、拉馬努金(Ramanujan)公式五、蒙特卡羅(MonteCarlo)法2021/6/27實(shí)驗指導(dǎo)

π是使人們最經(jīng)常使用的數(shù)學(xué)常數(shù)。人們對π的研究已經(jīng)持續(xù)了2500多年。在今天，這種探索還在繼續(xù)……2021/6/27世界上數(shù)學(xué)家們一致公認(rèn)：“歷史上一個國家計算圓周率的準(zhǔn)確度，可以作為衡量這個國家當(dāng)時數(shù)學(xué)水平的一個標(biāo)志?！睂?shí)驗指導(dǎo)2021/6/27π值——算法美的追求

π作為圓周率的符號，是由著名數(shù)學(xué)家歐勒于公元1737年首先使用的。古代的希伯來人，在描述所羅門廟宇中的“熔池”時曾經(jīng)這樣寫道：“池為圓形，對徑為十腕尺，池高為五腕尺，其周長為三十腕尺?！笨梢?，古希伯來人認(rèn)為圓周率等于3。不過，那時的建筑師們，似乎沒有人不明白，圓周長與直徑的比要比3大一些。公元前3世紀(jì)古希臘大數(shù)學(xué)家阿基米德求出了223/71<π

<22/7。2021/6/27“割圓術(shù)”中學(xué)問多

我國2000多年前的《周髀算經(jīng)》稱“周三徑一”，這是π的第一個近似值，叫做“古率”。據(jù)說，漢代大科學(xué)家、文學(xué)家張衡，有“圓周率一十之面”的推算。清代李潢考證這句話意思為π≈sqrt(10)。魏晉間劉徽由圓內(nèi)接正六邊形依次倍增到正192邊形，計算周長與值徑之比，得3.141024<π<3.142704實(shí)際應(yīng)用時取3.14，或分?jǐn)?shù)值157/50。2021/6/27“割圓術(shù)”中學(xué)問多

他的割圓術(shù)已含有無限逼近的極限思想，這是比求π值更可寶貴的。從方法上說，他得到了重要的“劉徽不等式”。

設(shè)圓內(nèi)接正n邊形的邊長為an，圓內(nèi)接正n邊形的面積為Sn。根據(jù)勾股定理，邊長有如下遞推公式：2021/6/27“割之彌細(xì)，失之彌少，割之又割，則與圓合體而無所失矣?！?/p>

面積與邊長有如下關(guān)系：

圓面積S與多邊形的面積Sn之間有如下關(guān)系：2021/6/27劉徽不等式借助于計算機(jī)來完成劉徽的工作：a(1)=sqrt(3);b(1)=3*sqrt(3)/2;fori=2:6a(i)=sqrt(2-sqrt(4-a(i-1)^2));b(i)=3*2^(i-2)*a(i);

c(i)=2*b(i)-b(i-1);endn=[3,6,12,24,48,96];size(b)result=[n;a;b]2021/6/27劉徽不等式result‘＝ans=3.00001.73212.598106.00001.00003.00003.401912.00000.51763.10583.211724.00000.26113.13263.159448.00000.13083.13943.146196.00000.06543.14103.14272021/6/27割圓術(shù)的意義劉徽創(chuàng)立的割圓術(shù)，其意義不僅在于計算出了Pi的近似值，而且還在于提供了一種研究數(shù)學(xué)的方法。這種方法相當(dāng)于今天的“求積分”，后者經(jīng)16世紀(jì)英國的牛頓和德國的萊布尼茨作系統(tǒng)總結(jié)而得名。鑒于劉徽的巨大貢獻(xiàn)，所以不少書上把他稱做“中國數(shù)學(xué)史上的牛頓”，并把他所創(chuàng)造的割圓術(shù)稱為“徽術(shù)”。2021/6/27韋達(dá)（VieTa）公式

1593年，韋達(dá)首次給出了計算Pi的精確表達(dá)式：韋達(dá)公式看起來有些神秘，其實(shí)它的導(dǎo)出過程所用的都是樸實(shí)簡潔的數(shù)學(xué)方法。2021/6/27韋達(dá)（VieTa）公式1、從sint開始2021/6/27韋達(dá)（VieTa）公式所以，對任意N，總有2021/6/27韋達(dá)（VieTa）公式2、從cos(pi/4)開始2021/6/27韋達(dá)（VieTa）公式3、使用VieTa公式計算Pi的近似值思考：如何利用韋達(dá)公式構(gòu)造出一種迭代算法？2021/6/27數(shù)值積分法計算Pi定積分計算出這個積分的數(shù)值，也就得到了Pi的值。2021/6/27數(shù)值積分法計算Pi1、梯形公式2021/6/27數(shù)值積分法計算Pi2、辛普森（Simpson）公式2021/6/27利用級數(shù)計算Pi1、萊布尼茨級數(shù)（1674年發(fā)現(xiàn)）2021/6/27利用級數(shù)計算Pi

1844年，數(shù)學(xué)家達(dá)什在沒有計算機(jī)的情況下利用此式算出了Pi的前200位小數(shù)。使用誤差估計式計算一下要精確到Pi的200位小數(shù)需要取級數(shù)的多少項？2021/6/27利用級數(shù)計算Pi2、歐拉的兩個級數(shù)（1748年發(fā)現(xiàn)）這兩個級數(shù)收斂也非常緩慢，計算時實(shí)用價值不大。2021/6/27利用級數(shù)計算Pi3、基于arctanx的級數(shù)對泰勒級數(shù)即為萊布尼茨級數(shù)，直接使用時收斂速度極慢，必須考慮加速算法。2021/6/27利用級數(shù)計算Pi觀察級數(shù)可知，x的值越接近于0，級數(shù)收斂越快。由此可以考慮令2021/6/27利用級數(shù)計算Pi因此，β＝4α－pi/4非常接近0。2021/6/27利用級數(shù)計算Pi加速效果非常明顯！2021/6/27蒙特卡羅（MonteCarlo）法單位圓的面積等于Pi，使用蒙特卡羅法，即用隨機(jī)投點(diǎn)的方法來求出這個面積Pi的近似值。具體方法如下：在平面直角坐標(biāo)系中，以O(shè)(0,0)，A(1,0)，C(1,1)，B(0,1)為四個頂點(diǎn)作一個正方形，其面積S＝1。以原點(diǎn)O為圓心的單位圓在這個正方形內(nèi)的部分是圓心角為直角的扇形，面積為S1＝Pi/4。2021/6/27蒙特卡羅（MonteCarlo）法在這個正方形內(nèi)隨機(jī)地投入n個點(diǎn)，設(shè)其中有m個點(diǎn)落在單位扇形內(nèi)。則

想隨機(jī)投點(diǎn)如何來實(shí)現(xiàn)？2021/6/27蒲豐（Buffon）擲針實(shí)驗另一種用蒙特卡羅法來計算Pi的方法是1777年法國數(shù)學(xué)家蒲豐（Buffon）提出的隨機(jī)擲針實(shí)驗。其步驟如下：（1）取一張白紙，在上面畫出許多間距為d的等距平行線。（2）取一根長度為l(l<d)的均勻直針，隨機(jī)地向畫有平行線的紙上擲去，一共擲n次。觀察針和直線相交的次數(shù)m。2021/6/27蒲豐（Buffon）擲針實(shí)驗（3）由幾何概率知道針和直線相交的概率為p=2L/πd，取m/n為p的近似值，則特別取針的長度L為d/2時，π=n/m。2021/6/27蒲豐（Buffon）擲針實(shí)驗（3）由幾何概率知道針和直線相交的概率為p=2L/πd，取m/n為p的近似值，則特別取針的長度L為d/2時，π=n/m。2021/6/27拉馬努金（Ramanujan)公式目前，計算pi的一個極其有效的公式為這個級數(shù)收斂得非?？?，級數(shù)每增加一項，可提高大約8位小數(shù)的精度。2021/6/27拉馬努金（Ramanujan)公式

1985年，數(shù)學(xué)家比爾.高斯帕依使用這個公式在計算機(jī)上算出了pi的1750萬位小數(shù)。這個神奇的公式歸功于印度年輕的傳奇數(shù)學(xué)家拉馬努金(Ramanujan,1887-1929).2021/6/27拉馬努金（Ramanujan)公式另一個經(jīng)過改進(jìn)的計算公式為：級數(shù)每增加一項，可提高14位小數(shù)的精度。2021/6/27迭代公式迭代公式1：1989年，BorWein發(fā)現(xiàn)了下列收斂于1/pi的迭代公式：2021/6/27迭代公式迭代誤差可以由下式估計迭代4次可精確到693位小數(shù)！8次后可以保證精確到小數(shù)點(diǎn)178814位?。?！2021/6/27迭代公式迭代公式2：

1996年，Baiey發(fā)現(xiàn)了另一個收斂于1/pi的迭代公式：2021/6/27迭代公式迭代誤差可以由下式估計2021/6/27結(jié)語隨著計算機(jī)技術(shù)的飛速發(fā)展計算方法的突破與創(chuàng)新，計算Pi的世界紀(jì)錄正在迅速地被刷新。目前，Pi的數(shù)值已計算到小數(shù)點(diǎn)后2061.5843億位。這一記錄是日本東京大學(xué)教授金田康