3. 节点导纳矩阵与网络方程构建
在电力系统稳态分析中，节点导纳矩阵 Ybus 是把网络拓扑、支路参数和变压器变比统一“压缩”成代数方程的核心数据结构。它描述的是：当网络中某一个节点注入单位电压、其余节点电压为零时，各节点注入电流的分布关系。换句话说，Ybus 的每一列对应一个节点的“激励—响应”关系，其对角元反映该节点与所有相连支路的自导纳之和，非对角元反映节点之间的互导纳，符号与连接关系直接相关。
从物理意义上看，Ybus 的对角元 Yii 等于连接到节点 i 的所有支路导纳之和（含对地支路），非对角元 Yij 等于节点 i 与 j 之间支路导纳的负值。若两节点之间没有直接支路，则对应非对角元为零。因此，Ybus 是一个对称、稀疏、对角占优的复数矩阵，这一结构特征直接决定了后续潮流计算和短路计算的求解效率。
支路追加法：从零开始构建 Ybus
实际工程中，网络规模大、支路多，直接按定义逐项填写容易出错。支路追加法是一种“逐条支路累加”的构建方式：先初始化一个全零的 n×n 矩阵，然后依次读入每条支路，按其两端节点编号和导纳值修改对应元素。对于一条连接节点 i 与 j、导纳为 yij 的普通支路，修改规则为：
Y[i][i] += y_ij
Y[j][j] += y_ij
Y[i][j] -= y_ij
Y[j][i] -= y_ij若支路还带有对地导纳 yi0 或 yj0，则只在对应对角元上追加。例如节点 i 侧有对地支路，则 Y[i][i] += y_i0。这种“逐条追加”的方式与网络实际接线顺序一致，便于程序实现，也便于人工核对。
以一个三节点系统为例：节点1—2之间支路阻抗 z12 = 0.02 + j0.04（标幺值），节点2—3之间支路阻抗 z23 = 0.01 + j0.03，节点1—3之间支路阻抗 z13 = 0.02 + j0.05，节点2还带有对地电容导纳 y20 = j0.02。按支路追加法，先计算各支路导纳：y12 = 1/(0.02+j0.04) = 10 − j20，y23 = 1/(0.01+j0.03) = 10 − j30，y13 = 1/(0.02+j0.05) ≈ 6.897 − j17.241。随后逐条累加，最终得到 3×3 的 Ybus。这一过程在程序中通常用稀疏矩阵存储，只记录非零元素及其行列索引。
变压器非标准变比处理
变压器是网络中唯一会“改变电压基准”的元件。当变压器变比不等于其两侧节点基准电压之比时，称为非标准变比。处理非标准变比的标准做法是采用 π 型等值电路：将变压器等效为串联导纳 yT 与两侧对地导纳的组合，其中非标准变比 k 以理想变压器形式串入。
设变压器阻抗归算到二次侧后的导纳为 yT，非标准变比 k = 实际变比 / 基准变比。则 π 型等值电路的三条支路导纳为：
串联支路：y_T / k
一次侧对地：y_T · (1 - k) / k²
二次侧对地：y_T · (k - 1) / k在 Ybus 中追加变压器支路时，不能简单套用普通支路公式，而应按上述三条支路分别累加。若变压器连接节点 p（一次侧）和 q（二次侧），则修改规则为：
Y[p][p] += y_T / k + y_T·(1-k)/k²
Y[q][q] += y_T / k + y_T·(k-1)/k
Y[p][q] -= y_T / k
Y[q][p] -= y_T / k注意此时 Ybus 不再对称，因为非标准变比破坏了互导纳的对称性。这是含变压器网络的重要特征，也是程序实现中必须单独处理的分支类型。
网络方程与节点电压方程的关系
网络方程的基本形式是节点电流注入方程：I = Ybus · U，其中 I 是节点注入电流列向量，U 是节点电压列向量。该方程由基尔霍夫电流定律（KCL）直接导出：每个节点的注入电流等于流出该节点所有支路电流之和。展开后，第 i 个节点的方程为：
I_i = Σ_j Y_ij · U_j在潮流计算中，节点注入电流并不直接已知，已知的是节点注入功率 Si = Pi + jQi。因此需要将电流方程转化为功率方程：Si = Ui · (Σj Yij Uj)*。这就是节点电压方程（潮流方程）的来源。可见，Ybus 是连接“网络结构”与“节点功率平衡”的桥梁：先由支路和变压器参数构建 Ybus，再代入节点功率约束，才能形成可迭代求解的非线性方程组。
从上述流程可以看出，Ybus 的构建不是简单的“填表”，而是对网络拓扑、支路参数和变压器变比的系统性编码。普通支路按“对角加、非对角减”的规则累加；变压器支路则必须先转换为 π 型等值电路，再按三条支路分别处理；对地支路只影响对角元。最终得到的 Ybus 代入节点电压方程，就形成了潮流计算所需的非线性方程组。理解这一构建过程，是后续牛顿-拉夫逊法、PQ 分解法能够正确实现的前提。