八 常微分方程初值问题的数值解2.ppt

* § 6.5 Runge-Kutta(龙格-库塔)方法 在高精度的单步法中,应用最广泛的是Runge-Kutta(龙格-库塔)方法 一、Runge-Kutta法的基本思想（1） Runge-Kutta法的基本思想（2） Runge-Kutta法的基本思想（3） 二、二阶龙格－库塔方法 三、三阶龙格－库塔方法 四、四阶龙格－库塔方法 解： 例2：用经典的Runge-Kutta方法 求解下列初值问题 。 经典的四阶Runge-Kutta公式： 1.3416 0.4 1.1832 0.2 1.0954 0.1 1.2649 0.3 1.4142 0.5 1.6733 0.9 1.5492 0.7 1.4832 0.6 1.6125 0.8 1.7321 1.0 同保留5位的精确值完全一致： 1.3416 0.4 1.1832 0.2 1.0954 0.1 1.2649 0.3 1.4142 0.5 1.6733 0.9 1.5492 0.7 1.4832 0.6 1.6125 0.8 1.7321 1.0 二、高阶和隐式Runge-Kutta方法 注：?对于显式N级R-K方法，最多只能得到N级方法； N-2 N-1 N 10,11,… 8,9 5,6,7 1,2,3,4 N ?已经证明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条件，即存在常数 则该单步法是收敛的，且 证明： 记 由截断误差的定义 因为单步法是 阶的： 满足 其中 二、绝对稳定性 /*Absolute Stibility*/ 计算过程中产生的舍入误差对计算结果的影响 首先以Euler公式为例，来讨论一下舍入误差的传播: 设实际计算得到的点 的近似函数值为 ， 其中 为精确值， 为误差 如果 ，则误差是不增的，故可认为是稳定的