大型九对角线性方程组的近似三角分解法算法实例“：窘＋-黄山学院.PDF

第ll卷第3期 黄山学院学报 、b1．11．N0．3 2009年6月 JournalofHuangshanUniversity Jun．2009 大型九对角线性方程组的近似三角分解法算法实例 刘文艳，龚学余，彭晓炜 (南华大学教理学院。湖南衡阳421001) 摘要：根据矩阵三角分解法的原理．对大型九对角线性方程组进行近似三角分解，再用追赶法得到线 性方程组的解。给出一个算法实例，表明该算法对大型九对角线性方程组的求解是快速和有效的。 关键词：线性方程组；椭圆型方程；九点差分；近似三角分解法 中图分类号：0175．2文献标识码：A 文章编号：1672-447X(2009)03-0001-04 对椭圆型方程(比如Poisson方程)进行求解时， 通常可以采用差分格式．最常见有五点差分格式和 2对Pisson方程的九点差分 九点差分格式。九点差分格式因其精度较大，因而 也常常被采用。UlPoisson方程在九点差分格式下得方程(1)，对D剖分网格如图l，在茗轴方向上将 到九对角线性方程组，根据方程组的特点，本文提 区间M一1等份，这样有Ⅳ1个节点，节点i=1，…，Ni； 出了近似因子分解法。并给出用近似因子分解法求 在Y方向上将区间／V2-3等份．这样有这样飓-2个 解大型九对角线性方程的一个算法实例，将结果与 节点，节点j=2，…，Ⅳrl。 解线性方程的高斯塞德尔算法进行比较。并进行误 为了保证在边界条件间和声上的对称性， 差分析，表明该算法是可行的。 加入点j=l和歹；Ⅳ2，并令 氓广uu，u讲Fu叫!-2， (2) 1 Poisson方程 差分后整理可得下面形式的方程组： ● ● ． ● ． ● ● Poisson方程 (3) (1) △“：窘+碧：似，)，)’(J，Y)毫D(1) △“2萨+萨2．，(工，)，)’(J)毫 这就是九 J-l J 少1 其中D是善7平面内一有界区域，其边界用彻 点格式，对‰， 来表示．假定它是分段光滑的曲线所组成，为方便 在等距网格t+l 中，它可以通 ． 起见，先取D为矩形区域D=f(茗，y)IOza，O郇 过周围的点表 脚 6}，其边界aD是由4条直线段所组成aD={(石，，，)l 扣1 示出来．如图 z=O，口，Dsys6，y．．-．-O，b，D≤菇s才o l所示。 田1九点差分格式 收稿日期：-2009-12-28 基金项目：田李白然科学基金童助10475039) 作者简介：埘K-艳．(19／9一)．：i工苏建湖人．南年夫擘教理擘院教舜．植科学技术学院博士研兜生．研完方向为收蒙受与等^子体物理．位置分方 租教值解。 万方数据 ·2·