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

对于预报-校正系统中预报与校正两种格式，他们的误差可以事后估计得到，将每次预报和校正值加上事后估计的误差，从而提高代数精度，因而基于这种误差的事后估计可以进一步优化预报校正系统.：*预报改进校正改进估算的误差估算的误差1、Euler及其改进方法2、Runge-Kutta方法3、Adams方法4、收敛性与稳定性5、方程组与高阶方程的情况6、边值问题**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格式的稳定性即所