chapter9常微分方程数值解法.ppt

2 3 32 3 31 32 2 , ( ) ( ) ( ) ( ) n n n n b h f x a h y h b b y x y x O h ? ? ? ?? ? ? ? ? ? ? ? ? ? ? 3 31 32 3 2 2 , ( ) [ ( , ) ( , ) ( ) ( , )] ( ) n n n n n x n n n y n n f x a h y hb y x hb f x y h h f x y y x f x y O h ? ? ? ? ? ? ? ? ? ? ? ? 取 3 31 32 1 1 , a b b ? ? ? 2 32 2 3 3 32 1 2 2 ( ) ( ) ( , ( ) ( ) ( )) ( ) n y n n n n h b y x f x h y hy x b h y x O h O h ?? ? ? ? ? ? ?? ? ? ? ?? 再利用 Taylor 展开 利用 Taylor 展开， 2 3 32 1 2 2 2 ( ) ( ) ( , ( )) ( ) n y n n n h h h b y x f x y y x O h ?? ? ? ? ? ? ? 3 4 32 2 12 2 2 ( ) ( ) ( , ( )) ( ) n y n n n h h h b y x f x y y x O h ?? ? ? ? ? ? ? 5 1 1 2 3 6 ( ) ( ) n n h h y x y O h ? ? ? ? ? ? ? ?? ? ? ( ) ( ) 代入 当 时， 32 2 b ? 4 1 1 ( ) ( ) n n y x y O h ? ? ? ? 第九章 常微分方程 数值解法 /*Numerical Method for Ordinary Differential Equations*/ 本章主要介绍 一阶 常微分方程初值问题的 数值解法 。 ? 初值问题及其 数值解 的概念 § 1 引言 常用的一些 解析解法 ： Lapalace 变换等 分离变量 法、 常数变易 法、 一阶 常微分方程初值问题： 0 ( , ), ( ) dy f x y a x b dx y a y ? ? ? ? ? ? ? ? ? ( ) ? 9 1 1 . . Th ( , ) f x y 对于初值问题 ，如果 在下列区域内连续： ( ) ? （解的 存在唯一 性） { ;| | } G a x b y ? ? ? ? ? 且关于 满足 Lipschitz 条件，即存在常数 ，使 y 0 L ? 1 2 1 2 | ( , ) ( , ) | | |; , f x y f x y L y y x y G ? ? ? ? ? 则初值问题 存在唯一解，且解是 连续可微 的。 ( ) ? 1 Def 所谓 数值解 是指：在解的 存在区间 上取一系列点 0 1 2 ... n a x x x x b ? ? ? ? ? ? 逐个求出 的近似值 1 2 ( , ,..., ) i y i n ? ( ) i y x 0 ; i x x ih ? ? 等距 节点： : h 步长 ? 初值问题 的解析解及其数值解的 几何 意义： ( ) ? o x y 初值问题 的解表示过点 的一条 曲线 ( ) ? 0 0 ( , ) x y n x ? ( , ) n n x y ? ( ) y y x ? 0 x ) , ( 0 0 y x ? 2 x ) , ( 2 2 y x ? ? ? 1 x ) , ( 1 1 y x ? 初值问题 的数值解表示一组 离散点列 ( ) ? ( , ) i i x y 可用 拟合 方法求该组数据 的 近似曲线 ( , ) i i x y 积分 曲线 § 2 Euler 方法 ? Euler 方法的导出 2 1 2 ( ) ( ) ( ) ( ) ! n n n n h y y x y x hy x ? ? ?? ? ? ? ? 将 在点 处进行 Taylor 展开 1 ( ) n y x ? n x 略去 项： 2 h 然后用 代替 ，即得 n y ( ) n y x 1 ( ) ( ) ( , ( )) n n n n y x y x hf x y x ? ? ? 1