资源描述：
《矩阵方程的解法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、两类矩阵方程的行对称矩阵解及AX二B的最佳逼近摘要本文首先介绍了行对称矩阵的定义及性质，利用矩阵的广义逆，奇异值分解，给出了矩阵方程AX二B有行对称解的充分必要条件及有解时通解的表达式；并给出了矩阵方程解集合中与给定矩阵的最佳逼近解的表达式。最后利用奇异值分解给出了矩阵方程AXAt=B有行对称解的充分必要条件及有解时通解的表达式。矩阵方程问题是指在满足一定条件的矩阵集合中求矩阵方程的解的问题。不同的约束条件，不同的矩阵方程，就导致了不同的约束矩阵方程问题。约束矩阵方程问题在结构设计，参数识别，主成分分析，勘测，遥感，生物学，电学，固体力学，结构动力学，分子光谱学，自动控制理论
2、，振动理论，循环理论等领域都有重要应用。约束矩阵方程问题的内容非常广泛.约束矩阵方程问题又分为线性约束矩阵方程问题和非线性约束矩阵方程问题.有关线性约束矩阵方程问题的研究成果相当丰富.其中最简单的矩阵方程力B是研究最透彻的一类问题.求解线性矩阵方程一般会遇到两种情况：一是当矩阵方程有解吋,如何求它的解及最佳逼近；二是当矩阵方程无解时，如何求它的最小二乘解。对于本文所研究的AX二B、AX"二B这两类简单矩阵方程，国内外学者已经作了大量研究。都在相应的文献中对其进行了大量的研究，解决了求此方程的一些约束解和最小二乘解的问题。自从针对工程应用领域提出了行对称矩阵概念之后，这方面研究
3、已经取得了一些成果，如对行对称矩阵的一些性质，行对称矩阵的QR分解。本文先对行对称矩阵进行介绍，再将行对称矩阵与约束矩阵方程结合起来，先研究了矩阵方程AX二B有行对称实矩阵解的充要条件,有解时，用奇异值分解及广义逆求出解及最佳逼近。再对矩阵方程axat=b有行对称实矩阵解的充要条件进行了研究，利用奇异值分解得出了有解时的充要条件及解的表达式。设肥初表示全体n*m阶实矩阵集合，rank(A)表示矩阵A的秩，人表示‘0…1、次对角线上元素全为1,其余元素全为0的方阵，即血二:；,显然有肓=人血成立。O严表示n阶正交矩阵全体。本文要讨论以下问题:问题1给定矩阵A,Bw口〃和，求实行
4、对称方阵X,使得AX二B。为问题1的解集。问题3给定矩阵A,Bw肿，求实行对称方阵X,使得AXAr=Bo定义ll2J设昇-(a.j)gRe,若力满足ai}=an_Mj,Z=1,2,…,zz;J=1,2,…，加,则称力为刀枷行对称矩阵.所有刀*〃7行对称矩阵的全体记为RSREo考查满足aijm’j的矩阵A,不难发现A是关于行具有某种对称性的矩阵，即当阶数n为奇数吋，以将凹行为对称线，矩阵A的2行关于该线对称；当阶数n为偶数时，在上行与出行间做一条直线，22则A的行关于该直线对称。或简单的说，将A进行上下翻转后矩阵不变，我们就称这种矩阵为行对称矩阵。为了更好的了解行对称矩阵，我们
5、介绍一下行对称矩阵的性质:(1)当n二2k时,RSRE二{A=(2)当n二2k+l时,RSRn*,H={A=*1aMA丿定义2⑴设A二⑺几恤，r(A)=r,A7A的大于零的特征值为…必。则瓦込，…皿称为A的奇异值。定义3山设矩阵AgR,rm,若矩阵XgR'^n满足如下四个Penrose方程：AXA二AXAX=X(AX)，二AX(XA)r=XA则称X为A的Penrose广义逆，记为X。设矩阵AgR^n,若矩阵Xw疋旳满足：从二厶⑷，XA二Pg，其中巴是子空间L上的止交投影矩阵，贝IJ称X为A的Moore广义逆矩阵。Moore广义逆矩阵与Penrose广义逆矩阵是等价的。因此X通
6、常称为Moore-Penrose广义逆。显然，当A为非奇异矩阵时,有£二屮。定义4山设A二(讥*亦严，令口A止送a泸二血(从4),则口•口称为心】尺泸“上的Frobenius范数。引理1⑵AeRSRn*m,当且仅当K=JnAo证明:血的第i行为(0,•••(),1,0.・・,0)«以的第i行j列位置的元素为勺*」•Oaij=%+ljo2RSRE设A二心馬，且扬，站,・・・,扬为A的奇异值分解，则A有如下分解：■■■A=UDVf,D=扬，0■其中U,V分别为III阶和n阶的正交矩阵。上式称为A的奇异值分解。。对任意Ag□"旳,A+存在并且唯一。给定矩阵A,BeD^,若矩阵A的奇
7、异值分解为A二U乙0I。0丿其中为二diag(di,。2,•••,/),at>0,(i=l,2,•••,r),r=rank(A),U=(t/i,t/2),V=(Vi.V2),U为m阶正交矩阵,V为n阶正交矩阵，t/ig□w*r,%g口"I则矩阵方程AX二B有解的充分必要条件是二0,且有解时的一般表达式为X二+其中Goe□(n-r)*rt是任意的。在尺沪"上，矩阵乘上一个正交矩阵后，它的Frobenius范数不变。2•问题1的解先对后面证明要用到的两个矩阵做奇异值分解：矩阵A的奇异值分解为A二UVT(1)
