矩阵的三角分解论文.docVIP

矩阵的三角分解 阿卜杜克力木.阿卜杜瓦尔斯 (伊犁师范学院数学与统计学院06（1）班 新疆 伊宁 835000) 摘要：把矩阵分解为一些矩阵的乘积或一些矩阵的和，无论是在矩阵的理论研究，还是在矩阵的应用中，都有重要的意义。在本文中我们实际上已经介绍了一些矩阵的分解，即对角化，若当标准型等，在本文中对一些常用的矩阵分解作进一步的介绍矩阵的三角分解。 关键词：矩阵分解，三角分解，线性方程组 中图分类号：0171 一、矩阵的三角分解 设A为阶方阵，线性方程组，如果能够分解为A=BC，且为三角矩阵，则线性方程组的求解问题转化为求方程组 的解，这为解线性方程组提供了一个简单的方法，也为解矩阵的其他问题给出了一个途径。 定理 1 设，则A可以惟一地分解为（）的充分必要条件是A的顺序主子式都不等于零。其中L为下三角阵,为单位上三角阵，（为单位上三角阵，U为下三角阵）。 证明: 只证即可,的证明与词此类似，读者可自证）。 必要性: 设A可以惟一地分解为，则 其中故，A顺序主子式=。因为A可逆及L与为可逆三角阵，故0所以，。 充分性 用归纳法，当时命题显然成立。假设命题在时成立，设，其中是下三角阵，为单位上三角阵， ，则 所以 下面证明惟一性：设，则因为，为单位上三角阵，下三角阵，所以是单位上三角阵，为下三角阵，从而，= 。以上的分解叫做矩阵的三角分解。 例1．求三阶方阵的LU分解。 解： 因为 所以A有惟一的LU分解。且 由可计算及如下： ， 于是 因此A的分解为 定理2 设，A的顺序主子式都不等于零，则A可以惟一地分解为,其中为单位下三角阵，D为对角阵，为单位上三角阵。 证明： 由定理1得A可以唯一地分解为A=L,设 L=则L==D 即A可以唯一地分解为A=D. 矩阵A=的Crout分解具体做法是：设下三角阵L=，单位三角阵，其中，i=1,2,….,n.令A=L,比较俩端的对应元素即可得到 , , 以上的分解叫做Crout分解。 例2 求A=的Crout分解和分解。 解 : 设A=L,即 = 比较两端的对应元素即得 , , 所以 A===L A== 定理3 （cholesky分解）正定的爱儿米特矩阵A可惟一地分解为A=L,其中L为正线下三角阵。 证明 设A,因为A正定，故A的顺序主子式都不等于零，由定理2，A可以惟一地分解为A=,因为A=,故,又由分解的惟一性得,所以.设，因为A正定且A与D合同，故记A=(),则D=,所以A可以惟一地分解为A===L,因为为单位下三角阵，A为正线对角阵，故L为正线下三角阵。以上的分解叫做矩阵的cholesky分解。 总结 矩阵是对数学学者从事数学研究工作的最灵活的一种工具。由于矩阵论既是一门发展完善，理论严谨，方法独特数学基础课，又广泛应用于工程科学的各个领域，尤其是在自然科学，工程技术，经济和管理科学，航空研究，国家边防，信息交流，新题研究等各个方面的研究不必缺少的科学工具之一。不过矩阵论课程的理论性强（其中矩阵分解令人最糊涂的知识之一），概念比较抽象，而且独特的数学思维方式和解题技巧。比如；随着科学技术的迅猛发展，特别是计算机的广泛应用，矩阵理论越来越受到数学工作者，科技和工程人员的重视。它不仅是一门重要的数学分支，而且在数值分析，最优化方法，微分方程，控制理论，数学模型等分支及各种工程学科有极其重要的应用。总而言之，这方面的研究还得继续开展。该文章介绍几种矩阵的分解方法，这几种方法能把有关矩阵的复杂问题简单化。 参考资料 [1] 魏丰，史荣昌.矩阵分析[M].北京理工大学出版社， 2005年6月,83-89 [2]罗家洪.矩阵分析引论[M].广州：华南理工大学出版社，1992年,129-132 [3]李乔.矩阵论八讲[M].上海：上海科学技术出版社， 1985年10月,1-5 [4]张凯院，徐仲陆全.矩阵论典型题解析及自测试题 [M].西安：西北工业大学出版社，2000年10月,149-162 [5]刘丁酉.矩阵分析[M].武汉大学出版社， 2003年8月,150-162 [6]苏育才，张跃辉.矩阵理论 [M].北京科学出版社，2003年8月,159-163 [7] 史荣昌，魏丰.矩阵分析学习指导[M].北京理工大学