《地球物理计算方法》第三章 微分方程差分法——Runge-Kutta方法.pptVIP

《地球物理计算方法》

第3章微分方程差分法——Runge-Kutta方法1、设计思想微分中值定理如果在上连续，在内可导，则在内至少存在一点，有于是有其中称作区间上的平均斜率。差商的定义：*问题：计算近似值的关键是如何选择算法确定平均斜率k*？对于欧拉格式精度低!改进的欧拉格式于是：其中预报-校正*Runge-Kutta算法的基本思想在区间上取几个点的斜率，然后计算它们的加权平均斜率k*，以此可以构造出更高精度的格式，这就是Runge-Kutta算法的基本思想。2、二阶Runge-Kutta有令和处的斜率值分别为和。其中：条件于是，有二阶Runge-Kutta格式yn+p如何确定？*采用预报-校正系统：（2）将近似值代入公式作迭代计算： （1）求得一个初步的近似值yn+p-得到斜率k2： 其中为待定参数希望合理的确定λ、p，以提高精度。*确定λ、p：（教材P103）*定理：当