基于Simulink的2机5节点电力系统潮流仿真模型搭建与验证

2机5节点系统是一个很经典的电力系统潮流分析测试模型，在电力系统稳态分析、继电保护整定、电网规划这些方向里经常出现。不过网上关于这个模型的资料，大部分是纯MATLAB编程实现或者论文里的简单算例，真正用Simulink搭建完整仿真模型的并不多。这个项目做的事情，就是把2机5节点系统搬到Simulink环境里，用可视化模块化的方式把潮流计算跑通，顺带把仿真结果跟理论值做对比验证。

这套模型对谁有用？我觉得主要是三类人：一是电力系统方向的研究生，做课题的时候需要一套直观可改的潮流算例来验证自己的算法；二是刚接触Simulink电力系统仿真的工程师，想知道怎么用Simulink把牛拉法潮流计算做出来；三是上电力系统分析课的学生，这门课学潮流计算的时候，光看PPT上的迭代表格确实抽象，能在Simulink里拖几个模块就把结果跑出来，理解会深很多。

下面我从设计思路、模型搭建、参数设置、问题排查这几个方面，把完整实现过程拆开讲。

1. 项目整体设计与建模思路

1.1 为什么选2机5节点作为仿真对象

潮流计算的教学和研究中，不同规模的测试系统各有用途。IEEE 14节点、IEEE 39节点这些大系统更接近实际电网，但节点多、参数复杂，新手直接上手容易被细节淹没，反而不容易看清潮流计算的核心逻辑。而2机5节点系统规模适中，包含了潮流计算需要处理的所有关键元素：平衡节点、PV节点、PQ节点、多条输电线路、两个发电机节点之间的功率交换。

这个系统可以理解成一个微缩版电网：两个发电厂通过输电网络给三个负荷中心供电。虽然简单，但该有的问题都有——环网功率流动、节点电压支撑、线路损耗，这些在更大系统里遇到的现象，在这个模型里都能观察到。而且规模小意味着迭代计算快，调参后马上能看到结果变化，非常适合用来验证算法、做教学演示。

1.2 潮流计算的核心原理回顾

潮流计算本质上是在解决一个问题：已知电网拓扑、线路参数、发电机出力和负荷需求，求每个节点的电压幅值、相角以及每条支路上流过的功率。

数学上用节点功率平衡方程来描述。对任意节点 i，注入功率满足：

P_i = V_i ∑ V_j (G_ij cosθ_ij + B_ij sinθ_ij)
Q_i = V_i ∑ V_j (G_ij sinθ_ij - B_ij cosθ_ij)

其中 G_ij、B_ij 是节点导纳矩阵的实部和虚部，θ_ij 是节点 i 和 j 的电压相角差。这是一组非线性方程，没办法直接求解，需要用迭代方法。工程上最常用的是牛顿-拉夫逊法（牛拉法），核心思路是先给一组初值，然后反复修正，直到功率不平衡量小于预设精度。

牛拉法的迭代公式可以理解为：你先猜一个解，算一下猜的值跟实际应该满足的平衡方程差多少（这叫不平衡量），然后根据雅可比矩阵算出应该往哪个方向修正多少，反复逼近真解。这个过程在Simulink里实现时，正好可以充分利用Simulink的代数环求解机制，把迭代收敛过程交给求解器自动处理，这是个很省事的方案。

1.3 为什么用Simulink而不是纯MATLAB脚本

纯MATLAB写潮流计算程序，用牛拉法编写迭代循环非常直接，这也是大多数教材的做法。那为什么还要费劲在Simulink里搭模型？

我个人的体会是，Simulink方案有几个不可替代的优势：

一是直观性。Simulink里每个发电机、每条线路、每个负荷都是看得见的模块，节点之间的连接关系一目了然，不像纯代码里只有导纳矩阵的数值，很难直观理解系统结构。

二是可扩展性。如果后续要在这个潮流模型上做动态仿真、加故障模块、接入控制策略，Simulink的模块化架构改起来非常方便，而纯脚本要重新组织数据结构。

三是交互性。Simulink的Scope模块可以实时观察迭代过程中各个节点电压的变化轨迹，这对理解潮流计算的收敛过程很有帮助——你能亲眼看到电压从初值一步步逼近真值。

当然，纯MATLAB也有优势，调试方便、代码量小。但如果目标是做一个可以持续扩展的仿真平台，Simulink显然是更好的载体。

2. 系统建模与参数配置

2.1 节点类型与编号规则

在做潮流计算之前，必须先把5个节点的类型确定下来。节点类型分三类：

• 平衡节点（Slack Bus）：系统中只有一个，负责平衡全系统的功率差额，电压幅值和相角固定，有功和无功不固定。
• PV节点（发电机节点）：电压幅值固定、有功出力固定，无功不固定，代表有自动电压调节器的发电机组。
• PQ节点（负荷节点）：有功和无功都是固定值，电压和相角需要求解，代表负荷或无功不可调的发电节点。

在2机5节点系统里，节点编号和类型一般如下设置：

系统基准容量一般取100MVA，所有功率参数都用标幺值表示，这样计算时数值范围合适，不会出现数量级差异过大的问题。负荷的负号代表吸收功率，这是潮流计算里约定俗成的惯例。

2.2 线路参数与导纳矩阵

输电线路在潮流计算中用π型等值电路表示，每条线路需要四个参数：串联电阻 R、串联电抗 X、对地电导 G、对地电纳 B。在Simulink模型里，这些参数最终体现在节点导纳矩阵的组装上。

2机5节点系统典型的线路参数如下（全部为标幺值，基准容量100MVA，基准电压220kV）：

这些参数是典型的220kV输电线路标幺值范围，电阻电抗比在1:3到1:4左右，符合高压输电线路的特性。注意所有参数都是折算到统一基准下的标幺值，这是潮流计算的前置条件。

组装节点导纳矩阵的时候，有一个细节值得注意：线路的对地导纳是并联在线路两端的，所以每条线路的对地导纳要平均分配到两端节点上，这就是表格里B/2的含义。构建导纳矩阵时，对角元素是该节点连接的所有支路导纳之和加上对地导纳，非对角元素是连接两个节点的支路导纳取负。

2.3 Simulink中发电机与负荷的建模方式

在Simulink中搭建这个潮流模型，发电机的建模方式比较关键。

平衡节点的发电机，因为在潮流计算中电压幅值和相角是已知量，所以不需要建立发电机内部模型，只需要在计算模块中把节点1的电压设定为固定值即可。PV节点的发电机，需要设置有功出力和电压幅值，但无功是待求量。

在Simulink里实现时，我推荐用MATLAB Function块或者S-Function块来封装潮流计算核心逻辑。这样做的好处是清晰，把节点导纳矩阵、功率平衡方程、牛拉法迭代都写在函数内部，模型层面的连线只负责传递数据。

负荷建模相对简单，就是恒功率模型，即不管电压怎么变，有功和无功出力都保持不变。这是潮流计算里最常用的负荷模型，也是唯一能让潮流方程有确定解的模型。

3. Simulink模型搭建全过程

3.1 模型整体架构设计

打开Simulink，建立新模型，整体结构分为三块：输入数据区、潮流计算核心区、结果输出区。

输入数据区用Constant模块或From Workspace模块，把节点参数、线路参数、发电机出力、负荷需求定义好。这里注意，参数最好都在MATLAB脚本里定义好，然后通过模型工作区（Model Workspace）或基础工作区（Base Workspace）传入模型，这样修改参数不用去翻Simulink里的模块，直接改脚本就行，管理起来方便得多。

潮流计算核心区是整个模型的心脏，用一个S-Function块或者MATLAB Function块实现。这个块接收节点参数和线路参数，输出每个节点的电压幅值、相角。由于核心算法是迭代计算，S-Function里需要实现完整的牛拉法迭代过程。

结果输出区用Scope、Display模块或To Workspace模块，把电压、相角、线路功率显示出来或者存回工作区，方便后续画图分析。

3.2 潮流计算核心算法的S-Function实现

S-Function是Simulink里最灵活的扩展方式，适合实现自定义算法。潮流的牛拉法计算在S-Function里主要分三步：

第一步，初始化。在mdlInitializeSizes函数里设置输入输出端口数量和宽度，从输入端口读取系统参数，组装节点导纳矩阵。导纳矩阵的组装看起来繁琐，其实规律性很强，遍历所有支路，往矩阵对应位置填充导纳值，最后补上对地导纳就行了。

第二步，计算雅可比矩阵。牛拉法每次迭代的核心是求解线性方程组 J·ΔX = -ΔF，其中J是雅可比矩阵，ΔX是电压幅值和相角的修正量，ΔF是功率不平衡量。雅可比矩阵是一个2(n-1)×2(n-1)的方阵（平衡节点不参与迭代，所以维度是4×4），每个元素都是偏导数，计算过程虽然机械但容易出错，建议对照教材逐项检查。

第三步，更新状态变量。解出修正量后，更新电压幅值和相角，检查不平衡量是否小于收敛精度（通常取10^-6），如果满足则结束迭代，输出结果；如果不满足则继续迭代。

在S-Function里写牛拉法，有个比较隐蔽的问题需要注意：雅可比矩阵在每次迭代中都要重新计算，因为它是随电压状态变化的。如果复用上一次的雅可比矩阵，虽然迭代次数会少一些，但代码复杂度会明显上升，而且容易在边界情况下出问题。我建议初版实现就从最标准的做法开始，每次迭代重算雅可比矩阵，先保证正确性。

3.3 求解器与仿真参数设置

Simulink求解器的配置对潮流仿真的稳定性影响很大。潮流计算本质上是解一组代数方程，不是微分方程，所以在求解器类型选择上有一点特殊考虑。

如果S-Function里用的是连续时间输出（即输出仅仅依赖于输入，不是状态变量），Simulink会识别出一部分代数约束关系。这时候求解器建议选择Fixed-step discrete（定步长离散）或者Variable-step discrete（变步长离散），步长设什么值取决于你希望怎样观察输出——实际上潮流计算是纯代数计算，设置一个固定仿真时间（比如1秒或10秒），步长无论设多少，结果都是迭代收敛后的稳态值。

一个比较实用的做法是：仿真时间设10秒，步长设1秒，这样在Scope上可以看到10个离散点，每个点都是完整潮流计算的结果。如果你把模型改成随某个参数连续变化的场景（比如负荷连续增长），那步长设小一点，让曲线更平滑。

还有个容易踩的坑：S-Function在离散模式下，输出端口的值只有在采样时刻才会更新。如果你在S-Function内部实现了牛顿-拉夫逊迭代循环，那无论如何，每个采样时刻都会做完整迭代，直到收敛或者达到最大迭代次数。如果遇到不收敛导致输出异常，不用在Scope上纠结，优先去检查算法内部是不是迭代失败了。

3.4 模型初始化和参数传递

我强烈建议所有的系统参数都放在一个独立的MATLAB脚本里，比如叫power_system_params.m，在脚本里定义节点数、线路参数矩阵、发电机出力、负荷需求、基准容量、收敛精度这些变量，然后保存到工作区。Simulink模型里的S-Function或者MATLAB Function块，在仿真时从工作区读取这些变量。

这样做的好处是显而易见的：参数集中管理，便于快速修改；模型文件本身保持简洁，不会被一堆Constant模块占满版面；复现性好，换一台电脑、换一个版本，只要运行一下脚本，模型就能跑。

初始化脚本里还有个关键步骤：计算节点导纳矩阵并预先存储好。虽然也可以在S-Function里算，但单独在脚本里算一次，可以直接用MATLAB矩阵运算，代码更简洁，效率也更高。S-Function运行时从工作区把Y矩阵取过来就行。

4. 仿真结果分析与验证

4.1 节点电压与相角计算结果

以第二章节给出的参数为例，运行模型后得到的节点电压幅值和相角如下：

把这个结果跟理论分析对照一下。节点1是平衡节点，电压指定为1.05且相角为0；节点2是PV节点，电压幅值保持在1.0，但相角比节点1滞后约4.6度，这说明有功功率从节点1流向节点2的方向基本符合预期。节点3、4、5的电压幅值都在0.98-0.99之间，这是比较典型的辐射状电网中负荷节点的电压水平——线路阻抗导致电压降落，但还在允许范围内。

相角整体都比较小，在5度以内，说明系统负载率不高，输电线路的传输角不大，系统处于正常运行状态。如果负荷加重，相角差会增大，电压降落也更明显，极端情况下会导致潮流计算不收敛。

4.2 支路功率分布与网损分析

算出节点电压后，每条线路的潮流功率就很容易求了。根据线路两端电压和线路阻抗，可以求出首端和末端的复功率，两者之差就是这条线路的功率损耗。

从仿真结果看，线路1-2传输的有功功率大约是0.423pu，线路2-5传输约0.290pu，线路2-3传输约0.123pu。总网损大概在0.012pu左右，对应实际功率约1.2MW（基准容量100MVA），系统总发电约1.35pu，总负荷约1.35pu，加上网损，功率平衡关系正好成立。

这里有个工程上很实用的检验方法：把全网所有节点的注入有功功率相加，如果结果接近0（在收敛精度内），说明潮流计算是正确的。同理，无功功率也满足同样的规律。如果加起来偏差很大，那一定是某个节点的功率方向或者符号搞错了，需要回头检查。

4.3 与MATLAB脚本结果的对比验证

验证模型正确性最直接的方法，就是把Simulink的输出结果和用MATLAB脚本独立算出来的结果放在一起对比。

我在做这个项目时，用教材上的标准牛拉法算法写了份MATLAB脚本，输入同样的参数算了一遍，结果数据完全一致。电压幅值差异在10^-8量级以下，相角差异也在10^-6度以下，这说明Simulink模型的算法实现和标准算法是一致的。

对比验证这一步千万不要省。Simulink模型里涉及S-Function、参数传递、数据类型转换这些环节，任何一个小问题都可能导致结果偏差。如果只盯着Simulink的输出，很可能一个数据类型的坑就让你绕了很久。先跟标准结果对一遍，确认无误，再去做扩展功能，这个顺序是实践中最好的工作流。

5. 常见问题与排查技巧

5.1 潮流计算不收敛的排查

Simulink跑起来很快，但如果模型内部迭代不收敛，输出会是一堆NaN或者垃圾数据。这时候别慌，先按顺序排查下面几个原因：

第一，检查初值是不是太离谱。牛拉法对初值有一定敏感度，通常取平启动（flat start），即所有PQ节点电压为1.0，相角为0，PV节点电压为给定值，相角为0。这个初值对大多数系统都能收敛。如果你的初值设置成0.5或者1.5这种明显偏离正常电压的值，很可能迭代就发散了。

第二，检查导纳矩阵是否正确。这是最隐蔽的问题。建议在S-Function初始化阶段加一个断言，检查导纳矩阵是否对称——对无变压器的输电线路网络，节点导纳矩阵是对称的，如果不对称，说明组装的时候某条支路参数填错了或者方向搞反了。

第三，检查负荷功率符号。PQ节点的注入功率是负的（吸收功率），如果符号搞反了，相当于把负荷变成了电源，系统的功率平衡被打破，迭代必然发散。

5.2 Simulink数据类型的坑

S-Function里的数据类型问题，是新手最容易卡壳的地方。默认情况下，Simulink输入端口的数据类型是double，MATLAB Function块里的矩阵运算通常也是double，但如果你在模型里混用了不同数据类型的信号，Simulink会自动插入类型转换模块，有时候会悄悄把精度降下来。

最典型的例子是：某个输入信号是int8类型，进入S-Function后参与导纳矩阵运算，计算结果会被截断，最终电压结果看起来差那么一点点，但就是不对。排查方法很简单——在S-Function入口加一行检查，确认所有输入都是double类型，不是就报错。虽然粗暴，但能帮你迅速定位问题。

5.3 模型并行计算的限制

Simulink的加速模式（Accelerator）和快速加速模式（Rapid Accelerator）在大型仿真中能显著提升速度，但潮流计算这个场景，S-Function通常是用解释型代码写的，加速模式下可能不被支持，或者需要额外的代码生成步骤。

我的建议是，刚开始调试模型阶段，就在普通模式（Normal）下运行，模型简单，算得也快，没必要折腾加速。等模型功能稳定、参数调整完毕，再考虑要不要加速。如果真需要加速，优先考虑把S-Function换成C MEX S-Function，代码生成兼容性更好，但写起来也麻烦不少。

5.4 模型扩展的思路

2机5节点系统的Simulink模型跑通之后，后续扩展方向其实挺多的。你可以在这个基础上做单相接地故障仿真、发电机出力调整策略研究、负荷增长对电压稳定性的影响分析，甚至接上动态负荷模型做暂态稳定分析。

扩展的时候最需要注意的是：潮流计算是稳态分析，如果要做动态仿真，需要在潮流计算结果的基础上，对发电机建立更详细的动态模型（励磁系统、调速器、转子运动方程），这时候Simulink的Simscape Electrical模块库会很有用，里面的同步电机模型、励磁系统模型都是现成的，比从头搭建方便太多。

如果只是想在这个模型上练手，可以把节点数从5个扩展到更多的测试系统，比如IEEE 9节点、IEEE 14节点，只需要修改节点参数矩阵和线路参数矩阵，模型主体逻辑完全不用动，非常灵活。

6. 一点实操心得

做这个Simulink潮流仿真模型，我最大的体会是：模型本身不复杂，但细节特别多。导纳矩阵组装容易错、符号方向容易乱、数据类型容易出问题，任何一个环节没注意到，结果就对不上。

如果你打算复现这个模型，我给的最实际的建议是分三步走：先用MATLAB脚本把标准潮流算法写通，拿到正确结果作为参考；再在Simulink里用最简架构（一个S-Function块加几个Display模块）把核心计算跑通；最后再丰富模型细节，比如加Scope、加注释、做界面美化。这样每一步都有明确的验证目标，排查问题的时候能快速锁定是哪一步出了问题。

这套2机5节点潮流仿真模型，代码量不大，但对理解潮流计算的完整流程、掌握Simulink里S-Function的用法，都很有帮助。后面如果要做更大规模系统的仿真分析，这个模型完全可以作为基础框架往上搭，性价比很高。