矩阵方程的解法.docxVIP

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