微分方程数值解法第二章.pptVIP

则我们可以令 ，连同格式（2.29）则可沿着 方向逐层把结点上的 值计算出来，以此作微分方程解在结点 处的值 的近似值。对初边值问题，同样使用格式（2.29）连同初、边值条件把求解区域Ω中的网络结点上的 值计算出来，作为微分方程初边值问题的近似解，由于格式（2.29)关于 明显解出来，所以称为显式差分格式。格式（2.29）称为解热传导方程（2.26）的古典显式格式。显然，对于求解来说，这种格式是最方便的，然而下面看到这种格式并不总是稳定的。 应该指出，为了推导古典显式差分格式（2.29），可以简单地应用导数的差商近似表达式 代入微分方程(2.26),并令差分方程解为 即可。虽然在边界结点上，差分方程和微分方程具有相同的初值解或者初边值条件，但是，一般而言，结点 上微分方程的精确解 和古典显式差分方程（2.29）的精确解 不相等，它们之间的差记为 ，且 假定 具有下面推导中所需要的有界偏导数，则有taylor展开，有 （2.30） 则 （2.31） 由式（2.26），（2.29），（2.30），（2.31）得 （2.32） 这就是 满足的方程。从式（2.31）我们有 或 （2.33） 从而，上式右边量描写了古典显式差分格式（2.29）在 点对微分方程的近似程度，将其定义为差分格式点 的截断误差，记为 ，即 （2.34） 因为假定 ， 在所考察的区域保持有界，所以当 时， 。古典显式差分格式的截断误差为 。从式（2.33）又可见到，如令 ，因 故截断误差 的阶可以提高。事实上，这时 ，为了提高截断误差的阶，我们也可用在式（2.27）中保留四阶中心差分项的办法达到，这时有差分格式 （2.35） 或者 相应的截断误差阶 。通常，格式可用图2.3表示。 图2.3 2.2.3 系数依赖于 的一维热传导方程的显式格式 考虑系数依赖于 的一维热传导方程 （2.36） 这时， 。方程（2.24）变成 此处“”表示对 求导。将微分算子 用差分算子代替，则可得到相应的差分格式，上式保留右边前两项，由 ，则差分方程 （2.37） 这一差分格式可用图2.4表示，其中 ，这是一个显式差分格式，其截断误差阶为 。 进一步，考虑下面热传导方程的差分近似 （2.38） 由方程右边 类似式（2.37）的推导，在上式中保留前两项，并且 和 分别用 和 代替，则得差分方程 （2.39） 此处 均在 处计算。 也可通过直接用中心差分算子 代替微分算子 的办法获得方程（2.38）的差分近似 （2.40） 这也是一个显式差分格式。 格式（2.39）和（2.40）的截断误差阶都是 。易见，由 代入格式（2.40）即为格式（2.39），差分格式（2.40）的推导方法，即在微分方程中直接用差分算子代替 ，正如前面已经指出的是推导差分格式的一个常用方法。 显然，微分方程（2.36），（2.38）中的 如果为 ，即其自变量包括空间变量和时间变量，这时差分格式（2.37），（2.39），（2.40）同样是微分方程的具有截断误差阶 的差分近似，这时格式（2.37），（2.39）中 和 ，格式（2.40）中 和 分别换成 ，