第9章微分方程初值问题的数值解法.ppt

* 第九章 微分方程初值问题的数值解法 内容提纲 引言 Euler法及其改进 Runge-Kutta方法 线性多步法 误差分析 数值解法的收敛性、相容性和稳定性 边值问题数值解法简介 引言 初值问题的数值解法:求初值问题的解在一系列节点的值 y ( xn )的近似值 yn 的方法.本章数值解法的特点:都是采用“步进式”,即求解过程顺着节点排列的次序一步步向前推进. 基本知识: (1) 定理1: 如果函数 f (x , y)在区域 上连续,且关于 y 满足Lipschitz条件 常微分方程初值问题: 求未知函数 y= y (x) . 此时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)是稳定的. 单步迭代: 计算 yn+1时仅用 yn ; 初值问题(1)与下列积分方程的解等价: 初值问题的数值解就是求一系列节点 上函数 y= y (x)的近似值 . 称为步长. 一般取等步长 h . 多步迭代: 计算 yn+1时除用 yn 外, 还要用到 yn-1, yn-2 ,…; k 步迭代要用到 yn-1, yn-2 ,…, yn-k+1 . 显式单步迭代: 隐式单步迭代: (2) 一、Euler方法及其改进 将[a , b ]n 等分, 记 微分法: 积分法: 积分项利用矩形公式计算 1. 显式Euler方法 (★) Taylor公式推导: Euler公式几何意义: P1 P2 Pk 也称折线法 P0 x y 2. 梯形法 称之为梯形公式.这是一个隐式公式,通常用迭代法求解.具体做法: 取 先用Euler法求出初值 ,即 ,将其代入梯形公式的右端,使之转化为显式公式,即 注: 当 f (x , y)关于y满足Lipschitz条件且步长h 满足 直至满足: 若采用梯形公式计算(★)中的积分项,则有 类似地,可得 (☆) 时,迭代格式 (☆) 收敛 . 3. 改进的Euler方法 把Euler法作为预报(称为预估公式),把隐式的梯形公式作为校正(称为校正公式 ),则得改进的Euler方法: 或 也称为预估-校正法. 有时为了方便,预估-校正格式也写成下面形式: 二、单步法的局部截断误差及精度 Def 1: 先假设 ,再估计误差 这种误差称为单步迭代法在 xk+1处的局部截断误差. Def 2: 若某种数值方法的局部截断误差为 ,则称该数值方法的精度为P 阶的. 注: 通常情况下,P 越大, h 越小,则截断误差越小,数值方法越精确. 设 10 .Euler方法是一阶方法. 所以Euler方法为一阶方法. 而 20 . 梯形法是二阶方法. Taylor展开 将 代入上式,得 而 代入上式得: 当h充分小时, 若 , 则可选取 h , 使得 故梯形法的精度为2 . 同样可以证明改进的Euler法也是二阶方法. 梯形法的局部截断误差为: 从而 例1: 取步长 h = 2/10, 2/20, 2/30, 2/40, 分别用欧拉法、改进的欧拉法和梯形法求解. 解: 记 f (x, y) = y－ x y2, xk= k h (k = 0, 1, 2,···, n ) (1) . Euler法: yk+1 =