第二节 龙格库塔方法.pptVIP

可将改进的欧拉格式改写成 的算术平均值作为平均斜率。 该公式可以看作是用 和 两个点处的斜率 和 由改进型欧拉公式我们可以猜想，如果在 内多预测几个点的斜率，再对他们进行加权平均， 可能得到精度更好的平均斜率！ 下面以2阶龙格-库塔方法为例来阐述这种思想 考察区间 上的一点 ， 用 和 的斜率 和 的加权平均作为平均 斜率 的近似值： 即取 其中 和 是待定常数。若取 ，则 问题在于如何确定 处的斜率 和常数 和 。 仿照改进的欧拉方法，用欧拉方法预测 的值， 并用它来估计斜率 ： 于是得到如下形式的算法： 通过适当选取参数 和 的值，使得公式具有 2阶精度！！ 由泰勒公式展开，要使公式具有 2 阶精度，只需 方程组有无穷多解：二级方法有无穷多种 常见的3种二级方法： ? 中点法（修正的Euler法） 取 ? 二阶龙格库塔方法 取 ?三级方法：N = 3 类似于N = 2的推导方法，可得到 常见的2种三阶方法： ? 库塔三阶方法 ? 四级方法：N = 4 局部截断误差 常见的2种四阶方法： ?经典龙格-库塔方法 解： 例2：用经典的龙格-库塔方法 求解下列初值问题 。 经典的四阶龙格-库塔公式： 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 二、高阶和隐式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阶 (1)一级二阶的隐式中点方法： (2)二级四阶的隐式R-K方法： 三、变步长方法 基本思想：根据精度自动地选择步长 对于经典Runge-Kutta方法： Step1:设从 出发，以 为步长，经过一步计算得到 Step2:取 为步长，再从 出发，经过两步计算得到 记 如果 ，则将步长折半进行计算，直到 为止 此时取 为最终结果； 如果 ，则将步长加倍进行计算，直到 为止 此时将步长折半一次计算，得到的为最终结果。 一、收敛性 /*Convergence*/ §3 单步法的收敛性、相容性和绝对稳定性 对于初值问题 的一种 单步法 产生的近似解，如果 对于任一固定的 ，均有 ， 则称该单步法是收敛的。 类似地可以定义隐式单步法、多步法（§4）的收敛性 设初值问题（*）对应的下列单步法是 阶的， 且函数 满足对 的Lipschitz条件，即存在常数 则该单步法是收敛的，且 证明： 记 由截断误差的定义 因为单步法是 阶的： 满足 其中 二、相容性 /*Consistency*/ 对于 阶方法： 若方法（**）的增量函数满足： 则称该方法与初值问题（*）相容。 设方法（**）与初值问题（*）相容，且 满足L-条件， 则该方法（**）是收敛的，即当 固定， 时 再由相容性得： 上式说明：当 时