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

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

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

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

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

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

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

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

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

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

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

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

圓面積S與多邊形的面積Sn之間有如下關(guān)系：2021/6/27劉徽不等式借助于計算機來完成劉徽的工作：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ù)學的方法。這種方法相當于今天的“求積分”，后者經(jīng)16世紀英國的牛頓和德國的萊布尼茨作系統(tǒng)總結(jié)而得名。鑒于劉徽的巨大貢獻，所以不少書上把他稱做“中國數(shù)學史上的牛頓”，并把他所創(chuàng)造的割圓術(shù)稱為“徽術(shù)”。2021/6/27韋達（VieTa）公式

1593年，韋達首次給出了計算Pi的精確表達式：韋達公式看起來有些神秘，其實它的導出過程所用的都是樸實簡潔的數(shù)學方法。2021/6/27韋達（VieTa）公式1、從sint開始2021/6/27韋達（VieTa）公式所以，對任意N，總有2021/6/27韋達（VieTa）公式2、從cos(pi/4)開始2021/6/27韋達（VieTa）公式3、使用VieTa公式計算Pi的近似值思考：如何利用韋達公式構(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ù)學家達什在沒有計算機的情況下利用此式算出了Pi的前200位小數(shù)。使用誤差估計式計算一下要精確到Pi的200位小數(shù)需要取級數(shù)的多少項？2021/6/27利用級數(shù)計算Pi2、歐拉的兩個級數(shù)（1748年發(fā)現(xiàn)）這兩個級數(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，使用蒙特卡羅法，即用隨機投點的方法來求出這個面積Pi的近似值。具體方法如下：在平面直角坐標系中，以O(shè)(0,0)，A(1,0)，C(1,1)，B(0,1)為四個頂點作一個正方形，其面積S＝1。以原點O為圓心的單位圓在這個正方形內(nèi)的部分是圓心角為直角的扇形，面積為S1＝Pi/4。2021/6/27蒙特卡羅（MonteCarlo）法在這個正方形內(nèi)隨機地投入n個點，設(shè)其中有m個點落在單位扇形內(nèi)。則

想隨機投點如何來實現(xiàn)？2021/6/27蒲豐（Buffon）擲針實驗另一種用蒙特卡羅法來計算Pi的方法是1777年法國數(shù)學家蒲豐（Buffon）提出的隨機擲針實驗。其步驟如下：（1）取一張白紙，在上面畫出許多間距為d的等距平行線。（2）取一根長度為l(l<d)的均勻直針，隨機地向畫有平行線的紙上擲去，一共擲n次。觀察針和直線相交的次數(shù)m。2021/6/27蒲豐（Buffon）擲針實驗（3）由幾何概率知道針和直線相交的概率為p=2L/πd，取m/n為p的近似值，則特別取針的長度L為d/2時，π=n/m。2021/6/27蒲豐（Buffon）擲針實驗（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ù)學家比爾.高斯帕依使用這個公式在計算機上算出了pi的1750萬位小數(shù)。這個神奇的公式歸功于印度年輕的傳奇數(shù)學家拉馬努金(Ramanujan,1887-1929).2021/6/27拉馬努金（Ramanujan)公式另一個經(jīng)過改進的計算公式為：級數(shù)每增加一項，可提高14位小數(shù)的精度。2021/6/27迭代公式迭代公式1：1989年，BorWein發(fā)現(xiàn)了下列收斂于1/pi的迭代公式：2021/6/27迭代公式迭代誤差可以由下式估計迭代4次可精確到693位小數(shù)！8次后可以保證精確到小數(shù)點178814位！??！2021/6/27迭代公式迭代公式2：

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