微分方程初值问题的数值解法.pptVIP

第1页，共41页，星期日，2025年，2月5日引言初值问题的数值解法:求初值问题的解在一系列节点的值y(xn)的近似值yn的方法.本章数值解法的特点:都是采用“步进式”,即求解过程顺着节点排列的次序一步步向前推进.基本知识:(1)定理1:如果函数f(x,y)在区域上连续,且关于y满足Lipschitz条件常微分方程初值问题:求未知函数y=y(x).第2页，共41页，星期日，2025年，2月5日此时Lipschitz条件显然成立.故常用在D上连续有界来代替f(x,y)关于y满足Lipschitz条件.注:如无特别说明,总假设(1)的解存在唯一且足够光滑.在

f(x,y)对变量y可微的情形下,若偏导数连续有界,则可取L为除了要保证(1)有唯一解外,还需保证微分方程本身是稳定的,即(1)的解连续依赖于初始值和函数f(x,y).也就是说,当初始值y0及函数f(x,y)有微小变化时,只能引起解的微小变化.(其中L称为Lipschitz常数),则对任何,初值问题(1)在[a,b]上存在唯一连续可微解y=y(x).定理2:如果函数f(x,y)在区域上关于y满足Lipschitz条件,则(1)是稳定的.第3页，共41页，星期日，2025年，2月5日单步迭代:计算yn+1时仅用yn;初值问题(1)与下列积分方程的解等价:初值问题的数值解就是求一系列节点上函数y=y(x)的近似值.称为步长.一般取等步长h.多步迭代:计算yn+1时除用yn外,还要用到yn-1,yn-2,…;k步迭代要用到yn-1,yn-2,…,yn-k+1.显式单步迭代:隐式单步迭代:(2)第4页，共41页，星期日，2025年，2月5日一、Euler方法及其改进将[a,b]n等分,记微分法:积分法:积分项利用矩形公式计算1.显式Euler方法(★)第5页，共41页，星期日，2025年，2月5日Taylor公式推导:Euler公式几何意义:P1P2Pk也称折线法P0xy第6页，共41页，星期日，2025年，2月5日2.梯形法称之为梯形公式.这是一个隐式公式,通常用迭代法求解.具体做法:取先用Euler法求出初值,即,将其代入梯形公式的右端,使之转化为显式公式,即注:当f(x,y)关于y满足Lipschitz条件且步长h满足直至满足:若采用梯形公式计算(★)中的积分项,则有类似地,可得(☆)第7页，共41页，星期日，2025年，2月5日时,迭代格式(☆)收敛.3.改进的Euler方法把Euler法作为预报(称为预估公式),把隐式的梯形公式作为校正(称为校正公式),则得改进的Euler方法:或也称为预估-校正法.有时为了方便,预估-校正格式也写成下面形式:第8页，共41页，星期日，2025年，2月5日二、单步法的局部截断误差及精度Def1:先假设,再估计误差这种误差称为单步迭代法在xk+1处的局部截断误差.Def2:若某种数值方法的局部截断误差为,则称该数值方法的精度为P阶的.注:通常情况下,P越大,h越小,则截断误差越小,数值方法越精确.设10.Euler方法是一阶方法.第9页，共41页，星期日，2025年，2月5日所以Euler方法为一阶方法.而20.梯形法是二阶方法.Taylor展开第10页，共41页，星期日，2025年，2月5日将代入上式,得而代入上式得:当h充分小时,若,则可选取h,使得第11页，共41页，星期日，2025年，2月5日故梯形法的精度为2.同样可以证明改进的Euler法也是二阶方法.梯形法的局部