毕业论文-广义逆矩阵的求法探讨

页,共21页0引言矩阵逆的概念只对非奇异方阵才有意义.但是，在实际问题中，我们碰到的矩阵并不都是方阵，即使是方阵，也不都是非奇异的。因此，有必要推广逆矩阵的概念.为此，本文给出了广义逆矩阵的定义，并利用广义逆的性质，给出其计算方法。1广义逆矩阵的概念与定理定义1.1设是的矩阵，若的矩阵满足如下四个方程的全部或者一部分，则称为的广义逆矩阵，简称广义逆.(1.1)(1.2)(1.3)(1.4)则称是的逆，记为.如果某个只满足（1.1）式，为的{1}广义逆，记为G{1}；如果另一个满足（1.1），（1.2）式，则称为的{1,2}广义逆，记为{1，2}；如果{1，2，3，4}，则是逆等.下面介绍常用的5种{1},{1,2}，{1,3},{1,4},{1,2,3,4}每一种广义逆矩阵又都包含着一类矩阵，分述如下：{1}中任意一个确定的广义逆，称作减号广义逆，或g逆，记为；{1,2}中任意一个确定的广义逆，称作自反减号逆，记为；{1,3}中任意一个确定的广义逆，称作最小范数广义逆，记为；{1,4}中任意一个确定的广义逆，称作最小二乘广义逆，记为；{1,2,3,4}：唯一一个，称作加号逆，或，记为.定义1.2设是的矩阵（,当时，可以讨论），若有一个的矩阵（记为）存在，使下式成立，则称为的减号广义逆或者逆：(1.5)当存在时，显然满足上式，可见减号广义逆是普通广义逆矩阵的推广；另外，由得可见，当为的一个减号广义逆时，就是的一个减号广义逆.定义1.3设的特征值为则称为矩阵的正奇异值，简称奇异值. 定义1.4设矩阵，如果时存在；或者当时，存在有，称这两种长方阵为最大秩方阵（满秩方阵），前者又称行最大秩矩阵（行满秩矩阵），后者又称为列最大秩矩阵（列满秩矩阵）.定义1.5设是矩阵,若有矩阵满足(或),则称为的右逆(或左逆),记为(或).定理1.1设是的矩阵，则的逆存在且唯一.证明先证的存在性.设的奇异值分解其中，是的非零奇异值，与是酉矩阵.令容易验证满足四个方程，因此存在.下面证的唯一性.假定也是满足4个方程，则因此,说明是唯一的，且若是非奇异矩阵，容易验证满足4个方程，此时.由此可见逆把逆推广到所有矩阵（甚至零矩阵）.定理1.2设，，存在阶的可逆矩阵及阶可逆矩阵，使则阶矩阵使得的充分必要条件是其中分别是阶任意矩阵.证明先证必要性，由条件有阶及阶可逆矩阵，使那么根据应满足的,有再令分块如题设要求，代入上式所以,于是有得到再证充分性，由于则引理1.1对于任意的矩阵，它的减号逆总存在，但不唯一，并且是的一个减号逆【1，2】.引理1.2对于任意的矩阵，它的极小范数总存在，但不唯一，并且是的一个极小范数逆【1‘2】.引理1.3对于任意矩阵,它的最小二乘逆总存在,但不唯一,并且它是的一个最小二乘逆【1，2】.引理1.4对于任意矩阵,它的加号逆总存在，并且唯一.其中这里是的满秩分解式【1，2，3】.定理1.3是矩阵,若是行满秩矩阵,则总有；是列满秩矩阵，则总有；，则总有，其中是的满秩分解式.定理1.4设则可将做满秩分解（或的最大秩分解）其中是阶矩阵，且.将一非列或非行满秩的非零矩阵表示为一列满秩和一行满秩的矩阵的积的分解称为满秩分解.在各种广义逆的直接计算方法中,几乎都要对矩阵进行满秩分解,例如分解等等.但当计算某些广义逆时,分解将带来大量非必要的计算,因而有必要对满秩分解的方法进行简化,为此,我们首先用构造性方法证明下述定理.定理1.5对任意矩阵,总存在着矩阵和矩阵，使得成立.证明设，则必有一个最大线性无关列，，…，，故令=[，，…，]于是有非奇异矩阵，使,亦即有(1.6)成立，其中为阶数适当的零矩阵，再另置换矩阵便有,于是由（1）知，=(1.7)其中,且显然有，.类似地可证存在着和，使有,成立，倘令(1.8)(1.9)同样有.特别，若A为行满秩或者列满秩，则与中之一为单位阵，定理依然成立.定理1.6对任何的矩阵，都有性质1.1(1)的充分必要条件是，此时,称为的一个左逆，记为.(2)的充分必要条件是，此时=称为的一个右逆，记为.证明(1)充分性，若则所以必要性，若，则存在阶及阶可逆矩阵,使或由定理1.2可得,则有即，于是有由于所以是可逆阵，那么所以，可取(2)同理可证性质(2)，可逆，有所以，可取2广义逆矩阵的计算方法2.1广义逆矩阵A+的奇异值分解法设矩阵，由定理1.1知存在并且唯一，当时，则有奇异值分解：其中，，为的奇异值，则具有如下形式：.例1用奇异值分解求，其中.解的奇异值分解为,所以=.例2设用奇异值分解法求.解因此特征值求出对应于所以=2.2广义逆矩阵的最大秩分解法的矩阵的秩，的最大秩分解为其中是阶矩阵，是阶矩阵，且，则(2.1)特别当时（行满秩阵）(2.2)当时（列满秩阵）(2.3)例3求矩阵的逆.解首先求得的满秩分解为，故==.2.3极限法求广义逆矩阵设是阶矩阵，则(2.4)证明因为由定理1.6得设用极限法求.解因为因此2.4广义逆矩阵的满秩分解法对任意矩阵，由定理1.5知，其中是阶矩阵，是阶矩阵，且，再由性质1.1可得

如果A是实矩阵，有设为矩阵的最大秩分解,则的广义逆矩阵的一般形式为.例5设，求其广义逆矩阵.解首先对进行最大秩分解,对作行初等变换如下:所以的最大秩分解为=由定理1.3知，这里为3阶可逆方阵，故为行满秩矩阵，故可取=从而=例6设矩阵=求.解有满秩分解为取=，从而=，得取,得得在依据性质1.1的(1.5)及(1.6)可分别求出于是得到2.4初等变换法求广义逆矩阵方法和步骤：经过一系列的初等行或初等列变换总可以将写成式的形式,这里分别是m和矩阵，由定理1.2，则的全部广义逆为这里、分别是任意的例7.解由上述定理,首先要将写成式的形式.为此,将作初等变换得=设,,,,则,从而，有=例8设,求广义逆.解=于是,.所以的减号广义逆为,其中.以上介绍了的初等变换法，那么我们现在给定一个矩阵,总有，有定理1.3知当时，有，当时，有,当时，有，其中是的满秩分解式.我们可以看出要求矩阵的任何一种广义逆矩阵,关键是求出一个.那么下给出了利用初等变换法求出的具体方法.设,(不必限制)则存在阶可逆矩阵使得则，令由于所以是的一个广义逆矩阵(.据此，我们对下面分块矩阵进行初等变换：=因此,.同理，对下面的分块矩阵施行初等变换：=因此，.这里、均指可逆矩阵.例9设=，求的最小二乘逆.解因为，所以.对下列矩阵施行初等行变换有所以=.例10设，求最小范数逆.解因为=，所以=,对下列矩阵施行初等变换有所以.上述例题给出的求广义逆矩阵和的方法,简便易行且使各种广义逆矩阵的计算得到了彻底解决.致谢本文是在的指导下完成的，在此衷心的感谢周教授的细心的指导，才能顺利完成本论文.参考文献[1]李宗铎.求逆矩阵的一个方法[J]数学通报，1983（11）：15—16.[2]南京大学数学系计算数学专业.线性代数[M].北京：科学出版社，1978：97.[3]任晓红.球广义逆矩阵Aˉ的初等变换法[J].西北轻工业学院学报，2000（2）：105—106.[4]周琳.介绍广义逆矩阵及其计算方法[J].本溪冶金高等专科学校学报，2001（2）：43—45.[5]北京大学数学力学系.高等代数[M].北京高等教育出版社，1978：187.[6]杨明，刘先忠.矩阵论[M].华中科技大学出版社，2005：95—98.[7]刘丁酉.矩阵分析[M].武汉大学出版社，2004：241—241.[8]苏育才，姜翠波等.矩阵理论[M].科学出版社，2003：192.[9]吴强.基于矩阵初等变换的矩阵分解法[J].数学理论与应用，2000，20（4）.[9]刘宣黄.广义逆矩阵的计算方法[J].江西电力职业技术学院学报，2008，21（1）：44—47.[10]Fuzhenzhang，MatrixTheory，Springer，1999.[11]HornRA，JohnsonCR.1989.MatrixAnalysis（矩阵分析）.杨奇.天津：天津大学出版社.[12]DJFieldWhatisthegoalofsensorycoding?1994(4)MHeiler.CSchnorrLearningsparse.[13]POHoyerNon-negativernatrixfactorizationwithsparsenessconstraints2004(9)基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统基于51单片机的嵌入式Internet系统的设计与实现单片机监测系统在挤压机上的应用MSP430单片机在智能水表系统上的研究与应用基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用单片机在高楼恒压供水系统中的应用基于ATmega16单片机的流量控制器的开发基于MSP430单片机的远程抄表系统及智能网络水表的设计基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计基于单片机