[第9章微分方程初值问题的数值解法-1.pptVIP

* 第九章 微分方程初值问题的数值解法 内容提纲 引言 Euler法及其改进 Runge-Kutta方法 线性多步法 误差分析 数值解法的收敛性、相容性和稳定性 边值问题数值解法简介 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 引言 初值问题的数值解法:求初值问题的解在一系列节点的值 y ( xn )的近似值 yn 的方法.本章数值解法的特点:都是采用“步进式”,即求解过程顺着节点排列的次序一步步向前推进. 基本知识: (1) 定理1: 如果函数 f (x , y)在区域 上连续,且关于 y 满足Lipschitz条件 常微分方程初值问题: 求未知函数 y= y (x) . Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 此时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)是稳定的. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 单步迭代: 计算 yn+1时仅用 yn ; 初值问题(1)与下列积分方程的解等价: 初值问题的数值解就是求一系列节点 上函数 y= y (x)的近似值 . 称为步长. 一般取等步长 h . 多步迭代: 计算 yn+1时除用 yn 外, 还要用到 yn-1, yn-2 ,…; k 步迭代要用到 yn-1, yn-2 ,…, yn-k+1 . 显式单步迭代: 隐式单步迭代: (2) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 一、Euler方法及其改进 将[a , b ]n 等分, 记 微分法: 积分法: 积分项利用矩形公式计算 1. 显式Euler方法 (★) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Taylor公式推导: Euler公式几何意义: P1 P2 Pk 也称折线法 P0 x y Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 2. 梯形法 称之为梯形公式.这是一个隐式公式,通常用迭代法求解.具体做法: 取 先用Euler法求出初值 ,即 ,将其代入梯形公式的右端,使之转化为显式公式,即 注: 当 f (x , y)关于y满足Lipschitz条件且步长h 满足 直