2019年最新-常微分方程数值解-精选文档.ppt

第九章 常微分方程数值解 ? 考虑 一阶 常微分方程的 初值问题 ? dy ? ? f ( x , y ) x ? [ a , b ] ? dx ? ? y ( a ) ? y 0 只要 f ( x , y ) 在 [ a , b ] ? R 1 上连续，且关于 y 满足 Lipschitz 条 f ( x , y ) ? f ( x , y ) | ? L | y ? y | 件 ，即存在与 x , y 无关的常数 L 使 | 1 2 1 2 对任意定义在 [ a , b ] 上的 y 1 ( x ) 和 y 2 ( x ) 都成立，则上述问题 存 在唯一解 。 要计算出解函数 y ( x ) 在一系列节点 a = x 0 x 1 … x n = b 处的近似值 y y ( x ) ( i ? 1 , ... , n ) i ? i h x ? x i ? 0 , ... , n ? 1 ) 节点间距 为步长，通常采用 等距节点 ， i ? i ? 1 i ( 即取 h i = h ( 常数 ) 。 在这些节点上采用离散化方法，（通常用数值积分、微分、 泰勒展开等）将上述初值问题化成关于离散变量的相应问题。 把这个相应问题的解 y n 作为 y ( x n ) 的近似值。这样求得的 y n 就是 上述初值问题在节点 x n 上的数值解。一般说来，不同的离散化 导致不同的方法。 9.1 欧拉 Euler 法与改进欧拉法 1. 欧拉法： y ( x ) ? y ( x ) 1 0 ? 向前差商近似导数 y ( x ) ? 0 h 亦称为 欧拉折线法 记为 ? y ( x ) ? y ( x ) ? h y ( x ) ? y ? h f ( x , y ) 1 0 0 0 0 0 y 1 x 0 x 1 y ? y ? h f ( x , y ) ( i ? 0 , ... , n ? 1 ) i ? 1 i i i 定义 在假设 y i = y ( x i ) ，即第 i 步计算是精确的前提下，考 虑的截断误差 R i = y ( x i +1 ) ? y i +1 称为 局部截断误差 定义 若某算法的局部截断误差为 O ( h p +1 ) ，则称该算法有 p 阶精度。 R i 的 主项 ? 欧拉法的局部截断误差： 2 h i 2 i 3 ? ? ? R ? y ( x ) ? y ? [ y ( x ) ? h y ( x ) ? y ( x ) ? O ( h )] ? [ y ? hf ( x , y )] i i ? 1 i ? 1 i i i i 3 ? ? ? y ( x ) ? O ( h ) 欧拉法具有 1 阶精度。 i 2 h 2 例 9.1 用欧拉法求初值问题 0 . 9 ? y ? y ? ? 1 ? 2 x ? ? y ( x 0 ) ? 1 x 0 ? 0 ? 当 h = 0.02 时在区间 [0, 0.10] 上的数值解。 ( x ) ? ( 1 ? 2 x ) 方程真解： y ? 0 . 45 n 0 1 2 3 4 5 x n 0 0.02 0.04 0.06 0.08 0.10 y n 1.0000 0.9820 0.9650 0.9489 0.9336 0.9192 y ( x n ) 1.0000 0.9825 0.9660 0.9503 0.9354 0.923 ? n = y ( x n ) - y n 0 0.0005 0.0005 0.0014 0.0018 0.0021 2 ． 改进欧拉法 一阶方程的初值问题与积分方程 y ( x ) ? y f ( t , y ( t )) dt 0 ? ? x 0 x 是等价的，当 x = x 1 时， y ( x ) ? y f ( t , y ( t )) dt 1 0 ? ? x 0 x 1 借助于数值积分，求 y ( x 1 ) 的值 用矩形公式