矩阵方程的解法

两类矩阵方程的行对称矩阵解

及AX=B的最佳逼近

摘要本文首先介绍了行对称矩阵的定义及性质，利用矩阵的广义逆，奇异值分解，给出了矩阵方程AX=B有行对称解的充分必要条件及有解时通解的表达式；并给出了矩阵方程解集合中与给定矩阵的最佳逼近解的表达式。最后利用奇异值分解给出了矩阵方程AXAT对称解的充分必要条件及有解时通解的表达式。

矩阵方程问题是指在满足一定条件的矩阵集合中求矩阵方程的解的问题。不同的约束条件，不同的矩阵方程，就导致了不同的约束矩阵方程问题。约束矩阵方程问题在结构设计，参数识别，主成分分析，勘测，遥感，生物学，电学，固体力学，结构动力学，分子光谱学，自动控制理论，振动理论，循环理论等领域都有重要应用。

约束矩阵方程问题的内容非常广泛. 约束矩阵方程问题又分为线性约束矩阵方程问题和非线性约束矩阵方程问题. 有关线性约束矩阵方程问题的研究成果相当丰富. 其中最简单的矩阵方程AX = B是研究最透彻的一类问题.

求解线性矩阵方程一般会遇到两种情况：一是当矩阵方程有解时，如何求它的解及最佳逼近；二是当矩阵方程无解时，如何求它的最小

B有行

二乘解。对于本文所研究的AX=B、AXAT B这两类简单矩阵方程，国

内外学者已经作了大量研究。都在相应的文献中对其进行了大量的研究，解决了求此方程的一些约束解和最小二乘解的问题。

自从针对工程应用领域提出了行对称矩阵概念之后，这方面研究已经取得了一些成果，如对行对称矩阵的一些性质，行对称矩阵的QR分解。

本文先对行对称矩阵进行介绍，再将行对称矩阵与约束矩阵方程结合起来，先研究了矩阵方程AX=B有行对称实矩阵解的充要条件，有解时，用奇异值分解及广义逆求出解及最佳逼近。再对矩阵方程

AXA B有行对称实矩阵解的充要条件进行了研究，利用奇异值分解

T

得出了有解时的充要条件及解的表达式。

设Rm*n表示全体n*m阶实矩阵集合，rank(A)表示矩阵A的秩,Jn表示

0

次对角线上元素全为1，其余元素全为0的方阵，即Jn=

1

1

,0 n*n

显然有J

1

n

Jn,Jn Jn成立。OR

T

n*n

表示n阶正交矩阵全体。

本文要讨论以下问题：

问题1 给定矩阵A,B m*n,求实行对称方阵X，使得AX=B。

问题 2给定X

1

R

n*n

，求X

SE

，使得 X X1 min X

X SE

X1 。其中SE

为问题1的解集。

问题 3给定矩阵

定义1[2]设A = (a) Rn*m，若A满足a

ij

ij

A,B R

mn

,求实行对称方阵X，使得AXAT=B。

an i 1,j,i 1,2, ,n;j 1,2, ,m

，

则称A为n *m行对称矩阵. 所有n *m行对称矩阵的全体记为

RSR

n*m

。

考查满足aij

an i 1,j的矩阵A，不难发现A是关于行具有某种

n 12

对称性的矩阵，即当阶数n为奇数时，以将

n2

行为对称线，矩阵A的

2

行关于该线对称；当阶数n为偶数时，在行与则A的行关于该直线对称。

n 2

行间做一条直线，

或简单的说，将A进行上下翻转后矩阵不变，我们就称这种矩阵为行对称矩阵。

为了更好的了解行对称矩阵，我们介绍一下行对称矩阵的性质：

（1）当n=2k时，RSR

n*m

=

A k*m

{A |A1 R}.

JkA1

A1

k*m1*m |A1 R, R} （2）当n=2k+1时，RSRn*m={A JA k1

定义2[1]设A=(a

ij

)n*m

，r(A)=r,ATA的大于零的特征值为 1, 2,

, r。则A的奇异值。

定义3[1]设矩阵A Rn*m，若矩阵X Rm*n满足如下四个Penrose方程： AXA=A XAX=X

(AX)(XA)

T

=AX =XA

则称X为A的Penrose广义逆，记为A 。 设矩阵A Rn*m，若矩阵X Rm*n满足： AX=P

R(A)

T

， ，

XA=P

R(X)

其中PL是子空间L上的正交投影矩阵，则称X为A的Moore广义逆矩阵。

Moore广义逆矩阵与Penrose广义逆矩阵是等价的。因此A 通常称为Moore-Penrose广义逆。

显然，当A为非奇异矩阵时，有A =A 1。 定义4

[1]

设A=(a

ij

)n*n R

n*n

，

令 A (

i,j 1

n

aij)

2/21

，则 称为R

n*n

上的Frobenius范数。

引理1[2]A RSRn*m，当且仅当A=JnA。

Jn的第i行为(0, 0,1,0 ,0)

n i 1

JnA

的第ｉ行j列位置的元素为a

n i 1,j

aij an i 1,j

A RSRn*m

设A=(a

ij

)m*n，

A的奇异值分解，则A有如下分解：

0 m*n

A=UDVT

，

D=

，

其中U，V分别为m阶和n阶的正交矩阵。 上式称为A的奇异值分解。。

对任意A m*n，A 存在并且唯一。

给定矩阵A,B m*n，若矩阵A的奇异值分解为

A=U

0

0 V0

T

其中 =diag(a1,a2, ,ar)，ai>0,(i=1,2,…,r),

r =rank(A),U=(U1,U2),V=(V1,V2),U为m阶正交矩阵，V为n阶正交矩阵，U1 m*r，V1

n*r

,

则矩阵方程AX=B有解的充分必要条件是U2TB=0,且有解时的一般表

X=A其中G0

在Rn*n上，矩阵乘上一个正交矩阵后，它的Frobenius范数不变。

2.问题1的解

先对后面证明要用到的两个矩阵做奇异值分解： 矩阵A的奇异值分解为 A=U

0

0 V0

T

B V2G0

(n r)*n

是任意的。

（1）

其中

=diag(a1,a2, ,ar)，ai>0,(i=1,2,…,r),

r =rank(A),U=(U1,U2),V=(V1,V2),U为m阶正交矩阵，V为n阶正交矩阵，U1 m*r，V1

n*r

。

矩阵W=V2 JnV2 Rn*(n r)的奇异值分解为

W=P

0

0 T

Q0

（2）

其中 =diag(b1b2 bt),bi>0,(i=1,2, ,t),t=rank(W),P=(P1,P2)

OR

n*n

,

n*t

1

1

Q=(Q1,Q2) OR(n r)*(n r),P R,Q

R

(n r)*t

.

给定矩阵A,B m*n,求实行对称方阵X，使得AX=B。

将A，W分别按（1），（2）进行分解，则问题1有解的充分必要条件是

U

2

T

B=0，PN=P（J

T2

T2

n

AB AB

）=0 （3）

且有解时的一般表达式为 X=A

B V2（WN Q2G

）（4）

其中G R(n r)*r是任意的。

证明：由引理4，

AX=B有解的充要条件为U2TB=0， 它的通解为X=A

B V2G0

。

由引理1，X为行对称矩阵的充要条件为X=JnX 即（A

B V2G0）=Jn（AB V2G0

）

（V2 JnV2）G0=J令 W=V2 JnV2 , N=J

n

n

AB AB

AB AB，

由引理4，则问题1有解的充要条件是…… 此处隐藏：2068字，全部文档内容请下载后查看。喜欢就下载吧 ……