2012数值分析 -lec25--Rung-Kutta方法.pptVIP

变步长Runge-Kutta方法 在单步法中每一积分步步长实际上是相互独立的，步长的 选择具有了灵活性。根据合理地选择每一积分步的步长， 既保证精度的要求，又可以减少计算量，从而减小舍入误 差。其方便的控制手段是基于误差的事后估计式。 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第二十五讲Runge-Kutta方法 第七章 常微分方程初值问题的数值解法 常微分方程及其求解概述 初值问题解法 1.Euler方法 2.线性多步法 3.Runge—Kutta方法 4.方程组及刚性问题的Gear方法 Runge--Kutta 方法 六个未知数，二个自由，故可取 故： N=4: 四级四阶R--K方法 ----最常用的古典Runge-Kutta方法。 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ? 由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 表 8 各种解法在例题中的结果比较 （2） 单步法，自动起步 （3） 易改为变步长 （4） 绝对稳定区域较同阶线性多步法大 （5） 计算工作量较大，有时大于隐式方法 绝对稳定性讨论 表 9 各级R-K方法的绝对稳定区域 R-K法的绝对稳定区域 表 10 不同步长对精度的影响 对于给定的精度 ,如果 ,反复将步长折半进行计算，直至 为止,这时取最终得到的作为结果； 如果 ，将反复将步长加倍，直到 为止，这时再将步长折半一次，就得到所要的结果。 这种通过加倍或折半处理步长的计算方法称为变步长方法。 注：推荐使用精度好计算量低的变步长方法。 用四阶显式R-K方法做变步长方法是实践中较好的方法！ 变步长Runge--Kutta方法框图 出发值的计算 作业 课后习题5、7、9、10、11 Thanks for your attention! Any questions?! Xm Euler 改进Euler R—K(4阶) 精确解 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.0000 1.2000 1.4420 1.7384 2.1041 2.5569 3.1183 3.8139 4.6747 5.7376 7.0472 1.0000 1.2210 1.4923 1.8284 2.2466 2.7680 3.4174 4.2257 5.2288 6.4704 8.0032 1.0000 1.2221 1.4977 1.8432 2.2783 2.8274 3.5201 4.3927 5.4894 6.8643 8.5834 1.0000 1.2221 1.4977 1.8432 2.2783 2.8274 3.5202 4.3928 5.4895 6.8645 8.5836 级 (1 D 1 2 3 4 (-2,0) (-2,0) (-2.51,0) (-2.78,0) x h1时的误差 h2时的误差 0.0 0.2 0.4 0.6 0.8 1.0 0 -0.092795 -0.012010 -0.001366 -0.000152 -0.000017 0 4.98 25.0 125.0 625.0 3125.0