抛物型偏微方程材料.ppt

* * 第七章　解抛物型偏微分方程的差分解法 稳定性讨论是差分法的重要问题。 常遇到的抛物型偏微分方程有： 与上一章的边值问题不同, 这些问题一般都是初值问题. 基本作法：按等间隔剖分，用各离散点上的差商来近似代替该点的偏微商，把偏微分方程转化为线性代数方程组，进而解出各离散点上的待求函数随时间的变化。 主要内容 §7.1简单差分法及其不稳定性 §7.2隐格式差分法与三对角方程 §7.3二维扩散方程的差分格式 §7.4含时间薛定格方程的差分格式 §7.1简单差分法及其不稳定性 一、差分格式 二、例 设 ，初条件和边界条件为： 容易检验，满足微分方程(7.4)和条件(7.8)的理论解为： 现在，我们利用递推式 作数值计算，并比较它与理论解(7.9)的误差。 取空间间隔，时间间隔，时间步数为20(t=0.015)和60(t=0.045)。 当取 时，数值解与理论解的误差很小。 但当取 时，结果就变得很糟。 DO 20 ITER=1,NITER POLD=0. D