- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《多种方法计算Pi并且精确度比较
多种方法计算圆周率并比较精确度
【摘要】本文介绍了多种方法求圆周率的近似值并对各种方法进行精确度的比较得出具体情况选择的方法,且通过mathematica编程模拟实验过程,得出各种方法的特点。
【关键字】圆周率 数值积分法 泰勒级数法 蒙特卡罗法 拉马努金公式法
0.引言
平面上圆的周长与直径之比是一个常数,称为圆周率,记作。在很长的一段时期,计算的值是数学上的一件重要的事情。有数学家甚至说:“历史上一个国家所得的圆周率的准确程度,可以作为衡量一个国家当时数学发展的一面旗帜。”足以见圆周率扮演的是角色是如此举足轻重。
作为经常使用的数学常数,它的计算已经持续了2500多年了,到今天都依然在进行着,中间涌现出许多的计算方法,它们都各有千秋,在此,我们选择几种较典型的方法,包括数值积分法,泰勒级数法,蒙特卡罗法,韦达公式法,拉马努金公式法以及迭代法来和大家一起体验的计算历程,同时利用mathematica通过对各种方法作精确度的比较得出选择的优先顺序,为相关的理论研究提供一定的参考价值。
数值积分法
以单位圆的圆心为原点建立直角坐标系,则单位圆在第一象限内的部分G是一个扇形,由曲线y=(x[0,1])及两条坐标轴围成,它的面积S=。算出了S的近似值,它的4倍就是的近似值。()
用一组平行于y轴的直线x=(1in-1,a=...=b)将曲边梯形T分成n个小曲边梯形,总面积S分成这些小曲边梯形的面积之和。如果取n很大,使每个小曲边梯形的宽度很小,可以将它上方的边界f(x)()近似地看作直线,将每个小曲边梯形近似当作梯形来求面积,就得到梯形公式,如果更准确一些将每个小曲边梯形的上边界近似地看作抛物线,就得到辛普森公式。
梯形公式:S
辛普森公式:S
Mathematica程序如下:
n=1000;y[x_]:=4/(1+x^2);
s1=(Sum[y[k/n],{k,1,n-1}]+(y[0]+y[1])/2)/n;
s2=(y[0]+y[1]+2*Sum[y[k/n],{k,1,n-1}]+4*Sum[y[(k-1/2)/n],{k,1,n}])/(6*n);
Print[{N[s1,20],N[s2,20],N[Pi,30]}]
注:以上s1,s2分别是用梯形共识和辛普森公式计算出的。最后一句中的N[s1,20]表示s1的前20为准确有效数字组成的近似值。N[Pi,30]是的前30位有效数字组成的近似值。Print[]语句表示将方括号内的数显示出来。
再取n=5000,10000时分别算出的近似,记录入表格如下:
梯形公式 辛普森公式 圆周率的准确值 3.14159265358979323846264338328 n=5000 3.1415926469231265718 3.1415926535897932385 3.14159265358979323846264338328 n=10000 3.1415926519231265718 3.1415926535897932385 3.14159265358979323846264338328 注:下划线标记的数字为用对应方法所求得的精确到的位数。
观察发现,当n的取值越大,即当[a,b]之间分的越细时,计算所得的的精确度越高。n 从1000到10000,梯形公式所求得的近似值精确度增加两位,而辛普森公式所求的的近似值有效位数没有变化,可见n的变化对梯形公式影响较大。但总的而言利用辛普森公式所求得的近似值精度比用梯形公式所求得的近似值很多。当要求精确度较高时,两种方法中可优先选择辛普森公式法。
泰勒级数法
利用反正切函数的泰勒级数:Arctan x=x-
?X=1代入上式得到莱布尼兹级数:
?令a=arctan,b=-a,则tan b=tan(-a)===
因此b=arctan,即-arctan=atctan,从而得到
=arctan+arctan 即=4(arctan+arctan)
?令a=arctan,tan a=,则容易算出tan2a=,tan4a=,
Tan(4a-)===
4a-=arctan
=4a-arctan=4arctan-arctan 即=16arctan-4arctan
利用mathematica编写???程序如下:
T[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/(2 k-1),{k,1,n}]
4N[T[1,1000],20]
N[4(T[1/2,1000]+T[1/3,1000]),50]
N[16T[1/5,1000]-4 T[1/239,1000],50]
N[Pi,50]
注:4N[T[1,1000],20
您可能关注的文档
- 《北京交通大学翻硕考研复试分数线未来趋势.doc
- 《北京交通大学考研情况详细解析.doc
- 《外文文献翻译译稿.doc
- 《北京十中分校已开发的校本课程名称.doc
- 《北京博飞港澳台全国联考中近史之新文化运动卷一 (含答案).doc
- 《外暖改施组.doc
- 《外来语.doc
- 《北京初二英语目标答案.doc
- 《外来语翻译.doc
- 《外文翻译要求and格式.doc
- 第十一章 电流和电路专题特训二 实物图与电路图的互画 教学设计 2024-2025学年鲁科版物理九年级上册.docx
- 人教版七年级上册信息技术6.3加工音频素材 教学设计.docx
- 5.1自然地理环境的整体性 说课教案 (1).docx
- 4.1 夯实法治基础 教学设计-2023-2024学年统编版九年级道德与法治上册.docx
- 3.1 光的色彩 颜色 电子教案 2023-2024学年苏科版为了八年级上学期.docx
- 小学体育与健康 四年级下册健康教育 教案.docx
- 2024-2025学年初中数学九年级下册北京课改版(2024)教学设计合集.docx
- 2024-2025学年初中科学七年级下册浙教版(2024)教学设计合集.docx
- 2024-2025学年小学信息技术(信息科技)六年级下册浙摄影版(2013)教学设计合集.docx
- 2024-2025学年小学美术二年级下册人美版(常锐伦、欧京海)教学设计合集.docx
文档评论(0)