第三章常微分方程差分方法--自控1401-05课件.ppt

* 补充：二元泰勒展开式 * 用二元泰勒公式展开 将 k1, k2 代入 中可得： * 二阶龙格-库塔方法 2 阶精度 * 四个未知变量，只有三个方程，有无穷多组解 每组解构成的龙格-库塔方法均为二阶 二阶龙格-库塔方法即为改进的欧拉方法 变形的欧拉法 中 点 方 法 * 三阶龙格-库塔方法 三阶龙格-库塔方法是用三个值 k1, k2, k3 的加权平均来近似取代 k* 要使三阶龙格-库塔方法具有三阶精度，必须使其局部截断误差为 O(h4) 将 k1, k2, k3 代入 yn+1 的表达式中，在 (xn, yn) 处用二元泰勒公式展开，与 y(xn+1) 在 xn 处的泰勒展开式比较 * 三阶龙格-库塔方法 类似二阶龙格-库塔方法的推导过程，8 个待定系数 c1, c2, c3, a2, a3, b21, b31, b32 应满足： 8 个未知参数，6 个方程，有无穷多组解 库塔公式 * 四阶龙格-库塔方法 类似可以推出四阶龙格-库塔公式，常用的有： 标准四阶龙格-库塔公式 * 四阶龙格-库塔方法 吉尔（Gill）公式 4 阶以上龙格-库塔公式的计算量太大，并且精度不一定提高，有时反而会降低，因此实际应用中一般选用四阶龙格-库塔已足可满足精度要求。 * 用经典四阶龙格-库塔方法求解前例的初值问题，并与改进欧拉法、梯形法在 x5 = 0.5 处比较其误差大小。 解：采用经典四阶龙格-库塔公式： 例： * 四阶R-K方法的精度比二阶方法高得多 精确解为： R-K方法的误差： 改进欧拉法的误差： 梯形法的误差： * * * * 第3章　常微分方程的差分方法 引言（常微分方程的差分方法概述） 显式欧拉法、隐式欧拉法、两步欧拉法 局部截断误差与精度 改进的欧拉方法 龙格-库塔方法 引言 一阶常微分方程初值问题： 微分方程 初始条件 定理：若 f (x, y) 在某闭区域 R ： 上连续，且在 R 域内满足李普希兹 (Lipschitz) 条件，即存在正数 L，使得对于 R 域内的任意两值 y1, y2，下列不等式成立： 则上述初值问题的连续可微的解 y(x) 存在并且唯一。 * * （显式）欧拉（Euler）方法 方法一　化导数为差商的方法 由于在逐步求解的过程中，y(xn) 的准确值无法求解出来，因此用其近似值代替。 为避免混淆，以下学习简记： y(xn)：待求函数 y(x) 在 xn 处的精确函数值 yn ：待求函数 y(x) 在 xn 处的近似函数值 * 代入初值问题表达式可得： 根据 y0 可以一步步计算出函数 y = y(x) 在 x1, x2, x3 x4, …上的近似值 y1, y2, y3, y4 , … 常微分方程数值解是一组离散的函数值数据，它的精确表达式很难求解得到，但可以进行插值计算后用插值函数逼近 y(x) * 方法二　数值积分法 同样以近似值 yn 代替精确值 y(xn) 可得： 将微分方程 y ? = f (x, y) 在区间 [xn, xn+1] 上积分： 左矩形公式 * 欧拉方法的几何意义 x y 0 * 隐式欧拉法 在数值积分法推导中，积分的近似值取为积分区间宽度与右端点处的函数值乘积，即： 这样便得到了隐式欧拉法： 含有未知的函数值 隐式欧拉法没有显式欧拉法方便 * 二步欧拉法 在数值积分法推导中，积分区间宽度选为两步步长，即积分区间为：[xn-1, xn+1]，则： 以 y(x) 在 xn -1, xn 上的近似值代替精确值可得： 需要前两步的计算结果 中矩形公式 * 梯形公式欧拉法 在数值积分法中，如果用梯形公式近似计算 f (x, y) 在区间 [xn, xn+1] 上的积分，即： 用近似值代替精确值可得梯形公式欧拉法： 上式右端出现了未知项，可见梯形法是隐式欧拉法的一种；实际上，梯形公式欧拉法是显式欧拉法与隐式欧拉法的算术平均。 * 例：用显式欧拉法、隐式欧拉法、梯形法求解初值问  题 取 h = 0.1，计算到 x = 0.5，并与精确解进行比较 解：由已知条件可得：h = 0.1，x0 = 0， y0 = 1， f (x, y) = - y + x + 1 显式欧拉法： * 隐式欧拉法： 化简得： 梯形公式欧拉法： * 计算结果： 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