6 广义逆

178 阅读1分钟

广义逆矩阵

motivation

不是方阵、不满秩,还不能求逆了?

definition

A和G是复数域矩阵(不一定方,但G的形状和A转置一样),如果满足下面4条件,则GGAA的广义逆矩阵。

  • AGA=AAGA = A
  • GAG=GGAG = G
  • (AG)H=AG(AG)^H = AG
  • (GA)H=GA(GA)^H = GA

性质定理

  • AA可逆,则A1=GA^{-1}=G
  • A=OA=O,则G=OG=O
  • G必然存在且唯一
  • (AHA)+=A+(AH)+(A^HA)^+=A^+(A^H)^+
  • (AAH)+=(AH)+A+(AA^H)^+=(A^H)^+A^+
  • A+=(AHA)+AH=AH(AAH)+A^+ = (A^HA)^+A^H = A^H(AA^H)^+
  • U,V酉矩阵,则(UAV)+=VHA+UH(UAV)^+=V^HA^+U^H
  • A+AB=A+ACAB=ACA^+AB=A^+AC \Leftrightarrow AB=AC

计算

先把AA满秩分解A=BCA=BC

之后A+=CH(CCH)1(BHB)1BHA^+=C^H(CC^H)^{-1}(B^HB)^{-1}B^H

应用

线性方程组Ax=bAx=b无解,如何求近似解?

方法是求最小二乘解,这个之前推过AHAx=AHbA^HAx=A^Hb

长度最小的称为极小最小二乘解

他的通解是x=A+b+(IA+A)y,yCnx = A^+b + (I-A^+A)y, \forall y \in C^n

其中A+bA^+b是唯一的极小最小二乘解