毕业论文矩阵分解方法的探讨.docx

毕业论文矩阵分解方法的探讨矩阵分解方法的探讨 The discussion about decomposition of Matrix 专 业: 数学与应用数学 作 者: 指导老师: 学校 二一 摘 要 矩阵是数学研究中一类重要的工具之一, 有着非常广泛的应用, 矩阵分解对矩阵理论及近代计算数学的发展起了关键作用. 本文从矩阵的LU分解、 矩阵的QR分解、 矩阵的满秩分解等几个方面对矩阵分解方法进行了论述: 给出了矩阵分解的几种方法. 关键词: 矩阵, 对称正定矩阵,矩阵的三角分解；矩阵的满秩分解；矩阵的QR分解. I Abstract The matrix is a important tool in class of mathematical research, and it has a very wide range of applications, matrix decomposition plays a key role in matrix theory and development of modern computational mathematics. This article begin at the discuss from the matrix of LU decomposition、Matrix of the QR Decomposition、Matrix decomposition of full rank and so on. given a matrix factorization method. Keywords: Matrix; Symmetric positive definite matrix, Triangular decomposition of matrix; matrix full rank decomposition; QRdecomposition of matrix. II 目 录 摘 要 . I Abstract . 错误！未定义书签。 0 引言 . 1 1 矩阵的三角(LU)分解 . 1 1.1 矩阵的三角分解基本概念与定理 . 1 1.2 常用的三角分解公式 . 7 1.2.1 杜利特分解 . 7 1.2.2 克劳特分解 . 7 1.2.3 乔累斯基分解 . 8 2 矩阵的满秩分解 . 15 2.1 矩阵的满秩分解基本概念与定理 . 15 3.1 矩阵的QR分解基本概念与定理 . 18 3.2 矩阵QR分解的常用方法 . 20 3.2.1利用Householder矩阵变换. 20 3.2.2利用QR分解公式 . 20 3.2.3利用列初等变换法 . 21 参考文献 . 243 矩阵的QR分解 . 18 0 引言 矩阵的三角分解、正交三角分解、满秩分解将矩阵分解为形式比较简单或性质比较熟悉的一些矩阵的乘积,这些分解式能够明显地反映出原矩阵的许多数值特征,如矩阵的秩、行列式、特征值及奇异值等. 另一方面, 构造分解式的方法和过程也能够为某些数值计算方法的建立提供了理论依据. 本文从矩阵的LU分解； 矩阵的QR分解； 矩阵的满秩分解等几个方面对矩阵分解方法进行论述: 探讨矩阵分解的方法. 1 矩阵的三角分解 1.1 矩阵的三角分解基本概念与定理 定义1.1 设AÎCm´n,如果存在下三角矩阵LÎCm´n和上三角矩阵UÎCn´m, 使5得A=LU, 则称A可作三角分解或LU分解. 定义1.2 设A为对称正定矩阵, D为行列式不为零的任意对角矩阵,则A=AT, U为一个单位上三角矩阵, 且有A=LDU成立： 1) 如果L是单位下三角矩阵, D是对角矩阵, U是单位上三角矩阵, 则称分解A=LDU为LDU分解. %为2) 如果L=LD是下三角矩阵, 而U是单位上三角矩阵, 则称三角分解A=LU克劳特(Crout)分解; %=DU是单位下三角矩阵, U%为%为上三角矩阵, 则称三角分解A=LU3) 如果U杜利特(Doolittle)分解; %-1U%, 称为不带平方根的乔累斯基(Cholesky)4) 如果A=LDU=LDD-1DU=LD分解; %, 则A=LDU=LDDU=LU%, 由于U%=L%T, 则L, DU=U5) 如果LD=%T, 称为带平方根的乔累斯基(Cholesky)分解. A=LL12121212第1页,共24页 定理1.1 n阶非奇异矩阵A可作三角分解的充要条件是Ak¹0(k=1,2,L,n-1),这里Ak为A的k阶顺序主子阵, 以下同. 证明 必要性. 设非奇异矩阵A有三角分解A=LU, 将其写成分块形式 æAA12öæLk çk÷=çèA21A22øèL210öæUk÷çL22øè0U12ö÷ U22ø这里Ak, Lk和Uk分别为A, L和U的k阶顺序主子阵. 首先由A¹0知L¹0, U¹0, 从而Lk¹0,Uk¹0; 因此Ak=LkgUk¹0(k=1,2,L,n-1). 充分性. 对阶数n作数学归纳法. 当n=1时, A1=,结论成立. 设对n=k结论成立, 即Ak=LkUk, 其中Lk和Uk分别是下三角矩阵和上三角矩阵. 若Ak¹0,则由Ak=LkUk易知Lk和Uk可逆. 现证当n=k+1时结论也成立, 事实上 æAkAk+1=çTèrk0öæUköæLk=÷÷çT-1çTak+1,k+1øçrU1èkkøè0cköL-1kck. ÷T-1-1÷ak+1,k+1-rkUkLkckø由归纳法原理知A可作三角分解. æ01ö定理 1.1 给出了非奇异矩阵可作三角分解的充要条件, 由于A=ç÷不满足定10èø理1.1的条件, 所以它不能作三角分解. 但 æ1æ00öæ00öæ11öæ00öçA=ç÷=ç÷ç÷=ç÷121101èøèøèøè12øç0è1ö÷1÷. 2ø上例表明对于奇异矩阵,它还能作三角分解未必要满足定理1.1的条件. 首先指出,一个方阵的三角分解不是唯一的, 从上面定义来看,杜利特分解与克劳特分解就是两种不同的三角分解,其实,方阵的三角分解有无穷多, 这是因为如果D是行列式不为零的任意对角矩阵, 有 %, A=LU(CD)(D-1U)=LU%也使A的一个三角分解. 因D的任意%,U%也分别是下、上三角矩阵, 从而A=LU其中L性, 所以三角分解不唯一. 这就是A的分解式不唯一性问题, 需规范化三角分解. 定理 1.2 设A为n阶方阵,则A可以唯一地分解为 第2页,共24页 A=LDU (1.1) 的充分必要条件是A的前n-1个顺序主子式Ak¹0(k=1,2,L,n-1). 其中L,U分别是单位下、上三角矩阵, D是对角矩阵D=diag(d1,d2,L,dn), dk=Ak(k=1,2,L,n), A0=1. Ak-1证明 充分性. 若Ak¹0(k=1,2,L,n-1), 则由定理1.1, 即实现一个杜利特分%, 其中L为单位下三角矩阵, U%为上三角矩阵,记 解A=LU(1)æu11u12Lu1nöæa11ç÷çu22Lu2n÷çU=ç=çOM÷çç÷çunnøçèè()a12La1(n)ö÷(2)(2)a22La2n÷(n)÷=A, OM÷(n)÷annø11(i)因为uiiºaii¹0(i=1,2,L,n-1). 下面分两种情况讨论： (1)(2)(n)(n)1) 若A非奇异,由式(1)有Dn=a11=A¹0, 所以anna22Lann=unn¹0, 这时 æ111ö(1)(2)(n),L,a22Lann令D=diaga11, 则D-1=diagç. (n)÷ça(1)a(2)÷annøè1122()于是有 %=LD(D-1U%)=LDU (1.2) A=LU是A的一个LDU分解. (i)(1)(2)(n-1)2)若A奇异,则uiiºaii=0,此时令D=diag(a11,a22,L,an-1,n-1,0), (1)(2)(n-1)Dn-1=diaga11,a22,L,an-1,n-1, a=(u1n,L,un-1,n), T()则 -1%æUaöæU0öæDnn-1n-1-1Un-1%UºçT÷=çT÷ç0øè00øè0Tè0-1öDn-1a÷=DU, 1ø因此不论哪种情况, 只要Ak¹0(k=1,2,L,n-1), 总存在一个LDU分解式(1.1), 第3页,共24页 kdk=akk=Ak(k=1,2,L,n-1),A0=1. Ak-1再证这个分解是唯一的, 仍分两种情况讨论： 1) 当A非奇异时,有A¹0, D¹0, U¹0, L¹0, 所以L、D、U均非奇异. 若还存在另一个LDU分解A=L1DU11, 这里L1, U1也非奇异, 于是有 1, D LDU=L1DU11 (1.3) -1上式两端左乘以L1以及右乘以U-1和D-1, 得 -1-1-1 L1L=DU11UD, (1.4) 但式(1.4)左端是单位下三角矩阵, 右端是单位上三角矩阵, 所以都应该是单位阵, 因此 -1-1L-1L=I, DU=I, 11UD即L1=L,U1U-1=D1-1D. 由后一个等式类似地可得U1U-1=I, D1-1D=I, 即有 U1=U, D1=D. 2) 若A奇异, 则式(1.3)可写成分块形式 %æL1çTèb1%0öæU%aöæL%0öæD111÷çT÷çT÷=çT1øè00øè01øèb%0öæU%0öæD÷ç÷ç1øè00øè0Taö÷, 1ø%, D%, U%, L%是n-1阶单位下三角阵; U%是n-1阶上三角阵; D%是n-1阶对其中L1111角阵; a, a1, b, b1是n-1维列向量. 由此得出 %D%æL11U1çT%èb1D1U1%D%aööæLDULLD11a1÷=çT%÷, T%T%b1D1a1øèbDUbDaø%, D%, U%均非奇异, 类似于前面的推理, 可得 %, D%, U%和L其中L111%=L%, D%=D%, U%=U%, a=a, bT=bT. L11111必要性. 假定A有一个唯一的LDU分解, 写成分块的形式便是 第4页,共24页 xöæLn-10öæDn-10öæUn-1aöæA çn-1, (1.5) ÷=çTç÷çT÷÷annøèbdnøè01øè01øèy其中Ln-1, Dn-1, Un-1, An-1分别是L, D, U, A的n-1阶顺序主子矩阵; x, y, a, b为n-1维列向量. 由式(1.5)有下面的矩阵方程: An-1=Ln-1Dn-1Un-1, (1.6) yT=bTDn-1Un-1, (1.7) x=Ln-1Dn-1a, (1.8) ann=bTDn-1a+dn. (1.9) 否则, 若An-1=0, 则由式(1.6)有An-1=Ln-1Dn-1Un-1=Dn-1=0. 于是有Ln-1Dn-1=Dn-1=0, 即Ln-1Dn-1奇异. 那么对于非其次线性方程组(1.8)有无穷多非零解, 不妨设有a¢, 使Ln-1Dn-1a¢=x, 而a¢=a. TT¢同理, 因Dn-1Un-1奇异, (Ln-1Dn-1)=Un-1Dn-1也奇异, 故有b¹b, 使TTTT¢T¢¢T¢Un-1Dn-1b=y, 或bDn-1Un-1=y. 取dn=ann-bDn-1a, 则有 æAn-1xöæLn-10öæDn-10öæUn-1a¢ö, ç÷=ç¢Tç÷÷ç÷T¢øè0dn1øè01øèyannøèb这与A的LDU分解的唯一性矛盾, 因此An-1¹0. 考察n-1阶顺序主子矩阵An-1由式(1.6)写成分块形式, 同样有An-2=Ln-2Dn-2Un-2. 由于Dn-1¹0, 所以Dn-2¹0, 可得An-2=Ln-2Dn-2Un-2=Dn-2¹0, 从而An-2¹0. 依此类推可得Ak¹0(k=1,2,L,n-1). 综上所述, 定理证明完毕. 推论 1 设A是n阶方阵, 则A可惟一进行杜利特分解的充分必要条件是A的前3n-1个顺序主子式 第5页,共24页 a11La1kAk=MM¹0,k=1,2,L,n-1, ak1Lakk其中L为单位上三角矩阵, 即有 æ1öç÷æu11u12Lu1nöl1÷ç21÷çuLu222n÷÷çA=çl31l32O ç÷OMç÷MMO1÷ç÷çunnøçl÷èlLl1n1n2n,n-1èø并且若A为非奇异矩阵, 则充要条件可换为: A的各阶顺序主子式全不为零, 即: Ak¹0,k=1,2,L,n. 推论 2 n阶方阵A可惟一地进行克劳特分解 3æl11öæ1u12Lu1nöç÷ç÷ll1Lu21222n÷ç÷ %=çA=LUçM÷çMOOM÷ç÷ç÷llLl1ønnøèèn1n2的充要条件为 a11La1kAk=MM¹0, k=1,2,L,n-1. ak1Lakk若A为奇异矩阵, 则lnn=0, 若A为非奇异矩阵, 则充要条件也可换为 Ak¹0, k=1,2,L,n. 定理 1.3 设A为对称正定矩阵, 则A可惟一地分解为A=LDLT, 3其中L为下三角矩阵, D为对角矩阵, 且对角元素是L对角线元素的倒数. 即 æ1çlç11öç÷÷, D=çç÷Oç÷Llnnøçççèö÷÷÷÷÷. ÷O÷1÷lnn÷øæl11çllL=ç2122çMMçèln1ln21l22第6页,共24页 j-1其中lij=aij-ålikljk/lkk, i=1,2,L,n, j=1,2,L,i. k=11.2 常用的三角分解公式 1.2.1 杜利特分解 设A为n阶方阵, 如何确定L和U这两个三角矩阵呢, 设A=LU, 其中 æçl11ö÷æu11u12Lu1nöL=çl21lç22ç÷u22Lu÷2nçMMO÷, U=çç÷ èln1ln2Ll÷çOM÷nnøèu÷nnø按矩阵的乘法, 有 min(åi,j)aij=lisusj,i,j=1,2,Ln s=1k-1由于lkk=1, 所以有akj=ukj+ålksusj, j=k,k+1,Ln. s=1故得 k-1ukj=akj-ålksusj,j=k,k+1,Ln. s=1同理 k-1læöik=çèaik-ålisusk÷øukk, i=k+1,k+2,Ln s=1即得到三角矩阵L和U. 1.2.2 克劳特分解 设A为n阶方阵, 有分解式A=LU, 即 æça11La1jLa1nöçMMM÷æçl11öæ1u12Lu1jLçça÷La=çMO÷çl÷ççOMi1ijLain÷i1Llii÷=ç1uj-1,jLçMMM÷ç÷çç÷çMO÷çOèan1LanjLann÷øçèln1LLLlnn÷øçè1当i³j时, 有 第7页,共24页 u1nöM÷u÷j-1,n÷M÷÷1÷ø aij=ålikukj=ålikukj+lij, k=1k=1ji-1得 lij=aij-ålikukj, i=1,2,Ln, j=1,2,L,i; k=1j-1当i<j时, 有 aij=ålikukj=ålikukj+liiuij, i=1,2,Ln, j=1,2,L,i; k=1k=1ji-1得 iæöuij=çaij-ålikukj÷/lii, i=1,2,Ln-1, j=i+1,i+2,Ln. k=1èø这样即可得到三角矩阵L和U. 1.2.3 乔累斯基(Cholesky)分解 设A为对称正定矩阵, 存在一个实的非奇异下三角矩阵L, 且L的对角元素为正时, 有惟一的分解式 A=LLT. 即 æa11öæl11öæl11Llj1Lln1öç÷ç÷ç÷MOMOOMMç÷ç÷ç÷çai1Laii÷=çli1Llii÷çljjLlnj÷, ç÷ç÷ç÷MMOMMOOMç÷ç÷ç÷çaLaLa÷çlLlLl÷çlnn÷ninnøninnøèèn1èn1ø当i³j时, 有 j-1æöaij=ålikljk=ålikljk+lijljj, 也即lij=çaij-ålikljk÷/ljj, i³j. k=1k=1k=1èøjj-1特别地, 当i=j时, 有 2lij=aii-ålik,i=1,2,L,n. k=1i-1例 1.1 求矩阵 第8页,共24页 æ52-4ç21-2A=çç-4-25çè010的LDU分解和Doolittle分解. 解 对A作矩阵 0ö÷1÷ ÷0÷2ø计算 对A(1)作矩阵 计算 对A(2)作矩阵 æç1öæ1ç2÷÷çç2L=çç51÷-1ç-511, 所以ç4ç-01÷L1=÷ç40ç5÷çç-5è0001÷øçè00æç52-40öç01-2÷L-1A=çç551÷÷1ç29=A(1) ç0-ç550÷÷÷è0102÷øæç1öæ1L01÷ç2=çç÷, L-1012=çç0-21÷ç021è0102÷øçè0-50æç52-40öç1÷L-1(1)05-251÷÷=A(2)2A=ç çç0012÷çè002-3÷÷ø第9页,共24页 ö÷÷÷1÷÷÷01÷øö÷÷÷ ÷ø 1 æ1öæ1öç÷ç÷0101÷, L-1=ç÷, L3=ç3ç001÷ç001÷ç÷ç÷è0001øè00-21ø计算 æ52-40öç÷12ç-1÷(3)-1(2)L3A=ç55÷=A ç12÷ç÷ç-7÷èø令 æ1öç÷2ç÷1ç5÷L=L1L2L3=ç÷ 4-21÷ç-ç5÷ç0521÷èø可得A的Doolittle分解为A=LA(3), A的LDU分解为 24æ5öæö1-0ç÷ç÷551ç÷ç÷1-25A=Lç5÷ç÷. ç÷ç112÷ç÷ç÷ç÷ç÷-71èøèø例 1.2 求三阶方阵 æ2-13öç÷A=ç121÷ ç243÷èø的LU分解与LDU分解. 解 因为 2-1=5¹0, A3=A=5¹0 A1=2¹0, A2=12所以A有唯一的LU分解和LDU分解. 且 第10页,共24页 3öæ1öæ2-1ç÷ç÷(1)-1L1=ç0.51LA=00.25-0.5=A, 1ç÷÷ç0ç101÷50÷èøèø由A(1)可计算L2及A(2)如下: 3öæ1öæ2-1ç÷ç÷(2)-1(1)L2=ç01LA=02.5-0.5=A, 2ç÷÷ç00ç021÷1÷èøèø于是 æ1öç÷L=L1L2=ç0.51÷ ç121÷èø因此A的LU分解为 3öæ1öæ2-1ç÷ç÷A=ç0.5100.25-0.5÷ç÷ ç121÷ç001÷èøèø且A的LDU分解为 æ1öæ2öæ1-0.51.5öç÷ç÷ç÷A=ç0.512.51-0.2÷ç÷ç÷ ç121÷ç1÷ç1÷èøèøèøæ1çça1ça12例 1.3 已知范德蒙(Vandermonde)矩阵V=ççMça1n-2ççan-1è1的三角分解. 11a2a22Ma3a32MLLLa2n-2a3n-2La2n-1a3n-1L1öan÷÷2an÷÷, 求VM÷ann-2÷÷n-1÷anø解 由于范德蒙矩阵V满足定理1.2的条件, 于是有唯一的三角分解: V=LU, 结合范德蒙矩阵的特点, 先对范德蒙矩阵V进行一系列初等行变换. 用n阶矩阵 第11页,共24页 æ1çç-a1ç-1L1=çççççè左乘范德蒙矩阵V得, 1-a11OO-a11-a1ö÷÷÷÷ ÷÷÷1÷ø11æ1ça2-a1a3-a1ç0ç0a2(a2-a1)a3(a3-a1)ç-1MML1V=çMn-4ç0a2(a2-a1)a3n-4(a3-a1)çn-3n-3ç0a2(a2-a1)a3(a3-a1)ç0an-2(a-a)an-2(a-a)221331è记 LLLLLLöan-a1÷÷an(an-a1)÷÷M÷ n-4an(an-a1)÷÷n-3an(an-a1)÷n-2an(a3-a1)÷ø1æ1çç01ç-a21-1L2=çOOçç-a2ççè则 1-a2ö÷÷÷÷ ÷÷÷1÷ø11æ1ça3-a1ç0a2-a1ç00(a3-a2)(a3-a1)-1-1L2L1V=çMMçMç00a3n-4(a3-a2)(a3-a1)çn-3ç00a(a3-a2)(a3-a1)3è一般地,记 LLLMLLö÷an-a1÷(an-a2)(an-a1)÷÷ M÷n-4an(an-a2)(an-a1)÷÷n-3an(an-a2)(an-a1)÷ø1第12页,共24页 æ1ççOç1ç-ak-1çLk=1ö÷÷÷÷÷,k=1,2,L,n-1. çç-ak1÷çOO÷ç÷ç-açk1÷è-ak1÷÷øL-1kd 左上角是k阶单位矩阵. 依次相乘有 æç111L1ça2-a1a3-a1Lan-1-a1çççÕ2(2a3-aj)LÕ(2an-1-aj)an-aj)j=1j=1Õ(j=1L-1LL-1-1n-12L1V=U=ççOLçççÕn-2(an-1-aj)j=1ççè从而 V=L1L2LLn-1U, 其中 æç1öçO÷ç÷ç1÷ Lak1÷k=çç÷a2çkak1÷, k=1,2,L,n-1. (1.10) çLOOO÷ç÷çLLLaçk1÷èan-kkan-k-1kLa2kak1÷÷ø在对U进行一系列初等列变换. 记 第13页,共24页 1öa÷n-a1÷÷÷÷L÷Õn-2÷(an-aj)÷j=1÷Õn-1÷(a)÷n-aj÷j=1øæ1-1-1L-1-1öç÷1ç÷ç÷1-1U1=ç÷ Oç÷ç÷1ç÷ç1÷èøD-1=diagæç1,11,11öèa2-a1a3-a,L,1a÷ n-a1ø有 æç100L00öç11L11÷çça3-a2Lan-1-a2a÷n-a2÷UU-1OLL÷÷1D-1ç1=çççÕn-2(an-1-aj)j=1Õn-2(a)÷n-aj÷j=1ç÷ççèÕn-1(a÷n-aj)÷j=1÷ø一般地,记 æç1öçO÷U-1çç1-11÷L-÷k=çO÷ ç÷ç1÷çè1÷÷ø7k个8D-1æ61ök=diagçç1,1,L1,a,11÷èk+1-aka,L,k+2-akan-ak÷ ø所以 æç1öçO÷ç111÷L÷Dç1÷kUk=çak+1-akç÷, k=1,2,L,n-1. (1.11) çak+2-ak÷çO÷ç÷èan-ak÷ø第14页,共24页 于是 V=L1L2LLn-1(Dn-1Un-1)L(DU11). 其中Lk由式(1.10)给出, 为下三角矩阵. 而DkUk由式(1.11)给出, 为稀疏上三角矩阵. 2 矩阵的满秩分解 2.1 矩阵的满秩分解基本概念与定理 定义2.1 若矩阵A的行(列)向量线性无关, 则称A为行(列)满秩矩阵. 4定义2.2 设A是秩为r的m´n矩阵, 若存在m´r列满秩矩阵F和r´n行2满秩矩阵G, 使得 A=FG 则称式为矩阵A的满秩分解. 定义2.32 设H是m´n的矩阵, rank(H)=r, 满足 1)H的前r行中每一行至少含有一个非零元素, 且每行第一个非零元素是1, 而后m-r行元素均为0; 2)设H中的第i行的第一个非零元素1位于第ji(i=1,2,L,r)列, 有j1<j2<L<jr 3)H的第j1, j2, L, jr列构成m阶单位矩阵I的前r列. 则称H为A的Hermite标准型. 定理2.1 设A为任一秩为r的m´n矩阵, 则A必有满秩分解式A=FG, 其中F为列满秩的, G为行满秩的. 证明 因为A的秩为r, 所以存在m阶可逆矩阵P和n阶可逆矩阵Q, 使得 æErPAQ=çè0若令 0ö÷ 0øæEöF=P-1çr÷,G=(Erè0ø第15页,共24页 0)Q-1 则F为m´n列满秩矩阵, G为r´n行满秩矩阵, 且有 æEA=P-1çrè0结论成立. 若记 0ö-1-1æErö÷Q=Pç÷(Er0øè0ø0)Q-1=FG F=(a1,a2,L,ar), GT=(b1,b2,L,br) 则有 T A=FG=a1b1T+a2b2+L+arbrT (2.2) 这里(2.2)式也是A的满秩分解的一种表示. 定理 2.2 任何非零矩阵AÎPm´n都存在满秩分解. 证明 设r(A)=r>0. 则可通过初等变换将A化为阶梯形矩阵B, 即 æGöA®B=ç÷, GÎPr´n且秩(G)=r. è0ø行于是存在有限个m阶初等矩阵的乘积P, 使得 æGöPA=B=ç÷ 或者 A=P-1B. è0ø于是 æGöA=P-1B=P-1ç÷ è0ø将P-1作相应的分块, P-1=(FS), FÎPm´n, SÎP则有 m´(n-r). A=PB=(F-1æGöS)ç÷=F×G+S×0=FG. è0ø其中F为列满秩矩阵,G为行满秩矩阵. 由于初等行变换有三种变换：1、调换两行;2、某一行乘以一个非零常数;3、某一行乘以一个非零常数加到另一行. 实际上只用第三种初等变换方法就可以将其化为阶第16页,共24页 梯形. 值得指出的是, A的满秩分解式为(2.1)与(2.2)并不是惟一的. 现对任一r阶可逆方阵H, 总有 % (2.3) A=FG=(FH)(H-1G)=FG%分别为m´r列满秩矩阵与r´n行满秩矩阵. 因而(2.3)式也是A的一个%, G成立, 且F满秩分解式. 定理2.33 设AÎCrr´r, 且 % A=BC=BC均为A的满秩分解, 则 1) 存在矩阵QÎCrr´r, 使得 %. %, C=Q-1CB=BQ2) CH(CC)(BB)H-1H-1%HCC%HB=CH()(-1%HB%B)-1%H. B定理2.42 设A是m´n的矩阵, rank(H)=r>0, 其Hermite标准型为H, 则在A的满秩分解中, 可取F为由A的j1, j2, L, jr列构成的m´r的矩阵, G为H的前r行构成的r´n的矩阵. 定理2.5 矩阵满秩分解的存在性定理 1´n1) 设AÎCm(r>0), 则使用初等行变换可将A化为Hermite标准型; r´n2) 设AÎCm(r>0), 则存在FÎCrm´r和GÎCrr´n, 使得A=FG. ræ12-3ö例 2.1 已知A是一个2´3矩阵ç÷, 则A的秩为1, 且它的满秩分解为 è24-6øæ1öA=FG=ç÷(12-3) è2ø显然, 分块矩阵 æErçè00öæErö÷=ç÷(Er0øè0ø0) 第17页,共24页 æ1332öç÷例 2.2 求矩阵A=ç2695÷的满秩分解. ç-1-330÷èø解 由题可知, rank(A)=2>0, 由定理2.4 可得 æ1332öæ1332öæ1332öç÷r3-2r2ç÷r1-r2æGör2-2r1A¾¾¾®0031¾¾¾®0031¾¾¾®，其中G=ç÷ç÷ r3+r1ç÷ç÷0èøè0031øç0062÷ç0000÷èøèøæ100öæ100öæ100öæ100öç÷ç÷ç÷ç÷P=E3E2E1=ç010÷ç010÷ç-210÷=ç-210÷ ç0-20÷ç101÷ç001÷ç1-21÷èøèøèøèøæ100öç÷则P-1=ç210÷对单位下三角矩阵求逆矩阵等于把严格下三角部分元素变号即可. ç-121÷èøæ10öæ10öç÷ç÷æ1332ö取P-1的前两列构成F=ç21÷, 则A=FG=ç21÷ç÷. ç-12÷ç-12÷è0031øèøèøæ1-111öç÷-11-1-1÷的满秩分解. 例2.3 求矩阵çç-1-111÷ç÷11-1-1èøæ1-111öæ1000öç÷ç÷-11-1-101-1-1行÷¾¾÷=B, rank(B)=2且B中的第1解 A=ç®çç-1-111÷ç0000÷ç÷ç÷11-1-10000èøèøæ1-1öç÷-11÷æ1000ö 列和第2列为单位矩阵的前两列, 故A=ç÷ç-1-1÷ç01-1-1èøç÷è11ø 3 矩阵的QR分解 3.1矩阵的QR分解基本概念与定理 第18页,共24页 定义3.15 设uÎRn是单位列向量,即uTu=1, 称矩阵 H=I-2uuH 为Householder矩阵. 由Householder矩阵确定的Rn上的现线性变换y=Hx称为Household变换er. 若u不是单位向量, 则定义H=I-2u22uuT为Householder矩阵, 对应的变换成为Householder变换. Householder矩阵具有如下性质: 1)HT=H; 2)HTH=E; 3)H2=E; 4)H-1=H; æE5)çrè00ö÷是n+1阶Householder矩阵; Hø6)H=-1. 定义3.2 如果实非奇异矩阵A能够转化成正交矩阵Q与实非5奇异上三角矩阵R的乘积, 即 A=QR, 则称上式为A的QR分解. 定理 3.1 任何实的非奇异n阶矩阵A可分解为正交矩阵Q和上三角矩阵R的乘5积,且除去相差一个对角线元素之绝对值全等于一定对角矩阵因子D外, 分解式A=QR是惟一的. 定理3.2 设A为m´n复矩阵(m´n), 且n个列向量线性无关, 则A有分解式 3A=UR