5、广义逆矩阵
逆,生来就是用于解方程组的。
-
逆:行列满秩。
-
单边逆:左逆列满秩;右逆行满秩。
-
求法:高斯消元。
-
列满秩矩阵:Ax=b
- 行初等变换,可以得到左逆
A
L
−
1
A_L^{-1}
- 有解的充要条件:
A
A
L
−
1
b
=
b
AA_L^{-1}b=b
- 唯一解:
x
=
(
A
H
A
)
−
1
A
H
b
=
A
L
−
1
b
Large x = (A^HA)^{-1}A^Hb=A_L^{-1}b
- 需要注意,左逆矩阵不唯一,
(
A
H
A
)
−
1
A
H
(A^HA)^{-1}A^H
- 行初等变换,可以得到左逆
-
行满秩矩阵:
- 列初等变换,可以得到右逆
A
R
−
1
A_R^{-1}
- 行满秩矩阵一定有解,且解不唯一。自由未知数的个数为n-m
-
A
H
(
A
A
H
)
−
1
A^H(AA^H)^{-1}
- 需要注意:
A
H
(
A
A
H
)
−
1
b
≠
A
R
−
1
b
large A^H(AA^H)^{-1}b {ne} A_R^{-1}b
- 行满秩矩阵的A+是A的一个右逆。
- 列初等变换,可以得到右逆
-
-
广义逆:任何矩阵都存在广义逆矩阵。
-
AGA=A 等价于 G是A的广义逆矩阵。
-
单边逆是广义逆的特殊情况。
-
Ax = b有解的充要条件:rank A = rank (A b)
-
当A列满秩,且rank A = rank (A b) ,则有唯一解。
-
当A非列满秩,且有解,则有无穷多解。
-
Gb = x , 且Ax = b,则称G是A的一个广义逆。
-
广义逆矩阵不唯一,零矩阵的广义逆矩阵是任意矩阵。
-
广义逆矩阵通常记作
A
−
A^-
A−,区别于
A
−
1
A^{-1}
A−1
-
广义逆矩阵的秩 大于等于 A的秩。取等号,等价于 G具有自反性。
-
广义逆矩阵是逆矩阵的推广,当A是可逆矩阵的时候,A–=A-1
-
广义逆矩阵有逆矩阵类似的性质:
性质 广义逆矩阵 逆矩阵 定义 对于矩阵A,存在矩阵B,使得ABA=A 对于方阵A,存在矩阵B,使得AB=BA=I 记号 A
−
A^-
A
−
1
A^{-1}
行为 对于任意向量b,满足 A
A
−
AA^-
对于任意向量b,满足 A
A
−
1
AA^{-1}
矩阵乘法 A
A
−
AA^-
AA −
1
^{-1}
矩阵的秩 rank( A
−
A
A^-A
A
A
−
AA^-
rank(AA −
1
^{-1}
逆的存在 广义逆存在于可逆和不可逆矩阵中 逆存在于方阵(可逆矩阵)中 唯一性 可以有多个不同的广义逆 逆是唯一的 幂等性 A
A
−
AA^-
A
−
A
A^-A
A
A
−
1
=
A
−
1
A
=
I
AA^{-1}=A^{-1}A=I
数乘 aA的广义逆为 1
a
A
−
,
a
≠
0
frac{1}{a} A^-,ane0
aA的逆为 1
a
A
−
1
,
a
≠
0
frac{1}{a} A^{-1},ane0
正交投影 若 (
A
A
−
)
H
=
A
A
−
(AA^-)^H=AA^-
A
A
−
=
P
R
(
A
)
AA^-=P_{R(A)}
A
A
−
1
=
I
AA^{-1}=I
-
-
自反广义逆
- 定义:AGA=A且GAG=G。存在且不唯一。
- 求法:
- 最大秩分解法:A=BD,
A
r
−
=
D
R
−
1
B
L
−
1
A_r^-=D_R^{-1}B_L^{-1}
- 构造法:X、Y是A的广义逆,则Z = XAY是自反广义逆,当然YAX也是自反广义逆。
- 公式法:X=(AHA)–AH,Y =AH(AAH)–都是自反广义逆
- R(AH)=R(AHA), N(A)=N(AHA)
- 存在D,使得AH=AHAD
- 代入可证明AXA=A
- rank(X)<=rank(AH)=rank(AHA) = rank(AHA(AHA)–AHA)=rank(AHAXA)<=rank(X)
- 所以X是自反广义逆
- 最大秩分解法:A=BD,
- A_是自反广义逆的充要条件是rank(A)=rank(A–)
- 必要性:AGA=A且GAG=G,则rank(A)=rank(AGA)<=rank(G)=rank(GAG)<=rank(A)
- 充分性:
- R(GA)属于R(G),R(GA)=R(A)=R(G),则R(G)=R(GA);
- GE=G,则存在X,GAX=G
- A=AGA=AGAXA=AXA
- 由构造法知G是自反广义逆
- 几何性质
-
R
(
A
)
⊕
N
(
A
H
)
=
C
m
R(A)oplus N(A^H)=C^m
-
R
(
A
H
)
⊕
N
(
A
)
=
C
n
R(A^H)oplus N(A)=C^n
-
R
(
A
)
⊕
N
(
A
r
−
)
=
C
m
R(A)oplus N(A^-_r)=C^m
-
R
(
A
r
−
)
⊕
N
(
A
)
=
C
n
R(A^-_r)oplus N(A)=C^n
-
R
(
A
r
−
)
=
R
(
A
H
)
R(A^-_r) = R(A^H)
-
N
(
A
r
−
)
=
N
(
A
H
)
N(A^-_r) = N(A^H)
-
-
MP广义逆
- 定义:AGA=A,GAG=G,(GA)H=GA, (AG)H=AG
- 存在且唯一
- 计算方法:
- 最大值分解法:
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
+
=
V
H
D
+
U
A^+=V^HD^+U
- 注意:最大值分解不唯一,然而最大值分解的这种乘积,即A+是唯一的。
- 最大值分解法:
- A+的性质:
- 自反性:(A+)+=A
- 唯一性:
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
- 几何性质:
R
(
A
+
)
=
R
(
A
H
)
R(A^+) = R(A^H)
- 正交投影性:
A
A
+
=
P
R
(
A
)
,
A
+
A
=
P
R
(
A
H
)
AA^+=P_{R(A)},A^+A=P_{R(A^H)}
- 行列子空间相等性:R(A)=R(AH)的充要条件是
A
+
A
=
A
+
A
A^+A=A^+A
-
(
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)^+
-
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是Hermite矩阵:
-
(
A
2
)
+
=
(
A
+
)
2
(A^2)^+=(A^+)^2
-
-
矩阵方程通解:
- AXB=D有解的充要条件:AA–DB–B=D
- 通解:X=A–DB–+Y-A–AYBB–
- Ax=b有解的充要条件:AA–b=b
- 通解:x = A–b +y-A-1Ay
-
相容方程组的最小范数解:
- 相容方程组:有解方程组
- AGA=A,(GA)H=GA
- Gb=x是最小范数解
-
不相容方程组的最小二乘解:
- 不相容方程组:无解方程组
- AGA=A,(AG)H=AG
- Gb=x是最小二乘解之一,最小二乘解的通解:
x
=
G
b
+
(
E
−
G
A
)
b
x = Gb+(E-GA)b
-
不相容方程组的最小二乘解有时候不够好,最小二乘解的最小范数解叫最佳逼近解:
- 最佳逼近解:A+b=x
- 如果方程组是相容方程组,则A+b是最小范数解
-
关于广义逆的运算规则
-
(
A
−
)
H
=
(
A
H
)
−
(A^-)^H=(A^H)^-
-
B
=
S
A
T
,
B
−
=
T
−
1
A
−
S
−
1
B=SAT,B^-=T^{-1}A^-S^{-1}
-
(
A
B
)
+
=
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)
-
(
A
B
)
+
=
B
+
A
+
(AB)^+=B^+A^+
-
原文地址:https://blog.csdn.net/weixin_45827703/article/details/135280176
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如若转载,请注明出处:http://www.7code.cn/show_52120.html
如若内容造成侵权/违法违规/事实不符,请联系代码007邮箱:suwngjj01@126.com进行投诉反馈,一经查实,立即删除!