第四章常微分方数值解法.pptxVIP

1 四　常微分方程数值解法 2 常微分方程数值解法 引言（常微分方程数值解法概述） 显式欧拉法、隐式欧拉法、二步欧拉法 局部截断误差与精度 改进的欧拉方法 龙格-库塔方法 收敛性与稳定性简述 一阶常微分方程组与高阶常微分方程 3 引言 一阶常微分方程初值问题： 上连续，且在 R 域内满足李普希兹 (Lipschitz) 条件，即存在正数 L，使得对于 R 域内的任意两值 y1, y2，下列不等式成立： 则上述初值问题的连续可微的解 y(x) 存在并且唯一。 4 引言（续） 实际生产与科研中，除少数简单情况能获得初值问题的初等解（用初等函数表示的解）外，绝大多数情况下是求不出初等解的。 有些初值问题即便有初等解，也往往由于形式过于复杂而不便处理。 实用的方法是在计算机上进行数值求解：即不直接求 y(x) 的显式解，而是在解所存在的区间上，求得一系列点 xn (n = 0, 1, 2, …) 上解的近似值。 5 欧拉（Euler）方法 方法一　化导数为差商的方法 由于在逐步求解的过程中，y(xn) 的准确值无法求解出来，因此用其近似值代替。 为避免混淆，以下学习简记： y(xn)：待求函数 y(x) 在 xn 处的精确函数值 yn ：待求函数 y(x) 在 xn 处的近似函数值 6 代入初值问题表达式可得： 根据 y0 可以一步步计算出函数 y = y(x) 在 x1