三对角矩阵的逆的算法及MATLAB实现—学士学位毕业论文.doc

2014届学士学位毕业论文三对角矩阵的逆的算法及MATLAB实现学 号：姓 名：班 级：指导教师：专 业：数学与应用数学系 别：数学系完成时间： 年 月学生诚信承诺书本人郑重声明：所呈交的论文《 》是我个人在导师 指导下进行的研究工作及取得的研究成果尽我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写的研究成果，也不包含为获得长治学院数学系或其他教育机构的学位或证书所使用过的材料所有合作者对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意签名： 日期： 论文使用授权说明本人完全了解长治学院数学系有关保留、使用学位论文的规定，即：学校有权保留送交论文的复印件，允许论文被查阅和借阅；学校可以公布论文的全部或部分内容，可以采用影印、缩印或其他复制手段保存论文签名： 日期： 指导教师声明书本人声明：该学位论文是本人指导学生完成的研究成果，已经审阅过论文的全部内容，并能够保证题目、关键词、摘要部分中英文内容的一致性和准确性。

学位论文指导教师签名： 时间 摘要三对角矩阵在现实生活中有很多的应用，因此三对角矩阵的计算近年来被广泛地研究分块周期三对角矩阵在科学和工程计算方面应用广泛，块三对角矩阵和分块带状矩阵在数学、物理和工程上的很多问题中都有重要的应用本文基于三对角矩阵的结构特点，给出了利用解线性方程组的方法、LU分解的方法求三对角矩阵逆矩阵的新算法，这些新算法运算量小，节省内存，在整个计算过程中，只需要进行较少次的乘除运算，新算法比传统算法的计算复杂度和计算时间要低 其次，通过算例来表示该算法的有效性和可行性 最后，利用MATLAB编程来实现三对角矩阵逆矩阵的新算法关键词：分块周期三对角矩阵；块三对角矩阵；分块带状三对角矩阵；解线性方程组；LU分解法；逆矩阵；MATLABTriple diagonal matrix inverse algorithm and MATLABAbstractTriple diagonal matrix in real life there are many applications, so the triple diagonal matrix calculation was widely studied in recent years. Block periodic triple diagonal matrix is applied widely in science and engineering calculation, and the block triple diagonal matrix block banded matrices in mathematics, physics and engineering has important applications in many of the problems, in this paper, based on the structure characteristics of triple diagonal matrices, is given by using the method of solving linear equations, the recursive method, LU decomposition of the new method to calculate the inverse matrix of triple diagonal matrix algorithm, the new algorithm computational complexity is small, save memory, in the whole computing process, only needs less arithmetic, a new algorithm than the traditional algorithm of computing complexity and computing time.Second by an example to show the feasibility and effectiveness of the algorithmFinally, using MATLAB to realize the triple diagonal matrix inverse matrix of the new algorithmKey words: Block periodic triple diagonal matrix; Block-triple diagonal matrix; Block banded triple diagonal matrix; Solution of linear equations; LU decomposition method; inverse matrix; MATLAB.目录1.引言 52.基础知识 62.1 定义1 62.2 定义2 62.3 定义3 73.分块周期三对角矩阵逆的新算法 73.1 分块三对角矩阵的一些性质 73.2 求分块周期三对角矩阵逆矩阵的新算法 104.块三对角矩阵的逆的算法 114.1 块三对角矩阵的一些性质 114.2 块三对角矩阵的逆 134.2.1 块三对角矩阵逆的性质 135.三对角矩阵逆元素的表示 145.1 一般三对角矩阵 145.2 用解线性方程组的方法求三对角矩阵的逆的算法 165.2.1 基本原理与算法 165.2.2三对角矩阵的逆矩阵的算法 186.三对角矩阵逆的算法的MATLAB实现 187.结束语 188.参考文献 18附录 19致谢 191.引言1.1 课题来源及选题意义三对角矩阵是计算数学的重要组成部分。

它是研究代数问题的三对角矩阵快速算法及有关理论的一门学科，它既涉及数学理论方面的研究，又涉及工程设计方面的研究随着科学技术的发展和计算机的普及，矩阵理论和方法得到了越来越广泛的应用在近代数学、工程技术、经济理论及管理科学中，大量地涉及到矩阵的理论，特别是一些具有特殊结构的三对角矩阵，相应的计算规模也越来越大近十几年来，国防科技和国民经济建设的许多领域中就不断地提出了大型或超大型科学计算问题由于矩阵在各个学术领域和重要应用课题中所起的不可替代的作用，故有必要对其进行细致的研究科学技术和工程应用中需要进行大量地矩阵计算，而这些矩阵自身往往具备一些特殊的结构，这既是本文所研究的一类重要而特殊的稀疏矩阵——三对角矩阵的求逆问题，该类矩阵经常出现在信号处理、图像处理和数值分析等学科的一些应用问题中在该类矩阵的有关研究中，求逆是一个重要的问题，且一直是人们的研究热点，目前已有一些研究三对角矩阵求逆的成果由于在许多科学技术与工程应用中，经常会出现大量的三对角矩阵的逆的算法进行计算，所以我们有必要对三对角矩阵的逆的算法进行研究1.2 研究现状对于三对角矩阵逆的算法及MATLAB实现，目前很多学者根据一些三对角矩阵的特殊结构，用不同的方法对三对角矩阵逆的算法及MATLAB实现做了很多研究，并取得一定的成就。

例如2012年杜永恩，陆全，徐仲利用LU和UL分解，并使用Sheman-Morrison-Woodbury 公式，得到一个求分块周期三对角矩阵逆矩阵的新算法(见[2])；冉瑞生和黄廷祝利用LU和UL分解给出了两个绞形块分解，建立了一个块三对角矩阵求逆的算法（见[3]）；刘长河,刘世祥,汪元伦用解线性方程组方法得到求逆的算法（见[5]）；余承依,陈跃辉,赵立群利用周期三对角矩阵的结构特点,借助矩阵的Crout分解的方法给出了一种求三对角矩阵逆矩阵的的算法（见[6]）；车毅,徐仲,雷小娜利用递归方法给出了求分块周期三对角矩阵的逆矩阵的一种新算法（见[7]）；冉瑞生,黄廷祝,刘兴平等研究了具有Doolittle分解的三对角矩阵的求逆，得到一个求逆的算法(见[8])不少学者研究了三对角矩阵的逆，并进一步给出了求三对角矩阵逆矩阵的新算法，而且新算法的计算量要比传统算法小，计算效率有显著提高，但其算法的实现还有待探究文中，为了讨论的方便，记三对角矩阵为（1.1）且定义n个数：=，=-（=2,…,n）为方便起见，我们约定若>, =1本文研究以求解分块周期三对角矩阵逆矩阵的新算法、块三对角矩阵逆矩阵的新算法、分块带状三对角矩阵求逆的算法、三对角矩阵逆元素的表示、稀疏矩阵的逆的算法，用算例来表示该算法的有效性和可行性。

最后用MATLAB编程来实现三对角矩阵求逆矩阵的算法2.基础知识2.1 定义1 阶矩阵 称为三对角矩阵. 如果，当. 2.2 定义2 设分块周期三对角矩阵有如下形式：（2.1）其中，的元素，，，,都是阶方阵若，则矩阵为分块三对角矩阵；若，矩阵中的元素，，，，都是实数，则矩阵为周期三对角矩阵，且若，且，则矩阵为对称周期三对角矩阵2.3 定义3 设块三对角矩阵具有如下形式（2.2）所有的块均是阶矩阵且非奇，负号仅是为了符号处理上的方便而添加的设是的顺序主子矩阵，其中的所有对角块矩阵即是的对角块矩阵，，…，假定的所有顺序主子块矩阵，，…，均非奇异为了讨论的方便，设，其中是矩阵3.分块周期三对角矩阵逆的新算法3.1 分块三对角矩阵的一些性质引理1 设是分块三对角矩阵，则可分解为：其中，，，可按：，，，，，计算证明(1)因为所以 (2)因为= 所以 .引理2(Sherman-Morrison公式) 设是阶可逆方阵，均是维列向量，则当且仅当时，是可逆的，且：引理3(Sherman-Morrison-Woodbury公式) 设是阶可逆方阵，,均为矩阵，则当且仅当可逆时，是可逆的，且：证明：令 则（3.1）令 则（3.1）式为 （3.2）在（3.2）式左右两端同时乘以 令，则（3.2）为 可得， 又所以 所以引理4 设是分块三对角矩阵，且是顺序主子阵可逆。

设 存在，则存在4个矩阵： ，，， 均为阶方阵，使得：，或其中对所有的都有且，，，可如下求得：给定，：（1），其中，，，有引理1得到3.2 求分块周期三对角矩阵逆矩阵的新算法给定阶可逆方阵, ,令, ,构造向量：, （3.3）则分块周期三对角矩阵可表示为：其中 （3.4） 由引理3的Sherman-Morrison-Woodbury公式可得：由此可得如下结论若是如（2.1）所示的分块周期三对角矩阵，，，如（3.3）和（3.4）中定义，设可逆，则可逆的充要条件是可逆，且的元素可由如下算法计算得到：任意选定，可逆，，给定，，， ，， ，， ， ，，, .以上所得即为分块周期三对角矩阵的逆矩阵，当取时，分块周期三对角矩阵子块都是1阶的实。