二-龙格库塔方法PPT课件.ppt

基本思想：利用 在某些特殊点上的函数值的线性组合来构造高阶单步法的平均斜率。第二节 龙格-库塔法什么叫平均斜率？对差商 应用微分中值定理，有，利用微分方程 ，有这里的 称为平均斜率。1.可将改进的欧拉格式改写成的算术平均值作为平均斜率。该公式可以看作是用 和 两个点处的斜率 和由改进型欧拉公式我们可以猜想，如果在内多预测几个点的斜率，再对他们进行加权平均，可能得到精度更好的平均斜率！2.下面以2阶龙格-库塔方法为例来阐述这种思想考察区间 上的一点 ，用 和 的斜率 和 的加权平均作为平均斜率 的近似值：即取其中 和 是待定常数。若取 ，则问题在于如何确定 处的斜率 和常数 和 。3.仿照改进的欧拉方法，用欧拉方法预测 的值，并用它来估计斜率 ：于是得到如下形式的算法：通过适当选取参数 和 的值，使得公式具有2阶精度！4.由泰勒公式展开，要使公式具有 2 阶精度，只需方程组有无穷多解：二级方法有无穷多种常见的3种二级方法：中点法（修正的Euler法）取 二阶龙格库塔方法取5.三级方法：N=3 类似于N=2的推导方法，可得到常见的2种三阶方法：库塔三阶方法6.四级方法：N=4局部截断误差常见的2种四阶方法：经典龙格-库塔方法7.解：例2：用经典的龙格-库塔方法求解下列初值问题 。经典的四阶龙格-库塔公式：8.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.73219.10.二、高阶和隐式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阶11.(1)(1)一级二阶的隐式中点方法：(2)(2)二级四阶的隐式R-K方法：12.三、变步长方法基本思想：根据精度自动地选择步长对于经典Runge-Kutta方法：Step1:设从 出发，以 为步长，经过一步计算得到Step2:取 为步长，再从 出发，经过两步计算得到13.记如果 ，则将步长折半进行计算，直到 为止此时取 为最终结果；如果 ，则将步长加倍进行计算，直到 为止此时将步长折半一次计算，得到的为最终结果。14.一、收敛性/*Convergence*/3 单步法的收敛性、相容性和绝对稳定性对于初值问题 的一种单步法 产生的近似解，如果 对于任一固定的 ，均有 ，则称该单步法是收敛的。类似地可以定义隐式单步法、多步法（4）的收敛性15.设初值问题（*）对应的下列单步法是 阶的，且函数 满足对 的Lipschitz条件，即存在常数则该单步法是收敛的，且证明：记由截断误差的定义16.因为单步法是 阶的：满足其中17.二、相容性/*Consistency*/对于 阶方法：若方法（*）的增量函数满足：则称该方法与初值问题（*）相容。18.设方法（*）与初值问题（*）相容，且 满足L-条件，则该方法（*）是收敛的，即当 固定，时再由相容性得：上式说明：当 时，方法（*）趋于原微分方程本章讨论的数值方法都是与原初值问题相容的 19.三、绝对稳定性 /*Absolute Stibility*/计算过程中产生的舍入误差对计算结果的影响首先以Euler公式为例，来讨论一下舍入误差的传播:设实际计算得到的点 的近似函数值为 ，其中 为精确值，为误差如果 ，则误差是不增的，故可认为是稳定的20.例如：对于初值问题精确解为而实际求解的初值问题为精确解为在 处的误差为可见误差随着 的增加呈指数函数增长如果初值问题为精确解为21.实际求解的初值问题为精确解为在 处的误差为可见误差随着 的增加呈指数函数递减当 时，微分方程是不稳定的；而 时，微分方程是稳定的。上面讨论的稳定性，与数值方法和方程中 有关22.实验方程：对单步法 应用实验方程，如果 ，当 时，则称该单步法是绝对稳定的，在复平面上复变量 满足的区域，称为该单步法的绝对稳定域，它与实轴的交集称为绝对稳定区间。若单步法是 阶的，则由实验方程可得：23.例3：分别求Euler法和经典的R-K法的绝对稳定区间。解：Euler公式：将其应用于实验方程绝对稳定域：当 时，绝对稳定区间：经典的R-K公式：24.25.当 时，绝对稳定区间：可以证明：存在唯一极小值点由 得26.例4：求梯形公式（隐式方法）的绝对稳定区间。解：梯形公式：将其应用于实验方程当 时，绝对稳定区间：27.