元件的组合降阶——节点消去

前面几篇文章里，我们已经把 R、L、C 支路一步步装配到了节点电压方程里。

如果把 RLC 串联支路完全展开，它会有外部节点$a,b$，也会有内部节点$c,d$。完整写法最清楚，因为每个元件的电压、电流和历史项都能看见；但它也会带来更多未知节点。对一个大型网络来说，如果每个组件内部节点都全部暴露给全局求解器，矩阵规模会很快变大。

所以工程程序里常常会做一件事：只保留组件对外可见的端口节点，把组件内部节点消去。

这件事在电网络里常叫 Kron reduction；在线性代数里，它就是 Schur complement；在有限元和数值计算里，也常被称为 static condensation。名字不同，本质都是同一件事：把一部分未知量从方程里代数消去，只留下外部关心的变量。

本文先不讨论大规模稀疏矩阵排序、代码生成或实时优化，只用我们熟悉的 RLC 串联支路说明一件最核心的事：节点消去会让外部矩阵变小，但也会把内部节点电压从主方程里拿掉。如果后面还想知道内部电压，就必须用恢复公式重新算回来。

从展开的 RLC 支路开始

仍然使用前文的 RLC 串联支路。节点从左到右为$a,c,d,b$，其中$a,b$是外部端口节点，$c,d$是内部节点。

三条支路电压分别是：

用 Dommel 形式写成节点矩阵时，支路方程可以整理为：

其中：

这里$G_R=1/R$，$G_L$和$G_C$分别来自电感、电容的 Dommel 等效导纳。这个矩阵是完整的：它既能求外部端口$a,b$，也能求内部节点$c,d$。

但如果这个 RLC 支路只是一个组件，对外只暴露$a,b$两个端口，那么全局网络求解器并不一定需要直接把$c,d$也作为全局未知量。我们可以把$c,d$在组件内部消去。

保留节点和内部节点

先从图上看这件事。展开时，R、L、C 各自保留自己的连接关系，中间有$c,d$两个内部节点；组合降阶以后，对外只剩$a,b$两个端口，内部结构被折算成一个等效导纳和一个历史电流源，也就是诺顿形式。

为了写成标准消去形式，先把节点按两类重新分组：

• retained nodes：保留节点，也就是对外端口节点$a,b$；
• internal nodes：内部节点，也就是组件内部的$c,d$。

于是：

这里$r$表示 retained，$k$表示 internal。下标用$k$只是为了和常见 Schur/Kron 记号一致；它不是时间步下标。

还要注意一个细节：前一节完整矩阵的节点顺序是$a,c,d,b$，这是沿着电路从左到右写的；但做消去时，更方便把保留节点放在前面、把要消去的内部节点放在后面，也就是把顺序重排成：

重排以后，原来的$G$矩阵会变成下面这个样子。虚线左上角是保留节点之间的$G_{rr}$，右下角是准备消去的内部节点块$G_{kk}$。

所以，分块公式不是凭空冒出来的，而是对“重排后的原始$G$矩阵”按虚线切开。按$\mathbf V_r,\mathbf V_k$分块，可以写成：

对这个 RLC 例子，各个分块就是：

历史项也随节点顺序分成两块：

这样写的好处是，哪些量要保留，哪些量准备消去，一眼就清楚了。

Schur 消去的标准形式

内部节点$c,d$没有额外外部电流注入，所以内部方程可以写成：

代入分块方程的第二行：

如果$G_{kk}$可逆，就可以解出内部节点：

为了写得短一点，记：

于是：

再把这个式子代回 retained 方程：

得到：

其中：

这就是 Schur complement / Kron reduction 的标准形式。消去以后，外部求解器只看见$\mathbf V_r=[u_a,u_b]^T$，也就是一个二端口支路。

代回 RLC，会得到什么

对 RLC 这个小例子，$G_{kk}$只是一个$2\times2$矩阵。它的行列式为：

也就是：

所以：

代入 Schur 公式，最终的 reduced matrix 会变成一个普通二端支路的形式：

其中：

这个结果很符合直觉。三个串联元件展开时有内部节点；消去内部节点以后，对外就只剩一个等效二端口。

历史项也会被组合到端口上，形式变成：

其中：

于是外部节点方程就可以写成：

这就是“元件组合降阶”的具体样子：原来 R、L、C 三个元件展开后有$a,c,d,b$四个节点；消去$c,d$以后，对外只剩$a,b$两个节点和一个等效二端支路。

消去不是免费午餐

节点消去的好处很明显：外部矩阵变小了。

原来这个支路如果完全展开，要把$c,d$作为未知节点一起放进全局方程；消去以后，全局方程只需要处理$a,b$。如果一个组件内部有几十、几百个内部节点，这种降阶对全局求解规模会很有帮助。

但代价也很明确：被消去的节点电压不再是主方程直接求出来的未知量。

在 RLC 例子里，消去后外部方程只能直接得到：

它不会直接给出：

这意味着，如果你只保存 reduced equation，那么下面这些量也不会自动出现：

• 电阻电压$u_R=u_a-u_c$；
• 电感电压$u_L=u_c-u_d$；
• 电容电压$u_C=u_d-u_b$；
• 内部节点监测量；
• 依赖内部节点电压的控制、保护或测量逻辑。

所以节点消去不是“什么信息都没少，只是矩阵变小了”。更准确地说，它是把内部节点从主求解器里拿掉了。外部端口行为可以保持等效，但内部变量需要额外恢复。

内部电压恢复公式

恢复公式其实前面已经推出来了：

对 RLC 例子来说，就是：

写成工程流程就是：

• 第一步，用 reduced equation 求出保留节点$u_a,u_b$；
• 第二步，用上一拍状态算出$I_{L,\mathrm{his}}$和$I_{C,\mathrm{his}}$；
• 第三步，用恢复公式算出内部节点$u_c,u_d$；
• 第四步，再由节点电压差回算$u_R,u_L,u_C$和支路电流。

注意顺序很重要。恢复内部节点时，不能凭空猜$u_c,u_d$，也不能只看 reduced matrix。恢复公式需要用到原来参与消去的那部分矩阵，也就是$G_{kk}$、$G_{kr}$和内部历史项$\mathbf I_{\mathrm{his},k}$。

这也是很多程序里容易出错的地方：主方程已经降阶了，但监测量、保护逻辑或内部状态更新仍然需要内部电压。如果只做了消去，没有配套恢复，就会出现“端口电流对了，但内部电压不知道从哪里来”的问题。

小结

节点消去做的是一件非常朴素的事：把内部节点方程解出来，再代回外部节点方程。

标准形式是：

对 RLC 串联支路来说，这一步会把展开的 R、L、C 三个元件组合成一个二端口等效支路。外部求解器只需要看$a,b$两个端口，矩阵规模更小。

但这个降阶会损失内部节点的直接可见性。$u_c,u_d$被消去以后，如果后续还要元件电压、内部监测量或状态更新，就必须用：

把内部节点电压恢复出来。

所以，节点消去和电压恢复最好放在一起理解：消去负责让外部矩阵变小，恢复负责把被拿掉的内部信息找回来。