矩阵方程的解法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
矩阵方程的解法
本文首先介绍了行对称矩阵的定义及性质,利用矩阵的广义逆,奇异值分解,给出了矩阵方程AX=B有行对称解的充分必要条件及有解时通解的表达式;并给出了矩阵方程解集合中与给定矩
阵的最佳逼近解的表达式。最后利用奇异值分解给出了矩阵方程
有行对称解的充分必要条件及有解时通解的表达式。矩阵方程问
题是指在满足一定条件的矩阵集合中求矩阵方程的解的问题。不
同的约束条件,不同的矩阵方程,就导致了不同的约束矩阵方程
问题。约束矩阵方程问题在结构设计,参数识别,主成分分析,
勘测,遥感,生物学,电学,固体力学,结构动力学,分子光谱学,自动控制理论,振动理论,循环理论等领域都有重要应用。
约束矩阵方程问题的内容非常广泛、约束矩阵方程问题又分为线性约束矩阵方程问题和非线性约束矩阵方程问题、有关线性约束矩阵方程问题的研究成果相当丰富、其中最简单的矩阵方程AX = B是研究最透彻的一类问题、求解线性矩阵方程一般会遇到两种情况:一是当矩阵方程有解时,如何求它的解及最佳逼近;二是
当矩阵方程无解时,如何求它的最小二乘解。对于本文所研究的AX=
B、这两类简单矩阵方程,国内外学者已经作了大量研究。都在相应的文献中对其进行了大量的研究,解决了求此方程的一些
约束解和最小二乘解的问题。自从针对工程应用领域提出了行对
称矩阵概念之后,这方面研究已经取得了一些成果,如对行对称矩阵的一些性质,行对称矩阵的QR分解。本文先对行对称矩阵进行介绍,再将行对称矩阵与约束矩阵方程结合起来,先研究了矩阵方程AX=B有行对称实矩阵解的充要条件,有解时,用奇异值分解及广义逆求出解及最佳逼近。再对矩阵方程有行对称实矩阵解的充要条件进行了研究,利用奇异值分解得出了有解时的充要条件及解的表达式。设表示全体n*m阶实矩阵集合,rank(A)表示矩阵A的秩,表示次对角线上元素全为1,其余元素全为0的方阵,即=,显然有成立。表示n阶正交矩阵全体。本文要讨论以下问题:问题1 给定矩阵A,B,求实行对称方阵X,使得AX=B。问题2 给定,求,使得。其中为问题1的解集。问题3 给定矩阵,求实行对称方阵X,使得=B。
定义设A = (),若A满足,则称A为n *m行对称矩阵、所有n *m行对称矩阵的全体记为。考查满足的矩阵A,不难发现A 是关于行具有某种对称性的矩阵,即当阶数n为奇数时,以将行为对称线,矩阵A的行关于该线对称;当阶数n为偶数时,在行与行间做一条直线,则A的行关于该直线对称。或简单的说,将A 进行上下翻转后矩阵不变,我们就称这种矩阵为行对称矩阵。为了更好的了解行对称矩阵,我们介绍一下行对称矩阵的性质:(1)当n=2k时,=、(2)当n=2k+1时,=定义设A=,r(A)=r,的大于零的特征值为。则称为A的奇异值。定义设矩阵A ,若矩阵X满足如下四个Penrose方程:
AXA=A XAX=X =AX =XA则称X为A的Penrose广义逆,记为。设矩阵A ,若矩阵X满足:AX= ,XA=,其中是子空间L上的正交投影矩阵,则称X为A的Moore广义逆矩阵。Moore广义逆矩阵与Penrose广义逆矩阵是等价的。因此通常称为Moore-Penrose广义逆。显然,当A为非奇异矩阵时,有=。定义设A=,令,则称为上的Frobenius范数。引理 A,当且仅当A=。证明:
的第i行为的第i行j列位置的元素为 A设A=,且为A的奇异值分解,则A有如下分解:
A=UD,D=,其中U,V分别为m阶和n阶的正交矩阵。上式称为A的奇异值分解。。对任意A,存在并且唯一。给定矩阵A,B,若矩阵A的奇异值分解为A=U其中=diag,>0,(i=1,2,…,r),r
=rank(A),U=,V=(),U为m阶正交矩阵,V为n阶正交矩阵,, ,则矩阵方程AX=B有解的充分必要条件是B=0,且有解时的一般表达式为 X= 其中是任意的。
在上,矩阵乘上一个正交矩阵后,它的Frobenius范数不变。
2、问题1的解先对后面证明要用到的两个矩阵做奇异值分解:矩阵A的奇异值分解为A=U (1)其中
=diag,>0,(i=1,2,…,r),r =rank(A),U=,V=(),U为m阶正交矩阵,V为n阶正交矩阵,,。矩阵W=的奇异值分解为W=P (2)其中=diag(),>0,(i=1,2,,t),t=rank(W),P=(),Q=(),,、给定矩阵A,B,求实行对称方阵X,使得AX=B。
将A,W分别按(1),(2)进行分解,则问题1有解的充分必要条件是B=0,N=()=0 (3)且有解时的一般表达式为 X=()(4)其中是任意的。证明:由引理4,AX=B有解的充要条件为B=0,它的通解为X=。由引理1,X为行对称矩阵的充要条件为X=X即()=()()=令 W= ,N= ,由引理4,则问题1有解的充要条件是B=0,N=()=0 当问题1有解时,可以解得=,是任意的。所以有解时,方程的行对称矩阵解为X=(),是任意的。问题2 给定,求,使得。其中为问题1的解集。给定,若问题1的解集合非空,则问题2在中存在唯一解,并且 = (3)其中
M=。证明:
因为问题1的解集合非空,则是Hilbert空间中一个非空闭凸锥。所以问题2有唯一解。
把(4)代入,有左乘,利用引理5,上式= = = 因为=0,所以上式 =+ =+ 左乘,上式 =++因此,要使得=min则 G=,其中
M=,所以=,其中M=。