广义逆矩阵ppt课件.ppt

《广义逆矩阵ppt课件.ppt》由会员分享，可在线阅读，更多相关《广义逆矩阵ppt课件.ppt（42页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第七章第七章 广义逆矩阵广义逆矩阵 广义逆矩阵是逆矩阵的推广，与线性方程组的求解有密切广义逆矩阵是逆矩阵的推广，与线性方程组的求解有密切联系。给定一个线性方程组联系。给定一个线性方程组 Ax=b，当矩阵当矩阵A可逆时，线性可逆时，线性方程组的解可表示为方程组的解可表示为x=A-1 b当矩阵当矩阵A是奇异矩阵或不是方阵时，线性方程组的解应如何表是奇异矩阵或不是方阵时，线性方程组的解应如何表示呢？当线性方程组是矛盾方程，或者说是不相容方程时，线示呢？当线性方程组是矛盾方程，或者说是不相容方程时，线性方程组能否有其它意义下的解，这种解又应当如何表示呢？性方程组能否有其它意义下的解，这种解又应当如何表

2、示呢？ 把逆矩阵推广到不可逆方阵或长方矩阵上，这就是所谓的把逆矩阵推广到不可逆方阵或长方矩阵上，这就是所谓的广义逆矩阵。广义逆矩阵。 广义逆矩阵具有通常逆矩阵的部分性质，并且在方阵可逆广义逆矩阵具有通常逆矩阵的部分性质，并且在方阵可逆时，它与通常的逆矩阵一致，而且广义逆矩阵可以给出线性时，它与通常的逆矩阵一致，而且广义逆矩阵可以给出线性方程组（包括相容的和矛盾方程组）各种解的统一形式。方程组（包括相容的和矛盾方程组）各种解的统一形式。主要内容主要内容：1广义逆矩阵及其分类广义逆矩阵及其分类2A+的计算的计算3几类弱逆几类弱逆4广义逆矩阵与线性方程组的解广义逆矩阵与线性方程组的解广义逆矩阵方程广

3、义逆矩阵方程 设设A A是是n n阶非奇异矩阵，则存在唯一的逆矩阵阶非奇异矩阵，则存在唯一的逆矩阵A A-1-1，它具有如下性质：它具有如下性质： AAAA1111 AAAAIAA1IAA1)()()()(4321PXAXAPAXAXPXXAXPAAXAHH或者说，或者说， A A-1-1是下述矩阵方程组的解是下述矩阵方程组的解-广义逆矩阵方程广义逆矩阵方程)()()()(4321PXAXAPAXAXPXXAXPAAXAHH设设,nmCA若矩阵若矩阵 满足如下四个（满足如下四个（Penrose) mnCX方程方程则称则称X为为A的的Moor Penrose逆，记为逆，记为A+例例：容易由定义直

4、接验算：容易由定义直接验算：若若,0011A则则021021A存在性证明存在性证明BCACCCBrnrrrmr使得使得在在由满秩分解定理知，存由满秩分解定理知，存设设, 0HHHHBBBCCCX11)()(令可以验证可以验证X满足广义逆矩阵方程满足广义逆矩阵方程mnCX设设 ，A+存在且唯一，即广义矩阵方程组存在且唯一，即广义矩阵方程组nmCA定理定理 有唯一解有唯一解设设, rrankA 若若, 0r则则A是是 阶零矩阵，可以阶零矩阵，可以 nm验证验证 阶零矩阵满足四个方程。阶零矩阵满足四个方程。mn 对于矩阵方程对于矩阵方程)()()()(4321PXAXAPAXAXPXXAXPAAXA

5、HH如果矩阵如果矩阵G仅满足其中的一个或几个时，可以定义仅满足其中的一个或几个时，可以定义不同的广义逆矩阵。不同的广义逆矩阵。因此，共可定义因此，共可定义1544342414CCCC类不同的广义逆。类不同的广义逆。由由A+的存在性可知，的存在性可知，15类广义逆都存在，除类广义逆都存在，除A+是是唯一确定的外，唯一确定的外，其余各类广义逆矩阵都不唯一确定。其余各类广义逆矩阵都不唯一确定。几几 类类 弱弱 逆逆Ai = |G满足第满足第i个个Penrose方程方程mnCG对于矩阵对于矩阵 ，记，记nmCAAi,j = |G满足第满足第i,j个个Penrose方程方程mnCGAi ,j ,k =

6、|G满足第满足第i,j,k个个Penrose方程方程mnCG广义逆集合广义逆集合各类广义逆的关系各类广义逆的关系几种常用的广义逆矩阵几种常用的广义逆矩阵,4, 3 , 2, 1ijikjiAAAAAA A11, ,它的形式记为它的形式记为A A1,21,2, ,它的形式记为它的形式记为A A1,31,3, ,它的形式记为它的形式记为A A1,41,4, ,它的形式记为它的形式记为rAAlAmA-最小二乘广义逆最小二乘广义逆-自反广义逆自反广义逆最小范数广义逆最小范数广义逆A1A1是指仅满足第一个是指仅满足第一个PenrosePenrose方程的广义逆，即若方程的广义逆，即若AAAA-1-1A=

7、A, A=A, 则记则记广义逆广义逆A A- -1AA 说明说明： 1)1)利用初等行变换，可以求得利用初等行变换，可以求得A A- -2 2）A A的减号逆的减号逆A A- -不唯一。不唯一。例例：设设,010101A容易验证容易验证100001,010001CB均均满足满足,AACAAABA故故B B，C C都是都是A A的减号逆的减号逆. .3)3)矩阵矩阵A A有唯一的有唯一的A A- -充分必要条件是充分必要条件是A A为非奇异矩阵，为非奇异矩阵，此时此时A-=A-1定理定理 A1的表示通式的表示通式AYAAAYAGACYAACAmnnmr的通式为则是任意矩阵是一个给定的广义逆，设1

8、,1此此定理表明：只要求出定理表明：只要求出 中的一个元素，就可得到中的一个元素，就可得到 中所有的元素。中所有的元素。1A1A广义逆矩阵广义逆矩阵A A+ +的计算：方法一的计算：方法一 利用利用满秩分解满秩分解如果矩阵如果矩阵A A有满秩分解有满秩分解A=BCA=BC，则有则有A A+ +的表达式，即的表达式，即 HHHHBBBCCCA11)()(1 AA因此广义逆因此广义逆A A+ +是通常逆矩阵概念的一种推广。是通常逆矩阵概念的一种推广。广义逆矩阵广义逆矩阵A A+ +与通常逆矩阵有许多类似的性质，但也有一些不与通常逆矩阵有许多类似的性质，但也有一些不同。同。 如果如果A A是非奇异矩

9、阵，则是非奇异矩阵，则 并且由上面的并且由上面的公式计算出公式计算出 ，从而，从而111BCA11BCA如果矩阵如果矩阵A A是行满秩的，是行满秩的，A A有满秩分解有满秩分解A = A = I Im m A A，则则A A+ +的表达式为的表达式为 1)(HHAAAA如果矩阵如果矩阵A A是列满秩的，是列满秩的，A A有满秩分解有满秩分解A = A I A = A I n n，则则A A+ +的表达式为的表达式为HHAAAA1)(特别地，设特别地，设 为为n维列向量，且维列向量，且 则则 , 0HH1)(1)(HH设设 为为n维行向量，且维行向量，且 则则 , 0例：例： 求广义逆求广义逆0

10、11001A1)(HHAAAAA是行满秩的，故由于解110101101100100101100110121110010111例：设例：设 求求,321AA由由A为列向量，即为列满秩，则为列向量，即为列满秩，则HHAAAA1)(从而从而321141A 若若A既不是行满秩也不是列满秩，则需首先对既不是行满秩也不是列满秩，则需首先对A进行满秩进行满秩分解，再求分解，再求A例：已知例：已知,122121211142A求求A矩阵矩阵A中分别有两行、两列对应成比例，因此中分别有两行、两列对应成比例，因此A既不是行既不是行满秩也不是列满秩满秩也不是列满秩首先利用初等行变换求出首先利用初等行变换求出A的的He

11、rmite标准型标准型H为：为：,000011001021H设设A的满秩分解为的满秩分解为 ，则，则BCA ,211112B11001021C于是于是HHHHBBBCCCA11)()(111112633621161110020111561651224112331广义逆矩阵广义逆矩阵A A+ +的计算：方法二的计算：方法二奇异值法奇异值法设矩阵设矩阵 的奇异值分解为的奇异值分解为A=UDVA=UDVH H其中其中U , V U , V 分别是分别是m m阶、阶、n n阶酉矩阵，阶酉矩阵，nmrCA),(,00021rdiagD则容易验证：则容易验证：.HUVDA其中其中),(,000112111

12、1rdiagD利用此方法，需首先对利用此方法，需首先对A进行奇异值分解。进行奇异值分解。例：设例：设010001000000A求求A先求先求A的奇异值分解。因为的奇异值分解。因为,000010001AAH为为. 0, 1, 1321对应的特征向量为：对应的特征向量为：.1 , 0 , 0,0 , 1 , 0,0 , 0 , 1321TTTxxx令令,21321VVxxxV其中其中,32211xVxxV设设111AVU100100100101000100000010010000 则则 的特征值的特征值AAH把把TTyy1 , 0 , 0 , 0,0 , 1 , 0 , 021扩充为扩充为 的一组

13、标准正交基得：的一组标准正交基得：4RTTyy0 , 0 , 1 , 0,0 , 0 , 0 , 143再令再令,4321yyyyU 则则,000000010001HVUA从而从而HUVA00010010000110000100000000100001100010001000010000100广义逆广义逆A A+ +的性质的性质 设设nmCAAA. 1HHAA. 3TTAA. 2HHAAAA4,. 5AA其中其中,C且且0,00,1AAAAHH5)()(7HHHHAAAAAAAAAA)()(8AAAAAAAAAAHHHH)()(. 6HHHHAAAAAAAnnmmCVCU,1111、 设设9

14、9、若有满秩分解式、若有满秩分解式A=BCA=BC，则则HHUAVUAV)()()()(10AArankAArankArankArankBCA都是酉矩阵，则都是酉矩阵，则12、当、当A 是是Hermite矩阵时，矩阵时，AAAAAAAAAAAA222222)()()()(举例说明广义逆不具有通常意义下逆矩阵的下列性质：举例说明广义逆不具有通常意义下逆矩阵的下列性质：（4 4）A A与与A A+ +的非零特征值并不互为倒数。的非零特征值并不互为倒数。 ABAB（1 1）kkAA（2 2）AAAA（3 3）例证例证1 10001A1011BAA 10111BB0001AB0101210011ABA

15、BAB由由0101A可得可得又又B为满秩矩阵，则为满秩矩阵，则例证例证2 20101A001121AkkAAAAAA可验证可验证由由0111A可得可得考虑非齐次线性方程组考虑非齐次线性方程组, bAx 其中其中 给定，而给定，而mnmCbCA,nCx为为待定向量待定向量。若若,rankAbArank方程组是相容方程组；否则，称为矛盾方程组或不相容方程组。方程组是相容方程组；否则，称为矛盾方程组或不相容方程组。则线性方程组则线性方程组 有解，则称该有解，则称该bAx 关于线性方程组关于线性方程组 的求解问题，常见的有以下几种情形：的求解问题，常见的有以下几种情形：bAx 1）在相容时，若系数矩阵

16、）在相容时，若系数矩阵 ，且非奇异，即，且非奇异，即 nnCA0detA则有则有唯一解唯一解;1bAx但当但当A是奇异方阵或长方矩阵时，它的解是奇异方阵或长方矩阵时，它的解不不唯一，我们可以利用减号逆给出方程组的通解。唯一，我们可以利用减号逆给出方程组的通解。 线性方程组线性方程组 求解求解bAx 2）如果方程组相容，且其解有无穷多个，可求出具有极小范数的如果方程组相容，且其解有无穷多个，可求出具有极小范数的解，即解，即,min xbAx其中其中 为欧氏范数，可以证明满足此条件为欧氏范数，可以证明满足此条件的解是的解是唯一的，称为唯一的，称为极小范数解极小范数解。3）若方程组不相容，则不存在通

17、常意义下的解，但在许多实际）若方程组不相容，则不存在通常意义下的解，但在许多实际问题中，需要求出这样的解：问题中，需要求出这样的解：bAxxnCxmin其中其中 为欧氏范数，称这个问题为求矛盾方程组的最小二乘问为欧氏范数，称这个问题为求矛盾方程组的最小二乘问题，相应的题，相应的x为矛盾方程组的为矛盾方程组的最小二乘解最小二乘解。4）一般说来，矛盾方程组的最小二乘解是不唯一的，但在最小）一般说来，矛盾方程组的最小二乘解是不唯一的，但在最小二乘解的集合中，具有最小范数的解二乘解的集合中，具有最小范数的解xbAxminmin是是唯一的，称之为唯一的，称之为极小范数最小二乘解，或最佳逼近解极小范数最小

18、二乘解，或最佳逼近解.（一）相容方程组的通解（一）相容方程组的通解为线性方程组为线性方程组 的解的充分必要条件是的解的充分必要条件是bAX mnnmCbCXCA,bAX )()()(ARbArankbArank)( ARb mnCGGbX 我们已知我们已知 相容相容bAX 1AG ，其中，其中定理定理 对于任意对于任意 ，都存在，都存在，使，使定理说明定理说明，对于任意的，对于任意的bAXAA,1是线性方程组是线性方程组 的一个特解。的一个特解。bAX 给定一个线性方程组给定一个线性方程组广义逆矩阵与线性方程组的求解有着密切关系。利用减号逆、广义逆矩阵与线性方程组的求解有着密切关系。利用减号逆

19、、最小范数广义逆、最小二乘广义逆以及加号逆可以给出上述诸最小范数广义逆、最小二乘广义逆以及加号逆可以给出上述诸问题的解。问题的解。定理定理 齐次线性方程组齐次线性方程组 的通解是的通解是0AXYAAIX)(证明证明: 对于任意向量对于任意向量 ，成立，成立nCY 其中其中 是任意向量。是任意向量。nCY 0)()(YAAAAYAAIAAX即即 是齐次线性方程组是齐次线性方程组 的解。的解。0AXYAAIX)(设设X0是齐次线性方程组是齐次线性方程组 的任一解，则的任一解，则0AX0000)()(XAAIXAAIAXAX因此，因此， 是齐次线性方程组是齐次线性方程组 的通解。的通解。YAAIX)

20、(0AX推论推论 相容线性方程组相容线性方程组 的通解为的通解为bAX YAAIbAX)(其中其中 是任意向量。是任意向量。nCY 例例1、求解、求解221232321xxxxx将将方程组改写为矩阵形式方程组改写为矩阵形式, bAx 其中其中210121A21b由于由于, 2bArankrankA所以该方程组是相容的。所以该方程组是相容的。首先求得首先求得A的一个减号的一个减号逆。 由由A是行满秩矩阵，则是行满秩矩阵，则1HHAAAA832645141从而从而原原方程组的通解为方程组的通解为YAAIbAX)(32132132123192461036913141yyyyyyyyy其中其中 为任意

21、向量。为任意向量。Tyyyy321,定义定义 相容线性方程组相容线性方程组 的所有解中的所有解中2范数最范数最小的解称小的解称 为方程组的最小范数解，记为为方程组的最小范数解，记为bAX（二）相容方程组的最小范数解（二）相容方程组的最小范数解bAXmX定理定理 相容线性方程组相容线性方程组 的最小范数解是唯一的最小范数解是唯一的，并且可表示为的，并且可表示为bAX 其中其中 是是A的最小范数广义逆。的最小范数广义逆。mA例例2、求方程组、求方程组221232321xxxxx的的最小范数解最小范数解由于由于A为行满秩矩阵，因此为行满秩矩阵，因此 为满秩方阵，则有为满秩方阵，则有TAA1TTAAA

22、A所以所以1THTHmAAAAAAA即即AAm832645141从而从而bAxm191013141此解此解即是即是32132132123192461036913141yyyyyyyyy中欧氏范数最小的一个中欧氏范数最小的一个 一个线性方程组一个线性方程组 是矛盾方程组或不是矛盾方程组或不相容方程组，它没有通常意义下的解，但可以寻求相容方程组，它没有通常意义下的解，但可以寻求该方程组在某种含义下的近似解。该方程组在某种含义下的近似解。bAX （三）不相容方程组的最小二乘解（三）不相容方程组的最小二乘解定义定义 不相容方程组不相容方程组 的最小二乘解的最小二乘解bAX X2minbAXbXAnCX

23、定义为满足下列条件的近似解定义为满足下列条件的近似解说明：和其它任何近似解相比较，说明：和其它任何近似解相比较， 所导致的误差平方和所导致的误差平方和X22bAX 最小。最小。矛盾方程组的最小二乘解导致的误差平方和矛盾方程组的最小二乘解导致的误差平方和 是是22bAX 唯一的，但最小二乘解不一定唯一。唯一的，但最小二乘解不一定唯一。定理定理：设：设 是一个最小二乘解，则矛盾方程组的最小二是一个最小二乘解，则矛盾方程组的最小二Gbx 乘解的通解为乘解的通解为yGAIGbx其中其中y为任意向量为任意向量定理定理 设设 ，则，则 是不相容方程组是不相容方程组mnCGGbX 3 , 1AGbAX 的最

24、小二乘解的充分必要条件是的最小二乘解的充分必要条件是例例3、求矛盾方程组、求矛盾方程组00212212121xxxxxx的的最小二乘解最小二乘解系数矩阵系数矩阵A和向量和向量b为为001,111221bA由由A为列满秩矩阵，则可求得为列满秩矩阵，则可求得A的一个最小二乘逆为：的一个最小二乘逆为：HHHHlAAAAAAA1147174111于是，求得一个最小二乘解为于是，求得一个最小二乘解为bAxl74111定理定理 不相容方程组不相容方程组 的最佳逼近解是唯一的最佳逼近解是唯一的，并且的，并且bAX 定义定义 不相容方程组不相容方程组 的最佳逼近解定义为满足下的最佳逼近解定义为满足下列条件的最

25、小二乘解，记为列条件的最小二乘解，记为bAX （四）不相容方程组的最佳逼近解（四）不相容方程组的最佳逼近解可以看出不相容方程组可以看出不相容方程组 的最佳逼近解是方程组的所的最佳逼近解是方程组的所有最小二乘解中范数最小的近似解。有最小二乘解中范数最小的近似解。 bAX *X22min*XXX 其中其中 是是方程组方程组 的最小二乘解的集合。的最小二乘解的集合。bAX bAX*说明说明：由于加号逆既是减号逆又是最小范数逆、最小二乘逆，由于加号逆既是减号逆又是最小范数逆、最小二乘逆，故对于方程组故对于方程组 ，不论其是否有解，均可用加号逆来，不论其是否有解，均可用加号逆来bAx 设设y为任意向量为

26、任意向量,则则：1 相容时，相容时，bAx yAAIbAx是是通解；通解；bAx是是最小范数解最小范数解2 不相容时，不相容时，bAx yAAIbAx是是最小二乘解的通解；最小二乘解的通解；bAx是是最小二乘解；最小二乘解；3bAx是是矛盾方程组矛盾方程组 的最佳逼近解解；的最佳逼近解解；bAx 例例4、求、求 的最佳逼近解。的最佳逼近解。202101A01bbAx 解解: 首先求首先求A的广义逆，对的广义逆，对A进行满秩分解进行满秩分解,BCA 其中其中101,21CB则由则由1HHCCCCHHBBBB1215110121210021101BCATbAx1, 0, 1101则最佳逼近解为则最佳逼近解为则则A的加号逆为：的加号逆为：