原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析求π
怎么计算π的值：
用莱布尼兹公式、韦达公式和斯特林公式分别迭代，当结果与π确切值之间误差小于0.005时结束迭代，通过分析三个公式的迭代次数得出最优的迭代公式。
莱布尼兹公式：
根据数列得出：
通过程序实现如下（laibulizi.m）：
for n=1:1000
n
a(n)=(-1).^(n-1)./(2*n-1);
%y=vpa(4*sum(a),30)
%error=vpa(y-pi,30)
y=4*sum(a)
error=y-pi
if(error=0.005error=-0.005)
break;
end
end;
运行结果为：
韦达公式：
通过分析公式可知：
则：
通过程序实现如下（weida.m）：
clear;
x=1;
for i=1:100;
i
x=x*cos(pi/2^(i+1));
pai=2/x
error=pai-pi
if(error=0.005error=-0.005)
break;
end
end
运行结果为：
斯特林公式：
程序实现如下（sitelin.m）：
for n=1:150
n
a=exp(n);
b=factorial(n);
c=sqrt(n);
d=(a*b)/((n^n)*c);
p=(d^2)/2
err=p-pi
if err=0.005
break;
end
end
分析：通过运行结果进行分析，当误差小于0.05时，使用莱布尼兹公式迭代了200次，得到的结果为3.1366，通过韦达公式迭代了5次，结果为3.1403，通过斯特林公式迭代了105次，结果为3.1466。通过比较得出韦达公式最快。
注：使用误差为0.005是因为在计算斯特林公式时计算到143项时结果就变成了无穷，所以选择了一个比较大的误差进行判断，其他两个公式（高斯公式和拉马努金公式）由于编程计算失败，故没有进行比较。
文档评论（0）