第二节 龙格-库塔方法[宣讲].pptVIP

由泰勒公式展开，要使公式具有 2 阶精度，只需 方程组有无穷多解：二级方法有无穷多种 常见的3种二级方法： ? 中点法（修正的Euler法） 取 ? 二阶龙格库塔方法 取 * 精品PPT | 实用可编辑 ?三级方法：N = 3 类似于N = 2的推导方法，可得到 常见的2种三阶方法： ? 库塔三阶方法 * 精品PPT | 实用可编辑 ? 四级方法：N = 4 局部截断误差 常见的2种四阶方法： ?经典龙格-库塔方法 * 精品PPT | 实用可编辑 解： 例2：用经典的龙格-库塔方法 求解下列初值问题 。 经典的四阶龙格-库塔公式： * 精品PPT | 实用可编辑 0.1 0.2 0.3 0.4 0.5 1.0954 1.1832 1.2649 1.3416 1.4142 0.6 0.7 0.8 0.9 1.0 1.4832 1.5492 1.6125 1.6733 1.7321 同保留5位的精确值完全一致： 0.1 0.2 0.3 0.4 0.5 1.0954 1.1832 1.2649 1.3416 1.4142 0.6 0.7 0.8 0.9 1.0 1.4832 1.5492 1.6125 1.6733 1.7321 * 精品PPT | 实用可编辑 * 精品PPT | 实用可编辑 二、高阶和隐式Runge-Kutta方法 注：?对于显式N级R-K方法，最多只能得到N级方法； N 1,2,3,4 5,6,7 8,9 10,11,… N N-1 N-2 ?已经证明N级R-K方法的阶 具有下列关系： ?若要得到N阶以上方法，则使用N级隐式R-K方法 ?N级隐式R-K方法的一般形式： N级隐式R-K法 可以达到2N阶 * 精品PPT | 实用可编辑 (1)一级二阶的隐式中点方法： (2)二级四阶的隐式R-K方法： * 精品PPT | 实用可编辑 三、变步长方法 基本思想：根据精度自动地选择步长 对于经典Runge-Kutta方法： Step1:设从 出发，以 为步长，经过一步计算得到 Step2:取 为步长，再从 出发，经过两步计算得到 * 精品PPT | 实用可编辑 记 如果 ，则将步长折半进行计算，直到 为止 此时取 为最终结果； 如果 ，则将步长加倍进行计算，直到 为止 此时将步长折半一次计算，得到的为最终结果。 * 精品PPT | 实用可编辑 一、收敛性 /*Convergence*/ §3 单步法的收敛性、相容性和绝对稳定性 对于初值问题 的一种 单步法 产生的近似解，如果 对于任一固定的 ，均有 ， 则称该单步法是收敛的。 类似地可以定义隐式单步法、多步法（§4）的收敛性 * 精品PPT | 实用可编辑 设初值问题（*）对应的下列单步法是 阶的， 且函数 满足对 的Lipschitz条件，即存在常数 则该单步法是收敛的，且 证明： 记 由截断误差的定义 * 精品PPT | 实用可编辑 因为单步法是 阶的： 满足 其中 * 精品PPT | 实用可编辑 二、相容性 /*Consistency*/ 对于 阶方法： 若方法（**）的增量函数满足： 则称该方法与初值问题（*）相容。 * 精品PPT | 实用可编辑 设方法（**）与初值问题（*）相容，且 满足L-条件， 则该方法（**）是收敛的，即当 固定， 时 再由相容性得： 上式说明：当 时，方法（**）趋于原微分方程 本章讨论的数值方法都是与原初值问题相容的 * 精品PPT | 实用可编辑 三、绝对稳定性 /*Absolute Stibility*/ 计算过程中产生的舍入误差对计算结果的影响 首先以Euler公式为例，来讨论一下舍入误差的传播: 设实际计算得到的点 的近似函数值为 ， 其中 为精确值， 为误差 如果 ，则误差