数学实验(算圆周率).docVIP

下载本文档

64
0
约1.59万字
约 3页
2017-06-16 发布于河南
举报
版权申诉

数学实验(算圆周率).doc

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数学实验(算圆周率)

数学实验报告实验问题人们很早就发现圆的周长和直径的比是个常数，也就是我们都熟悉的圆周率π，那么这个常数的值究竟是多少哪？我国的伟大数学家祖冲之将π精确到3.1415926-3.1415927之间，西方的数学研究者也从来没有停止过对π的计算，那么在今天计算机技术如此发达的情况下，我们如何将π的近似值精确到更多位数哪？二、问题的分析 π的精确度问题困扰大家已久，所以我们不妨在MATLAB中采用循环来实现对π的更高位数的精确。不妨以精确到小数点后9位为例。我们知道因此可以利用三、实验内容在MATLAB中键入以下程序 clear; n=0; r=1; p=0; k=-1; a=1; b=1; while r=1.0e-9 n=n+1; k=k*(-1); a=4*a;b=9*b; pl=p+k/(2*n-1)*(2/a+3/b); r=abs(4*(pl-p)); fprintf(n=%.0f,p=%.10f\n,n,4*pl); p=pl; end 四、问题求解结果和结论 n=1,p=3.3333333333 n=2,p=3.1172839506 n=3,p=3.1455761317 n=4,p=3.1408505618 n=5,p=3.1417411974 n=6,p=3.1415615879 n=7,p=3.1415993410 n=8,p=3.1415911844 n=9,p=3.1415929813 n=10,p=3.1415925796 n=11,p=3.1415926705 n=12,p=3.1415926497 n=13,p=3.1415926545 n=14,p=3.1415926534 n=15,p=3.1415926536 对泰勒级数，当｜x｜越小时级数收敛速度越快，这启示我们另外构造相关级数来逼近π。事实上，容易证明： =4artan-arctan =2artan+arctan 五、问题的进一步拓展和实验用复化梯形法计算π的值在MATLAB中键入以下程序： clear; f=inline(4./(1+x.*x)); a=0;b=1;n=1; h=(b-a)/n; t1=h/2*(f(a)+f(b)); er=1;k=1; while er1.0e-5; s=0; for i=1:n s=s+f(a+(i-1/2)*h); end t2=(t1+h*s)/2; er=abs(t2-t1); fprintf(n=%.0f,p=%.6f,r=%.6f\n,k,t2,er); n=2*n;h=h/2;t1=t2; k=k+1; end 六、实验的总结和体会通过用MATLAB对π的近似值的计算，加深了对MATLAB循环结构的认识，方便了以后对数学软件的进一步使用。