原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验9π的近似计算309宿舍
数学实验报告
实验序号： 日期：2011年 10月 25日
班级 姓名 姓名 （） 姓名 姓名 姓名 实验
名称 π的近似计算 问题背景描述：π是人们经常使用的数学常数，对π的研究已经持续了2500多年，今天，这种探索还在继续中。
1.割圆术。2.韦达（VieTa）公式。3.利用级数计算π。4.拉马努金（Ranmaunujan）公式。5.迭代方法。6.π的两百位近似值。 实验目的：1.了解圆周率π的计算历程。
2.了解计算π的割圆术、韦达公式、级数法、拉马努金公式、迭代法。
3.学习、掌握MATLAB软件有关的命令。 实验所用软件及版本：
主要内容（要点）：利用韦达（VieTa）公式
计算π的近似值。 实验过程记录：练习1 用刘徽的迭代公式
计算π的近似值。
相应的MATLAB代码为
clear;
x=1;
for i=1:30
x=vpa (sqrt(2-sqrt(4-x^2)),15)%计算精度为15位有效数字
S=vpa(3*2^i*x,10)
end
计算可得
x =.517638090205042 S =3.105828541
x =.261052384440103 S =3.132628613
x =.130806258460286 S = 3.139350203
x =.654381656435522e-1 S =3.141031951
x =.327234632529735e-1 S =3.141452472
x =.163622792078742e-1 S =3.141557608
x =.818120805246955e-2 S =3.141583892
x =.409061258232818e-2 S =3.141590463
x =.204530736067660e-2 S =3.141592106
x =.102265381402739e-2 S =3.141592517
x =.511326923724832e-3 S =3.141592619
x =.255663463951308e-3 S =3.141592645
x =.127831732236766e-3 S =3.141592651
x =.639158661510219e-4 S =3.141592653
x =.319579330795908e-4 S =3.141592653
x =.159789665403054e-4 S =3.141592654
x =.798948327021645e-5 S =3.141592654
x =.399474163511619e-5 S =3.141592654
x =.199737081755909e-5 S =3.141592654
x =.998685408779670e-6 S =3.141592654
x =.499342704389851e-6 S = 3.141592654
x =.249671352194927e-6 S =3.141592654
x =.124835676097464e-6 S =3.141592654
x =.624178380487320e-7 S =3.141592654
x =.312089190243660e-7 S =3.141592654
x =.156044595121830e-7 S =3.141592654
x =.780222975609150e-8 S =3.141592654
x =.390111487804575e-8 S =3.141592654
x =.195055743902288e-8 S =3.141592654
x =.975278719511453e-9 S =3.141592654
练习2 用韦达公式计算π的近似值。
相应的MATLAB代码为
clear;
x=1;
for i=1:20;
i
x=vpa(x*cos(pi/2^(i+1)),30);%计算精度为30位有效数字
pai=vpa(2/x,22)
error=vpa(pai-2/x,22)
end
计算可得
i =1 pai =2.828427124746190097603 error =.3774484188e-21
i =2 pai = 3.061467458920718232298 error =.2223605969e-21
i =3 pai =3.121445
文档评论（0）