数值分析5---常微分方程初值问题.pptVIP

第五章常微分方程初值问题引言根本概念Euler方法及其改进

§1引言1.常微分方程的定解问题与应用应用：自然科学领域，如物理；工程技术问题，如石油勘探。为的，该区域人口的自然增长率为。人口的增长与人口的常微分方程的定解问题主要有初值问题和边值问题两大类，我们仅考虑初值问题。实例：马尔萨斯人口模型：假设某特定区域在时刻的人口总数成正比，所以t时刻的人口总数满足如下的微分方程：一般地，称这样的方程为模型方程。

2.常微分方程的解法：(1)解析法：给出精确解析解。只适合少数简单情况。(2)近似解法：给出解的近似表达式。如级数法,逐步逼近法。(3)数值方法：给出方程在离散点上的近似解。它适合计算机求解,应用广泛,具有理论应用价值。3.常微分方程初值问题的数值方法〔理论和计算方法〕单步法Euler方法Taylor方法和Runge-Kutta方法多步法Adams方法和一般线性多步法线性多步法的收敛性与稳定性

§2根本概念一、常微分方程初值问题的一般提法问题:求函数满足其中为已知函数,是已知值.(可能是观察值或实验值)根本条件:设(2)f(x,y)在D上关于变量y满足Lipschitz连续条件：(1)f(x,y)在D上连续；其中L为Lipschitz常数。定理1假设f(x,y)在D上满足根本条件,一阶常微分方程初值问题式(1),(2)对任意给定的存在唯一解且在[a,b]上连续可微.(2)(1)(3)

关于解y(x)的适定性：定义1方程(1),(2)的解y(x)称为适定的,若存在常数对任意满足条件常微分方程(1),(2)上各加一个摄动(扰动)项.存在唯一解z(x),且有初值问题：摄动(扰动)误差定理2假设f(x,y)在D上满足根本条件,那么微分方程(1),(2)的解y(x)是适定的.(1)适定问题的解y(x)连续依赖于(1)式右端的f(x,y)和初值。或者说解y(x)关于(1)式右端的f(x,y)和初值稳定.注:(2)假设f(x,y)在D上满足根本条件,从而方程(1),(2)的解y(x)存在且适定.(4)

二、初值问题数值解的根本概念因为初值问题的数值解法是通过微分方程离散化而给出解在某些离散点上(节点上)的近似值，为了讨论问题方便,引入以下概念。在上引入节点常用等步长:则有(1)，(2)的准确解记为y(x),求初值问题数值解的方法是步进法,即逐个节点计算,由称为步长。的近似解记为步进法单步法:多步法:仅仅由计算计算共用到l个值.即称为l步法。单步法与多步法的区别:(1)计算方面:l步方法只用于的计算,的计算要用其它方法。

(2)理论分析:单步法比的多步法容易分析(稳定性).(3)选步长方面:单步法容易改变步长.(4)精度:多步法精度高一些.单步法与多步法又都有显式方法和隐式方法之分.计算公式依次可写成:显式单步法:隐式单步法:该式右端项含有因此若求需要解方程。注:显式多步法:隐式多步法:线性多步法:注:(9)关于都是线性的。其中是独立于k和f的常数。(9)是显式(右端不含有)；那么(9)是隐式的。(5)(6)(7)(8)(9)

§3Euler方法考虑问题特点：简单，精度低.一、显式Euler方法(折线法)1.显式欧拉公式设节点为:那么(1),(2)的Euler方法为步长其中三、微分方程初值问题的数值解法讨论的问题:1.方法构造2.误差分析3.稳定性(1)(2)(10)

推导公式:(1)Taylor展开法将在点进行Taylor展开忽略这一高阶项,由于因此分别用近似得结合初值条件,并(2)向前差分近似微分法向前差商近似得(11)于是，即得(10)式：(12)

(3)左矩数值积分法将(1)式两端从到积分,得数值积分采用左矩形公式,即由初始条件亦得(10)式.o用近似得将近似号改为等号，近似并结合初始条件即得(10)式。

2.几何意义方程(1),(2)的解曲线过点具有斜率从出发以为斜率作直线段,交于点此时从出发,以为斜率于点此时以次类推.最终到达点这样得到了一条折线用折线作为(1),(2)解曲线的近似曲线，折线法。o过的解曲线如图示,具有斜率作直线段,交过的解曲线如图.它在点的右侧具有斜率与过的解曲线相切.折线称为欧拉折线，所以欧拉方法又称为

二、隐式Euler方法和梯形方法(Euler方法的改进)将在点进行Taylor展开，忽略项,用分别近似可得隐式欧拉方法:1.隐式Euler方法或说明:(1)隐式Euler方法也可用近似微分或者用右矩数值求积公式来建立.向后差商,即用(2)隐式Euler方法(13)是关于的方程，若求需要解方程(13).(13)

2.梯形方法取平均,得忽略项,用得梯形方法:与说明:梯形方法也是隐式方法,要求,需解方程(14).由分