抛物方程差分方法初步.ppt

若希望得到M层 的值，由于计算格式要用到M-1层 处的值，故初始时刻应给出 的初始值。 2。第三边值—初值问题 (1)边界条件处理 法一 局部截断误差O(h) 法二 为提高精度, 采用中心差商 同理 局部截断误差 五、二维热传导方程 以下列初边值问题为例 选取网格 记 1。显式格式 稳定性条件 对n维热传导问题，稳定性条件为 维数越高，时间步长越小，计算量越大！ 2。隐式格式 每一层上要求解一个(N-1)(M-1)阶的 五对角方程组。 交替方向法 将二维问题化为两个一维问题来求解。 引进新的差分格式，满足 （1）无条件稳定； （2）有合理的精确度； （3）相应的代数方程组易解。 3。P-R(Peaceman-Rachford)方法 M=N 时，每层上要求解N-1个含N-1个未知量的三对角方程组。 计算量比隐格式低得多！无条件稳定！ 缺点：不能推广到三维问题！ 4。D-R(Douglas-Rachford)方法 每层上求解三对角方程组。引入过渡层，需要较大的存储量！ 可推广到三维！无条件稳定！ 小结： 1。二层方法 显格式：条件稳定 隐格式：无条件稳定 2。初值问题和第三边值条件的处理 3。二维问题 显格式：条件稳定