跳到正文
格致开物MATHWIKI

广义逆

内容维护:格致开物 · 编写与核验 · 修订记录 · 纠错

广义逆(generalized inverse)将普通矩阵逆推广到非方阵和奇异矩阵。最基本的一种约定是:给定 A∈ℝm×n,若 G∈ℝn×m 满足 AGA=A, 就称 G 为 A 的一个广义逆,也称 {1}-逆或内逆。这个条件要求:对于确实来自 A 的输出,经过 G 再经 A 后,输出保持不变。不同资料可能附加其他条件,使用“广义逆”时应说明采用哪组条件。摩尔—彭若斯伪逆是由四个条件唯一选定的特殊广义逆。[1]

本文着重说明基本条件 AGA=A 的非唯一性,以及它怎样表达线性方程的相容性和通解。

一个矩阵可以有许多个广义逆

取 A=(1000),设 G=(abcd)。两次相乘得到 AGA=(a000). 所以广义逆条件只要求 a=1,其余三个参数任意。可见基本广义逆通常不唯一,也不一定满足 AG=GA=I。

若右端为 b0=(2,0)𝖳,方程 Ax=b0 的全部解是 (2,t)𝖳。广义逆给出的 Gb0=(2,2c)𝖳 是其中一个解,参数 c 决定它在不可识别方向上走多远;基本条件并未要求选择长度最小的那个。若右端改成 (2,1)𝖳,原方程无解,任意 G 也不能凭空补出缺失的第二输出坐标。

相容性与通解公式

定理 固定一个满足 AGA=A 的 G。方程 Ax=b 有解,当且仅当 AGb=b。相容时全部解为 x=Gb+(I−GA)z,z∈ℝn.

先证相容性。若 b=Ax0,则 AGb=AGAx0=Ax0=b。反过来,若 AGb=b,直接取 x=Gb 就满足方程。

对通解公式,注意 A(I−GA)=A−AGA=0. 所以任意 z 都给出 Ax=AGb=b。反向若 x 是任一解,在公式中取 z=x,右端成为 Gb+x−GAx=Gb+x−Gb=x。因此没有漏掉任何解。公式把一个特解与核空间中的自由部分分开,而不是要求每个参数 z 都对应不同的解。

平面点二一被不同投影送到横轴上的三零和二零,两个投影均保持横轴的每个点。
同一个输出子空间可以有不同的投影方向;基本广义逆没有规定正交性。

图采用 G=(1100)。此时 AG 把 (u,v) 送到 (u+v,0),所以图中 (2,1) 被送到 (3,0);取伪逆 A+=A 时则送到 (2,0)。两个映射都保持横轴不动,但前者沿斜方向投影,后者沿垂直方向投影。图里的三零不是原来不相容方程的精确解,而是经过 AG 改写后的右端。

基本广义逆总存在

设 A 的秩为 r。通过选取输入、输出空间的基,可将它写成 A=P(Ir000)Q, 其中 P,Q 分别是可逆的 m×m 和 n×n 矩阵,中间块为 m×n。令 J 为这个中间矩阵,令 J− 为尺寸反向、左上块同为 Ir、其余为零的 n×m 矩阵。块乘法给出 JJ−J=J。于是 G=Q−1J−P−1 满足 AGA=PJJ−JQ=A,证明每个有限矩阵都至少有一个基本广义逆。这个构造中的基可以多种方式选择,因此也与非唯一性一致。

幂等投影与伪逆的附加条件

由 AGA=A 得到 (AG)2=AGAG=AG,(GA)2=GAGA=GA. 因此 AG 与 GA 均为幂等投影。AG 的像空间恰为 A 的列空间:它的像显然包含于列空间,而对每个 y=Ax,AGy=AGAx=Ax=y,故列空间中的点都被保留。它不一定是正交投影,图中的斜投影就是反例。

摩尔—彭若斯伪逆除第一条件外,还满足 GAG=G,(AG)𝖳=AG,(GA)𝖳=GA. 后两条将投影限制为正交投影,四条合在一起唯一确定 A+。复数情形使用共轭转置。最小二乘与最小范数性质依赖这些额外结构,不能仅从 AGA=A 推出。[1]

变式取 A=(1  1)。其广义逆为所有列向量 G=(t,1−t)𝖳,因为 AGA=(t+1−t)A=A。对右端二,得到特解 (2t,2−2t);只有 t=1/2 给出最小欧几里得范数解 (1,1)。若 A 本来是可逆方阵,左右乘普通逆则强迫 G=A−1,非唯一性随之消失。

参考文献

  1. ↑ 1.0 1.1 Roger Penrose,A generalized inverse for matrices,Mathematical Proceedings of the Cambridge Philosophical Society,51(3),1955,406—413,原文扫描,定理1。