矩阵求逆引理(Sherman-Morrison-Woodbury公式)详解
1. 前言在电力系统分析、控制理论、数值计算中我们经常会遇到形如( I A B ) − 1 (IAB)^{-1}(IAB)−1或A ( I B A ) − 1 A(IBA)^{-1}A(IBA)−1的矩阵运算尤其是在节点阻抗矩阵推导、系统等效建模、状态估计等场景中。矩阵求逆引理Sherman-Morrison-Woodbury公式简称SMW公式是解决这类问题的核心工具它可以将高维矩阵求逆转化为低维矩阵求逆大幅降低计算复杂度同时也为电力系统中等效阻抗的推导提供了理论支撑。本文将从基础定义、核心恒等式、详细推导、电力系统应用实例、数值实现等方面带你彻底掌握矩阵求逆引理。2. 矩阵求逆引理的核心形式2.1 通用形式Sherman-Morrison-Woodbury公式设A ∈ R n × n A \in \mathbb{R}^{n \times n}A∈Rn×nU ∈ R n × k U \in \mathbb{R}^{n \times k}U∈Rn×kC ∈ R k × k C \in \mathbb{R}^{k \times k}C∈Rk×kV ∈ R k × n V \in \mathbb{R}^{k \times n}V∈Rk×n且A AA和C − 1 V A − 1 U C^{-1}VA^{-1}UC−1VA−1U均可逆则( A U C V ) − 1 A − 1 − A − 1 U ( C − 1 V A − 1 U ) − 1 V A − 1 (AUCV)^{-1}A^{-1}-A^{-1}U(C^{-1}VA^{-1}U)^{-1}VA^{-1}(AUCV)−1A−1−A−1U(C−1VA−1U)−1VA−1这是最通用的形式适用于任意矩阵维度的情况。2.2 常用特例单位矩阵形式在电力系统和控制理论中最常用的是A I AIAI单位矩阵的情况此时公式简化为( I U C V ) − 1 I − U ( C − 1 V U ) − 1 V (IUCV)^{-1}I-U(C^{-1}VU)^{-1}V(IUCV)−1I−U(C−1VU)−1V进一步当C I CICI时得到Sherman-Morrison-Woodbury公式的简化形式( I U V ) − 1 I − U ( I V U ) − 1 V (IUV)^{-1}I-U(IVU)^{-1}V(IUV)−1I−U(IVU)−1V其中U ∈ R n × m U \in \mathbb{R}^{n \times m}U∈Rn×mV ∈ R m × n V \in \mathbb{R}^{m \times n}V∈Rm×n。此外还有一个与矩阵求逆引理密切相关的矩阵恒等式常称为push-through恒等式( I A B ) − 1 A A ( I B A ) − 1 (IAB)^{-1}AA(IBA)^{-1}(IAB)−1AA(IBA)−1其中A ∈ R n × m A \in \mathbb{R}^{n \times m}A∈Rn×mB ∈ R m × n B \in \mathbb{R}^{m \times n}B∈Rm×n且I A B IABIAB和I B A IBAIBA均可逆。下面我们给出该恒等式从SMW公式出发的详细推导。从SMW公式到push-through恒等式的推导已知SMW简化形式取C I CICI( I U V ) − 1 I − U ( I V U ) − 1 V (IUV)^{-1}I-U(IVU)^{-1}V(IUV)−1I−U(IVU)−1V设U A , V B UA,\ VBUA,VB代入得( I A B ) − 1 I − A ( I B A ) − 1 B (IAB)^{-1}I-A(IBA)^{-1}B(IAB)−1I−A(IBA)−1B将等式两边同时右乘矩阵A AA( I A B ) − 1 A [ I − A ( I B A ) − 1 B ] A (IAB)^{-1}A\big[I-A(IBA)^{-1}B\big]A(IAB)−1A[I−A(IBA)−1B]A右侧展开( I A B ) − 1 A I A − A ( I B A ) − 1 B A (IAB)^{-1}AIA-A(IBA)^{-1}BA(IAB)−1AIA−A(IBA)−1BA利用单位矩阵性质I A A IAAIAA整理得( I A B ) − 1 A A − A ( I B A ) − 1 B A (IAB)^{-1}AA-A(IBA)^{-1}BA(IAB)−1AA−A(IBA)−1BA提取公因子A AA( I A B ) − 1 A A [ I − ( I B A ) − 1 B A ] (IAB)^{-1}AA\Big[I-(IBA)^{-1}BA\Big](IAB)−1AA[I−(IBA)−1BA]对括号内做恒等变形凑出( I B A ) − 1 (IBA)^{-1}(IBA)−1I − ( I B A ) − 1 B A ( I B A ) − 1 ( I B A ) − ( I B A ) − 1 B A ( I B A ) − 1 ( I B A − B A ) ( I B A ) − 1 I ( I B A ) − 1 \begin{aligned} I-(IBA)^{-1}BA (IBA)^{-1}(IBA)-(IBA)^{-1}BA \\ (IBA)^{-1}\big(IBA-BA\big) \\ (IBA)^{-1}I \\ (IBA)^{-1} \end{aligned}I−(IBA)−1BA(IBA)−1(IBA)−(IBA)−1BA(IBA)−1(IBA−BA)(IBA)−1I(IBA)−1代回上式最终得到( I A B ) − 1 A A ( I B A ) − 1 \boldsymbol{(IAB)^{-1}AA(IBA)^{-1}}(IAB)−1AA(IBA)−1至此由Sherman-Morrison-Woodbury公式推导出push-through恒等式。3. 总结矩阵求逆引理Sherman-Morrison-Woodbury公式是电力系统分析、控制理论和数值计算中的重要工具它的核心价值在于实现了( I A B ) − 1 A (IAB)^{-1}A(IAB)−1A与A ( I B A ) − 1 A(IBA)^{-1}A(IBA)−1的等价转换大幅降低高维矩阵求逆的计算复杂度尤其是设备数远小于网络节点数的场景为电力系统中等效阻抗推导、控制器设计提供了简洁的理论支撑。