《地球物理计算方法》第3章常微分方程的数值求法1.pptx

1;2;3;4;5;6;7;1、Euler及其改进方法;1、Euler方法;求解过程从边界条件y|x0=y0即y(x0)=y0开始：;11;12;13;精度分析;15;16;17;;19;;21;22;23;;;26;Euler方法分为：显式格式和隐式格式。 显式格式： Euler格式， Euler两步格式 特点：计算量小，稳定性差； 隐式格式：隐式Euler格式， 梯形格式 特点：迭代法求解，计算量大，稳定性好； 综合应用： 改进的Euler格式（预报校正系统） 即：显式格式得到预报值，隐式格式迭代一次得到校正值。 27;;1、设计思想;30;改进的欧拉格式;32;33;34;35;2、二阶Runge-Kutta;?y n ?1 ? y n ? h[(1? ? )k1 ? ?k2 ];38;39;40;41;42;43;44;;;例： 求解以下初值问题;48;;;51;52;53;Runge-Kutta方法的实质： 就是通过几个点上的预报斜率值(xn+p，xn+q)估算需求点xn+1的斜 率值。而没有用到已给点上的斜率值，从而造成了重复计算。;;;;58;59;60;61;62;63;64;65;66;67;68;;yn?1 ? (1? ?) pn?1 ? ?cn?1 考察其相应的近似关系式: y(x ) ? y(x ) ? (1? ?) h (3y?(x ) ? y?(x