求解大规模稀疏线性方程组的算法.PDF

第27卷第3卿 北耵 [虫，大学学报 v【、117、h)j 二00】年 9，j J叭IRNALo}BFuIN(；r)oL、rEcHhIfLIN】v}瞄lnl 船”100l 求解大规模稀疏线性方程组的算法 陈 志，高旅端 (北京【1L大学应用数理学院，北京】000二2) 摘 要缔出了求解大规模稀疏线性方程组的一个窭用算法，遮算法能够保持娟阵的稀疏|生和减少存∞：量，并 且能够求解一些大规模的问题．而这些问题所对应的系数矩阵可能不是稀疏的． 关键词：大规模，稀疏；Lu分解 中围分类号o22l2 文献标识码：A 04 文章编号2001)03一0262 0引 言 一般说来，太规模问题的系数矩阵具有某些稀疏的结拘，称为大规模稀疏问题．对大规模稀疏问题， 如果能够通过调整方程及未知量的顺序使得方程组的系数矩阵威带状结构，则对系数矩阵使用通常的Lu 分解，可以保障单位下三角矩阵￡及上三角矩阵U仍为带状结构．尽管如此，若在分解中使用了选主元的 技术，将会使分解后的矩阵的带宽增大．例如，设方程组系数矩阵的上带宽为，，下带宽为，，在分解中使用 选取列主元的技术，则上三角矩阵D的带宽可能是f+r．这样，算法的计算量和存贮量将会增大．因此， 相对系数矩阵的阶数H，当上、下带宽，、r都很小时、可以使用通常的Lu分解方法求解方程组；对于下 H：shcnhe佗阵，，=1，r；”一1．则以其为系数矩阵的方程组不宜使用通常的Lu分解方法求解． 艾献【1]和【2】中提出了一种新的Lu分解方法．本文将其改造为一种求解线性代数方程组的实用方 法二对系数矩阵的稀疏性的要求更为宽松，仅要求上、下带宽，、r中有一个比一小得多即可，例如下 Hessenber叠阵．新方法的存贮量比通常的Lu分解方法要小． 1基本算法 考虑线性代数方程组 口jr=6(j=l，2，…，月) (I) 其中工～．∈只”“=l，2，…，Ⅳ)设√=(口．，口：，…，n。)，则(|)可以写为 爿1x=6 (2) 这里：占=伯。，6，，…，6。)7∈尺I 不失一般性，设一1的上带宽7较小国为对较小的下带宽，总可以通过对矩阵行列的初等变换转化成 较小上带宽的形式． 新力’法是一种迭代方法．从初始点』．出发，第≈次迭代求出J。+。，满足(1)的前A个方程 口．。J=6． (j=【，2．…，t) (3) 算法l ①取初始电J．∈月”；置u．=，。，皮=I． ②计算口．=一Ⅱ．∈∥小。 收稿日期．二0f肌D9一{5． 基金项H：北京ff】敷委基金资助项甘(98080lO!)． 作者简介：陈志(1944一)，异，教授，硕L；高旅端(1944一)．男，教授．硕上 万方数据 第3期 陈志等：求解犬规模稀疏线性方程组的算法 @计算^。=q毛，≮的选取使得口一。=a!≈≠o @计算J。+，=x。+^。^。，其中A。=(6。—n》。)／a■。 ⑤若＆=，I，则停止；否则，转@． ⑥校正q：峨+．=％吼，这里蛾∈R‘4““…一“，并且满足咖。=o．置女=k+1，转2． 定理l 由算法1产生的J。+．陋=l，2，…，n)满足式(3)，并且 口 口 = f庶=I，2，-一，”一I) (4j 迕 对 峨^ —一=nj￡‘‘=njq．吩一。％=njq缉目+。…哆。≈=a阻置+，…弓一，：，=o 从而对汪1，2．一，t，