《地球物理计算方法》第三章 微分方程差分法——收敛性与边值问题.pptVIP

《地球物理计算方法》

第3章微分方程差分法—收敛性与边值问题1、收敛性与稳定性2、方程组与高阶方程的情况3、边值问题**1、收敛性定义：对于任何固定的，当步长，有，则称此方法收敛；即截断误差随步长减小趋于0；收敛性讨论的是截断误差的性质。*以显式欧拉公式为例：*李普希兹（Lipshitz）条件（解存在且唯一条件）对于边值问题：如果存在实数L0，使得f关于L满足Lips条件，L为f的Lips常数。**记：en=|y(xn)-yn|,L是f(x)关于y的李普希兹常数，经反复递推，可得其中C,T为常数。若初始y0是准确的，即e0=0，则当h-0时，有en-0。欧拉方法是收敛的。来讨论.（稳定性讨论的是舍入误差δ的性质）2、稳定性*设模型方程Euler公式：设节点值上有大小为的扰动，此误差的传播使节点值产生大小为的扰动值。1、显式Euler格式的稳定性即所以要使必条件满足误差方程模型方程的计算公式为设节点值上有大小为的扰动，此误差的传播使节点值产生大小为的扰动值，若计算过程不再引进新的误差，则即2、隐式Euler格式的稳定性*由于