矩阵三角分解探讨.docVIP

矩阵三角分解探讨 　　摘 要：矩阵是数学中最重要的基本概念之一，是代数的一个主要研究对象，也是数学研究及应用的一个重要工具。对矩阵三角分解的相关内容进行系统的介绍，给出分解的思想、分解方法、分解的存在性及分解的条件等相关内容，然后就矩阵三角分解思想研究它在解决具体问题中的应用。 　　关键词：矩阵；三角分解；LU分解 　　在近代数学、工程技术、经济理论管理科学中，大量涉及矩阵理论的知识，很多问题都可以归结为矩阵并最终通过矩阵来解决。经查阅发现，目前关于矩阵三角分解的应用研究不少，但对三角分解缺乏系统的研究。 　　矩阵三角分解法是指高斯消去法解线性方程组的变形解法。其实质就是将系数矩阵A分解为两个三角形矩阵L和U相乘，即A=LU。 　　一、矩阵的直接三角分解 　　矩阵的直角三角分解即可以不经过消元步骤，直接将矩阵进行分解。 　　定义1 设A∈Rn×n，若A能分解为一个下三角矩阵L与一个上三角矩阵U的乘积，即A=LU，则称这种分解为矩阵A的三角分解。 　　（1）如果A可分解为A=LDU，其中L是单位下三角矩阵，D是对角矩阵，U是单位上三角矩阵，则称A可作LDU分解； 　　（2）如果在A=LU中，L是单位下三角矩阵，U为上三角矩阵，则称此三角分解为杜利特（Doolittle）分解； 　　（3）如果在A=LU中，L是下三角矩阵，U是单位上三角矩阵，则称此三角矩阵为克劳特（Crout）分解。 　　定理1 n阶方阵A非奇异的充要条件为（或A经行、列变换后）存在LDU分解。其中L为n阶单位下三角矩阵，D为n阶非奇异对角阵，U为n阶单位上三角矩阵。 　　推论1 奇异矩阵不能进行LDU分解。 　　 推论2 若矩阵A有奇异主子矩阵，则A不能直接进行LDU分解。 　　定理2 方阵A的LDU分解唯一的充要条件为A的各主子矩阵非奇异。 　　推论1 设n阶矩阵A的各阶顺序主子式皆不为零，则A存在唯一Crout的分解。 　　推论2 设n阶矩阵A的各顺序主子式皆不为零，则A存在唯一的分解A=LDU 　　其中L为单位下三角矩阵，D为对角矩阵，U为单位上三角矩阵。 　　二、常用的三角分解公式及其应用 　　1.杜利特（Doolittle）分解 　　下面直接用矩阵乘法求L及U的元素，由 　　A=a11 a12…a1na21 a22…a2n■ ■ ■an1 an2…ann=1a21 1■ ■■1an1 an2…lnn-1 1u11 u12…u1nu22…u2n■■unn 　　步骤如下：步骤1 u1i=a1i，i=2，3，…n l1i=a1i／u11，i=2，3，…n 　　计算U第r行，L的第r列元素，r=2，3，…n 　　步骤2 uri=ari-■lrkuki（i=r，r+1，…n） 　　步骤3 lir=（ari-■likukr）／urr（i=r，r+1，…n，r≠n） 　　求解Ly=b，Ux=y计算公式。 　　步骤4 y1=b1y1=b1-■likyk（i=2，3，…n） 　　步骤5 xn=yn=unnxi=（yi-■uikxk）／uii（i=n-1，n-2，…，1） 　　【例】 用直接三角分解法解2 1 11 3 21 2 2x1x2x3=465 　　解：由计算公式，得A=11／2 11／2 3／5 12 1 1 5／2 3／2 3／5=LU. 　　求解Ly=（4，6，5）T，得y=（4，4，3／5）T， 　　求解Ux=（4，4，3／5）T，得x=（1，1，1）T. 　　2.克劳特（Crout）分解 　　直接三角分解法Ax=b（要求A所有顺序主子式都不为零）克劳特分解的计算公式如下： 　　 步骤1 u1i=a1i（i=1，2，…n），lil=ail／u11（i=2，3，…n） 　　计算U第r行，L的第r列元素，i=2，3，…n 　　步骤2 lli=aij-■likukj，i=1，2，…n，j=1，2，…i， 　　步骤3 uij=aij-■likuiii=1，2…n-1，j=i+1，…，n 　　【例】试将下列矩阵进行克劳特分解：A=4 8 42 7 21 2 3 　　解 首先我们约定，当上限小于下限时，和号“■”取0，根据上式有l11=a11=4， 　　 u12=a12／l11=2，u13=a13／l11=1， 　　 l21=a21=2， 　　 l22=a22-l21u12=3，u23=（a23-l21u13）／l22=0， 　　 l31=a31=1， 　　 l32=a32-l31u12=0， 　　 l33=a33-l31u13-l32u23=2， 　　则有A=LU=4 2 32 3 01 0 21 2 1 1 0 1 　　3.三对角方程组的追赶法 　　在许多科学计算问