1. 程式人生 > >024矩陣求逆引理

024矩陣求逆引理

  記得以前在學序貫平差的時候就用過矩陣求逆引理,但是當時只是死記,當然早就忘了。在此作個筆記記錄一下引理的推導過程。 (A+BCD)1=A1+X \begin{aligned} (A+BCD)^{-1} &= A^{-1} + X \end{aligned}

移項 \Downarrow \text{移項}

(A+BCD)(A1+X)=E \begin{aligned} (A+BCD) (A^{-1} + X) &= E \end{aligned}

展開 \Downarrow \text{展開}

E+AX+BCDA1+BCDX=E \begin{aligned} E + AX + BCDA^{-1} + BCDX &= E \end{aligned}

求解 \Downarrow \text{求解}

X=(A+BCD)1BCDA1=[BC(C1B1A+D)]1BCDA1=(C1B1A+D)1(BC)1BCDA1=(C1B1A+D)1DA1=[(C1+DA1B)B1A]1DA1=A1B(C1+DA1B)1DA1 \begin{aligned} X &= -(A+BCD)^{-1}BCDA^{-1} \\ &= -[BC(C^{-1}B^{-1}A+D)]^{-1}BCDA^{-1} \\ &= -(C^{-1}B^{-1}A+D)^{-1} (BC)^{-1} BCDA^{-1} \\ &= -(C^{-1}B^{-1}A+D)^{-1} DA^{-1} \\ &= -[(C^{-1}+DA^{-1}B)B^{-1}A]^{-1} DA^{-1} \\ &= -A^{-1}B(C^{-1}+DA^{-1}B)^{-1} DA^{-1} \\ \end{aligned}

代入 \Downarrow \text{代入}

(A+BCD)1=A1A1B(C1+DA1B)1DA1 \begin{aligned} (A+BCD)^{-1} &= A^{-1} - A^{-1}B(C^{-1}+DA^{-1}B)^{-1} DA^{-1} \end{aligned}

  記住這種表示方式:(A+BCD)1=A1+X\begin{aligned} (A+BCD)^{-1} &= A^{-1} + X \end{aligned},這樣只需求解XX便可推匯出結論。   提取相同矩陣時,為了便於理解也可以一個一個往外提取。   另外可以這樣記憶第二項:

A1BA^{-1}B (C1C^{-1} +DA1B)1DA1+DA^{-1}B)^{-1} DA^{-1}

A1BA^{-1}B (C1+D(C^{-1}+D A1BA^{-1}B )1DA1)^{-1} DA^{-1}

A1B(C1+A^{-1}B (C^{-1} + DA1DA^{-1} B)1B)^{-1} DA1DA^{-1}

  如果對類似的形式推導求逆引理,例如(ABC1D)1(A-BC^{-1}D)^{-1},可以類比得到:

(ABC1D)1=A1+A1B(CDA1B)1DA1 (A-BC^{-1}D)^{-1} = A^{-1} + A^{-1}B(C-DA^{-1}B)^{-1} DA^{-1}

  用同樣的方法也可求得:

(A+BCT)1=A1A1B(I+CTA1B)1CTA1 (A+BC^T)^{-1} = A^{-1} - A^{-1}B(I+C^TA^{-1}B)^{-1} C^TA^{-1}