In mathematics, specifically linear algebra, the Woodbury matrix identity – named after Max A. Woodbury – says that the inverse of a rank-k correction of some matrix can be computed by doing a rank-k correction to the inverse of the original matrix. Alternative names for this formula are the matrix inversion lemma, Sherman–Morrison–Woodbury formula or just Woodbury formula. However, the identity appeared in several papers before the Woodbury report. The Woodbury matrix identity is
( A + U C V ) − 1 = A − 1 − A − 1 U ( C − 1 + V A − 1 U ) − 1 V A − 1 , {\displaystyle \left(A+UCV\right)^{-1}=A^{-1}-A^{-1}U\left(C^{-1}+VA^{-1}U\right)^{-1}VA^{-1},}
where A, U, C and V are conformable matrices: A is n×n, C is k×k, U is n×k, and V is k×n. This can be derived using blockwise matrix inversion. While the identity is primarily used on matrices, it holds in a general ring or in an Ab-category. The Woodbury matrix identity allows cheap computation of inverses and solutions to linear equations. However, little is known about the numerical stability of the formula. There are no published results concerning its error bounds. Anecdotal evidence suggests that it may diverge even for seemingly benign examples (when both the original and modified matrices are well-conditioned).
Discussion To prove this result, we will start by proving a simpler one. Replacing A and C with the identity matrix I, we obtain another identity which is a bit simpler:
( I + U V ) − 1 = I − U ( I + V U ) − 1 V . {\displaystyle \left(I+UV\right)^{-1}=I-U\left(I+VU\right)^{-1}V.}
To recover the original equation from this reduced identity, replace U {\displaystyle U} by A − 1 U {\displaystyle A^{-1}U} and V {\displaystyle V} by C V {\displaystyle CV} . This identity itself can be viewed as the combination of two simpler identities. We obtain the first identity from
I = ( I + P ) − 1 ( I + P ) = ( I + P ) − 1 + ( I + P ) − 1 P , {\displaystyle I=(I+P)^{-1}(I+P)=(I+P)^{-1}+(I+P)^{-1}P,}
thus,
( I + P ) − 1 = I − ( I + P ) − 1 P , {\displaystyle (I+P)^{-1}=I-(I+P)^{-1}P,}
and similarly
( I + P ) − 1 = I − P ( I + P ) − 1 . {\displaystyle (I+P)^{-1}=I-P(I+P)^{-1}.}
The second identity is the so-called push-through identity
( I + U V ) − 1 U = U ( I + V U ) − 1 {\displaystyle (I+UV)^{-1}U=U(I+VU)^{-1}}
that we obtain from
U ( I + V U ) = ( I + U V ) U {\displaystyle U(I+VU)=(I+UV)U}
after multiplying by ( I + V U ) − 1 {\displaystyle (I+VU)^{-1}} on the right and by ( I + U V ) − 1 {\displaystyle (I+UV)^{-1}} on the left. Putting all together,
… excerpt ends here. Continue reading the full article.
