广义逆矩阵及其应用

题目广义逆矩阵及其应用学院专业通信与信息系统学生学号目录第一章前言………………………………………………………………………1第二章广义逆矩阵………………………………………………………………2§2.1广义逆矩阵的定义………………………………………………………2§2.2广义逆矩阵的性质………………………………………………………3第三章广义逆矩阵的计算………………………………………………………12§3.1一般广义逆求解…………………………………………………………12§3.2Moore-Penrose广义逆…………………………………………………16结论………………………………………………………………………………19-1-第一章前言线性方程组的逆矩阵求解方法只适用于系数矩阵为可逆方阵,但是对于一般线性方程组，其系数矩阵可能不是方阵或是不可逆的方阵，这种利用逆矩阵求解线性方程组的方法将不适用。为解决这种系数矩阵不是可逆矩阵或不是方阵的线性方程组，我们对逆矩阵进行推广，研究广义逆矩阵，利用广义逆矩阵求解线性方程组。广义逆矩阵在数据分析、多元分析、信号处理、系统理论、现代控制理论、网络理论等许多领域中有着重要的应用，本文针对广义逆矩阵的定义、性质、计算及其在线性方程组中的应用进行研究，利用广义逆矩阵求解线性方程组的通解及极小范数解。逆矩阵的概念只对非奇异矩阵才有意义，但在实际问题中，遇到的矩阵不一定是方阵，即使是方阵也不一定非奇异，这就需要将逆矩阵的概念进行推广。为此，人们提出了下述关于逆矩阵的推广：（1）该矩阵对于奇异矩阵甚至长方矩阵都存在；（2）它具有通常逆矩阵的一些性质；（3）当矩阵非奇异时，它即为原来的逆矩阵。满足上面三点的矩阵称之为广义逆矩阵。1903年，瑞典数学家弗雷德霍姆开始了对广义逆矩阵的研究，他讨论了关于积分算子的一种广义逆。1904年，德国数学家希尔伯特在广义格林函数的讨论中，含蓄地提出了微分算子的广义逆。美国芝加哥的穆尔(Moore)教授在1920年提出了任意矩阵广义逆的定义，他以抽象的形式发表在美国数学会会刊上。我国数学家曾远荣和美籍匈牙利数学家冯·诺伊曼及其弟子默里分别在1933年和1936年对希尔伯特空间中线性算子的广义逆也作过讨论和研究。1951年瑞典人布耶尔哈梅尔重新给出了穆尔(Moore)广义逆矩阵的定义，并注意到广义逆矩阵与线性方程组的关系。1955年，英国数学物理学家彭罗斯(Penrose)以更明确的形式给出了与穆尔(Moore)等价的广义逆矩阵定义，因此通称为Moore-Penrose广义逆矩阵，从此广义逆矩阵的研究进入了一个新阶段。现如今，Moore-Penrose广义逆矩阵在数据分析、多元分析、信号处理、系统理论、现代控制理论、网络理论等许多领域中有着重要的应用，使这一学科得到迅速发展，并成为矩阵论的一个重要分支。-2-第二章广义逆矩阵§2.1广义逆矩阵的定义一、Penrose广义逆矩阵的定义为了推广逆矩阵的概念，我们引进了广义逆矩阵的定义，下面给出广义逆矩阵的Moore-Penrose定义。定义2.1设矩阵nmCA，若矩阵mnCX满足如下四个Penrose方程AAXA（ⅰ）XXAX（ⅱ）AXAXH)(（ⅲ）XAXAH)(（ⅳ）中的一部分或全部方程，则称X为A的一个广义逆矩阵。若X只满足（ⅰ）式，则X成为A的一个}1{-逆，可记为1A，所有满足}1{-逆的X构成的集合记为1A。若X满足四个方程中的第kji,,,个方程，则称X为A的一个kji,,,-逆，记为kjiA,,,，所有满足kji,,,-逆的X构成的集合记为kjiA,,,。二、常见广义逆定义按照广义逆定义，分别满足一个、两个、三个和四个方程的广义逆矩阵一共有44342414CCCC=15类，其中常见的有1A，2,1A，3,1A，4,1A，4,3,2,1A。定义2.2设有复矩阵nmCA。若有一个mn复矩阵X存在，使下式成立，则称X为A的减号逆：AAXA(2.1)当1A存在时，显然1A满足上式，可见减号逆X是普通逆矩阵1A的推广；另外，由AAXA得HHAAXA)(，即-3-HHHHAAXA可见，当X为A的一个减号逆时，HX就是HA的一个减号逆。定义2.3设复矩阵nmCA，若有一个mn矩阵X，满足：AAXA且XXAX称X为A的一个自反逆矩阵，记作为rA，rA满足Penrose方程的（ⅰ），（ⅱ）式，所以}2,1{AAr。显然，自反广义逆为减号逆的子集。对矩阵X是矩阵A的1-逆，即1AX,若矩阵A也是矩阵X的1-逆，即1XA，则X为A的一个自反逆矩阵。定义2.4设复矩阵nmCA，若有一个mn矩阵X，满足：AAXA及AXAXH)(，则称X为A的最小二乘广义逆，记作lA，lA满足Penrose方程的（ⅰ），（ⅲ）式，所以}3,1{AAm。最小二乘广义逆是用条件AXAXH)(对减号逆进行约束后所得到的子集。定义2.5设复矩阵nmCA，若有一个mn矩阵X，满足：AAXA及XAXAH)(，则称X为A的最小范数广义逆，记作mA，mA满足Penrose方程的（ⅰ），（ⅳ）式，所以}4,1{AAl。显然，最小范数广义逆也是减号逆的子集。若X满足全部四个方程，则称X为A的Moore-Penrose广义逆矩阵，记为A。§2.2广义逆矩阵的性质将一个非零矩阵分解为一个列满秩矩阵与一个行满秩矩阵的乘积，是矩阵分解理论中的常见问题。特别是在广义逆矩阵的计算与研究中有着重要的应用。定义2.6设矩阵nmrCA（r＞0），如果存在一个列满秩矩阵rmrCF与一个行-4-满秩矩阵nrrCG使得FGA，则称上式为A的一个满秩分解。定理2.1对任意矩阵nmrCA（r＞0），必存在着矩阵rmrCF和nrrCG使FGA。证明：由rrankA，对A进行若干次初等行变换后，可将A化为行阶梯矩阵B，0GB，其中rrankG。故存在若干个m阶初等矩阵的乘积P，使得BPA，即BPA1，将1P分块为MFP,1,rmrCF，)(rmmCM，便有FGGMFA0,。因F是可逆矩阵1P的前r列，所以F是一个rm列满秩矩阵，G是nr行满秩矩阵，故FGA是A的一个满秩分解。上式FGA是A的一个满秩分解，但是A的满秩分解并不是唯一的。任意取一个r阶非奇异矩阵B,若FGA是一个满秩分解，则显然GBFBA1也是A的一个满秩分解。一、{1}-逆的性质定理2.2设nmCA，则A的Moore-Penrose逆存在且唯一。证设rrankA.若r=0，则A是nm零矩阵，可以验证mn零矩阵满足四个Penrose方程。若r0，则A有满秩分解分解FGA,-5-取HHHHFFFGGGX11，则X满足4个Penrose方程，所以，X是Moore-Penrose广义逆矩阵。设X，Y均满足四个Penrose方程，则YYYAYYAYYAXAXAYAYAXXAYAXXAYAXXAXXAXXXHHHHHHHHHHHHHHHH综上所诉，A存在且唯一。A满足四个Penrose方程的所有方程，所以，A属于15类广义逆矩阵中的任意一类。上面我们证明了A的存在性，所以，任意的类广义逆矩阵都是存在的。对任意的C，定义为00,0,1（2.4）下面给出{1}-逆的一些性质。定理2.3设nmCA，nmCB，C，则（1）}1{)()1(HHAA；（2）}1){()1(AA；（3）若S和T非奇异，则}1){(1)1(1SATSAT；（4）rankArankA1；（5）1AA和AA1均为幂等矩阵且与A同秩；（6）；)())((),()(),()()1()1()1(HHARAARANAANARAAR（7）nIAA1的充要条件是nrankA，mIAA1的充要条件是mrankA；（8）AAABAB1的充要条件是rankAABrank)(，BABABB1的充要条件是rankBABrank)(。证（1）由}1{1AA，有AAAA1，两边同时求共轭转置得-6-HHAAAA1，即HHHHAAAA)(1，由定义知}1{1HHAA。（2）AAAAAAA11,由{1}-逆定义得，}1{1AA。（3）SATSATSASATTSATSATSAT111111,由{1}-逆定义得，}1{111SATSAT。（4）rankAAAArankAArankrankA111,故rankArankA1.。（5）11121AAAAAAAA，故1AA为幂等矩阵，又由AAAAAAAA11121，故AA1为幂等矩阵，所以rankAAArankAAArankrankA)()()1()1(，也即rankAAArank)()1(。同理，rankAAArank)()1(。（6）由)()()()()1()1(ARAAARAARAR,得)()(1ARAAR，类似的，由)()()()()1()1(ANAAANAANAN，得)()1(ANAAN。又因为，)())(())(())(()()1()1()1(HHHHHHHHARAAARAARAARAR，所以AARAARHH1。（7）充分性：nrankA，所以，nAArank1，由AA1为幂等矩阵且非奇异，易知nIAA1。必要性：由nIAA1，nAArank1，故nrankA。另一式同理可证明。-7-（8）充分性：)()(ARABR，rankAABrank)(，所以，)()(ARABR。所以存在矩阵X，使ABXA，从而AABXABXABABAABAB)1()1()()(。必要性：rankAabrankAABABrankrankA)(])([1，故rankAABrank)(。另一式同理可证明。性质（5）逆命题仍然成立，即定理2.4设nm复矩阵A，若存在mn矩阵X，使AX为幂等矩阵，且rankAAXrank)(，则矩阵}1{AX。证明：AX幂等，则AXAXAX，而)()(ARAXR,又rankAAXrank)(，所以，)()(ARAXR，存在矩阵Y,使得AXYA，有AAXYAXAXYAXA，即1AX。二、2,1-逆的性质因为在Penrose方程（1）（2）中，A和X的位置是对称的，所以A{1,2}X与{1,2}AX是等价的，即A和X总是互为2,1-逆。这与通常矩阵A的逆的逆是A本身是一样的。定理2.5设矩阵1,AZY,又设YAZX，则2,1AX。证明：1,AZY，则AAYA，AAZA，AAZAZAAYAAYAZAAXA)(，-8-XYAZYAYAZYAZAZAYYAZAYAZXAX)(，由上2式得，2,1AX。定理2.6给定矩阵A，若1AX，则2,1AX的充要条件是rankArankX。证明：充分性：若1AX，则AAXA，且AX和XA幂等，ArankXArankAXrank)()(，又rankArankX，所以，rankXrankAXArank。由定理2.3得1XA，所以，2,1AX。必要性：1AX，则rankArankX，又2,1AX，根据X为自反广义逆，有X{1}A，则rankXrankA所