《地球物理计算方法》习题课（二）.pptVIP

*预报改进改进校正*预报改进校正改进1、Euler及其改进方法2、Runge-Kutta方法3、Adams方法4、收敛性与稳定性5、方程组与高阶方程的情况6、边值问题**1、收敛性定义：对于任何固定的，当步长，有，则称此方法收敛；即截断误差随步长减小趋于0；收敛性讨论的是截断误差的性质。*记：en=|y(xn)-yn|,L是f(x)关于y的李普希兹常数，经反复递推，可得其中C,T为常数。若初始y0是准确的，即e0=0，则当h-0时，有en-0。欧拉方法是收敛的。来讨论.（稳定性讨论的是舍入误差δ的性质）2、稳定性*设模型方程Euler公式：设节点值上有大小为的扰动，此误差的传播使节点值产生大小为的扰动值。1、显式Euler格式的稳定性即所以要使必须满足条件误差方程模型方程的计算公式为设节点值上有大小为的扰动，此误差的传播使节点值产生大小为的扰动值，若计算过程不再引进新的误差，则即2、隐式Euler格式的稳定性*由于,则恒有,故恒有因此，隐式Euler格式是绝对稳定的（无条件稳定的）（对任何h0