[工学]第四章常微分方程数值解法.ppt

1 四　常微分方程数值解法 2 常微分方程数值解法 引言（常微分方程数值解法概述） 显式欧拉法、隐式欧拉法、二步欧拉法 局部截断误差与精度 改进的欧拉方法 龙格-库塔方法 收敛性与稳定性简述 一阶常微分方程组与高阶常微分方程 3 引言 一阶常微分方程初值问题： 上连续，且在 R 域内满足李普希兹 (Lipschitz) 条件，即存在正数 L，使得对于 R 域内的任意两值 y1, y2，下列不等式成立： 则上述初值问题的连续可微的解 y(x) 存在并且唯一。 4 引言（续） 实际生产与科研中，除少数简单情况能获得初值问题的初等解（用初等函数表示的解）外，绝大多数情况下是求不出初等解的。 有些初值问题即便有初等解，也往往由于形式过于复杂而不便处理。 实用的方法是在计算机上进行数值求解：即不直接求 y(x) 的显式解，而是在解所存在的区间上，求得一系列点 xn (n = 0, 1, 2, …) 上解的近似值。 5 欧拉（Euler）方法 方法一　化导数为差商的方法 由于在逐步求解的过程中，y(xn) 的准确值无法求解出来，因此用其近似值代替。 为避免混淆，以下学习简记： y(xn)：待求函数 y(x) 在 xn 处的精确函数值 yn ：待求函数 y(x) 在 xn 处的近似函数值 6 代入初值问题表达式可得： 根据 y0 可以一步步计算出函数 y = y(x) 在 x1, x2, x3 x4, …上的近似值 y1, y2, y3, y4 , … 常微分方程数值解是一组离散的函数值数据，它的精确表达式很难求解得到，但可以进行插值计算后用插值函数逼近 y(x) 7 欧拉方法（续） 方法二　数值积分法 同样以近似值 yn 代替精确值 y(xn) 可得： 将微分方程 y  = f (x, y) 在区间 [xn, xn+1] 上积分： 8 欧拉方法的几何意义 x y 0 9 隐式欧拉法 在数值积分法推导中，积分的近似值取为积分区间宽度与右端点处的函数值乘积，即： 这样便得到了隐式欧拉法： 隐式欧拉法没有显式欧拉法方便 10 二步欧拉法 在数值积分法推导中，积分区间宽度选为两步步长，即积分区间为：[xn-1, xn+1]，则： 以 y(x) 在 xn -1, xn 上的近似值代替精确值可得： 中矩形公式 11 梯形公式欧拉法 在数值积分法中，如果用梯形公式近似计算 f (x, y) 在区间 [xn, xn+1] 上的积分，即： 用近似值代替精确值可得梯形公式欧拉法： 上式右端出现了未知项，可见梯形法是隐式欧拉法的一种；实际上，梯形公式欧拉法是显式欧拉法与隐式欧拉法的算术平均。 12 例 用显式欧拉法、隐式欧拉法、梯形法求解初值问题： 取 h = 0.1，计算到 x = 0.5，并与精确解进行比较 解：由已知条件可得：h = 0.1，x0 = 0， y0 = 1， f (x, y) = - y + x + 1 13 例：（续） 14 计算结果： xn 显式法 yn 隐式法 yn 梯形法 yn 精确解 y (xn) 0.0 1 1 1 1 0.1 1.000000 1.009091 1.004762 1.004837 0.2 1.010000 1.026446 1.018594 1.019731 0.3 1.029000 1.051315 1.040633 1.040818 0.4 1.056100 1.083014 1.070097 1.070320 0.5 1.090490 1.120922 1.106278 1.106531 15 局部截断误差 为了简化分析某常微分方程数值算法的误差，现假设 yn = y(xn)，即在前一步 yn 准确的前提下，估计： 称上述误差 Tn+1 为该常微分方程数值算法的局部截断误差 如果某个常微分方程数值算法的局部截断误差可表示为 O(h p+1)，则称该数值算法的精度是 p 阶 欧拉法的精度为一阶；二步欧拉法的精度为二阶；梯形公式欧拉法的精度为二阶。 16 泰勒展开法 如果初值问题中的 f (x, y) 充分可微，则可将 y(xn+1) 在点 xn 处展开： 如果只保留线性项，忽略 h2 及以上各项，则： 17 局部截断误差的分析 利用泰勒公式展开，比较各算法与展开式的前几项 将 y(xn+1) 在 xn 点处用泰勒公式展开： 显式欧拉法的局部截断误差： 1 阶精度 18 补充：二元函数微分中值定理 19 1 阶精度 20 分别将 y(xn+1), y(xn-1) 在 xn 点处用泰勒公式展开： 二步欧拉法的局部截断误差： 2 阶精度 21 2 阶精度 22 各种欧拉法的比较 方法 精度 评述 显式欧拉法 1 最简单，精度低 隐式欧拉法 1 不便计算，稳定性好 二步欧拉法 2 需要