第七章常微分程数值解1.docVIP

PAGE 第七章 常微分方程数值解 本章介绍求解微分方程数值解的基本思想和方法. 常微分方程 含有自变量、未知函数和它的一阶导数和高阶导数的方程. 它是描述运动、变化规律的重要数学方法之一，分为两类： 1.初值问题，即给出未知函数及导数在初始点的值 2.边值问题，即给出未知函数及（或）它的某些导数在区间两个端点的值 常微分方程初值问题 SKIPIF 1 0 其中 SKIPIF 1 0 为 SKIPIF 1 0 的已知函数， SKIPIF 1 0 为给定的初值. 这里仅讨论一阶标量微分方程初值问题的数值解法. 而高阶微分方程通常可化为一阶微分方程组来研究. 数值解法 寻求微分方程初值问题之解 SKIPIF 1 0 在一系列离散点 SKIPIF 1 0 上的近似值: SKIPIF 1 0 SKIPIF 1 0 的方法. SKIPIF 1 0 : 问题的数值解 数值解所满足的离散方程统称为差分格式. 步长： SKIPIF 1 0 ，一般取定步长 SKIPIF 1 0 初值问题的适定性(其解是否唯一存在) 记(带形)区域： SKIPIF 1 0 为 SKIPIF 1 0 ，即 SKIPIF 1 0 = SKIPIF 1 0 . 设 SKIPIF 1 0 为连续映射，若存在常数 SKIPIF 1 0 使得不等式 SKIPIF 1 0 对一切 SKIPIF 1 0 都成立，则称 SKIPIF 1 0 在 SKIPIF 1 0 上关于 SKIPIF 1 0 满足Lipschitz条件，而式中的常数L称为Lipschitz常数. 定理 初值问题 SKIPIF 1 0 ，当 SKIPIF 1 0 在 SKIPIF 1 0 上连续，且关于 SKIPIF 1 0 满足Lipschitz条件，则其解存在且唯一. § 7.1 Euler方法 Euler公式 将初值问题 SKIPIF 1 0 的求解区间 SKIPIF 1 0 SKIPIF 1 0 等分,分点: SKIPIF 1 0 ,其中 SKIPIF 1 0 将 SKIPIF 1 0 写成等价的积分方程形式： SKIPIF 1 0 , 在上式中令 SKIPIF 1 0 ，并用左矩形公式计算右端积分，得到 SKIPIF 1 0 , (*) SKIPIF 1 0 ——余项 将(*)中的余项 SKIPIF 1 0 截去，可得 SKIPIF 1 0 则有 SKIPIF 1 0 的近似值 SKIPIF 1 0 的递推公式 SKIPIF 1 0 (*1) Euler公式 SKIPIF 1 0 表示当 SKIPIF 1 0 为精确值时，利用(*1)即Euler公式计算 SKIPIF 1 0 时的误差. SKIPIF 1 0 为Euler方法的局部截断误差. 例1 用Euler公式解初值问题 SKIPIF 1 0 解：取 SKIPIF 1 0 ，Euler公式的具体形式为 SKIPIF 1 0 其中 SKIPIF 1 0 已知 SKIPIF 1 0 ，则有 SKIPIF 1 0 SKIPIF 1 0 … … … 依次计算可得 SKIPIF 1 0 其部分结果见下表 SKIPIF 1 0 数值解 SKIPIF 1 0 准确解 SKIPIF 1 0 局部截断误差 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 可见Euler方法的计算结果精度