




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
§4.3变步长求积公式复化求积方法对提高精度是行之有效的,但是使用复化求积公式之前必须给出适宜的步长,步长取得太大精度难以保证,步长太小那么会导致计算量的增加,而事先给出一个恰当的步又往往是困难的。实际计算时通常采用变步长的求积方案,即在步长逐次折半〔或称步长二分〕的过程中,反复利用复化的求积公式进行计算,直到二分前后的两次积分近似值相当符合为止。变步长求积公式一、变步长的梯形法那么探讨梯形法那么的变步长法那么计算规律:设将求积区间(a,b)分成n等分,那么一共有n+1个分点,按梯形公式(5)计算积分Tn〔对f调用n+1次〕。如果将各求积区间再二分一次,那么分点增至2n+1个,假设仍直接用梯形公式计算二分后的积分值T2n,将需要对f调用2n+1次。一、变步长的梯形法那么一、变步长的梯形法那么注意到T2n的全局部点当中,有一半n+1个是二分前的原有分点,重复计算这些“老分点”上的函数值显然是个浪费。变步长的梯形法那么为了防止浪费,将二分前后的两个积分值联系起来加以考察。注意到每个子区间(xk-1,xk),经过二分再增加一个新分点xk-1/2后,用复化梯形公式求得该子区间上的积分值为因此有整理得:变步长的梯形法那么再利用(5)式,得:这个式子的前一项Tn是二分前的积分值,在求T2n时可作为值使用,而它的后一项只涉及二分时新增加的分点xk-1/2,所要调用f的次数为n。可见递推公式(9)由于防止了老结点的重复计算,而使计算量节约了一半。变步长梯形法那么的计算流程变步长梯形法那么的程序框图变步长梯形法那么的算法框图:其中T1和T2分别代表二分前后的积分值。各框的含义是:[框1]准备初值。[框2]按递推公式(9)求二分后的积分值。从第一个分点x=a+h/2出发,取h为步长逐步向右跨,即可依次确定公式(9)中的各个分点。图中将所得到的分点暂存于单元x中。[框3]控制精度。[框4]修改步长。变步长梯形法那么举例[例4-3-1]用变步长的梯形法那么计算积分值[解]先对整个区间(0,1)使用梯形公式。计算函数在端点的值f(0)=1,f(1)=0.8414710那么T1=[f(0)+f(1)]/2=0.9207355将区间二等分,求中点的函数值f(1/2)=0.9588511[例4-3-1]按递推公式(9),得T2=T1/2+f(1/2)/2=0.9397933进一步二分求积区间,并计算新的分点上的函数f(1/4)=0.9896158f(3/4)=0.9088517再利用(9)式,得这样不断二分下去,计算结果列于表下表中。[例4-3-1]积分I*的实际值是0.9460831用变步长梯形法那么二分10次得到了这个结果。复化梯形公式和变步长法的比较复化梯形法变步长法计算T1时,b-a=1-0=1b-a=1-0=1T1=[f(0)+f(1)]/2=0.9207355计算T2时,h=(b-a)/2=1/2h=b-a=1新增节点f
(1/2)=0.9588511
T2=h/2×[f(1)+f(0)+2f(1/2)]
T2=T1/2+f(1/2)/2=0.9397933
=T1/2+f(1/2)/2=0.9397933复化梯形公式和变步长法的比较复化梯形法变步长法计算T4时,h=(b-a)/2/2=1/4h=(b-a)/2=1/2新增节点f(1/2)=0.9588511
f(1/4)=0.9896158
f(1/4)=0.9896158f(3/4)=0.9088517
f(3/4)=0.9088517T4=h/2×[f(1)+f(0)+2f(1/2)+2f(1/4)+2f(3/4)]T4=T2/2+h/2×[f(1/4)+f(3/4)]复化梯形公式和变步长法的比较复化梯形法变步长法计算T8时,h=(b-a)/2/2/2=1/8h=(b-a)/2/2=1/4新增节点f(1/2)=0.9588511
f(1/8)
f(1/4)=0.9896158f(3/4)=0.9088517
f(3/8)
f(1/8)
f(3/8)f(5/4)f(7/8)
f(5/8)
f(7/8)T8=h/2×[f(1)+f(0)+2f(1/2)+2f(1/4)+2f(1/8)+2f(3/8)+2f(5/8)+2f(7/8)]T8=h/2×[f(1)+f(0)+2f(1/2)+2f(1/4)
T8=T4/2+h/2×[f(1/8)+2f(3/4)]+f(3/8)+f(5/8)+
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 架桥协议书模板
- 2025年陕西申论c类试题及答案
- 2025年数控加工专业测试题及答案
- 2025年前端面试题及答案js
- 2025年销售人才综合测试题及答案
- 2025年EQ产品知识考试题目及答案
- 2025年邮政人员考试试题及答案
- 基于2025年物联网的无土栽培蔬菜产业链发展研究报告
- 校区施工协议书
- 校园整洁协议书
- GB/T 20485.33-2018振动与冲击传感器校准方法第33部分:磁灵敏度测试
- GA/T 1393-2017信息安全技术主机安全加固系统安全技术要求
- GA 304-2012塑料管道阻火圈
- 螺旋藻片的10大销售话术(附4个联合用药方案)
- 2023年燕山大学辅导员招聘考试笔试题库及答案解析
- 食堂从业人员健康情况晨检记录表
- 冷却系统橡胶软管设计基础规范
- 地源热泵埋管冬夏季换热平衡计算
- (中职中专)财经法规与会计职业道德课件完整版电子教案
- 常用呼吸机报警及故障的排除课件
- 微生物的代谢调节 ppt课件
评论
0/150
提交评论