差分格式.doc

§3. 热传导方程 上一节曾指出，由于定解问题中的每一个偏导数都有多种差分近似，所以一个定解问题可以有多个不同的差分格式。下面以热传导方程为例，对此展开讨论。为简单起见，先不给出定解条件。 考虑热传导方程 方程中出现了一阶时间导数 和二间空间导数 。 对于一阶时间导数 ，常用的差分近似就有三种，记 向前差分近似 向后差分近似 中心差分近似 这里，我们用“”代表上一节推导差分近似的过程。 前面已经看到，二阶导数通常用中心差分近似。对这里的二阶空间导数的中心差分近似为 但是也可以考虑其他可能的方案。由泰勒展开，有 所以，如果用 代替 ，虽然会引入新的误差，但这种误差与差分近似已有的误差为同一量级的，因而还是可以接受的。这样一来，二阶空间导数的差分近似又有了两种新的方案 将前面讨论过的一阶时间导数的三种差分近似与这里给出的二阶空间导数的三种差分近似方案 一阶时间导数的三种差分近似 向前差分近似 向后差分近似 中心差分近似 二阶空间导数的三种差分近似方案 方案一 模板 ( ( ( 方案二 模板 ( ( ( 方案三 模板 ( ( ( （注）二阶空间导数差分近似的三个方案 的取值固定在 、 或 三个不同的时刻。 在上面的表中，除了给出差分格式，还给出了这个格式所用到的网格点的图示，称为格式的模板。另外，表中打 表示这个格式在网格点的使用上明显不合理，下面将不予考虑。在余下的格式中，我们挑出三个作为基本格式（打(者，后面将会解释这三个名称的含义），它们是 两层显式格式 两层隐式格式 三层显式格式 对表中标有 (、(、(、( 的这四个格式，有以下的分析。 对于 ( 和 ( ，令 ，则 ，于是它们可写成 显然它们与前两个基本格式是类似的。 而对于格式 ( 和 ( ，仍令 并取 ，则由于 ，所以如果用 作为时间步长， 就是 的下一个时刻，可记作 ，于是这两个格式也可写成 可见它们也与前两个基本格式类似。 综合以上的分析，这四个格式与基本格式重复，也不必再考虑。 以下将集中分析上面选出的三个基本格式。前两个基本格式，只用到 和 两个时间层的网格点，故称为两层格式。而第三个格式用到了 、 、 三个时间层的网格点，所以称为三层格式。 下面为热传导方程补上定解条件，考虑热传导方程定解问题 对这个定解问题，第一个基本格式为 将它改写成便于计算的形式 式中 。这个格式的求解步骤如下： 时刻（） ： 时刻（） ： ，() ， 时刻（） ： ，() ， 一般地，如果已经求出 时刻的近似解，则 时刻： ，() ， 对 反复进行这一步骤，直到求出 时刻的近似解。 从这个计算过程可以看出，只要用初始条件给出第0层网格点上的解，就可以用这个格式一步一步地计算出各个时间层的近似解。这样的格式称为显式格式。 再来考虑第二个基本格式，它可以改写成 显然，第 层的三个近似解 、、 是第 层的 方程组中的 ， 都是已知的（由边界条件事先确定的），所以可以移到方程右边。 如果将第 层的近似解 ， 再将方程组的系数矩阵记作 则方程组可写成 这个格式的求解步骤如下： 时刻（） ：由初始条件 计算 。 时刻（） ：利用 和边界条件计算右端项 ， 解方程组 ，求出 。 时刻（） ：利用 和边界条件计算右端项 ， 解方程组 ，求出 。 一般地，如果已经求出 时刻的近似解 ，则 时刻 ：利用 和边界条件计算右端项 ， 解方程组 ，求出 。 对 反复求解方程组，直到求出 时刻的近似解。 从上述计算过程可以看出，这个格式一个时间步的计算都需要求解一个方程组，才能得到这个时间层的近似解。这样的格式称为隐式格式。 上面对显式格式和隐式格式的分析表明，从计算量和编程工作量的角度来说，显式格式由于不需要求解方程组，计算量小，便于编程，所以优于隐式格式。 最后，简单分析一下第三个基本格式，即三层格式，它可以改写成 可以看出，第 层的近似解第 层第 层第层第层第层 ，它只能给出第层第层 ，下面给出的格式也是如此。由此可见， 是一个重要的量，称为（热传导方程差分格式的）网格比。 除了本节已经讨论过的三个基本格式之外，热传导方程还有没有其它形式的差分格式？答案显然是肯定的。这里再给出一个比较常用的格式，作为这一节的结束。 在三个基本格式中有两个是两层格式，它们的区别在于：用中心差分近似二阶空间导数时，时间自变量 的取值固定