7常微分方程初值问题解法_第1页
7常微分方程初值问题解法_第2页
7常微分方程初值问题解法_第3页
7常微分方程初值问题解法_第4页
7常微分方程初值问题解法_第5页
已阅读5页，还剩31页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1,第7章常微分方程初值问题数值解法,7.1引言,2,本章研究的问题:,3,7.1欧拉方法及改进的欧拉法7.1.1欧拉公式,图7.1欧拉折线法,4,1欧拉公式,(2),(3),5,6,2欧拉公式的截断误差,(4),7,3单步法的局部截断误差与阶,局部截断误差可以理解为计算一步的误差.,8,则称该方法具有P阶精度.,9,4后退的欧拉方法,(5),10,(6),(6)式称为后退的欧拉方法,它是隐式的,欧拉公式(2)是显式的,11,(7),12,13,后退的欧拉方法的局部截断误差:,14,5梯形方法,(8)式称为梯形方法.,(8),15,梯形方法的局部截断误差:,16,6.改进欧拉法及局部截断误差,预测步,校正步,或者写成,（1）改进的欧拉公式:,17,（2）改进的欧拉方法的局部截断误差,18,考虑改进Euler法,如果将其改成,-(1),7.2Runge-Kutta法,19,改进Euler法是由梯形公式和Euler公式复合而成,梯形公式具有2阶精度,(1)式为一种二阶Runge-Kutta法,同样可以证明,改进Euler法也具有2阶精度。,20,Runge-Kutta方法的推导,21,Runge-Kutta方法的一般形式：,确定了阶数之后，再通过Taylor展开、比较两边系数的方法，确定各待定系数：,22,二阶显式Runge-Kutta方法,23,24,25,26,例,27,结果及比较,28,三阶显式Runge-Kutta方法,在推导二阶显式方法的过程中，注意到局部截断误差表达式中h3项包含了以下表达式：,因此若要在局部截断误差中消去h3项，必须增加包含了以上各项的多个方程，同时我们注意到r=2时，只有等四个待定系数，少于方程的数目，所以这样的系数不存在。故：r=2时Runge-Kutta方法只能是二阶的。要得到三阶的方法，则必须有r=3。,29,三阶显式Runge-Kutta方法,30,四阶显式Runge-Kutta方法,31,例,32,结果及比较,33,关于Runge-Kutta方法,34,提高Runge-Kutta方法的精度的方法,提高精度最简单的方法是缩短步长，但要以牺牲计算速度和积累舍入误差为代价。,35,变步长的Runge-Kutta方法,作为妥协，如果能在计算过程中实时控制步
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论