袁暉坪,呂福起,何 靜,江維瓊,易 強(qiáng),呂希元

(重慶財(cái)經(jīng)學(xué)院 軟件學(xué)院,重慶 401320)

1 引言與預(yù)備知識(shí)

矩陣的極分解在大數(shù)據(jù)、人工智能和數(shù)值分析等領(lǐng)域應(yīng)用廣泛[1-8],如在對(duì)寬帶信號(hào)測(cè)向研究中,通過(guò)對(duì)方向矩陣進(jìn)行極分解構(gòu)造聚焦矩陣,不需譜峰搜索便可確定來(lái)波方向的估計(jì)值,從而極大提高其計(jì)算精度和分辨率[6]. 矩陣的廣義逆在數(shù)理統(tǒng)計(jì)、軟件工程和隨機(jī)規(guī)劃等領(lǐng)域具有重要作用[9],許多實(shí)際問(wèn)題中關(guān)于行(列)或?qū)蔷€的對(duì)稱圖像(矩陣),若用計(jì)算機(jī)直接對(duì)高維數(shù)據(jù)矩陣進(jìn)行分解,則計(jì)算量大、效率低. 若能發(fā)現(xiàn)矩陣具有某種行或列的對(duì)稱性,則問(wèn)題即易解決,故尋找矩陣中某一塊與其他塊之間的結(jié)構(gòu)關(guān)系非常重要[7-15]. 文獻(xiàn)[7-8]探討了行(列)對(duì)稱矩陣及酉對(duì)稱矩陣的極分解; 文獻(xiàn)[12-13] 研究了泛延拓矩陣的QR分解和奇異值分解. 本文進(jìn)一步研究泛延拓矩陣的極分解與廣義逆,給出泛延拓矩陣的極分解與廣義逆的公式和快速算法,并給出泛延拓矩陣極分解的若干擾動(dòng)界. 本文用AH和A+分別表示矩陣A的共軛轉(zhuǎn)置陣與Moore-Penrose逆,m×n表示m×n復(fù)陣集,表示秩為r的m×n復(fù)陣集,‖‖F(xiàn)表示Frobenius范數(shù).

定義1[12]設(shè)A∈m×n,Q1,Q2,…,Qk-1均為m階正交矩陣,稱

為A的k次泛行延拓矩陣,A稱為其母矩陣,其中Ai=QiA,i=1,2,…,k-1. 當(dāng)Q1=Q2=…=Qk-1=Q時(shí),記

R(A;Q1,…,Qk-1)=Rk(A;Q).

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

C(A;Q1,…,Qk-1)=(A,A2,…,Ak-1)

為A的k次泛列延拓矩陣,A稱為其母矩陣,其中Ai=AQi,i=1,2,…,k-1. 特別地,當(dāng)Q1=Q2=…=Qk-1=Q時(shí),簡(jiǎn)記為

C(A;Q1,…,Qk-1)=Ck(A;Q).

顯然,當(dāng)Q1=Q2=…=Qk-1=I(單位矩陣)時(shí),R(A;Q1,…,Qk-1)=Rk(A)即為A的第一類k次行延拓,C(A;Q1,…,Qk-1)=Ck(A)=(A,A,…,A)即為A的第一類k次列延拓[10]; 當(dāng)Q1=Q2=…=Qk-1=J(單位反對(duì)角矩陣)時(shí),R(A;Q1,…,Qk-1)=Rk(A;J)即為A的k次行周期對(duì)稱陣,C(A;Q1,…,Qk-1)=Ck(A;J)即為A的k次列周期對(duì)稱陣[11]；當(dāng)Q1,Q2,…,Qk-1為實(shí)酉變換矩陣時(shí),R(A;Q1,…,Qk-1)即為A的k次行酉對(duì)稱矩陣,C(A;Q1,…,Qk-1)即為A的k次列酉對(duì)稱矩陣[7].

2 泛行(列)延拓矩陣的極分解和廣義逆公式

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

均為kn階酉矩陣.

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

同理可證(P1(U))HP1(U)=Ikn,故P1(U)為kn階酉矩陣. 同理可證P2(U)為kn階酉矩陣. 證畢.

以下若無(wú)特殊說(shuō)明,所有酉陣P1(U),P2(U)均與引理1相同.

引理2[15]設(shè)A∈m×n,則對(duì)任何酉矩陣U∈m×m,V∈n×n均有UAV的Moore-Penrose逆：

(UAV)+=VHA+UH.

證明：1) 由引理1知P1(U)為酉矩陣,因?yàn)?/p>

又由引理1知P2(U)為酉矩陣,因?yàn)?/p>

2) 由1)、引理2及文獻(xiàn)[15]知,

定理2設(shè)正規(guī)陣A∈n×n的極分解為A=HU=UH,其中U為酉陣,H為正定Hermite陣,且AAH=H2,則存在兩個(gè)酉陣P1(U),P2(U)∈kn×kn,使得：

證明：與定理1的證明類似,故略.

即PPH=Ikn,使得

證明：易驗(yàn)證PPH=Ikn,且

即PPH=Ikm,使得

證明：易驗(yàn)證PPH=Ikm,且

3 泛行(列)延拓矩陣極分解的擾動(dòng)分析

引理3設(shè)μ1,μ2,…,μn,ν1,ν2,…,νn均為復(fù)數(shù),則

證明：由復(fù)數(shù)的性質(zhì)及Cauchy-Schwarz不等式,有

引理41) 設(shè)A∈m×n,Bij∈n×s,i,j=1,2,…,k,則

2)

證明：由矩陣范數(shù)的定義及引理3 可知結(jié)論成立.

證明：由定理2、引理4和引理5,知

證畢.

證明：類似于定理5的證明,故略.

泛行延拓陣R(A;Q1,…,Qk-1)的極分解也有類似于定理5和定理6的擾動(dòng)界.

4 泛延拓矩陣的極分解與廣義逆算法

根據(jù)上述討論,可得下列算法.

步驟1) 求矩陣A的的極分解A=UH2；

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

步驟1) 求矩陣A的極分解A=H1U；

步驟2) 計(jì)算定理1中的酉矩陣P2(U)；

類似可得與定理2、定理4、定理5的相應(yīng)分解算法.

5 數(shù)值實(shí)例

使得

則

綜上所述,本文給出了泛行(列)延拓矩陣與母矩陣的極分解、廣義逆與擾動(dòng)界之間的結(jié)構(gòu)關(guān)系及一些計(jì)算公式和算法,結(jié)果表明,用母矩陣代替泛行(列)延拓矩陣計(jì)算極分解、廣義逆與擾動(dòng)界,既能簡(jiǎn)化計(jì)算,又不會(huì)降低數(shù)值精度.