电力系统潮流计算入门：从数学模型到MATLAB实战（附完整代码）

电力系统潮流计算：从核心原理到MATLAB工程实践

对于每一位踏入电力系统分析领域的研究者或工程师而言，潮流计算都是一个绕不开的基石。它不仅是理解电网稳态运行的“钥匙”，更是进行系统规划、安全评估和市场运营等高级分析的起点。然而，面对复杂的数学模型和迭代算法，许多初学者往往感到无从下手，停留在理论理解的层面，难以将知识转化为解决实际问题的能力。本文旨在打破这一壁垒，我们将抛开繁复的纯理论推导，聚焦于如何将抽象的数学模型转化为可运行的代码。我将以一个实践者的视角，带你深入理解潮流计算的核心骨架，并手把手地构建一个完整的MATLAB计算程序。我们的目标不仅是“看懂”，更是“亲手实现”，最终能够解读母线电压、线路潮流等计算结果背后的物理意义，为后续的科研与工程应用打下坚实基础。

1. 潮流计算的数学内核：模型构建与节点分类

要编写一个潮流计算程序，首要任务是清晰地定义我们所要解决的数学问题。这就像盖房子前必须先有精确的蓝图。潮流计算的核心，是求解一组描述电力网络功率平衡的非线性代数方程。

1.1 潮流方程的两种坐标系表达

电力网络中的每个节点（或称母线）的电气状态，通常由四个变量描述：有功注入功率 (P)、无功注入功率 (Q)、电压幅值 (V) 和电压相角 (\theta)。对于一个具有 (N) 个节点的系统，我们就有 (4N) 个变量。网络本身的物理约束——基尔霍夫电流定律，为我们提供了 (2N) 个实数方程（每个节点对应一个有功平衡方程和一个无功平衡方程）。因此，我们必须指定其中 (2N) 个变量，才能求解剩余的 (2N) 个变量。这些指定方式，就引出了节点的分类。

在建立方程时，我们有两种常用的坐标系选择：极坐标和直角坐标。它们各有优劣，适用于不同的算法。

极坐标形式 更符合我们的物理直觉。对于节点 (i)，其注入功率与全网节点电压的关系为： [ P_i = V_i \sum_{j=1}^{N} V_j (G_{ij}\cos\theta_{ij} + B_{ij}\sin\theta_{ij}) ] [ Q_i = V_i \sum_{j=1}^{N} V_j (G_{ij}\sin\theta_{ij} - B_{ij}\cos\theta_{ij}) ] 其中，(\theta_{ij} = \theta_i - \theta_j)，(G_{ij} + jB_{ij}) 是节点导纳矩阵 (Y_{bus}) 的第 (i) 行第 (j) 列元素。这种形式下，待求变量是电压幅值 (V) 和相角 (\theta)，方程数量与PQ节点数相关。

直角坐标形式 则将电压表示为实部和虚部：(V_i = e_i + jf_i)。相应的潮流方程为： [ P_i = e_i \sum_{j=1}^{N}(G_{ij}e_j - B_{ij}f_j) + f_i \sum_{j=1}^{N}(G_{ij}f_j + B_{ij}e_j) ] [ Q_i = f_i \sum_{j=1}^{N}(G_{ij}e_j - B_{ij}f_j) - e_i \sum_{j=1}^{N}(G_{ij}f_j + B_{ij}e_j) ] 直角坐标形式在处理PV节点时，方程形式相对统一，但变量数目加倍。

提示：牛顿-拉夫逊法在极坐标和直角坐标下均可实现。极坐标形式方程数少，雅可比矩阵计算略复杂；直角坐标形式方程数多，但雅可比矩阵元素公式更规整。初学者可从极坐标形式入手，其物理意义更清晰。

1.2 节点类型的工程意义与处理

如何指定那 (2N) 个已知变量？这取决于节点的物理特性，由此产生了三种经典节点类型：

在程序实现中，我们需要用数组明确记录每个节点的类型和对应的已知量、未知量。例如，可以定义一个结构体数组 bus，每个元素包含 type, P_sch, Q_sch, V_sch, theta_sch 等字段。平衡节点的处理是关键，它在迭代过程中不参与有功和电压幅值的修正方程，但其功率在迭代收敛后由线路潮流总和反推计算得出。

2. 算法核心：牛顿-拉夫逊法详解与实现

在众多潮流算法中，牛顿-拉夫逊法因其二阶收敛速度和良好的收敛可靠性，成为实际工程和学术研究中最主流的方法。它的核心思想是局部线性化：在某个近似解附近，用泰勒展开忽略高阶项，将非线性方程组转化为线性方程组进行迭代求解。

2.1 算法原理与迭代流程</