矩阵方程的解法.docx
- 文档编号:18270826
- 上传时间:2023-08-14
- 格式:DOCX
- 页数:7
- 大小:37.44KB
矩阵方程的解法.docx
《矩阵方程的解法.docx》由会员分享,可在线阅读,更多相关《矩阵方程的解法.docx(7页珍藏版)》请在冰点文库上搜索。
矩阵方程的解法
两类矩阵方程的行对称矩阵解
及AX=B的最佳逼近
摘要本文首先介绍了行对称矩阵的定义及性质,利用矩阵的广义逆,奇异值分解,给出了矩阵方程AX=B有行对称解的充分必要条件及有解时通解的表达式;并给出了矩阵方程解集合中与给定矩阵的最佳逼近解的表达式。
最后利用奇异值分解给出了矩阵方程AXATB有行对称解的充分必要条件及有解时通解的表达式。
矩阵方程问题是指在满足一定条件的矩阵集合中求矩阵方程的解的问题。
不同的约束条件,不同的矩阵方程,就导致了不同的约束矩阵方程问题。
约束矩阵方程问题在结构设计,参数识别,主成分分析,勘测,遥感,生物学,电学,固体力学,结构动力学,分子光谱学,自动控制理论,振动理论,循环理论等领域都有重要应用。
约束矩阵方程问题的内容非常广泛.约束矩阵方程问题又分为线性约束矩阵方程问题和非线性约束矩阵方程问题.有关线性约束矩阵方程问题的研究成果相当丰富.其中最简单的矩阵方程AX=B是研究最透彻的一类问题.
求解线性矩阵方程一般会遇到两种情况:
一是当矩阵方程有解时,如何求它的解及最佳逼近;二是当矩阵方程无解时,如何求它的最小
二乘解。
对于本文所研究的AX=B、AXATB这两类简单矩阵方程,国内外学者已经作了大量研究。
都在相应的文献中对其进行了大量的研究,解决了求此方程的一些约束解和最小二乘解的问题。
自从针对工程应用领域提出了行对称矩阵概念之后,这方面研究已经取得了一些成果,如对行对称矩阵的一些性质,行对称矩阵的QR分解。
本文先对行对称矩阵进行介绍,再将行对称矩阵与约束矩阵方程结合起来,先研究了矩阵方程AX=B有行对称实矩阵解的充要条件,有解时,用奇异值分解及广义逆求出解及最佳逼近。
再对矩阵方程AXATB有行对称实矩阵解的充要条件进行了研究,利用奇异值分解得出了有解时的充要条件及解的表达式。
设Rm*n表示全体n*m阶实矩阵集合,rank(A)表示矩阵A的秩,Jn表示
1
0n*n
n*n
0次对角线上元素全为1,其余元素全为0的方阵,即Jn=
1
显然有Jn1Jn,JnTJn成立。
ORn*n表示n阶正交矩阵全体。
本文要讨论以下问题:
问题1给定矩阵A,Bm*n,求实行对称方阵X,使得AX=B。
问题2给定X1Rn*n,求XSE,使得XX1minXX1。
其中SEXSE
为问题1的解集
问题3给定矩阵A,BRmn,求实行对称方阵X,使得AXAT=B
定义1[2]设A=(aij)Rn*m,若A满足aijani1,j,i1,2,,n;j1,2,,m,
则称A为n*m行对称矩阵.所有n*m行对称矩阵的全体记为
考查满足aijani1,j的矩阵A,不难发现A是关于行具有某种对称性的矩阵,即当阶数n为奇数时,以将n1行为对称线,矩阵A的2行关于该线对称;当阶数n为偶数时,在n行与n2行间做一条直线,
22则A的行关于该直线对称。
或简单的说,将A进行上下翻转后矩阵不变,我们就称这种矩阵为行对称矩阵。
为了更好的了解行对称矩阵,我们介绍一下行对称矩阵的性质:
(1)当n=2k时,RSRn*m={A|A1Rk*m}.
JkA1
A1
2)当n=2k+1时,RSRn*m={A|A1Rk*m,R1*m}
定义2[1]设A=(aij)n*m,r(A)=r,ATA的大于零的特征值为1,2,,r。
则
1,2,,r称为A的奇异值
定义3[1]设矩阵ARn*m,若矩阵XRm*n满足如下四个Penrose方程:
AXA=A
XAX=X
(AX)T=AX
(XA)T=XA
则称X为A的Penrose广义逆,记为A设矩阵ARn*m,若矩阵XRm*n满足:
XA=PR(X),
其中PL是子空间L上的正交投影矩阵,则称X为A的Moore广义逆矩阵。
Moore广义逆矩阵与Penrose广义逆矩阵是等价的。
因此A通常称为Moore-Penrose广义逆。
显然,当A为非奇异矩阵时,有A=A1。
n
称为R
定义4[1]设A=(aij)n*nRn*n,令A(a2i/j21)tr(ATA),则
i,j1
上的Frobenius范数。
证明:
ni1
Jn的第i行为(0,0,1,0,0)
JnA的第i行j列位置的元素为ani1,j
aijani1,j
设A=(aij)m*n,且1,2,,r为A的奇异值分解,则A有如下分解:
1
A=UDVT,D=r,
0
0m*n
其中U,V分别为m阶和n阶的正交矩阵。
上式称为A的奇异值分解。
。
对任意Am*n,A存在并且唯
给定矩阵A,Bm*n,若矩阵A的奇异值分解为
其中=diag(a1,a2,,ar),ai>0,(i=1,2,⋯,r),
r=rank(A),U=(U1,U2),V=(V1,V2),U为m阶正交矩阵,V为n阶正交矩
则矩阵方程AX=B有解的充分必要条件是U2TB=0,且有解时的一般表
达式为
X=ABV2G0
其中G0(nr)*n是任意的
在Rn*n上,矩阵乘上一个正交矩阵后,它的Frobenius范数不变
2.问题1的解先对后面证明要用到的两个矩阵做奇异值分解:
矩阵A的奇异值分解为
0T
A=U0VT
(1)
00
其中
=diag(a1,a2,,ar),ai>0,(i=1,2,⋯,r),
r=rank(A),U=(U1,U2),V=(V1,V2),U为m阶正交矩阵,V为n阶正交矩矩阵W=V2JnV2Rn*(nr)的奇异值分解为
W=P000QT
(2)
其中=diag(b1b2bt),bi>0,(i=1,2,,t),t=rank(W),P=(P1,P2)
Q=(Q1,Q2)OR(nr)*(nr),P1Rn*t,Q1R(nr)*t.给定矩阵A,Bm*n,求实行对称方阵X,使得AX=B。
将A,W分别按
(1),
(2)进行分解,则问题1有解的充分必要条件
且有解时的一般表达式为
X=ABV2(WNQ2G)(4)
其中GR(nr)*r是任意的。
证明:
由引理4,
AX=B有解的充要条件为U2TB=0,它的通解为X=ABV2G0。
由引理1,X为行对称矩阵的充要条件为X=JnX
即(ABV2G0)=Jn(ABV2G0)
(V2JnV2)G0=JnABAB
令W=V2JnV2,
N=JnABAB,
由引理4,则问题1有解的充要条件是
U2TB=0,
P2TN=P2T(JnABAB)=0当问题1有解时,可以解得
G0=WNQ2G,GR(nr)*r是任意的。
所以有解时,方程的行对称矩阵解为X=ABV2(WNQ2G),GR(n*rr)是任意的。
问题2给定X1Rn*n,求XSE,使得XX1minXX1。
其中SEXSE
为问题1的解集。
给定X1Rn*n,若问题1的解集合SE非空,则问题2
在SE中存在唯一解X,并且
X=ABV2[Q2T(MWN)](3)
其中M=V2T(X1AB)。
证明:
因为问题1的解集合SE非空,则SE是Hilbert空间Rn*n中一个非空闭凸锥。
所以问题2有唯一解XSE[4,5,6]。
把(4)代入,有
ABV2(WNQ2G)X1
V1T
T
V2T
T2
=VT(ABV2(WNQ2G)X1)2
2
V2(WNQ2G)(X1AB)2
2
(ABV2(WNQ2G)X1)2
WNQ2GV2T(X1AB)2+V1T(X1AB)2
2T2
Q2G(MWN)2+V1T(X1AB)2
左乘QT,上式
T2T2T2
=GQ2T(MWN)2+Q1T(MWN)2+V1T(X1AB)2
因此,要使得XX12=min
则G=Q2T(MWN),其中M=V2T(X1AB),
所以X=ABV2[Q2T(MWN)],其中M=V2T(X1AB)。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 矩阵 方程 解法