常微分方程数值解第六章_第1页
常微分方程数值解第六章_第2页
常微分方程数值解第六章_第3页
常微分方程数值解第六章_第4页
常微分方程数值解第六章_第5页
已阅读5页，还剩53页未读， 继续免费阅读

下载本文档

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

文档简介

计算方法,江西理工大学,第六章常微分方程数值解法,6.1引言1、一阶常微分方程初值问题的一般形式是,其中:,2、方程6.1解存在定理,3、数值解的分类,6.2.1Euler公式假设初值问题(6.1)式(6.2)式的解y=y(x)存在且足够光滑，对求解区域a,b分成n+1个节点:,6.2Euler方法,比较6.4式和6.5式，为求得,只需用到,这种方法称为单步法，而6.6式需要,。这种方法称为多步法。,6.4式和6.6式中的,被显式的表示出来了，故,被称为显式公式，而6.5式的两边都含有,项，因,而被称为隐式公式。,从计算结果可见,步长h越小,数值解的精度越高.,近似值,准确值,6.2.2改进的Euler方法,从计算结果可见,改进的Euler方法明显地改善了精度.,6.2.3Euler公式的误差分析,y(xn)表示精确值,6.2.4*Taylor展开方法,6.3Runge-Kutta方法,由Euler方法:,比较例6.2与例6.3的计算结果,显然四阶R-K方法的精度高.尽管四阶R-K方法的计算量比改进的Euler方法大,但由于放大了步长,在求相同节点上的近似值时,所需的计算量几乎相同.以上讨论的是显式R-K方法,同样也可以构造隐式R-K方法,其一般形式,四阶Runge-Kutta算法,6.4单步方法的收敛性和稳定性初值问题的数值解法是经过某种离散化过程导出的,因此需要对数值解法进行定性分析.本节主要讨论单步方法的收敛性与稳定性.,6.4.1单步方法的收敛性,定义6.1,定理6.1,证明,Lipschiz条件：|f(x,y2)-f(x,y1)|L|y2-y1|,等比数列公比为,在收敛性的讨论中,我们已假定差分方程是精确求解的,但实际情况并非如此.例如,初始数据可能存在误差,计算过程中也不可避免地产生计算舍入误差,这些误差的传播和积累都会影响到数值解.那么实际计算得出的数值解能否作为精确解的近似呢?这取决于计算误差是否可控制,这就是数值方法稳定性的问题.,定义6.2,6.4.2单步方法的稳定性,综上所述，收敛性是反映差分公式本身的截断误差对数值解的影响；稳定性是反映计算过程中舍入误差对数值解的影响。单步显式方法的稳定性与步长密切相关，在一种步长下是稳定的差分公式，取大一点步长就可能是不稳定的，只有既收敛又稳定的差分公式才有实用价值。6.5线性多步方法,6.5.1利用待定参数法构造线性多步方法,6.5.2利用数值积分构造线性多步方法,Adams显式公式：,Adams隐式公式：,预校算式每一步只需重新计算f(x,y)的函数值二次，因此比四阶标准R-K公式的计算量小。其缺点是要用其他方法计算起始值，计算过程中改变步长困难。,6.6常微分方程组与高阶微分方程的数值解法一阶常微分方程初值问题的数值方法，原则上都可推广到一阶方程组和高阶方程情
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

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

最新文档

评论

0/150

提交评论