擬行(列)對稱矩陣的Schur分解及正交對角分解

袁 暉 坪

(重慶工商大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,重慶 400067；經(jīng)濟(jì)社會應(yīng)用統(tǒng)計(jì)重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400067)

矩陣的Schur分解、正交對角分解與Hermite矩陣分解是矩陣的基本分解方法之一,在因子分析、數(shù)理統(tǒng)計(jì)、系統(tǒng)論、控制論、信息論、商務(wù)智能、最優(yōu)化及各種工程問題中應(yīng)用廣泛[1-9].廣義逆矩陣在測量學(xué)、經(jīng)濟(jì)學(xué)、數(shù)值分析、優(yōu)化理論、系統(tǒng)論和控制論等領(lǐng)域應(yīng)用廣泛[10].當(dāng)數(shù)據(jù)矩陣維數(shù)較大時,用計(jì)算機(jī)對其進(jìn)行直接分解,計(jì)算量極大,效率較低.但若能找到矩陣中某一部分與其他部分之間的某種定量關(guān)系,會極大提高計(jì)算效率,因此尋找矩陣的特殊結(jié)構(gòu)關(guān)系有一定的理論意義[11-16].尤其當(dāng)矩陣具有某種行或列對稱性時,矩陣的Schur分解、正交對角分解與Hermite矩陣分解很容易求得,從而可以節(jié)省大量存儲量空間和計(jì)算量.文獻(xiàn)[11-14]研究了行(列)對稱矩陣的QR分解與奇異值分解,文獻(xiàn)[15-16]研究了擬行(列)對稱矩陣的極分解,獲得了一些有意義的結(jié)果.本文進(jìn)一步探討擬行(列)對稱矩陣的Schur分解、正交對角分解與Hermite矩陣分解,給出擬行(列)對稱矩陣的Schur分解、正交對角分解、Hermite矩陣分解及其Moore-Penrose逆的公式,推廣了文獻(xiàn)[7-8]的結(jié)果,擴(kuò)充了文獻(xiàn)[15-16]的結(jié)果.拓寬了應(yīng)用范圍.本文用AH表示矩陣A的共軛轉(zhuǎn)置矩陣,m×n與m×n分別表示m×n實(shí)矩陣集與復(fù)矩陣集,A+表示A的Moore-Penrose逆.

定義1[15]設(shè)A∈m×n,Q1,Q2,…,Qk-1均為m階置換矩陣,則稱

為A的k次擬行對稱矩陣,A稱為其母矩陣.特別地,當(dāng)Q1=Q2=…=Qk-1=Q時,簡記R(A;Q1,…,Qk-1)=Rk(A;Q).

定義2[15]設(shè)A∈m×n,Q1,Q2,…,Qk-1均為n階置換矩陣,則稱

C(A;Q1,…,Qk-1)=(A1,A2,…,Ak-1)(其中Ai=AQi,i=1,2,…,k-1)

為A的k次擬列對稱矩陣,A稱為其母矩陣.特別地,當(dāng)Q1=Q2=…=Qk-1=Q時,簡記C(A;Q1,…,Qk-1)=Ck(A;Q).

由定義1和定義2可知,擬行(列)對稱矩陣是行(列)對稱矩陣[12]、行(列)反對稱矩陣[8]、行(列)延拓矩陣[13]與行(列)酉對稱矩陣[14]的進(jìn)一步推廣.

1 擬對稱矩陣的3種分解及其Moore-Penrose逆

引理1設(shè)Q1,Q2,…,Qk-1均為n階置換矩陣,U為n階酉矩陣,則

為kn階酉矩陣.

證明：因?yàn)閁UH=UHU=I,QQH=QHQ=I,所以容易驗(yàn)證

若無特殊說明,下面出現(xiàn)的酉矩陣P1均與引理1相同.

1.1 擬對稱矩陣的Schur分解及其Moore-Penrose逆

定理1(Schur分解) 設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知A∈n×n的k次擬行對稱矩陣為R(A;Q1,…,Qk-1)∈kn×n,A的Schur分解為A=UHLU,其中U為酉矩陣,L為上(下)三角陣且其對角元為A的特征值.則存在酉矩陣P1, 使得：

證明：1) 由引理1知,P1為酉矩陣,故

2) 由1)及文獻(xiàn)[5]知,

定理1推廣了文獻(xiàn)[7]的定理2及文獻(xiàn)[8]的定理1.

定理2(Schur分解) 設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知A∈Cn×n的k次擬列對稱矩陣為C(A;Q1,…,Qk-1)∈n×kn,A的Schur分解為A=UHLU,其中U為酉矩陣,L為上(下)三角陣且其對角元為A的特征值.則存在酉矩陣P1,使得：

證明：1) 由引理1知,P1為酉矩陣,故

2) 由1)及文獻(xiàn)[5]知,

定理2推廣了文獻(xiàn)[7]的定理3和文獻(xiàn)[8]的定理2.

1.2 擬對稱矩陣的正交對角分解及其Moore-Penrose逆

定理3(正交對角分解) 設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知可逆矩陣A∈n×n的k次擬行對稱矩陣R(A;Q1,…,Qk-1)∈kn×n,且UHAV=D1=diag(σ1,σ2,…,σn),其中σi>0(i=1,2,…,n),U,V為正交矩陣，則存在正交矩陣V和P1,使得:

2) 由1)及文獻(xiàn)[5]知,

定理4(正交對角分解) 設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知可逆矩陣A∈n×n的k次擬列對稱矩陣C(A;Q1,…,Qk-1)∈n×kn,且UHAV=D1=diag(σ1,σ2,…,σn),其中σi>0(i=1,2,…,n),U,V為正交矩陣.則存在正交矩陣V和P1,使得：

2) 由1)及文獻(xiàn)[5]知,

1.3 擬對稱矩陣的Hermite矩陣分解及其Moore-Penrose逆

定理5設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知n階Hermite矩陣A的k次擬行對稱矩陣為R(A;Q1,…,Qk-1)∈kn×n,且A=UHDU,其中D為實(shí)對角矩陣,對角元為A的特征值,U為酉陣.則存在酉陣P1,使得:

證明：1) 由引理1知P1為酉矩陣,故

2) 由1)及文獻(xiàn)[5]知,

定理6設(shè)Q1,Q2,…,Qk-1均為m階置換矩陣,已知n階Hermite矩陣A的k次擬列對稱矩陣為C(A;Q1,…,Qk-1)∈n×kn,且A=UHDU,其中D為實(shí)對角矩陣,對角元為A的特征值,U為酉陣.則存在酉陣P1,使得：

證明：1) 由引理1知P1為酉矩陣,故

2) 由1)及文獻(xiàn)[5]知,

2 擬對稱矩陣Schur分解的算法

算法1擬行對稱矩陣R(A;Q1,…,Qk-1)的Schur分解算法.

步驟1) 求矩陣A的Schur分解A=UHLU；

步驟2) 計(jì)算引理1中的酉矩陣P1；

算法2擬列對稱矩陣C(A;Q1,…,Qk-1)的Schur分解算法.

步驟1) 求矩陣A的Schur分解A=UHLU；

類似地,可得出其他幾種分解的算法,故略.

例1在彩色圖像數(shù)字水印算法中,假設(shè)原始像素塊為

則由定理2知,存在酉矩陣

使得

綜上,本文研究了擬行(列)對稱矩陣的性質(zhì),給出了它們的Schur分解、正交對角分解、Hermite矩陣分解和Moore-Penrose逆的公式及快速算法,導(dǎo)出了擬行(列)酉對稱矩陣與母矩陣的Schur分解、正交對角分解、Hermite矩陣分解和Moore-Penrose逆之間的定量關(guān)系.結(jié)果表明,用母矩陣代替擬行(列)對稱矩陣進(jìn)行Schur分解、正交對角分解、Hermite矩陣分解及Moore-Penrose逆,既能極大減少計(jì)算量和儲存量,又不會喪失數(shù)值精度.因此,本文不僅推廣了文獻(xiàn)[7-8]的相關(guān)結(jié)果,也豐富了文獻(xiàn)[15-16]的研究結(jié)果,拓寬了實(shí)際應(yīng)用領(lǐng)域的范圍.