基本定义
1 矩阵的左右逆
从广义的角度讲,任何一个矩阵 G 都可以称为矩阵 A 的逆矩阵,若它与矩阵 A 的乘积等于单位矩阵 I,即满足 GA=I。根据矩阵 A 本身的特点,满足这一定义的矩阵 G 存在以下三种可能的情况:
- G 存在,并且唯一;
- G 存在,但不唯一;
- G 不存在。
当rankA<minn,m时不存在左右逆,但可以有广义逆
例:
A1=211−210−1−2−1,A2=45−28−73,A3=[123511]
A1 有唯一逆矩阵
A2 存在多个 2×3 矩阵 L 使得 LA2=I2×2(有左逆)
例如
A2=45−28−73,L=[76821700],L=[003275],⋯
对矩阵 A3,没有任何 3×2 矩阵使得 G3A3=I3×3,但存在多个 3×2 矩阵 R,使得 A3R=I2×2(有右逆)
例如
A3=[123511],R=1−1310−1,R=−10210−1,⋯
1.1 左右伪逆矩阵
满足 LA=I,但不满足 AL=I 的矩阵 L 称为矩阵 A 的左逆矩阵
满足 AR=I,但不满足 RA=I 的矩阵 R 称为矩阵 A 的右逆矩阵
一个矩阵的左逆矩阵或右逆矩阵往往非唯一,下面考虑它们的一个特殊解:
- 考查 n>m 且 A 具有满列秩 (rankA=m) 的情况。此时,m×m 矩阵 ATA 是可逆的,容易验证 L=(ATA)−1AT 满足 LA=I,这个左逆矩阵是唯一确定的,常称为 左伪逆矩阵。
设 A 为 n×m 矩阵,且列满秩。则存在左逆矩阵 L 满足 LA=Im。一个特解为 L0=(ATA)−1AT。所有左逆可表示为 L=L0+Z,其中 Z 是任意 m×n 矩阵且满足 ZA=0。Z 有 m×n 个未知数,但 ZA=0 只给出了 m×m 个方程。于是每个行向量有 n−m 个自由参数,m 行共有 m(n−m)≥1 个自由参数。自由参数可以取无穷多个值,所以 Z 有无穷多个解。
- 考查 n<m 且 A 具有满行秩 (rankA=n) 的情况。此时,n×n 矩阵 AAT 是可逆的,容易验证 R=AT(AAT)−1 满足 AR=I,这个右逆矩阵是唯一确定的,常称为 右伪逆矩阵。
2 矩阵的广义逆
考虑一个 n×m 维的秩亏缺矩阵 A(即rank(A)<min{m,n})。n×m 维的秩亏缺矩阵的逆矩阵称为 广义逆矩阵。令 A− 表示 A 的广义逆矩阵。
无论 A−A=Im×m 还是 AA−=In×n 都不可能成立。有必要使用三个矩阵的乘积定义一个秩亏缺矩阵 A 的逆矩阵。
考虑线性矩阵方程 Ax=y 的求解。两边左乘 AA−,则有 AA−Ax=AA−y。若 A− 是 A 的广义逆矩阵,则 Ax=y⇒x=A−y,代入 AA−Ax=AA−y 可得 AA−Ax=Ax,我们希望该式对任意非零向量 x 均应成立,则必须要求下列约束条件满足:AA−A=A。满足 AA−A=A 的矩阵 A− 称为 A 的 广义逆矩阵。
2.1 Moore-Penrose逆矩阵
满足 AA−1A=A 的矩阵 A−1 不是唯一的,存在明显的缺陷。只能保证 A−1 是 A 的广义逆矩阵,并不能保证 A 是 A−1 的广义逆矩阵
考虑原矩阵方程 Ax=y 的解方程 x=A−1y 的求解。两边左乘 A−1A,则有 A−1Ax=A−1AA−1y。由于矩阵 A 是 A−1 的广义逆矩阵,则 x=A−1y⇒Ax=y,代入 Ax=y 可得 A−1y=A−1AA−1y,因此
A−1AA−1=A−1也需要满足
我们还希望当 A 是列满秩或者行满秩时,广义逆矩阵 A−1 能够包括左伪逆矩阵和右伪逆矩阵在内的特例:左伪逆矩阵 L=(ATA)−1AT 满足 LA=I 但不存在 AL=I,但 AL=A(ATA)−1AT=(AL)T,同样,右伪逆矩阵满足 RA=(RA)T。所以对称性也要满足
- Moore-Penrose逆矩阵
令 A 是任意 n×m 矩阵,称 A+ 是 A 的广义逆矩阵,若 A+ 满足以下四个条件(常称 Moore-Penrose 条件)则称为Moore-Penrose逆矩阵
- AA+A=A
- A+AA+=A+
- AA+=(AA+)T
- A+A=(A+A)T
与逆矩阵、左伪逆矩阵、右伪逆矩阵一样,Moore-Penrose逆矩阵也是唯一的
逆矩阵和上面提到的两种伪逆矩阵都是 Moore-Penrose 逆矩阵的特例
求解方法
-
Moore-Penrose 逆矩阵 A+ 具有的性质
- Moore-Penrose 逆矩阵 A+ 是唯一的
- Moore-Penrose 逆矩阵的广义逆矩阵等于原矩阵:(A+)+=A
- 若 D=diag(d11,⋯,dnn) 为 n×n 对角矩阵,则 D+=diag(d11+,⋯,dnn+),其中 dii+=dii−1(若 dii=0)或者 dii+=0(若 dii=0)
如此定义后D满足Moore-Penrose条件
-
利用 SVD 分解计算 Moore-Penrose 逆矩阵 A+
若 n×m 矩阵 A 的奇异值分解为 A=UΣVT,其中 U 是 n×n 的正交矩阵,V 是 m×m 的正交矩阵,Σ 是 n×m 的对角矩阵,则 A+=VΣ+UT
可带入Moore-Penrose条件验证