矩阵方程AXB=C

矩阵方程AXB?C的定秩解及其最佳逼近问题

第1章 绪论

对于矩阵方程AXBT?C,刘瑞娟[2]利用矩阵对的商奇异值分解,得到了解非空的充要条件和解的最大(小)秩.1955年Penrose[5]得到了AXB?C有解的充要条件和通解的表达式;1970年Lnacaster[6]利用Kornecker乘积和拉直映射也得到了它的一般解的条件和显式解;1976年，S.K.Mitra[7]研究了它的Hermitian解及非负定解的条件,并给出了两种解的表达式;1990年,戴华[8]研究了此方程在对称矩阵集合中的求解问题，利用矩阵对的广义奇异值分解得到了问题有对称解的条件和解的表达式；邓远北、彭向阳、雷渊[9?12]等系统地研究了此方程在对称化矩阵、（反）对称矩阵、半正定矩阵、正交（反）对称矩阵、（反）自反矩阵集合中的等式约束解以及最小二乘解；2003年,廖安平[13]利用广义奇异值分解研究了它在对称半正定矩阵集合中的最小二乘解;2004年,彭亚新[14]用迭代法系统地研究了矩阵方程AXB?C的一般解 对称解 (反)中心对称解 (反)自反解 双对称解与对称次反对称解等问题.

对于矩阵方程的定秩解问题及其最佳逼近,1972年，S．K．Mitra[16]提出了线性矩阵方程(组)的定秩求解问题；于1984年，Sujit Kumar Mitra[17]利用空间有关理论及秩的相关不等式，给出了矩阵方程组AX?C,XB?D 的极小秩及其它定秩的通解;Uhlig[18]于1987年给出了矩阵方程 AX?B的可能秩的解；于1990年Mitra[15]研究了矩阵方程组A1XB1?C1 ,A2XB2?C2 的公共解的最小秩；Gross

[19]使用了广义逆给出了矩阵方程 AXA*?B的最大秩和最小秩的

Hermitian非负定解;Xiao Q F,Hu X Y,Zhang L[35]于2009年研究了矩阵方程

AX?B的对称最小秩解和最佳逼近解;2007年,雷渊[12]利用矩阵对的广义奇异,AXAT?BYBT?D , ?ATXA,BTXB???C,D? , ?AXB,GXH???C,D?

值分解和标准相关分解研究了以下不相容矩阵方程和矩阵方组

AXB?D的最小二乘问题等价转换为相容矩阵方程的求解问题,并得到了相应的最小二乘解的表达式;于2009年，肖庆丰[4]对矩阵对利用广义奇异值分解，研究了矩阵方程AX?B的自反、反自反、中心对称、反中心对称矩阵的定秩求解的问题及最小秩解的最佳逼近问题，还研究了矩阵方程 AX?B 的对称、反对称矩阵反问题的定秩求解问题，并得到了相应的成果.

第2章 秩约束下矩阵方程AXB2.1引言

?C的一般解及其最佳逼近

秩约束下矩阵的求解问题是线性代数的重要研究课题之一.关于矩阵方程AXB?C的在秩约束下的一般解、最小秩和最小秩解的通式及其最佳逼近解、最大秩和最大秩解的通式.本章采用了RSVD分解,对矩阵方程AXB?C的一般解的秩的情况进行了详细的分析,成功地获取了解的最大(小)秩及其定秩解的表达式.利用相应的结果,也获得了最小秩解的表达式和唯一最佳逼近解的表达式.

2.2 提出问题

问题I 给定矩阵A?Rm?n,B?Rn?m,C?Rm?m，非负整数s，

(i)求X?Rn?n,使得AXB?C，且rank(X)?s; (ii)若有解,设SX?{X?R~X?SXn?n|AXB?C},求m,M~使得 ,

m?minrank(X),M?maxrank(X)X?SX以及S?{X|rank(X)?m,X?SX}.

~~m问题II 给定X?R?n?n,求X?S使得

~~m ||X?X?||?min||X?X?||

X?Sm~2.3 解决问题 2.3.1问题Ⅰ的解

给定矩阵A?Rm?n,B?Rn?m,C?Rm?m，记:

ra?rank(A),rb?rank(B),rc?rank(C)

rca?rank?CA?,rcb?C??C?rank??,rcab?rank??B??BA??0?

k1?rcab?ra?rb,k2?rca?rb?rcab

k3?rcb?ra?rcab,k4?rc?rcab?rca?rcb

引理 2.3.1 (RSVD定理[1])给定矩阵A?Rm?n,B?Rn?m,C?Rm?m,则存在非奇异

矩阵M?Rm?m,N?Rm?m及正交矩阵U?ORn?n,V?ORn?n，使得

A?M?AU,B?V?BN,C?MT?CN (1.1)

其中

??0?0??0??0??0?0?n?r?a????????????k1?00I000k3000I00k4?A?0?k10?k2?0?k30?k4?I?rca?rc?0m?rca?rca?rc?，??????B?????0000k1k20I000000k300I0k4??m?r00b?00?k2， ?k004?I0?rcb?rcrcb?rcm?rcb??I000000I0000k200I000k3k4000SCAB00000000rcb?rc?C?0?k1?0?k20?k3?0?k4?r?rc0?ca0?m?rcam?rcb??4，

其中SCAB?diag{?1,?,?k},?1??2????k?0.

4于是有如下定理:

定理2.3.1 设矩阵A?Rm?n,B?Rn?m,C?Rm?m分解如引理1，非负整数s，则矩阵方程AXB?C有解X?Rn?n的充要条件是:

?k1?0??k2?0,?k?0?3?rcab?ra?rb?0?即?rca?rb?rcab?0， (1.2) ?r?r?r?0acab?cb若矩阵方程AXB?C有解,则有

(1) 矩阵方程AXB?C的通解表达式为:

??Y11 X?U?Y31??Y41?Y13SCAB0?Y14?T0?V, (1.3) ?0??其中Y11,Y13,Y14,Y31,Y41是具有相应维数的任意矩阵,U,V,SCAB见引理1.

(2) 矩阵方程AXB?C的解矩阵中最小秩及最小秩解为:

最小秩为m?rc?rcab?rca?rcb, (1.4) 且最小秩解为:

??1YSY1331CAB? X?U?Y31?0??Y13SCAB0?0?0?V, (1.5) ?0??T~其中Y13,Y31是具有相应维数的任意矩阵,SCAB,U,V见引理1.

(3) 矩阵方程AXB?C的解矩阵中最大秩及最大秩解为: 最大秩为,要考虑两种情况：

第一，当min{n?rc,m?rc}?min{m?rc,rca?rc}?min{n?rc,rcb?rb}时， M1?maxrank(X)?rc?rcab?rca?rcb?min{n?rc,m?rc}， (1.6)

X?SX且相应的最大秩解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?0?V?0??T, (1.7)

?1Y13为行满其中,Y14,Y41是具有相应维数的任意矩阵，我们总可以选取Y11?Y31SCAB秩或列满秩,SCAB,U,V见引理1.

第二，当min{n?rc,m?rc}?min{m?rc,rca?rc}?min{n?rc,rcb?rb}，

M2?maxrank(X)?rc?rcab?rca?rcb?min{n?rc,rcb?rb}?min{m?rc,rca?rc},

X?SX(1.8)

且相应的最大秩解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?0?V?0??T, (1.9)

其中,Y11,Y13,Y31是具有相应维数的任意矩阵，我们总可以选取Y14,Y41为行满秩或列满秩,SCAB,U,V见引理1.

(3) 也要分两种情况:

~?s?M,矩阵方程AXB?C的解矩阵中具有秩为s的解为: 第一,对于m1??Y11X?U?Y31??Y41?Y13SCAB0?Y14?T0?V, (1.10) ?0??其中Y41,Y14是具有相应维数的任意矩阵,SC,U,V见引理1,我们可以选Y11,Y31,Y13

?1Y13)?s?rc?rcab?rca?rcb. 使得rank(Y11?Y31SCAB~?s?M,矩阵方程AXB?C的解矩阵中具有秩为s的解为: 第二,对于m2??Y11X?U?Y31??Y41?Y13SCAB0?Y14?T0?V, (1.11) ?0??其中Y11,Y31,Y13是具有相应维数的任意矩阵,SC,U,V见引理1,我们可以选Y41,Y14 使得rank(Y14)?rank(Y41)?s?rc?rcab?rca?rcb.

证明 给定A?Rm?n,B?Rn?m,C?Rm?m,RSVD分解由引理1给出，则

C?AXB?M?CN?M?AUTXV ?M(?令Y?UTXV,则有

rank(C?AXB)?rank(M(? ?rank(?CCC?BBN

??AUXVT?)N

??AUXVAT?B)N)

??Y?B)

将Y作如下分块:

?Y11?Y21 Y???Y31??Y41Y12Y22Y32Y42Y13Y23Y33Y43Y14??Y24?, Y34??Y44?则

??I0?I?0?0?Y22???0?Y32?0?Y42??00?k1k2?00I000k3?C??Y?AB?000?k1?000?k2?Y23?Y240?k3?SCAB?Y33?Y340?k4??Y43?Y440rca?rc?000?m?rcak4rcb?rcm?rcb??，

在上式中，由Yij(i?2,3,4,j?2,3,4)的任意性，得到

minrank(C?AXB)?minrank(?C??AY?B)

?k1?k2?k3 ?rca?rcb?rcab

因此，矩阵方程AXB?C是可解的?minrank(C?AXB)?0

X?k1?0?rcab?ra?rb?0?? ??k2?0,即?rca?rb?rcab?0

?k?0?r?r?r?0acab?3?cb(1)易知若矩阵方程AXB?C是有解,即k1?k2?k3?0时,则解的一般表达式为

X?UZVT,其中

??YYYn?rc1314??110?k4, Z??Y31SCAB??Y00?41?rca?rc?m?rck4rcb?rb?其中Y11,Y13,Y14,Y31,Y41是具有相应维数的任意矩阵,SCAB,U,V见引理1. (2)记SX为相容矩阵AXB?C的解集合,即 SX?{X|X?UZVT};

由高斯块变换，得

??Y11 ?Y31??Y41?Y13SCAB0?Y14??1?Y11?Y31SCABY13?0??T??0??0?Y41??0SCAB0Y14??0?, 0??则有

X?SXminrank(X)?minrank(Z)?minrank(T)?k4，

Z?1因此有，Y11?Y31SCABY13,Y14?0,Y41?0

容易得知:

~?minrank(X)?minrank(Z)?k?r?r?r?r,且最小秩解为: mccabcacb4X?SXZ??1YS?31CABY13 X?U?Y31?0??Y13SCAB0?0?0?V, ?0??T其中Y13,Y31是具有相应维数的任意矩阵,SCAB,U,V见引理1. 记Sm~为相容矩阵AXB?C的最小秩解集合,即

~ Sm~?{X|rank(X)?m,X?SX}

(3)由(2)中的矩阵T可得, 矩阵方程AXB?C的解矩阵中最大秩为: 要考虑两种情况：

第一，当min{n?rc,m?rc}?min{m?rc,rca?rc}?min{n?rc,rcb?rb}时， M1?maxrank(X)?rc?rcab?rca?rcb?min{n?rc,m?rc}，

X?SX且相应的最大秩解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?0?V?0??T,

?1Y13为行满其中,Y14,Y41是具有相应维数的任意矩阵，我们总可以选取Y11?Y31SCAB秩或列满秩,SCAB,U,V见引理1.

第二，当min{n?rc,m?rc}?min{m?rc,rca?rc}?min{n?rc,rcb?rb}，

M2?maxrank(X)?rc?rcab?rca?rcb?min{n?rc,rcb?rb}?min{m?rc,rca?rc},

X?SX

且相应的最大秩解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?0?V?0??T,

其中,Y11,Y13,Y31是具有相应维数的任意矩阵，我们总可以选取Y14,Y41为行满秩或列满秩,SCAB,U,V见引理1.

(4) 也要分两种情况:

~?s?M,我们可以选Y,Y,Y使得 第一,对于m1131131rank(Y11?Y31SCABY13)?s?rc?rcab?rca?rcb,

?1则矩阵方程AXB?C的解矩阵中具有秩为s的解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?T0?V, ?0??其中Y41,Y14是具有相应维数的任意矩阵,SC,U,V见引理1.

~?s?M,我们可以选Y,Y使得 第二,对于m24114rank(Y14)?rank(Y41)?s?rc?rcab?rca?rcb,

则矩阵方程AXB?C的解矩阵中具有秩为s的解为:

??Y11X?U?Y31??Y41?Y13SCAB0?Y14?T0?V, ?0??其中Y11,Y31,Y13是具有相应维数的任意矩阵,SC,U,V见引理1 .证毕.

第三章 矩阵方程AXB3.1 引言

?C的自反解及其最佳逼近

秩约束下矩阵的求解问题是线性代数的重要研究课题之一.关于矩阵方程

的在秩约束下的对称解、最小秩和最小秩对称解的通式及其最佳逼近

解、最大秩和最大秩对称解的通式.本章采用了RSVD分解,对矩阵方程AXB?C的对称解的秩的情况进行了详细的分析,成功地获取了对称解的最大(小)秩及其定秩解的表达式.利用相应的结果,也获得了最小秩对称解的表达式和唯一最佳逼近解的表达式.

AXB?C3.2 提出问题

n?n问题Ⅰ给定矩阵A?Cm?n,B?Cn?m,C?Cm?m， 求X?CR(P),使得AXB?C，

若有解,记SX?{X?Cn?n|AXB?C}

问题II 给定X?C?n?n,求X?SX使得

~X?SX~ ||X?X?||?min||X?X?|| 3.3问题Ⅰ的解

设P?C义反射矩阵.

设X?Cn?n,若X满足X?PXP,则称X为关于P的自反矩阵.所有关于P的

n?nn?n自反矩阵的全体记为CR(P)?{X?C|X?PXP}.

n?n，且PH?P,P2?I，则称P为广义反射矩阵.本文中的P均为广

引理3.3.1[33] 矩阵X?CRn?n(P)的充要条件是X可以表示为:

(3.1)

?X1 X?U??0?0??UX2??H其中X1?Cp?p,X2?C(n?p)?(n?p),r?rank(In?P),U为酉矩阵且由P唯一确定.

引理3.3.2[34] (广义奇异值分解(GSVD)) 给定A?Cm?n,C?Cm?l，则存在酉矩

阵U1?OCn?n,V1?OCl?l以及非奇异矩阵W1?Cm?m使得

A?W1?AU1,C?W1?CV1 (3.2) 这里?A?Cm?n,?C?Cm?l，并且k?r?AC?,r?r(A),t?r(A)?r(C)?r?AC?,

?IA?0??A??0??0r?t0SA00t

0?r?0?0A?k?0?mn?r?tt?r?k，?C?0C?0???0??00SC00t0?r?0?IC?k?0?mk?r?tt?r?k,

l?r?t?k其中,IA和IC为单位矩阵,0A和0C为零矩阵，并且 SA?diag(?1,?,?t),SC?diag(?1,?,?t),

其中1??1????t?0,0??1????t?1，而?i2??i2?1,i?1,?,t. 3.3.1问题Ⅰ的解

现在考虑问题Ⅰ，取U如（3.1）式，记

AU?(A1,A2),UH?B1?B???, (3.3)

?B2?其中A1?Cm?r,A2?Cm?(n?r),B1?Cr?m,B2?C(n?r)?m,r?rank(In?P)

而矩阵对(A1,A2),(B1H,B2H)的广义奇异值分解如下:

A1?M1?A1U1,A2?M1?A2V1,

(3.4)

B1?M2?HB1HU2,r?rB2?M2?,V1?OCHB2HV, (3.5)

,U2?OCr?r其中酉矩阵U1?OCM1?Cm?m(n?r)?(n?r),V2?OC(n?r)?(n?r),非奇异矩阵

,M2?Cm?rm?m,

?Cm?(n?r)这里?A1?C,?A2,并且k1?rank(A1,A2),r1?rank(A1),t1?rank(A1)?

m?rrank(A2)?rank(A1,A2);?HHB1H?C,?B2H?CHm?(n?r),并且k2?rank(BH,BH),r2?r

12ank(B1),t2?rank(B1)?rank(B2)?rank(B1,B2),

?IA1?0???0??0r1?t1HH0SA00t11?A10?r1?t1?t10?0A?k1?r11?0?m?k1r?r1， ?A2?0A2?0???0??00SA00t120?r1?t1?t10?IA?k1?r12?0?m?k1k1?r1,

n?r?r1?t1?k1 ?BH1?IBH1?0???0??00SBH100t20?r2?t2?t20?0BH?k2?r21?0?m?k2r?r2, ?BH2?0BH2?0???0??00SBH200t20?r2?t2?t20?IBH?k2?r22?0?m?k2k2?r2,

r2?t2n?r?r2?t2?k2于是有以下定理:

定理3.3.1给定A?Cm?n,B?Cn?m,C?Cm?m，非负整数s和广义反射矩阵P?Cn?n

AU,UHB按(3.3)式进行分块,矩阵对(A1,A2),(B1H,B2H)的广义奇异值分解由

(3.4)(3.5)式给出,则矩阵方程AXB?C有自反矩阵解的充要条件是:

C13?0,C31?0,C4i?0(i?1,2,3,4),Cj4?0(j?1,2,3), (3.6)

并且若(3.6)式成立,则矩阵方程AXB?C有自反矩阵解的一般表达式为:

?U1HZUX?U?0??C11??1C21其中Z??SA1?Z31?C12SBH12??UHV1YV2?0H, (3.7)

?1SA1(C22?SA2Y22SBH)SBH21?1?1Z32r2?t2t2r?t2Z13?r1?t1?Z23?t1Z33?r?r1?,

?Y11?Y??Y21?Y?31Y12Y22C32SBH2?1?n?r?r?t?k111??1SA2C23?t1C33?k1?r1?Y13t2k2?r2,Zi3,Y1i(i?1,2,3),Z31,Z32,Y21,Y22,Y31,是

n?r?r2?t2?k2具有相应维数的任意矩阵,U1,U2,V1,V2,SA,SA,SB,SB见(3.4)(3.5)式,Cij(i?1

12H1H2,2,j?1,2),C23,C32,C33见①式.

n?n证明 给定A?Cm?n,B?Cn?m,C?Cm?m,由引理3.3.1，对X?CR(P),有

X?U???X1?00??UX2??H,

其中X1?Cp?p,X2?C(n?p)?(n?p),由Ai,Bi(i?1,2)的定义,矩阵方程AXB?C有自反矩阵解等价于矩阵方程：

A1X1B1?A2X2B2?C, 有一般解.

下面来解矩阵方程A1X1B1?A2X2B2?C,

矩阵对(A1,A2),(B1H,B2H)分解如(3.4),(3.5)式 令F(X1,X2)?A1X1B1?A2X2B2?C,有

F(X1,X2)?M1??M1?A1U1X1(M2?U1X1U2HHB1U2)HH?M1?A2A2V1X2(M2?HHB2V2)HH?C

A1?HHB1M2?M1?V1X2V2?HBH2M2?C

?M1(?A1U1X1U2H?HHB1??A2V1X2V2H?HHB2?M1CM?1?H2)MH2

令Z?U1X1U2H,Y?V1X2V2H,则有 rank(F(X1,X2))?rank(M1(?

将Z,Y,M1?1CM?Z11?Z?Z21???Z31r2?t2?12A1Z?HHB1HHB1??A2Y?HHB2HBH2?M1CM?1?H2?1?H2)MH2)

?rank(?A1Z???A2Y??M1CM)

作如下分块：

?Y11?Y?21??Y31Y12Y22Y32t2Z12Z22Z32t2Z13?r1?t1?Z23t?1Z33??r?r1r?t2， Y?Y13?n?r?r1?t1?k1?Y23t1?Y33?k1?r1?k2?r2，

n?r?r2?t2?k2M1CM?1?12?C11?C21???C31??C41r2?t2C12C22C32C42t2C13C23C33C43k2?r2C14?r1?t1?C24t1?, ① C34?k1?r1?C44?m?k1m?k2则有

?A1Z?HHB1??A2Y?HBH2?M1CM?1?H2

?C13?C14?r1?t1??C24t1?， ?C34?k1?r1??C44?m?k1m?k2?Z11?C11?SA1Z21?C21????C31??C41?r2?t2Z12SBH?C121SA1Z22SBH?SA2Y22SBH?C2212SA2Y23?C23Y33?C33?C43k2?r2Y32SBH?C322?C42t2在上式中，由于Zij(i?1,2,j?1,2),Ykl(k?2,3,l?2,3)的任意性，可知

minrank(F(X1,X2))?minrank(?A1Z?HHB1??A2Y?HBH2?M1CM?1?H2)

则矩阵方程F(X1,X2)?0有解的充要条件是：

C13?0,C31?0,C4i?0(i?1,2,3,4),Cj4?0(j?1,2,3).

易知,若阵方程F(X1,X2)?0有解,解的一般表达式为:

?U1HZUX?U?0?2??UHV1YV2?0H,

?C11??1其中Z??SAC211?Z31?C12SBH1?1SA1(C22?SA2Y22SBH)SBH21?1?1Z32r2?t2t2r?t2Z13?r1?t1?Z23?t1Z33?r?r1?,

?Y11?Y??Y21?Y?31Y12Y22C32SBH2?1?n?r?r?t?k111??1SA2C23?t1C33?k1?r1?Y13t2k2?r2,Zi3,Y1i(i?1,2,3),Z31,Z32,Y21,Y22,Y31,是

n?r?r2?t2?k2具有相应维数的任意矩阵,U1,U2,V1,V2,SA,SA,SB,SB见(3.4)(3.5)式,Cij(i?1

12H1H2,2,j?1,2),C23,C32,C33见①式.

3.3.2问题Ⅱ的解

由问题Ⅰ的解,易知矩阵方程AXB?C的自反矩阵解集合是一个闭凸集,因此问题Ⅱ一定存在唯一的最佳逼近.

n?n?n?nn?nn?n因为CR的一个子空间，令(CR(P))表示CR(P)的正交补空间,(P)是C则对任意的X??Cn?n,有

X??X1??X2?

?n?n?n?n?n?n其中X1??CR(P),X2?(CR(P)).因为X1?CR(P),由引理3.3.1知:

X?1??X11?U??0?0??U??X22?H

利用矩阵的广义奇异值分解表达式中的酉矩阵U1,U2,V1,V2和引理(3.1)式中

?H?H的U,对给定的U1X11U2,V1X22V2进行如下分块:

U1X11U2?H???W11???W21???W31?W12?????W13?W22W32W23?,V1X22V2??W33????H???W44???W54???W64?W45W55W65?????W46?W56?, ??W66???(3.8)

于是有如下定理:

定理2.3.2 设矩阵A?Cm?n,B?Cn?m,C?Cm?m分解式为(3.4)(3.5)式,且(3.6)式成立.给定X?C示为:

?n?nn?n(P),且它可表,则问题Ⅱ存在唯一的最佳逼近解X?CR~?U1HZU~X?U?0??C11??1其中,Z??SAC211?W?31?C12SBH12??UHV1YV2?0?H, (3.8)

?1S?1A1(C22?SA2W55SBH)SBH21??1W32r2?t2t2r?t2?W13?r1?t1??W23?t1?r?r1W33??

?W?44??Y??W54?W??64W45W55C32SBH2???1?W46?n?r?r1?t1?k1??1SA2C23?t1C33?k1?r1?t2k2?r2,Zi3,Y1i(i?1,2,3),Z31,Z32,Y21,Y22,Y31,是

n?r?r2?t2?k2具有相应维数的任意矩阵,U1,U2,V1,V2,SA,SA,SB,SB见(3.4)(3.5)式,Cij(i?1

12H1H2,2,j?1,2),C23,C32,C33见①式.

证明 对于给定的X??Cn?n,有如下分解:

X??X1?X2?? ,

?n?n?n?n其中X1??CR(P),X2?(CR(P)).

利用酉矩阵对Frobenius范数的不变性和范数的基本性质，对于X?S,有

~m2222X?X??X?X1?X22???X?X1???X2?

因此,minX?X?~X?Sm?minX?X1~X?Sm,而

2?X?X1?2?UZU?U?10??U1ZUH2H2??UHV1YV2?0?2HH?X1?UZU??0?H12??X???HV1YV2??00?110???X22?2

?X11?H2?V1YV2?X22?H2?2

?Z?U1X11U2?112?11?Y?V1X22V2?122

?SC21?W?2?1A1?212?C11?W?11?C12SBH?W?11?Z13?W?222?132?Z23?W2?232

2?SA(C22?SA2Y22SBH)SBH?W2?Z31?W31?2?Z32?W32?2??Z33?W33?2?1?

?2?Y11?W44?2?Y12?W45?12?2?Y13?W46?652?Y21?W54?662?Y22?W55?SA2C23?W56

?Y31?W?642?C32SBH?W?C33?W

X?SXminX?X??minX?X1~X?Sm? ?

??1?11?Z13?W13,Z23?W23,SA1(C22?SA2Y22SBH)SBH?W22,Z31?W312???

Z32?W32,Z33?W33,Y11?W44,Y12?W45,Y13?W46,Y21?W54,Y22?W55

???????将上式代入(3.7)式,得(3.8)式.证毕.

参考文献：

[1] D.L.Chu,B.D.Moor.On a variational formulation of QSVD and the RSVD

L.A.A.,2000(311)：61-78

[2] 刘瑞娟.几类约束矩阵方程的定秩解的问题.长沙理工大学,2008,1-43 [3] 肖庆峰.秩约束下几类矩阵方程问题及其最佳逼近问题.湖南大学数学与计

量经济学院,2009,1-106

[5] Penrose R.A geralized inverse for matrices.Proc Cambridge Philos

Soc,1955,5:406-413

[6] Lancaster P.Explicit solutions of Linear matrix equations.SIAM

Review,1970,72(4):544-566

[7] Khatai C G,Mitra S K.Hermitian and nongative definite solutions of

linear matrix equations. SIAM Journal on Applied Mathematics,1976, 31(4):576-585

[8] Dai H.On the symmetric solutions of linear matrix equations.Linear

Algebre and Its Application.1990,131:1-7

[9] 邓远北.几类线形矩阵方程的解与PROCRUSTES问题:[湖南大学博士论文].长

沙:湖南大学数学与计量经济学院,2003,1-99

[10]Peng X Y,Hu X Y,Zhang L.The reflexive and anti-reflexive solutions of the matrix equation AHXB?C. Journal of Computational and Applied Mathematics,2007,200:749-760

[11]彭向阳.几类特殊约束矩阵方程问题及其最佳逼近问题:[湖南大学博士文].

长沙:湖南大学数学与计量经济学院,2006,1-138

[12]雷渊.求解一类矩阵最佳逼近问题的理论和算法:[湖南大学博士论文].长沙:

湖南大学数学与计量经济学院,2007,1-134

[13]Liao A P,Bai Z Z.Least squares solutions of AXB?Dover symmetric

positive semidefinte mareices. Journal of ComputationalMathematics ,2003,21(2):175-182

[14]彭亚新.求解约束矩阵方程及其最佳逼近的迭代法研究:[湖南大学博士文].

长沙:湖南大学数学与计量经济学院,2004,1-142

[15]Mitra S K.A pair of simultaneous linear matrix equations A1XB1?C1,

A2XB2?C2 and a matrix programming problem. Linear Algebre and Its Application.1990,131:107-123

[16]Mitra S K.Fixed rank solutions of linear matrix equations.Sankhya

Ser.A.,1972, 35: 387-392

[17]Mitra S K.The matrix equation AX?C,XB?D.Linear Algebre and Its

Application,1984,59:171-181 [18]Uhlig F .On the matrix equation AX?B with applications to the

generators of controllability matrix.Linear Algebre and Its Application,1987,85:203-209

[19]Gross J.Nonnegtive-definite and positive-definite solution to the

matrix equation AXA??B revisited. Linear Algebre and Its Applica- tion.2000,321:123-129

[33]Peng Z Y,Hu X Y.The reflexive and anti-reflexive solutions of the

matrix equation AX?B. Linear Algebre and Its Application,2003, 375:147-155

[34]Paige C C,Saunders M A.Towards a generalized singular value decom- Position.SIAM Journal on Numerical Anaysis,1981,18:398-405

[35]Qing-Feng Xiao,Xi-Yan Hu,Lei Zhang.The Symmetric Minimal Rank

Solution of the Matrix Equation AX?B and the Optimal Approxim- ation.Electronic Journal of LINEAR Algebra.2009,(18):264-273(SCI)

Application,1984,59:171-181 [18]Uhlig F .On the matrix equation AX?B with applications to the

generators of controllability matrix.Linear Algebre and Its Application,1987,85:203-209

[19]Gross J.Nonnegtive-definite and positive-definite solution to the

matrix equation AXA??B revisited. Linear Algebre and Its Applica- tion.2000,321:123-129

[33]Peng Z Y,Hu X Y.The reflexive and anti-reflexive solutions of the

matrix equation AX?B. Linear Algebre and Its Application,2003, 375:147-155

[34]Paige C C,Saunders M A.Towards a generalized singular value decom- Position.SIAM Journal on Numerical Anaysis,1981,18:398-405

[35]Qing-Feng Xiao,Xi-Yan Hu,Lei Zhang.The Symmetric Minimal Rank

Solution of the Matrix Equation AX?B and the Optimal Approxim- ation.Electronic Journal of LINEAR Algebra.2009,(18):264-273(SCI)

本文来源：https://www.bwwdw.com/article/6oqp.html