矩阵分解及应用1 引言矩阵是研究图形（向量）变换的基本工具许多数学模型都可以用矩阵表示,矩阵理论既是学习数学的基础,又是一门最有实用价值的数学理论.它不仅是数学的一个重要分支,而且业已成为现代各科领域处理大量有限维空间形式与数量关系强有力的工具.矩阵在代数学习课程中占有重要的地位,而矩阵的分解在矩阵理论研究及其应用中有着重要意义,是其他一些研究课题解决问题的工具.本文介绍了矩阵的几种分解方法:三角分解、正交分解、满秩分解、奇异值分解以及各种分解方法的应用.三角分解在求线性方程组的过程中占有十分重要的作用;正交三角)(QR 分解在计算数学中扮演十分重要的角色,尤其是以QR 分解所建立的QR 方法,已对数值线性代数理论的近代发展起了关键的作用;矩阵的满秩分解和奇异值分解是近几十年来求各类最小二乘问题和最优化问题的重要数学工具. 2 矩阵的三角分解及应用 2.1 杜利特尔分解法定义 2.1]1[ 对于n 阶矩阵A =)(ij a ,n j i ,,2,1, =.如果LU A =,其中L 为单位下三角矩阵,U 为上三角矩阵,则称LU A =为矩阵A 的杜利特尔分解. 确定三角矩阵L 和U 的方法:设LU A =,其中L =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1112121 n n l l l ,U =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n u u u u u u 22211211 按矩阵的乘法有ij a =∑=),m in(1j i s sj isu l,n j i ,,2,1, =由于kk l =1所以有kj a =+kj u ∑-=11k s sj ksu l,n k k j ,,1, +=所以kj u =-kj a ∑-=11k s sj ksu l,n k k j ,,1, +=同理ik l =kkskk s is ik u u l a ∑-=-11,n k k i ,,2,1 ++=这样便可以得到三角矩阵L 和U . 2.2 克劳特分解法定义 2.2]1[ 对于n 阶矩阵A =)(ij a ,n j i ,,2,1, =,如果LU A =,其中L 为下三角矩阵,U 为单位上三角矩阵,称LU A =为矩阵A 的克劳特分解.确定三角矩阵L 和U 的方法:设LU A =,其中L =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡nn n n l l l l l l 21212111,U =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1112112 n n u u u 按矩阵的乘法有ij a =∑=),m in(1j i s sj isu l,n j i ,,2,1, =由于kk u =1所以有ik a =+ik l ∑-=11k s sk isu l,n k k i ,,1, +=所以ik l =-ik a ∑-=11k s sk is u l .n k k i ,,1, +=.同理=kj u kksjk s ks kj l u l a ∑-=-11,1,2,,j k k n =++这样便可以得到三角矩阵L 和U . 2.3 矩阵三角分解的应用例2.1 用LU 分解法求解下列方程组123121325282117x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解 （1）杜利特尔分解法:原方程组的系数矩阵为A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡1102021211 令 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡101001323121l l l 111213222333000u u u uu u ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡1102021211 则有=j u 1j a 1;=1i l 111u a i 1i a = =22u -22a 21l 12u 1=;=23u -23a 21l =13u 2- =32l 22123132a u u l -2-=;=33u 33a 31(l -13u 32l +)23u 3=所以=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡101001323121l l l ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-122011001 , ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡232322131211000u u uu u u ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=300210211 这样原方程组就化为依次求下列两个三角方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-122011001⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321y y y ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=785 （2-1） ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-300210211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321x x x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=321y y y （2-2）解方程组（2-1）可得=1y 5,=2y 3,=3y 3 代入（2-2）可得=1x 2-,=2x 5,=3x 1.（2）克劳特分解法:令 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡333231222111000l l l l l l ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡10010*******u u u ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1102021211 则有=1i l 1i a ;=i u 1111l a ii a 1= =22l 22a 21l -12u 1= ;=32l 31l 12u 2-= =23u 22132123l u l a -2-=;=33l -33a 31(l +13u 32l 3)23=u所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡322011001000333231222111l l l l l l ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡10021021110010*******u u u 这样原方程组就化为依次求下列两个三角方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-322011001⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321y y y ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=785 （2-3） ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-100210211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321x x x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=321y y y （2-4） 则由方程组（2-3）可得51=y ,=2y 3,=3y 1 代入（2-4）可得21-=x ,52=x ,13=x .由上述解题过程可以看出,为求方程组的解,可先把方程组的系数矩阵进行分解,方程组化为b LUX =,则求解此方程组要化为依次求方程组Ly bUx y =⎧⎨=⎩,由Ly b =可求出y ,再将y代入Ux y =，求出x 即得方程组的解.这两种方法在解更多元的方程组中比其它的消元法更加方便.3 矩阵的正交三角)(QR 分解及应用 3.1 矩阵的正交三角)(QR 分解定义3.1]3[ 如果实（复）非奇异矩阵A 能够化成正交（酉）矩阵Q 与实（复）非奇异上三角矩阵R 的乘积,即QR A =,则称式QR A =为矩阵A 的QR 分解.把一个矩阵A 进行正交分解的方法:设n 阶实（复）非奇异矩阵A 的n 个列向量依次为n 21α,,α,α ,由于A 非奇异,所以这n 个向量线性无关,将它们正交化,可得到n 个标准正交列向量n q q q ,,,21 ,将n21α,,α,α 正交化,可得⎪⎪⎩⎪⎪⎨⎧---=-==--1111,1212211ββαββαβαβn n n n n n k k k 其中,(,)()(,)i j ij j j k j i =<αβββ,将上式改写为⎪⎪⎩⎪⎪⎨⎧+++=+==--nn n n n n k k k βββαββαβα11,112121211 用矩阵形式表示为)(n 21α,,α,α =C n ),,,(21βββ ,其中C =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1112121 n n k k k 再将n βββ,,,21 单位化,可得1i i i=q ββ,(1,2,)i n =…, 于是有()=12n α,α,,α()n C 12=β,β,,β()n 12q ,q ,,q ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡n βββ21C 令⎩⎨⎧⋅==C diag R Q n n ),,,(),,,(2121βββq q q （3-1） 则Q 是正交（酉）矩阵,R 是上三角矩阵,且有QR A =.例3.1 求矩阵011110101A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦的QR 分解. 解 令 1230111,1,0101⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ααα 将它们进行正交化可得1122133212130111121,,2232311223⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥===-==--=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥-⎢⎥⎣⎦⎣⎦βαβαββαββ 根据（3-1）构造矩阵0Q ⎡⎤⎢⎥⎢⎥=⎥⎥⎦111221131R ⎤⎡⎤⎥⎢⎥⎥⎢⎢⎥⎢⎥⎢⎢⎥==⎢⎥⎢⎢⎥⎢⎥⎢⎢⎥⎢⎢⎢⎥⎢⎢⎢⎥⎣⎦⎣⎣则有QR A =.3.2 矩阵正交三角()QR 分解的应用通过对实系数n 次方程求根问题及上Hessenberg 阵的QR 分解理论的分析和研究,将实系数n 次方程的求根问题转化为求上Hessenberg 阵的特征值问题,利用上Hessenberg 阵的QR 分解法求实系数n 次方程11100n n n n a x a x a x a --++⋅⋅⋅++=的全部实根.对于实系数n 次方程1110()0n n n n n P x a x a x a x a --=++⋅⋅⋅++=,令ii na b a =,(0,1,2,,1)i n =-则()0n P x =化为首项系数为1的代数方程1110()0n n n n Q x x b x b x b --=++⋅⋅⋅++=由线性代数可知,()0n Q x =可以看成是实矩阵A 的特征多项方程0E A λ-=, 即1110()0n n n n Q b b b λλλλ--=++⋅⋅⋅++=,其中⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=---01001000000100000101321b b b b b A n n n 因此,求方程()0n Q x =的全部实根就转化为求上述矩阵A 的全部实特征值.可以直接用QR 分解法求出矩阵A 的全部实特征值. 4 矩阵的满秩分解及应用 4.1 矩阵的满秩分解定义 4.1[5] 设矩阵(0)m n r A C r ⨯∈>,如果存在矩阵m r r n r r F C G C ⨯⨯∈∈和,使得A FG =,则称式A FG =为矩阵A 的满秩分解.对矩阵A 作满秩分解的方法:设矩阵(0)m n r A C r ⨯∈>,rankA r =,根据矩阵的初等变换理论,对矩阵A 进行初等行变换,可将矩阵A 化为阶梯形矩阵B ,即0G A B ⎡⎤→=⎢⎥⎣⎦行,r n r G C ⨯∈,于是存在有限个m 阶初等矩阵的乘积,记作P ,使得PA B = 或者 A P B -1=将P -1分块为[]P FS -1=,m r m r n rF C S C ⨯⨯-∈∈（n-r ）， 则有[]10G A P B FS FG -⎡⎤===⎢⎥⎣⎦其中F 是列满秩矩阵,G 是行满秩矩阵.例4.1 求矩阵02042100036330211441j j jA j j +⎡⎤⎢⎥=----⎢⎥⎢⎥-⎣⎦（j =. 解 根据上述的理论提供的算法,需要求出阶梯形矩阵B 及诸初等矩阵的乘积P ,为此,对矩阵[]A I 进行初等行变换,当A 所在的位置成为阶梯形矩阵B 时,I 所在的位置就是进行初等变换对应的初等矩阵的乘积P .[]=I A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----+10014411200103363000001124020j j jj j ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-+---−→−13100000003101210000021212102110j j j j j 行所以11100010122221000121003000000113j j j B j j ⎡⎤-⎡⎤⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥=+-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦，P=可求得1200030211j P -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦所以201101012032200012121j j j A j ⎡⎤⎡⎤--⎢⎥⎢⎥=-⎢⎥⎢⎥+⎢⎥⎣⎦⎣⎦.4.2 矩阵满秩分解的应用利用矩阵的满秩分解处理一些问题时,有时会十分方便,下面介绍运用矩阵的满秩分解,给出线性方程组最小二乘解和极小范数最小二乘解.定理 4.1[3] 设矩阵{13}m n m A C b C A A ⨯∈∈∈（1，3），，，,则x A b =（1，3）是方程组b Ax =的最小二乘解.定理 4.2[6] 设矩阵m n m A C b C ⨯∈∈，,则x A b +=是方程组Ax b =的极小范数最小二乘解.例4.2 取例4.1中的矩阵A 和11j b -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,求方程组Ax b =的极小范数最小二乘解. 解 由例4.1知矩阵A 的一个满秩分解式为201101012032200012121j j j A FG j ⎡⎤⎡⎤--⎢⎥⎢⎥==-⎢⎥⎢⎥+⎢⎥⎣⎦⎣⎦ 于是有82210HF F ⎡⎤=⎢⎥⎣⎦,H 1339222GG 39722j j ⎡⎤-⎢⎥=⎢⎥⎢⎥+⎢⎥⎣⎦从而H H H H F F F GG G A 11)()(--+=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡--+-⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-+=130202411538113939314461121221100210100j j j j j j⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--+-++-------+-+---+---=j j j j j j j j jj j jj j j 2091177138137399682636709481121655690960543973918120108781461800017481所以方程组Ax b =的极小范数最小二乘解为⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+-------==+j j j j j b A x 59462129551813362608741.5 矩阵的奇异值分解及应用矩阵的奇异值分解不仅是矩阵理论和矩阵计算的最基本和最重要的工具之一,而且在控制论、系统辨识、信号处理、最优化问题、特征值问题、最小二乘问题及统计学等方面都有直接而重要的应用,以下介绍矩阵奇异值分解及在广义逆中的应用. 5.1矩阵的奇异值分解定义5.1[6] 设矩阵(0)m n r A C r ⨯∈>,H A A 的特征值为12r r 1n 0λλλλλ+≥≥≥>===则称i σ=n i ,,2,1 =）为矩阵A 的奇异值.定义5.2]6[ 设矩阵(0)m n r A C r ⨯∈>,则存在m 阶酉矩阵U 和n 阶酉矩阵V ,使得000H U AV ⎡⎤=⎢⎥⎣⎦∑（5-1）其中∑=),,,(21r diag σσσ ,而i σ),,2,1(r i =为矩阵A 的全部非零奇异值,则000H A U V ⎡⎤=⎢⎥⎣⎦∑为矩阵A 的奇异值分解.对一个矩阵A 作奇异值分解的方法:设矩阵(0)m n r A C r ⨯∈>,记Hermite 矩阵H A A 的特征值为12r r 1n 0λλλλλ+≥≥≥>===则存在n 阶酉矩阵V ,使得120()00H H n V A A V λλ⎡⎤⎡⎤⎢⎥==⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦∑ （5-2）将V 分块为[]()1212,,n r n n r r n rV V V V C V C ⨯⨯--=∈∈ 并改写式（5-2）为2000HA AV V ⎡⎤=⎢⎥⎢⎥⎣⎦∑则有2112,0H H A AV V A AV ==∑（5-3） 由式（5-3）的第一式可得211H H V A AV =∑ 或者 1111)()Hr AV AV I --=∑∑(由式（5-3）的第二式可得22()()0H AV AV = 或者 2AV =0令111U AV -=∑,则11H r U U I =,即1U 的r 个列是两两正交的单位向量,记作1(,,,)r U 12=υυυ.所以可将,,,r 12υυυ,扩充为m C 的标准正交基,记增添的向量为,,r m +1υυ,并构造矩阵2,,r m U +1=(υυ),则[]121,,,,,r r m U U U 12+==(υυυυυ)，是m 阶酉矩阵,且有11210H H r U U I U U ==， 于是可得 []111121212000H H HHH H U U U U AV UAV AV U U UU ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦∑∑∑=== 000⎡⎤⎢⎥⎣⎦∑= 改写式（5-1）为000HA U V ⎡⎤=⎢⎥⎣⎦∑所以矩阵A 可以进行奇异值分解,其分解式为 HV U A ⎥⎦⎤⎢⎣⎡=∑000. 例5.1 求矩阵100111A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦的奇异值分解. 解 2112TB A A ⎡⎤==⎢⎥⎣⎦的特征值是13λ=,21λ=,对应的特征向量依次为 111ξ⎡⎤=⎢⎥⎣⎦, 211ξ-⎡⎤=⎢⎥⎣⎦于是可得2=rankA,001⎤=⎥⎣⎦∑所以存在正交矩阵V =⎥⎥⎦此时 V V =1,计算1110U AV -⎥==⎥⎥⎥⎥⎦∑构造 []120U U U ==⎥⎥⎦ 则A 的奇异值分解为00100H A U V ⎤⎥=⎢⎥⎢⎥⎣⎦.5.2 奇异值分解在广义逆中的应用定义5.2[10] 矩阵m n A C ⨯∈,满足方程AXA A = (1) XAX X = (2) ()H AX AX = (3)()H XA XA = (4)中的)(i )(,j ,, )(l 的n m X C ⨯∈称为A 的一个},,,{l j i -逆,记作(,,,)i j l A ,另外还可以记作{,,,}X A i j l ∈,其中{,,,}A i j l 表示(,,,)i j l A 的全体,特别记(1,2,3,4)A A +=,称A +为A 的Penrose Moore -广义逆,有了奇异值分解,我们就能很容易地利用矩阵A 的奇异值分解表示出矩阵A 的各种广义逆(,,,)i j l A .矩阵A 的奇异值分解的应用:设m n A C ⨯∈有奇异值分解HV U A ⎥⎦⎤⎢⎣⎡=∑000（5-4） 其中U 和V 均为酉矩阵,01>⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=∑r σσ . )(A rank r =,对于 X ＝(,,,)i j l A ,记⎥⎦⎤⎢⎣⎡==22211211^X X X X XU V X H 则)1HU X X X V A⎥⎥⎦⎤⎢⎢⎣⎡=-∑2221121)1( （5-5）)2HU X X X X V A⎥⎥⎦⎤⎢⎢⎣⎡=∑∑-122121121)2,1( （5-6） )3HU X X V A⎥⎥⎦⎤⎢⎢⎣⎡=-∑22211)3,1(0 （5-7）)4HU X X V A⎥⎥⎦⎤⎢⎢⎣⎡=-∑22121)4,1(0 （5-8）)5HU X V A⎥⎥⎦⎤⎢⎢⎣⎡=-∑00211)3,2,1( （5-9）)6HU X V A⎥⎥⎦⎤⎢⎢⎣⎡=-∑00121)4,2,1(（5-10） )7HU V A ⎥⎥⎦⎤⎢⎢⎣⎡=-+∑0001（5-11）证明 )1 由式AXA A =可知,HHHV U V XU V U ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡∑∑∑000000000 即⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡∑∑∑00000000022211211X X X X 故可知∑∑∑=11X ,即111-∑=X ,则式（5-5）得证.)2 把（5-4）式和（5-5）式代入式XAX X =可得 HHU X X X V U X X X X X X V ⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡---∑∑∑∑222112122211212221121000 即⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡--∑∑∑2221121122121121X X X X X X X 亦即122122X X X ∑=,则（5-6）式得证.)3 由（5-4）式和（5-5）式可知 HHU X I U U X X X U AX ⎥⎦⎤⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=∑∑∑-00000122221121但由式H AX AX )(=,故必有12X ＝0,则（5-7）式得证.)4 由（5-4）式和（5-5）式可知 HHH V X I V V U U X X X V XA ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=∑∑∑-0000021*******但由式()HXA XA =,则21210()000HHHI I X XA V V V V X⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎣⎦∑∑ 故必有021=X ,则（5-8）式得证.)5 由（5-4）式和（5-6）式可知HH H U X IU U X X X X V V U AX ⎥⎦⎤⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=∑∑∑∑-0000021211221121但由式()H AX AX =,则有HH H U X I U U X IU AX ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=∑∑0000)(1212 所以有012=X .由（5-4）式和（5-7）式可知HHH H U X V U X X V V U U X X V XAX ⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=---∑∑∑∑00000002112221122211但由式XAX X =,则有HHU X X V U X V ⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡--∑∑22211211000 所以有022=X ,则（5-9）式得证.)6 由（5-4）式和（5-6）式可知HHH V X I V V U U X X X X V XA ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=∑∑∑∑-0000021211221121但由式()H XA XA =,则有HHV X I V V X IV ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∑∑00002121 所以有021=X .由（5-4）式和（5-8）式可知HHH H U X V U X X V V U U X X V XAX ⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=---∑∑∑∑00000001212212122121但由式XAX X =,则有HHU X X V U X V ⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡--∑∑2212112100所以有022=X ,则（5-10）式得证.)7 由（5-4）式和（5-9）式可知HHH V X I V V U U X V XA ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=∑∑∑-000000021211但由式()H XA XA =,则有HH V X IV X I V ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∑∑00002121 所以有021=X .由（5-4）式和（5-10）式可知HHH U X I U U X V V U AX ⎥⎦⎤⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=∑∑∑-000000012121但由式()H AX AX =,则有HH U X I U U X IU ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∑∑00001212 所以有012=X ,由)5的证法可知022=X .则（5-11）式得证. 6 结束语本文论述了矩阵的分解及这些分解在其它学科中的应用，在论述过程中着重讨论了三角分解、QR 分解、满秩分解及奇异值分解的具体应用.对于矩阵的分解可以借助计算机和一些数学软件进行分解，利用这些工具可以很方便的求出矩阵的某一种分解，这样在计算矩阵分解时就可以节省很多的时间.但矩阵的分解并不是只有这四种分解，还存在有其它的分解及其相应的应用，在这就不再做深入的研究了，要进行更深入的分析，使其更加完善，有待于今后的进一步研究.参考文献[1] 石东洋.数值计算方法[M].郑州:郑州大学出版社,2007.[2] 关红钧,苏艳华.关于n阶矩阵的三角分解[J].沈阳航空工业学院学报,2001,18(4)：38-39.[3] 程云鹏,张凯院,徐仲.矩阵论[M].西安:西北工业大学出版社,1999.[4] 禹利萍.用QR分解法求实系数代数方程的全部实根[J].承德民族师专学报,2001,21(2)：14-15.[5] 刘轩黄.矩阵的满秩分解及其应用[J].江西电力职工大学学报,1999,12(4)：5-7.[6] 戴华.矩阵论[M].北京:科学出版社,2001.[7] 黄有度,狄成恩,朱士信.矩阵论及其应用[M].合肥:中国科学科技大学出版社,1995.[8] 张凯院,徐仲等.矩阵论典型试题解析及自我测试题[M].西安:西北工业大学出社,2001.[9] 燕列雅.初变换进行矩阵的QR分解[J].数学通讯报,1998,12(1):13-15.[10] 骈俊生.奇异值分解在广义逆中的若干应用[J].阜阳师范学院学报,2005,22(1)：16-17.[11] 刘轩黄.满秩分解的两种简单实现方法及应用[J].华中师范大学学报,1986,20（4）:425-426.[12] 北京大学数学系.高等代数[M].北京:高等教育出版社,1978.。
![矩阵分解的研究[开题报告]](https://img.taocdn.com/s1/m/b0eb5e464028915f814dc253.png)
