广义逆矩阵

2026/6/4 · 离散数学

基本定义

1 矩阵的左右逆

从广义的角度讲,任何一个矩阵 GG 都可以称为矩阵 AA 的逆矩阵,若它与矩阵 AA 的乘积等于单位矩阵 II,即满足 GA=IGA = I。根据矩阵 AA 本身的特点,满足这一定义的矩阵 GG 存在以下三种可能的情况:

  • GG 存在,并且唯一;
  • GG 存在,但不唯一;
  • GG 不存在。

rankA<minn,m\operatorname{rank} A < min{n,m}时不存在左右逆,但可以有广义逆

例:

A1=[221112101],A2=[485723],A3=[131251]A_1 = \begin{bmatrix} 2 & -2 & -1 \\ 1 & 1 & -2 \\ 1 & 0 & -1 \end{bmatrix}, \quad A_2 = \begin{bmatrix} 4 & 8 \\ 5 & -7 \\ -2 & 3 \end{bmatrix}, \quad A_3 = \begin{bmatrix} 1 & 3 & 1 \\ 2 & 5 & 1 \end{bmatrix}

A1A_1 有唯一逆矩阵

A2A_2 存在多个 2×32 \times 3 矩阵 LL 使得 LA2=I2×2LA_2 = I_{2 \times 2}(有左逆) 例如

A2=[485723],L=[72068170],L=[037025],A_2 = \begin{bmatrix} 4 & 8 \\ 5 & -7 \\ -2 & 3 \end{bmatrix}, \quad L = \begin{bmatrix} 7 & 2 & 0 \\ 68 & 17 & 0 \end{bmatrix}, \quad L = \begin{bmatrix} 0 & 3 & 7 \\ 0 & 2 & 5 \end{bmatrix}, \cdots

对矩阵 A3A_3,没有任何 3×23 \times 2 矩阵使得 G3A3=I3×3G_3 A_3 = I_{3 \times 3},但存在多个 3×23 \times 2 矩阵 RR,使得 A3R=I2×2A_3 R = I_{2 \times 2}(有右逆) 例如

A3=[131251],R=[111031],R=[110021],A_3 = \begin{bmatrix} 1 & 3 & 1 \\ 2 & 5 & 1 \end{bmatrix}, \quad R = \begin{bmatrix} 1 & 1 \\ -1 & 0 \\ 3 & -1 \end{bmatrix}, \quad R = \begin{bmatrix} -1 & 1 \\ 0 & 0 \\ 2 & -1 \end{bmatrix}, \cdots

1.1 左右伪逆矩阵

满足 LA=ILA = I,但不满足 AL=IAL = I 的矩阵 LL 称为矩阵 AA 的左逆矩阵 满足 AR=IAR = I,但不满足 RA=IRA = I 的矩阵 RR 称为矩阵 AA 的右逆矩阵

  • 仅当 nmn \geq m 时,矩阵 ARn×mA \in \mathbb{R}^{n \times m} 可能有左逆矩阵;

  • 仅当 nmn \leq m 时,矩阵 ARn×mA \in \mathbb{R}^{n \times m} 可能有右逆矩阵。

一个矩阵的左逆矩阵或右逆矩阵往往非唯一,下面考虑它们的一个特殊解:

  • 考查 n>mn > mAA 具有满列秩 (rankA=m\operatorname{rank} A = m) 的情况。此时,m×mm \times m 矩阵 ATAA^T A 是可逆的,容易验证 L=(ATA)1ATL = (A^T A)^{-1} A^T 满足 LA=IL A = I,这个左逆矩阵是唯一确定的,常称为 左伪逆矩阵

AAn×mn \times m 矩阵,且列满秩。则存在左逆矩阵 LL 满足 LA=ImLA = I_m。一个特解为 L0=(ATA)1ATL_0 = (A^T A)^{-1} A^T。所有左逆可表示为 L=L0+ZL = L_0 + Z,其中 ZZ 是任意 m×nm \times n 矩阵且满足 ZA=0ZA = 0ZZm×nm \times n 个未知数,但 ZA=0ZA = 0 只给出了 m×mm \times m 个方程。于是每个行向量有 nmn - m 个自由参数,mm 行共有 m(nm)1m(n - m) \ge 1 个自由参数。自由参数可以取无穷多个值,所以 ZZ 有无穷多个解。

  • 考查 n<mn < mAA 具有满行秩 (rankA=n\operatorname{rank} A = n) 的情况。此时,n×nn \times n 矩阵 AATA A^T 是可逆的,容易验证 R=AT(AAT)1R = A^T (A A^T)^{-1} 满足 AR=IA R = I,这个右逆矩阵是唯一确定的,常称为 右伪逆矩阵

2 矩阵的广义逆

考虑一个 n×mn \times m 维的秩亏缺矩阵 AA(即rank(A)<min{m,n}\operatorname{rank}(A) < \min\{m, n\})。n×mn \times m 维的秩亏缺矩阵的逆矩阵称为 广义逆矩阵。令 AA^- 表示 AA 的广义逆矩阵。

无论 AA=Im×mA^-A = I_{m \times m} 还是 AA=In×nAA^- = I_{n \times n} 都不可能成立。有必要使用三个矩阵的乘积定义一个秩亏缺矩阵 AA 的逆矩阵。

考虑线性矩阵方程 Ax=yAx = y 的求解。两边左乘 AAAA^-,则有 AAAx=AAyAA^-Ax = AA^-y。若 AA^-AA 的广义逆矩阵,则 Ax=yx=AyAx = y \Rightarrow x = A^-y,代入 AAAx=AAyAA^-Ax = AA^-y 可得 AAAx=AxAA^-Ax = Ax,我们希望该式对任意非零向量 xx 均应成立,则必须要求下列约束条件满足:AAA=AAA^-A = A。满足 AAA=AAA^-A = A 的矩阵 AA^- 称为 AA广义逆矩阵

2.1 Moore-Penrose逆矩阵

满足 AA1A=AAA^{-1}A = A 的矩阵 A1A^{-1} 不是唯一的,存在明显的缺陷。只能保证 A1A^{-1}AA 的广义逆矩阵,并不能保证 AAA1A^{-1} 的广义逆矩阵

考虑原矩阵方程 Ax=yAx = y 的解方程 x=A1yx = A^{-1}y 的求解。两边左乘 A1AA^{-1}A,则有 A1Ax=A1AA1yA^{-1}Ax = A^{-1}AA^{-1}y。由于矩阵 AAA1A^{-1} 的广义逆矩阵,则 x=A1yAx=yx = A^{-1}y \Rightarrow Ax = y,代入 Ax=yAx = y 可得 A1y=A1AA1yA^{-1}y = A^{-1}AA^{-1}y,因此

A1AA1=A1A^{-1}AA^{-1} = A^{-1}也需要满足

我们还希望当 AA 是列满秩或者行满秩时,广义逆矩阵 A1A^{-1} 能够包括左伪逆矩阵和右伪逆矩阵在内的特例:左伪逆矩阵 L=(ATA)1ATL = (A^TA)^{-1}A^T 满足 LA=ILA = I 但不存在 AL=IAL = I,但 AL=A(ATA)1AT=(AL)TAL = A (A^TA)^{-1}A^T = (AL)^T,同样,右伪逆矩阵满足 RA=(RA)TRA = (RA)^T。所以对称性也要满足

  • Moore-Penrose逆矩阵 令 AA 是任意 n×mn \times m 矩阵,称 A+A^+AA 的广义逆矩阵,若 A+A^+ 满足以下四个条件(常称 Moore-Penrose 条件)则称为Moore-Penrose逆矩阵
    1. AA+A=AAA^+A = A
    2. A+AA+=A+A^+AA^+ = A^+
    3. AA+=(AA+)TAA^+ = (AA^+)^T
    4. A+A=(A+A)TA^+A = (A^+A)^T

与逆矩阵、左伪逆矩阵、右伪逆矩阵一样,Moore-Penrose逆矩阵也是唯一的

逆矩阵和上面提到的两种伪逆矩阵都是 Moore-Penrose 逆矩阵的特例


求解方法

  • Moore-Penrose 逆矩阵 A+A^+ 具有的性质

    1. Moore-Penrose 逆矩阵 A+A^+ 是唯一的
    2. Moore-Penrose 逆矩阵的广义逆矩阵等于原矩阵:(A+)+=A(A^+)^+ = A
    3. D=diag(d11,,dnn)D = \operatorname{diag}(d_{11}, \cdots, d_{nn})n×nn \times n 对角矩阵,则 D+=diag(d11+,,dnn+)D^+ = \operatorname{diag}(d_{11}^+, \cdots, d_{nn}^+),其中 dii+=dii1d_{ii}^+ = d_{ii}^{-1}(若 dii0d_{ii} \neq 0)或者 dii+=0d_{ii}^+ = 0(若 dii=0d_{ii} = 0

    如此定义后DD满足Moore-Penrose条件

  • 利用 SVD 分解计算 Moore-Penrose 逆矩阵 A+A^+
    n×mn \times m 矩阵 AA 的奇异值分解为 A=UΣVTA = U \Sigma V^T,其中 UUn×nn \times n 的正交矩阵,VVm×mm \times m 的正交矩阵,Σ\Sigman×mn \times m 的对角矩阵,则 A+=VΣ+UTA^+ = V \Sigma^+ U^T

    可带入Moore-Penrose条件验证