偏微分方程的交替分组显式迭代方法.pdfVIP

2008年第 2 期 新疆财经大学学报 No. 2. 2008 (总第 32 期 ) Journal of X inj iang Un iversity of F inance and Econom ic s General No. 32 偏微分方程的交替分组显式迭代方法 1, 2 2 1 张春梅 ,程 　蓓 ,李 　静 ( 1. 新疆大学 ,新疆 乌鲁木齐 830046; 2. 西安交通大学 ,陕西 西安 7 10049) 摘 　要 :求解复杂的偏微分方程或方程组时 ,对方程构造的差分格式可分为显式和隐式两大 类 。显式格式虽适用于并行计算 ,但其稳定性条件有严格限制 ; 隐式格式稳定性虽好 ,但不能 直接用于并行计算 。本文利用交替分组显式算法的思想来设计隐式差分方程组的迭代解法 , ( ) 得到交替分组显式迭代法 A GE I 。这种方法可用于主对角 占优的一般三对角方程组的迭代 求解 ,不仅格式容易实现而且可直接进行并行计算 。 关键词 :交替分组显式迭代方法 ;收敛性 ;绝对稳定 ;并行计算 中图分类号 : O24 1. 82　　　文献标识码 : A 　　　文章编号 : 167 1 - 9840 (2008) 02 - 0096 - 03 一 、引言 随着并行机和向量机的问世与发展 ,传统的有限差分方法逐渐暴露出一些缺陷 :古典显式虽适用于 并行计算 ,但其稳定性条件有严格限制 ;而古典隐式和 Crank - N ico lson 格式稳定性虽好 ,但不便于直接 在并行机上应用 。20世纪 70年代初 ,M iranker指出 ,用有限差分逼近偏微分方程时 ,主要是组织传统差 分方法的并行实现 。自 20世纪 80 年代以来 , Evan s和 A bdu llah 设计的分组显式方法 [ 1 ] 不但保证了数 值计算的稳定性 ,并且由于可以显式求解而使该方法具有很好的并行性质 。1957 年 , Sau l ’yev提出了 一个半隐式格式 ,它具有一般隐式格式稳定性好的特点 ,可以从右向左 、从左向右显式地交替求解 ,这样 就不需要求解线性代数方程组 ,因而受到重视 。而并行计算机和向量计算机的出现 ,使 Sau l ’yev格式 ( ) [ 2 ] 不易使用 。1983年 Evan s等人在 Sau l’yev格式的基础上 ,构造了分组显式 GE方法或 GE格式 。 　　二 、交替分组显式迭代方法 ( ) [ 3 ] 交替分组显式 A GE 计算的思想可以用来设计迭代求解隐式差分方程组 ,这就是 A GE迭代法 。 A GE迭代法可用于具有主对角优势的一般三对角方程组的迭代求解 。假设差分格式的矩阵形式为 : k +1 a b u1 f 1 k +2 c a b