圆周率的计算培训资料.ppt

《圆周率的计算培训资料.ppt》由会员分享，可在线阅读，更多相关《圆周率的计算培训资料.ppt（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、圆周率的计算1 1、实验时期、实验时期 通过实验进行估算通过实验进行估算,这是计算圆周率这是计算圆周率的的第一阶段的的第一阶段中国中国:(1)“圆径一而周三圆径一而周三”-周髀算经周髀算经(2)“周三径一周三径一,方五斜七方五斜七”-木工口诀木工口诀古埃及：数谷粒与称重量：古埃及：数谷粒与称重量：1579年年,韦达证明韦达证明1630年年,德国人鲁道夫德国人鲁道夫:小数点以后小数点以后35位位1150年年,印度数学家婆什迦罗印度数学家婆什迦罗:3.141614241424年年,中亚细亚数学家卡西中亚细亚数学家卡西:3.14159265358979325.3.141592653589794238

2、462643383279502883、分析方法、分析方法从十七世纪中叶起从十七世纪中叶起,人们开始用分析方法来求人们开始用分析方法来求的近似值的近似值,其中其中应用的主要工具是收敛的无穷乘积和无穷级数应用的主要工具是收敛的无穷乘积和无穷级数.麦琴麦琴(Machin)给出给出(Machin公式公式)1656年年,沃里斯沃里斯(Wallis)证明证明取取k=10取取k=20欧拉证明了欧拉证明了(1735)注：注：BaselProblem,1644年提出年提出1、Buffon投针投针1.在白纸上画上许多条间距为在白纸上画上许多条间距为d的平行直线的平行直线2.取长为取长为l(ld)的针的针,随机地投

3、掷在白纸上随机地投掷在白纸上,投掷投掷n次次,观察与直线相观察与直线相交的次数交的次数,记为记为m.4、概率方法、概率方法针线相交的概率针线相交的概率随机投针的概率含义随机投针的概率含义(1)针的中点针的中点M与平行线的距离与平行线的距离x均匀分布于区间均匀分布于区间0,d/2(2)针与平行线的交角均匀分布于区间针与平行线的交角均匀分布于区间记针的中点为记针的中点为M,x表示点表示点M与较近的平行线间距离与较近的平行线间距离,表示针与表示针与平行线间夹角平行线间夹角在间隔为在间隔为d的平行线间随机投掷长度为的平行线间随机投掷长度为l的针的针0,d/2中随机选取中随机选取x,0,中随机产生中随机

4、产生,构成平面中点构成平面中点x,计算针和直线相交的概率是计算针和直线相交的概率是设投掷设投掷n次次,相交相交m次次,则针与线相交的频率为则针与线相交的频率为m/n针与平行线相交的条件针与平行线相交的条件n=2212,Buffon:=3.142;n=5000,Wolf:=3.1593d=45;l=36;n=20000;x=;y=;P=;Q=;fori=1:n,x1=rand*d*0.5;y1=rand*pi;if2*x1=l*sin(y1),x=x,x1;y=y,y1;elseP=P,x1;Q=Q,y1;endendm=length(x),p=2*l*n/(d*m)s=0:0.01:pi;pl

5、ot(s,l*sin(s)/2,k,Linewidth,2)holdonplot(y,x,r.),plot(Q,P,b.),plot(0,pi,0,0,k,Linewidth,2),plot(0,pi,22.5,22.5,k,Linewidth,2)plot(0,0,0,22.5,k,Linewidth,2),plot(pi,pi,0,22.5,k,Linewidth,2)axis(-0.2,3.3,-2,24)利用蒙特卡洛算法计算圆周率利用蒙特卡洛算法计算圆周率利用蒙特卡洛算法计算定积分利用蒙特卡洛算法计算定积分2、互素数分布、互素数分布取一个大整数取一个大整数N,在在1到到N之间随机取一对

6、整数之间随机取一对整数a,b,它们互素的概率它们互素的概率注：注：随机整数随机整数randint(1,1,)求求a,b最大公约数最大公约数gcd(a,b)5、数值积分方法、数值积分方法数值积分方法数值积分方法6、代数迭代、代数迭代对正数对正数a0,b0,定义算术均值数列和几何均值数列定义算术均值数列和几何均值数列若两数列极限相等若两数列极限相等,则称此极限为它们的算术几何均值则称此极限为它们的算术几何均值,记为记为AGM(a0,b0)取取则有则有记记年代年代19491973198919992011精确位数精确位数2035100万万10亿亿2061亿亿2000万亿万亿“十位小数就足以使地球周界准

7、确到一英寸以内十位小数就足以使地球周界准确到一英寸以内,三十位小数便三十位小数便能使整个可见宇宙的四周准确到连最强大的显微镜都不能分辨能使整个可见宇宙的四周准确到连最强大的显微镜都不能分辨的一个量。的一个量。”西蒙西蒙纽克纽克姆姆 实验一小结实验一小结1、几何算法、几何算法:割圆法割圆法,正多边形逼近圆正多边形逼近圆2、分析算法、分析算法:无穷级数无穷级数,有限项逼近无穷项有限项逼近无穷项3、数值分析、数值分析:定积分计算定积分计算,矩形逼近曲边梯形矩形逼近曲边梯形4、概率统计、概率统计:蒲丰投针、整数互素蒲丰投针、整数互素提示提示:思考题思考题(1)能设计其它计算圆周率的算法吗能设计其它计算圆周率的算法吗?(2)如何计算无理数如何计算无理数e?结束语结束语谢谢大家聆听！谢谢大家聆听！23