矩陣分解的應(yīng)用淺述

金少華　徐勇　程俊明

【摘要】矩陣論既是數(shù)學(xué)的一個(gè)重要分支，又是一門(mén)具有實(shí)用價(jià)值的數(shù)學(xué)理論.本文給出了矩陣分解的若干應(yīng)用.

【關(guān)鍵詞】矩陣論;正規(guī)矩陣分解;滿(mǎn)秩分解;奇異值分解;特征值

【基金項(xiàng)目】河北工業(yè)大學(xué)研究生《矩陣論》示范課程建設(shè)項(xiàng)目（121031）;2019年河北省研究生示范課程立項(xiàng)建設(shè)項(xiàng)目（KCJSX2019011）

《矩陣論》是工科研究生的一門(mén)重要的學(xué)位課（參見(jiàn)文獻(xiàn)[1]-[4]）.本文給出了矩陣分解的若干應(yīng)用.

一、利用正規(guī)矩陣的分解計(jì)算矩陣函數(shù)

例1 設(shè)矩陣A=100-101-100-110-1001，試求eA，cos A.

解 對(duì)于正規(guī)矩陣A，可求得P-1AP=Λ，其中

P=1001011001-10100-1，? Λ=0022，? P-1=12P，于是

eA[ZK（]=P11e2e2P-1=121+e21-e21+e21-e21-e21+e21-e21+e2.[ZK）]

cos A=P11cos 2cos 2P-1=121+cos 21-cos 21+cos 21-cos 21-cos 21+cos 21-cos 21+cos 2.

二、利用矩陣的滿(mǎn)秩分解求矩陣的特征根

設(shè)n階方陣A的滿(mǎn)秩分解為A=F·G，則由文獻(xiàn)[5]，有

|λEn-A|=|λEn-F·G|=λn-rλEr-G·F.（1）

因此，若A的秩rn，則可通過(guò)（1）式右邊的行列式來(lái)求A的特征根，這將大大簡(jiǎn)化計(jì)算.

例2 求n階矩陣A=22…222…222…2 的特征值.

解 容易看到，矩陣A的秩為1，其滿(mǎn)秩分解為A=22211…1，由等式（1），得 |λEn-A|=λn-1（λ-2n），所以矩陣A的特征值為λ1=λ2=…=λn-1=0，λn=2n.

三、利用矩陣A的滿(mǎn)秩分解求解齊次線(xiàn)性方程組AX=0

結(jié)論1[6] 設(shè)矩陣A的滿(mǎn)秩分解為A=F·G，則 GX=0AX=0.

例3 求解方程組AX=0，其中A=14-1562000-14-1-2-4011-2-1-1-6.

解 先做A的滿(mǎn)秩分解.對(duì)A進(jìn)行行的初等變換，有

A→1000-70100-900106000111，

所以，矩陣A的滿(mǎn)秩分解為A=F·G=14-152000-1-2-401-2-1-11000-70100-900106000111.解方程組GX=0，得方程組AX=0的通解為X=k79-6-111，k∈R.

四、利用矩陣A的滿(mǎn)秩分解求其廣義逆A+

例4 已知A=1i-i0i-1111i-i0，求其廣義逆A+.

解 對(duì)A進(jìn)行行的初等變換，有A→1i-i000010000，

得A的滿(mǎn)秩分解為A=10i1101i-i00001=F·G.

于是A+[ZK（]=G+F+=GHGGH-1FHF-1FH=16101-i0-ii0i-3i6-3i.[ZK）]

五、利用矩陣A的奇異值分解求其廣義逆A+

設(shè)矩陣A的奇異值分解為A=Uσ1σr00VH，則A+=Vσ-11σ-1r00n×mUH.

例5 設(shè)A=111111，求A 的廣義逆A+.

解 首先注意到A的秩為1，AAT的特征值是λ1=6，λ2=0，由兩兩正交的單位特征向量構(gòu)成的正交矩陣可取為U=121-111，又ATA的特征值是λ1=6，λ2=0，λ3=0，由兩兩正交的單位特征向量構(gòu)成的正交矩陣可取為V=33-12-163312-1633026，則A的奇異值分解為A=U600000VT，

從而A+=V1600000UT=16111111.

六、應(yīng)用奇異值分解求解方程組AX=0

設(shè)矩陣A∈Rm×nr的奇異值分解為A=UDVT，D=Σ000，其中矩陣U是m階正交矩陣，矩陣V是n階正交矩陣，Σ=diag（σ1，σ2，…，σr）.于是方程組Ax=0可寫(xiě)為UDVTx=0.該方程兩邊左乘UT，得DVTx=0.令VTx=y，則有Dy=0，從而該方程組的通解為

y=k1er+1+…+kn-ren（k1，…，kn-r∈R）.

設(shè)正交矩陣V的第j個(gè)列向量為vj（j=1，2，…，n），則AX=0的通解為：

x=Vy=k1vr+1+…+kn-rvn（k1，…，kn-r∈R）.

即求出矩陣ATA的屬于特征根0的線(xiàn)性無(wú)關(guān)的特征向量，則這些特征向量的任意線(xiàn)性組合即為Ax=0的通解[6].

例6 求解方程組AX=0，其中A=-112-221-1-2011-2.

解 ATA的屬于特征根0的線(xiàn)性無(wú)關(guān)的特征向量為：ξ1=1-110，ξ2=0201，則方程組AX=0的通解為x=k1ξ1+k2ξ2（k1，k2∈R）.

七、利用矩陣A的廣義逆A+求解方程組Ax=b

對(duì)于給定的方程組Ax=b，只要求出A+，則A+b便給出了該方程組的各種意義下的解，即當(dāng)Ax=b相容時(shí)，A+b或?yàn)槠湮ㄒ唤猓驗(yàn)槠湮ㄒ坏臉O小范數(shù)解，而當(dāng)Ax=b不相容時(shí)，A+b或?yàn)槠湮ㄒ蛔钚《私猓驗(yàn)槠湮ㄒ坏臉O小最小二乘解.從而廣義逆A+將方程組Ax=b的求解問(wèn)題從理論上及方法上都圓滿(mǎn)解決了.

例7 已知A=101011112， b=111，求方程組Ax=b的極小范數(shù)解或極小范數(shù)最小二乘解x0，并指出x0的類(lèi)型.

解 對(duì)A進(jìn)行行的初等變換，有A→101011000，

得A的滿(mǎn)秩分解為A=100111101011=F·G.

于是A+[ZK（]=G+F+=GTGGT-1FTF-1FT=195-41-451112.[ZK）]

所求x0=A+b=19224， 因AA+b=23112≠b， 故x0是極小范數(shù)最小二乘解.

例8 已知A=110101101211，b=314，求方程組Ax=b的極小范數(shù)解或極小范數(shù)最小二乘解x0，并指出x0 的類(lèi)型.

解 先求廣義逆A+，對(duì)A進(jìn)行行的初等變換，有故Ax=b有解，其極小范數(shù)解為x0=A+b=（1 1 0 1）T.

【參考文獻(xiàn)】

[1]徐仲，張凱院，陸全，等.矩陣論簡(jiǎn)明教程[M].北京：科學(xué)出版社，2014.

[2]張紹飛，姚慕生.矩陣論教程[M].北京：機(jī)械出版社，2012.

[3]靳全勤.初等變換的一個(gè)應(yīng)用：矩陣的滿(mǎn)秩分解[J].大學(xué)數(shù)學(xué)，2009，25（5）：195-197.

[4]金少華，徐勇，金大永.矩陣論教學(xué)的幾點(diǎn)注記[J].數(shù)學(xué)學(xué)習(xí)與研究，2019（23）：13.

[5]同濟(jì)大學(xué)應(yīng)用數(shù)學(xué)系.高等代數(shù)與解析幾何[M].北京：高等教育出版社，2005.

[6]金少華，金大永，徐勇.矩陣分解在求解線(xiàn)性方程組中的應(yīng)用[J].高師理科學(xué)刊，2016，36（4）：61.