CH02抛物方程差分法CH2.1-2.4,2.5-2.8.pptVIP

在上式中保留前二项，并且 和 分别用 和 代替，则得差分方程 (2.39) 也可通过直接用中心差分算子 代替微分算子 的办法获得方程(2.38)的差分近似 (2.40) 这也是一个显式差分格式。 格式(2.39)和(2.40)的截断误差阶都是 。易见，由 注： 均在 处计算。 Delta 代入格式(2.40)即为格式(2.39)，差分格式(2.40)的推导方法,即在微分方程中直接用差分算子代替 正如前面已经指出的是推导差分格式的一个常用方法。 显然，微分方程(2.36)，(2.38)中的 如果为 ，即其自变量包括空间变量和时间变量，这时差分格式(2.37)，(2.39)，(2.40)同样是微分方程的具有截断误差阶 的差分近似，这时格式(2.37)，(2.39)中 和 ,格式(2.40)中 和 分别换成 , 。 2.3 隐式差分格式 隐式差分格式特点： 1. 具有二个或二个以上结点处的值未知； 2. 计算工作量较大； 3. 稳定性较好。 得 由 推导其最简单的隐式差分逼近─古典隐式格式。 现在对热传导方程 2.3.1 古典隐式格式 17 15 格式用图2.5表示，其截断误差阶为 ，与古典显式差分格式相同。 或者 (2.41) 保留二阶导数项，且以 替代 ，则得差分格式 我们也可通过直接用差分算子代替 的方法，即 代入微分方程，得到格式(2.41)。 古典隐式差分格式 m,n+1 m-1,n+1 m+1,n+1 m,n 图2.5 隐式差分格式是解热传导方程(2.26)的常用的差分格式，由式(2.24)，有 2.3.2 隐式格式 由 得 (2.42) 42 两边仅保留二项，用 代替 ，则得差分格式 (2.43) 这是一个隐式差分格式，称为 差分格式，截断误差阶为 。 (2.44) 由于格式(2.44)中包括六个结点，故也称为六点格式(如图2.6所示)。 m,n+1 m-1,n+1 m+1,n+1 m,n 图2.6 m-1,n m+1,n 44 也可将 代入微分方程(2.26)，得到 格式。 由式(2.19.3)，可令 则可得 另一精度较高的六点差分格式，如前在式(2.42)中仅保留直到 的项，即有 13 代入上式，则有如下差分格式： (2.45) 称为 差分格式。 38 截断误差阶 23 23 因为 48 前面，我们已经推导了热传导方程(2.26)的古典显式格式。古典隐式格式及 格式等。实际上，它们都可以作为本节推导的加权六点隐式格式的特殊情形。 2.3.3 加权六点隐式格式 由 得到 即 用 代替 ，则得差分格式 或者 (2.46) 这是一个六点差分格式(如图2.7所示)，称为加权六点差分格式。 40 时, 为古典显式格式； 时, 为古典隐式格式； 时, 为 格式； 加权六点格式亦可直接由差商代替导数得到 2.3.4 系数依赖于 的一维热传导方程的一个隐式格式的推导 由其 展开式可得 (2.47) 的差分逼近。 考虑方程 已知 11 令 代入式(2.48)，则 (2.48) 因此 43 格式(2.49.1)具有截断误差阶 。 这是一个隐式差分格式(如图2.8所示)。 (2.49.1) 因此得差分方程 (2.49.2) 可写成形式 m,n+1 m-1,n+1 m+1,n+1 m,n 图2.8 m-1,n m+1,n 前节引进的隐式差分方程，在要求解未知函数值的时间层 上包括三个未知函数值 。因此，这些隐式差分格式仅仅适合于解如图 中所示的边值问题。在每一时间层，需要求解的隐式差分方程形成了一个线性代数方程组，它的系数矩阵是三对角形矩阵，即仅在主对角线及其相邻二条对角线上有非零元素。方程组写成一般