《地球物理计算方法》第三章 微分方程差分法——Euler及其改进方法、Runge-Kutta方和法Adams方法.pptVIP

***采用预报-校正系统：（2）将近似值代入公式作迭代计算： （1）求得一个初步的近似值yn+p-得到斜率k2： 其中为待定参数希望合理的确定λ、p，以提高精度。*确定λ、p：（教材P103）*定理：当时，Runge-Kuta格式具有二阶精度。特别地，当时，有改进欧拉格式采用精度分析方法：对y=x2能够精确成立。*特别地，当时，有变形的欧拉格式*3、三阶Runge-Kutta有xn+p,xn+q。令xn，xn+p，xn+q处的斜率值分别为k1,k2和k3。其中：k1=f(xn,yn);条件问题如何确定k2,k3，使其加权平均得到近似的K*，使得所构造的格式具有三阶精度?*加权平均的k*，计算格式为：其中为待定参数。合理的确定，以提高精度。*采用逐步预报方式：首先k1,k2预报平均yn+q斜率(k1,k2用前面的方法计算）:得到采用泰勒展开方法（或代数精度方法）选择参数p,q,α,λ,μ，保证上述格式具有三阶精度。*根据预报校正，得到下列差分格式：三阶Runge-Kutta格式4、四阶Runge-Kutta