已阅读5页,还剩13页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
Pi的计算,1.Pi计算(古典方法与级数逼近) 2.数值积分梯形方法 3.数值积分抛物线方法 4.数值微分,1、Pi的计算古典方法,(17世纪以前) 通过圆内接正多边形和圆外切正多边形来夹逼. 早在1500年前的南朝宋齐时期,我国古代数学家祖冲之就用 古典方法给出了应介于3.1415926和3.1415927之间的结果. 由于古典方法计算十分繁琐,费时费劲,改进速度慢,因 此有不少致力于计算的人,获得的小数点后十几位最多几十 位的精度(最好记录39位,1630 德国格林伯格)。,2、Pi的计算级数逼近,从17世纪中叶开始,微积分方法开始诞生,人们逐步掌握了 解决一些实际问题的分析方法和有用工具。,clear; k=0;r=1;p=0;K=-1; while r=1.0e-5 k=k+1; K=K*(-1); p1=p+K/(2*k-1); r=abs(4*(p1-p); fprintf(k=%.0f,p=%.6fn,k,4*p); p=p1; end,MATLAB程序,改进算法:,clear; k=0;r=1;p=0;K=-1; m=1;n=1; while r=1.0e-5 k=k+1; K=K*(-1); m=4*m;n=9*n; p1=p+K/(2*k-1)*(2/m+3/n); r=abs(4*(p1-p); fprintf(k=%.0f,p=%.8fn,k,4*p); p=p1; end,MATLAB程序,改进算法:,麦琴(Machin)公式:,clear; k=0;r=1;p=0;K=-1; m=1;n=1; while r=1.0e-5 k=k+1; K=K*(-1); m=25*m;n=239*239*n; p1=p+K/(2*k-1)*(20/m-239/n); r=abs(4*(p1-p); fprintf(k=%.0f,p=%.8fn,k,4*p); p=p1; end,MATLAB程序,3、Pi的定积分计算,,,定积分的近似计算,在高等数学中有一类积不出的积分,如,(概率积分) (椭圆积分),输入: x=0:0.1:1; y=inline(exp(-x.2); y= exp(-x.2); quad(y,0,1) trapz(x,y) 输出: ans = 0.746211 ans = 0.746826,MATLAB命令 梯形方法 trapz(x,y) 抛物线方法quad(f,a,b),如求积分的近似值,梯形方法,分析:,条件:,当区间划分为n等分时,trapz(x,y),复化梯形方法,当区间等距划分为n个子区间时,现对每个子区间再二等分,得到2n个子区间,举例,求Pi的近似值.,取,运用复化梯形算法,输入初值:,输出结果:STOP,No,Yes,clear;a=0;b=1; f=inline(4/(1+x*x); t1=(b-a)/2*(f(a)+f(b); er=1;n=1; while er1.0e-6 h=(b-a)/n; s=0; for i=1:n s=s+f(a+i*h-h/2); end t2=(t1+h*s)/2; er=abs(t2-t1); fprintf(t=%.6f,r=%.6fn,t2,er); n=2*n; t1=t2; end,抛物线方法,分析:,条件:,梯形法是对每个子区间用梯形面积近似曲线下面积累加而成;而抛物线法是对每个子区间考虑其中点,用三点决定的抛物线下面积来近似.,当区间划分为n(n1)等分时,辛浦生(Simpson)方法,clear; a=0;b=1; f=inline(4/(1+x*x); t1=(b-a)/2*(f(a)+f(b); er=1;s1=0; n=1;,while er1.0e-6 h=(b-a)/n; s=0; for i=1:n s=s+f(a+i*h-h/2); end t2=(t1+h*s)/2;
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 卫星行业十四五竞争格局分析及投资前景与战略规划研究报告2024-2034版
- 中国锌氧化物行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 中国都市景观行业发展分析及发展前景与趋势预测研究报告2024-2034版
- 中国蓝鳍金枪鱼行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 中国碘伏棉球行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 中国测量杆行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 中国智能节电器行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 中国护理模型行业发展前景及发展策略与投资风险研究报告2024-2034版
- 中国套趾拖行业发展趋势及发展前景研究报告2024-2034版
- 中国唇膏瓶行业市场现状分析及竞争格局与投资发展研究报告2024-2034版
- 2024届高考英语二轮复习完形填空 夹叙夹议文课件(共72张)
- 2024年北京市各法院聘用制审判辅助人员招聘笔试参考题库附带答案详解
- 2024年广东省2024届高三二模生物试卷(含标准答案)
- 期中综合检测试题-2023-2024学年五年级下册语文统编版
- 运营计划岗位认知
- 2024中智集团财务岗位招聘笔试参考题库附带答案详解
- 2024年01月天津医科大学招考聘用46人笔试历年(2016-2023年)真题荟萃带答案解析
- 2024年郑州市高三二模(高中毕业年级第二次质量预测)语文试卷
- 环卫保洁日常巡检方案
- 小学心理健康家长会课件
- 预防甲状腺癌指南
评论
0/150
提交评论