龙格-库塔课件.pptVIP

* * * * * 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 * * * * 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 暮春三月，草长莺飞 2011-3-16 * 暮春三月，草长莺飞 3.3 龙格-库塔讲义 ——《数值分析简明教程》 主讲：王老师 黄淮学院数学科学系 2011-3-16 暮春三月，草长莺飞 * 考察差商 ，根据微分中值定理存在点 ，使得 从而利用所给方程 得 我们称 为区间 上的平均斜率K*，这样只要对平均斜率K*提供一种算法，相应地我们便导出一种计算格式。 龙格－库塔（Runge－Kutta）方法设计思想就是设法在 内多预报几个点的斜率值，然后把它们加权平均作为平均斜率，以期望构造出更高精度的计算格式。 2011-3-16 暮春三月，草长莺飞 * 这里运用了3.2节学过的改进的欧拉公式： 他用Xn与Xn+1 两点的斜率值K1和K2 的算术平均值作为平均斜率K* ，而Xn 处的斜率值 则利用已知信息 通过欧拉方法来预报。 K* 2011-3-16 暮春三月，草长莺飞 * 随意考察区间 内 一点 用两个点 的斜率的 的加权平均代替平均斜率 ,于是我们就得到如下计算格式: (类似于改进的欧拉公式) 其中有两个待定参数λ，p 。 2.二阶龙格-库塔方法 2011-3-16 暮春三月，草长莺飞 * 如下图 即 则 我们同样可以求出 2011-3-16 暮春三月，草长莺飞 * 适当选取它们的值，就可使上述格式有较高的精度。若 ，该格式是二阶的 ，故统称满足这一条件的一族格式为二阶龙格－库塔格式。特别地，当 时，上述格式即为改进的欧拉格式，如果 ，则上述格式称为变形的欧拉格式，亦称为中点格式。 2011-3-16 暮春三月，草长莺飞 * 假定 ，分别将 和 进行Taylor展开： 2011-3-16 暮春三月，草长莺飞 * 3.三阶龙格-库塔方法 为了进一步提高精度，我们可以考虑用三个点 的斜率值 加权平均得出平均斜率 的近似值，其中， 于是就可以构造所谓的三阶龙格－库塔格式 ，下列库塔格式是其中的一种： 2011-3-16 暮春三月，草长莺飞 * 4.四阶龙格-库塔方法 继续上述过程，我们可以导出四阶龙格－库塔格式，下列经典格式是其中的一种： 值得注意的是，龙格－库塔法的推导基于泰勒展开法，因而它要求解具有较好的光滑性。如果解的光滑性差，则该方法得到的解反而不好。 2011-3-16 暮春三月，草长莺飞 * 变步长的龙格-库塔方法 同积分的数值计算一样，微分方程的数值解法也需要选择步长。同样，我们可以采取步长加倍或折半的办法选择步长，即通过检查步长折半前后的两种计算结果的偏差： 来判断选取的步长是否合适，具体可以分为两种情况来处理：对于给定精度 ，若 ，则反复将步长折半进行计算直到 为止，取步长折半后的“新值”作为结果；相反的，反复将步长加倍直到 ，取步长加倍前的“老值”作为结果。 这种通过步长加倍或折半的手续处理步长的方法称为变步长方法。 * 例1. 使用高阶R-K方法计算初值问题 解: (1) 使用三阶R-K方法 RK.m * 其余结果如下: (2) 如果使用四阶R-K方法 n xn k1 k2 k3 yn 1.0000 0.1000 1.0000 1.1025 1.2555 1.1111 2.0000 0.2000 1.2345