数值分析第8章常微分方程初值问题的数值解法.pdfVIP

第八章第八章 第八章第八章 常微分方程初值问题的数值解法常微分方程初值问题的数值解法 常微分方程初值问题的数值解法常微分方程初值问题的数值解法 §8.1 引言 §8.2 几种简单的单步法 §8.3 Runge –Kutta 公式 §8.4 单步法的收敛性、 相容性与稳定性 §8.5 线性多步法 数值算例 2009-09-26 1 §8.1 引言 一、问题 本章着重讨论一阶常微分方程初值问题  ′ dy y = = f （x ，y ）   dx 在区间[a,b]上的数值解法。  y （x ）= y  0 0 这些问题多数情况下求不出解析解，只能用近似 的方法求解。常用的近似方法有两类。一类称为 近似解析法，如级数解法，逐次逼近法等。另一 类称为数值解法，它可以给出解在一些离散点上 的近似值。 2009-09-26 2 二、解的存在唯一性 { } 若f （x ,y ）在区域D= a ≤ x ≤ b，y ∈R 上连续，且 关于y 满足 李普希兹 （Lipschitz ）条件,即存在常 数L ,使 f ( x , y 1 ) − f ( x , y 2 ) ≤ L y 1 − y 2 ， y ,y 对 内任两个 1 2 均成立，其中L 是与x ，y 无 关 G 的常数，则上面的初值问题存在唯一解，且解是 连续可微的。 Remark ：在f (x ,y )对y 可微的情况下，若偏导数有 ∂f ( x , y ) 界，则可取 L = max ，此时有 ( x , y )∈ D ∂y ∂f (x,ξ ) f (x,y ) − f (x,y ) = (y − y ) ≤ L y − y ,ξ 介于y 与y 之间。 1 2 ∂y 1 2 1 2 1 2 此时Lipschitz连续条件显然成立。这是验证该条件 的最简便的方法。 2009-09-26 3 三、解的适定性 解的适定性是指解的存在唯一性以及数值稳定性。 此处主要是指解对于右端项以及初值扰动的适应 性。关于适定性有如下的结论： { } 定理：若f （x ,y ）在区域D= a ≤ x ≤ b，y ∈R 上满 足Lipschitz连续条件，则初值问题是适定的。 四、等价的积分方程 若y (x)是初值问题的解，对方程两边同时积 分，利用初始条件可得： x y (x) =