kmr隐式差分方程课稿fmj.pptVIP

2.3 隐式差分格式 与显式差分格式不同，隐式差分格式中包括了（n+1）时间层上二个或二个以上结点处的未知值（例如 ），使用隐式差分格式和使用显式差分格式求解完全不同。相对而言，使用隐式差分格式求解，每时间层包含有较多的计算工作量。从后面对差分格式的稳定性分析可知，隐式格式的优点在于，其稳定性要求对步长比的限制大为放宽，而这正是我们所期望的。 2.3.1 古典隐式格式 现在对热传导方程 推导其最简单的隐式差分逼近——古典隐式格式。由 故 式中左边如果仅保留二阶导数项，且以 替代 ，则得差分格式 或者 （2.41） 格式用图2.5表示，其截断误差阶为 ,与古典差分格式相同。 图2.5： 为了求得第（n+1）时间层上的 的值，必须通过解线性代数方程组。这是一个隐式差分格式，必须联合其初边值条件求解。格式（2.41）通常称为古典隐式格式。 我们也可以通过直接用差分算子代替 的方法，即 代入微分方程，得到格式（2