记得以前在学序贯平差的时候就用过矩阵求逆引理,但是当时只是死记,当然早就忘了。在此作个笔记记录一下引理的推导过程。
(A+BCD)−1=A−1+X
⇓移项
(A+BCD)(A−1+X)=E
⇓展开
E+AX+BCDA−1+BCDX=E
⇓求解
X=−(A+BCD)−1BCDA−1=−[BC(C−1B−1A+D)]−1BCDA−1=−(C−1B−1A+D)−1(BC)−1BCDA−1=−(C−1B−1A+D)−1DA−1=−[(C−1+DA−1B)B−1A]−1DA−1=−A−1B(C−1+DA−1B)−1DA−1
⇓代入
(A+BCD)−1=A−1−A−1B(C−1+DA−1B)−1DA−1
记住这种表示方式:
(A+BCD)−1=A−1+X,这样只需求解
X便可推导出结论。
提取相同矩阵时,为了便于理解也可以一个一个往外提取。
另外可以这样记忆第二项:
A−1B (
C−1
+DA−1B)−1DA−1
扫描二维码关注公众号,回复:
3476486 查看本文章
A−1B
(C−1+D
A−1B
)−1DA−1
A−1B(C−1+
DA−1
B)−1
DA−1
如果对类似的形式推导求逆引理,例如
(A−BC−1D)−1,可以类比得到:
(A−BC−1D)−1=A−1+A−1B(C−DA−1B)−1DA−1
用同样的方法也可求得:
(A+BCT)−1=A−1−A−1B(I+CTA−1B)−1CTA−1
参考:
矩阵求逆引理(matrix inversion lemma)