本文介绍: 矩阵理论知识点总结

5、广义逆矩阵

逆,生来就是用于解方程组的。

  1. 逆:行列满秩。

  2. 单边逆:左逆列满秩;右逆行满秩。

    1. 求法:高斯消元。

    2. 列满秩矩阵:Ax=b

      1. 行初等变换,可以得到左逆

        A

        L

        1

        A_L^{-1}

        AL1

      2. 有解的充要条件:

        A

        A

        L

        1

        b

        =

        b

        AA_L^{-1}b=b

        AAL1b=b (可拓展为AGb=b)

      3. 唯一解:

        x

        =

        (

        A

        H

        A

        )

        1

        A

        H

        b

        =

        A

        L

        1

        b

        Large x = (A^HA)^{-1}A^Hb=A_L^{-1}b

        x=(AHA)1AHb=AL1b

      4. 需要注意,左逆矩阵不唯一,

        (

        A

        H

        A

        )

        1

        A

        H

        (A^HA)^{-1}A^H

        (AHA)1AH也是列满秩矩阵A的左逆。

    3. 行满秩矩阵:

      1. 列初等变换,可以得到右逆

        A

        R

        1

        A_R^{-1}

        AR1

      2. 行满秩矩阵一定有解,且解不唯一。自由未知数的个数为n-m
      3. A

        H

        (

        A

        A

        H

        )

        1

        A^H(AA^H)^{-1}

        AH(AAH)1是A的一个右逆。右逆矩阵不唯一。

      4. 需要注意:

        A

        H

        (

        A

        A

        H

        )

        1

        b

        A

        R

        1

        b

        large A^H(AA^H)^{-1}b {ne} A_R^{-1}b

        AH(AAH)1b=AR1b 通常情况求出来的是两个不同的解,均满足Ax=b。这是因为行满秩矩阵Ax=b的解不唯一。

      5. 行满秩矩阵的A+是A的一个右逆。
  3. 广义逆:任何矩阵都存在广义逆矩阵。

    1. AGA=A 等价于 G是A的广义逆矩阵。

    2. 单边逆是广义逆的特殊情况。

    3. Ax = b有解的充要条件:rank A = rank (A b)

    4. 当A列满秩,且rank A = rank (A b) ,则有唯一解。

    5. 当A非列满秩,且有解,则有无穷多解。

    6. Gb = x , 且Ax = b,则称G是A的一个广义逆。

    7. 广义逆矩阵不唯一,零矩阵的广义逆矩阵是任意矩阵。

    8. 广义逆矩阵通常记作

      A

      A^-

      A,区别于

      A

      1

      A^{-1}

      A1

    9. 广义逆矩阵的秩 大于等于 A的秩。取等号,等价于 G具有自反性。

    10. 广义逆矩阵是逆矩阵的推广,当A是可逆矩阵的时候,A=A-1

    11. 广义逆矩阵有逆矩阵类似的性质:

      性质 广义逆矩阵 逆矩阵
      定义 对于矩阵A,存在矩阵B,使得ABA=A 对于方阵A,存在矩阵B,使得AB=BA=I
      记号

      A

      A^-

      A

      A

      1

      A^{-1}

      A1

      行为 对于任意向量b,满足

      A

      A

      AA^-

      AAb=b

      对于任意向量b,满足

      A

      A

      1

      AA^{-1}

      AA1b=b

      矩阵乘法

      A

      A

      AA^-

      AAA=A

      AA

      1

      ^{-1}

      1A=A

      矩阵的秩 rank(

      A

      A

      A^-A

      AA)=rank(

      A

      A

      AA^-

      AA)=rank(A)

      rank(AA

      1

      ^{-1}

      1)=rank(A)=n(满秩)

      逆的存在 广义逆存在于可逆和不可逆矩阵中 逆存在于方阵(可逆矩阵)中
      唯一性 可以有多个不同的广义逆 逆是唯一的
      幂等性

      A

      A

      AA^-

      AA

      A

      A

      A^-A

      AA是幂等矩阵

      A

      A

      1

      =

      A

      1

      A

      =

      I

      AA^{-1}=A^{-1}A=I

      AA1=A1A=I是幂等矩阵

      数乘 aA的广义逆为

      1

      a

      A

      a

      0

      frac{1}{a} A^-,ane0

      a1Aa=0

      aA的逆为

      1

      a

      A

      1

      a

      0

      frac{1}{a} A^{-1},ane0

      a1A1a=0

      正交投影

      (

      A

      A

      )

      H

      =

      A

      A

      (AA^-)^H=AA^-

      (AA)H=AA,则

      A

      A

      =

      P

      R

      (

      A

      )

      AA^-=P_{R(A)}

      AA=PR(A)

      A

      A

      1

      =

      I

      AA^{-1}=I

      AA1=I,I是从Cn到R(A)的正交投影

  4. 自反广义逆

    1. 定义:AGA=A且GAG=G。存在且不唯一。
    2. 求法:
      1. 最大秩分解法:A=BD,

        A

        r

        =

        D

        R

        1

        B

        L

        1

        A_r^-=D_R^{-1}B_L^{-1}

        Ar=DR1BL1

      2. 构造法:X、Y是A的广义逆,则Z = XAY是自反广义逆,当然YAX也是自反广义逆。
      3. 公式法:X=(AHA)AH,Y =AH(AAH)都是自反广义逆
        1. R(AH)=R(AHA), N(A)=N(AHA)
        2. 存在D,使得AH=AHAD
        3. 代入可证明AXA=A
        4. rank(X)<=rank(AH)=rank(AHA) = rank(AHA(AHA)AHA)=rank(AHAXA)<=rank(X)
        5. 所以X是自反广义逆
    3. A_是自反广义逆的充要条件是rank(A)=rank(A)
      1. 必要性:AGA=A且GAG=G,则rank(A)=rank(AGA)<=rank(G)=rank(GAG)<=rank(A)
      2. 充分性:
        1. R(GA)属于R(G),R(GA)=R(A)=R(G),则R(G)=R(GA);
        2. GE=G,则存在X,GAX=G
        3. A=AGA=AGAXA=AXA
        4. 由构造法知G是自反广义逆
    4. 几何性质
      1. R

        (

        A

        )

        N

        (

        A

        H

        )

        =

        C

        m

        R(A)oplus N(A^H)=C^m

        R(A)N(AH)=Cm

      2. R

        (

        A

        H

        )

        N

        (

        A

        )

        =

        C

        n

        R(A^H)oplus N(A)=C^n

        R(AH)N(A)=Cn

      3. R

        (

        A

        )

        N

        (

        A

        r

        )

        =

        C

        m

        R(A)oplus N(A^-_r)=C^m

        R(A)N(Ar)=Cm

      4. R

        (

        A

        r

        )

        N

        (

        A

        )

        =

        C

        n

        R(A^-_r)oplus N(A)=C^n

        R(Ar)N(A)=Cn

      5. R

        (

        A

        r

        )

        =

        R

        (

        A

        H

        )

        R(A^-_r) = R(A^H)

        R(Ar)=R(AH)

      6. N

        (

        A

        r

        )

        =

        N

        (

        A

        H

        )

        N(A^-_r) = N(A^H)

        N(Ar)=N(AH)

  5. MP广义逆

    1. 定义:AGA=A,GAG=G,(GA)H=GA, (AG)H=AG
    2. 存在且唯一
    3. 计算方法:
      1. 最大值分解法:

        A

        +

        =

        D

        H

        (

        D

        D

        H

        )

        1

        (

        B

        H

        B

        )

        1

        B

        H

        A^+=D^H(DD^H)^{-1}(B^HB)^{-1}B^H

        A+=DH(DDH)1(BHB)1BH

      2. 奇异值分解法:

        A

        +

        =

        V

        H

        D

        +

        U

        A^+=V^HD^+U

        A+=VHD+U (A = UDVH

      3. 注意:最大值分解不唯一,然而最大值分解的这种乘积,即A+是唯一的。
    4. A+的性质:
      1. 自反性:(A+)+=A
      2. 唯一性:

        A

        +

        =

        (

        A

        H

        A

        )

        +

        A

        H

        =

        A

        H

        (

        A

        A

        H

        )

        +

        =

        D

        H

        (

        D

        D

        H

        )

        1

        (

        B

        H

        B

        )

        1

        B

        H

        A^+=(A^HA)^+A^H=A^H(AA^H)^+=D^H(DD^H)^{-1}(B^HB)^{-1}B^H

        A+=(AHA)+AH=AH(AAH)+=DH(DDH)1(BHB)1BH

      3. 几何性质:

        R

        (

        A

        +

        )

        =

        R

        (

        A

        H

        )

        R(A^+) = R(A^H)

        R(A+)=R(AH) (这是自反广义逆具有的性质)

      4. 正交投影性:

        A

        A

        +

        =

        P

        R

        (

        A

        )

        ,

        A

        +

        A

        =

        P

        R

        (

        A

        H

        )

        AA^+=P_{R(A)},A^+A=P_{R(A^H)}

        AA+=PR(A),A+A=PR(AH)

      5. 行列子空间相等性:R(A)=R(AH)的充要条件是

        A

        +

        A

        =

        A

        +

        A

        A^+A=A^+A

        A+A=A+A

      6. (

        A

        H

        A

        )

        +

        =

        A

        +

        (

        A

        H

        )

        +

        =

        A

        +

        (

        A

        A

        H

        )

        +

        A

        =

        A

        H

        (

        A

        A

        H

        )

        +

        (

        A

        H

        )

        +

        (A^HA)^+=A^+(A^H)^+=A^+(AA^H)^+A=A^H(AA^H)^+(A^H)^+

        (AHA)+=A+(AH)+=A+(AAH)+A=AH(AAH)+(AH)+

      7. A

        +

        A

        =

        (

        A

        H

        A

        )

        +

        (

        A

        H

        A

        )

        =

        (

        A

        H

        A

        )

        (

        A

        H

        A

        )

        +

        A^+A=(A^HA)^+(A^HA)=(A^HA)(A^HA)^+

        A+A=(AHA)+(AHA)=(AHA)(AHA)+

    5. 若A是Hermite矩阵:
      1. (

        A

        2

        )

        +

        =

        (

        A

        +

        )

        2

        (A^2)^+=(A^+)^2

        (A2)+=(A+)2

  6. 矩阵方程通解:

    1. AXB=D有解的充要条件:AADBB=D
    2. 通解:X=ADB+Y-AAYBB
    3. Ax=b有解的充要条件:AAb=b
    4. 通解:x = Ab +y-A-1Ay
  7. 相容方程组的最小范数解:

    1. 相容方程组:有解方程组
    2. AGA=A,(GA)H=GA
    3. Gb=x是最小范数解
  8. 不相容方程组的最小二乘解:

    1. 不相容方程组:无解方程组
    2. AGA=A,(AG)H=AG
    3. Gb=x是最小二乘解之一,最小二乘解的通解:

      x

      =

      G

      b

      +

      (

      E

      G

      A

      )

      b

      x = Gb+(E-GA)b

      x=Gb+(EGA)b

  9. 不相容方程组的最小二乘解有时候不够好,最小二乘解的最小范数解叫最佳逼近解:

    1. 最佳逼近解:A+b=x
    2. 如果方程组是相容方程组,则A+b是最小范数解
  10. 关于广义逆的运算规则

    1. (

      A

      )

      H

      =

      (

      A

      H

      )

      (A^-)^H=(A^H)^-

      (A)H=(AH)

    2. B

      =

      S

      A

      T

      ,

      B

      =

      T

      1

      A

      S

      1

      B=SAT,B^-=T^{-1}A^-S^{-1}

      B=SAT,B=T1AS1

    3. (

      A

      B

      )

      +

      =

      B

      +

      A

      +

      (AB)^+=B^+A^+

      (AB)+=B+A+的充要条件:

      R

      (

      A

      H

      A

      B

      )

      R

      (

      B

      )

      ,

      R

      (

      B

      B

      H

      A

      H

      )

      R

      (

      A

      H

      )

      R(A^HAB)sub R(B),R(BB^HA^H)sub R(A^H)

      R(AHAB)R(B),R(BBHAH)R(AH)

    4. (

      A

      B

      )

      +

      =

      B

      +

      A

      +

      (AB)^+=B^+A^+

      (AB)+=B+A+的充分条件:A列满秩(AHA满秩),B行满秩(BBH满秩)

原文地址:https://blog.csdn.net/weixin_45827703/article/details/135280176

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。

如若转载,请注明出处:http://www.7code.cn/show_52120.html

如若内容造成侵权/违法违规/事实不符,请联系代码007邮箱:suwngjj01@126.com进行投诉反馈,一经查实,立即删除!

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注