求三对角和周期三对角矩阵逆矩阵的一种新算法.pdfVIP

长江大学学报 (自然科学版) 2010年3月第7卷第 1 ：理工 · 126 · JournalofYangtzeUniversity(NatSciEdit) Mar．2010．Vo1．7No．1：Sci Eng 求三对角和周期三对角矩阵逆矩阵的一种新算法 余承依，陈跃辉，赵立群 (1NJI币范学院数学与信息科学系，福建漳州363ooo) [摘要]利用周期三对角矩阵的结构特点，通过适当的矩阵分块，将周期三对角矩阵的求逆转化为三对角 矩阵的求逆问题，同时借助矩阵的Crout分解方法给出了一种求三对角矩悻逆矩阵的的简单算法，并将其 应用到求解周期三对角矩阵逆矩阵中。数值试验表明此算法是有效的。 [关键词]三对角矩阵；周期三对角矩阵；逆矩阵 [中图分类号]O241．6 [文献标识码]A [文章编号]1673—1409(2010)01一N126一O3 在许多科学与工程计算中经常会出现三对角矩阵和周期三对角矩阵，如构造具有周期型边界插值条 件的三次样条函数、热传导方程等都会涉及到此类矩阵。不少学者对三对角矩阵 引和周期三对角矩 阵[。 进行了研究，目前已有一些研究成果。笔者在此基础上进一步研究此类矩阵逆矩阵的求解，讨论 了周期三对角矩阵逆矩阵与三对角矩阵逆矩阵的关系，将周期三对角矩阵求逆转化为三对角矩阵的求 逆，并进一步给出了一种三对角矩阵求逆的快速算法。 1 周期三对角矩阵求逆 设周期三对角矩阵： C1 n I b2 C2 T == 口 b C n，1 6”+1 这里假定r非奇异，对 T做分块 ： b1 C a2 b2 f2 T 一 [On{：] — ∈ n，r1 6一 C 一1 n b I口1 l0 一 ∈ R 一 1 10 对 矩阵的逆矩阵做类似的分块 ，记 ： r === [ 定理1设T为周期三对角矩阵，对T做分块L-『 ]，记r的分块矩阵为f_ ]，则有： a )，-I Lm； d _J d ===1／(7~一aT1 ) m 一一 TT ===一 C 一 1 +p ／ 证明 由： [收稿日期]2009—12—06 [基金项目]福建省教育厅资助科技项 目 (JB07147)；漳州师范学院科研项目 (SK07008)。 [作者简介]余承依 (1976一)，男，1998年大学毕业，硕士，讲师，现主要从事矩阵计算和不确定性信息处理方面的教学与研究工作。 第 7卷 第 1期：理工 余承依等：求三对角和周期三对角矩阵逆矩阵的一种新算法 1厂G 旷 一 y 儿mT耋]I~一 一[ C + — L 户 + 一 (洲 口T + ===1 m 1 + 柑 — O1× (1) 由式(1)第2式得P 一一d 。将P代人式(1)第3式得d 一 1／(一a 卢)。再由式 (1)第4 式易得 一一 。将 mT代人式 (1)第 1式可得： G 一 (L一 )一 +PmT／d