Simulink实现2机5节点系统潮流仿真模型-代码聚汇网

Simulink实现2机5节点系统潮流仿真模型

胖厨胡学斌·2026-08-29·阅读 21 GEO检测（免费）

Simulink潮流计算2机5节点

|难度：普通|阅读约 35 分钟

做仿真这些年，最常被学生和刚入行的同事问到的，就是“有没有一个简单又完整的系统模型可以让我先跑通流程”。今天分享的这个 2机5节点系统潮流仿真模型（Simulink仿真实现） ，正好是这类需求里性价比最高的一套。它不像几十节点的电网模型那样复杂到无从下手，也比单机无穷大系统更能反映多电源、多负荷之间的功率分配关系。用Simulink搭一套能跑通的潮流模型，你基本就摸清了电力系统稳态分析从建模到求解的全流程。

这个模型能做什么？一句话：在给定发电机出力、负荷大小和网络参数后，计算出每个节点的电压幅值、相角，以及各条支路流过的有功和无功功率。它既是电力系统分析课程里的必修内容，也是你做暂态稳定、短路计算、无功优化等后续研究时的“前置条件”。适合正在学《电力系统分析》的本科生、准备做毕业设计的研究生，以及刚接触电力系统仿真、想快速上手Simulink的工程师。

1. 项目概述与仿真目标拆解

1.1 什么是“2机5节点系统”

所谓2机5节点系统，指这个仿真模型里有2台发电机（电源）、5条母线（节点），以及连接这些母线的变压器和输电线路。节点之间通过线路和变压器构成一个小型输电网，负荷接在部分母线上。

节点编号是有讲究的。我在建模时一般把平衡节点放在1号母线，两台发电机中一台做平衡机，另一台做PV节点，其余带负荷的母线作为PQ节点。这种节点类型划分不是随便定的，潮流计算里每个节点的已知量不同，直接决定了方程能不能解出来。

拿我常用的一个典型算例来说：

节点编号	节点类型	已知量	待求量
1	平衡节点（Slack）	电压幅值、相角	有功、无功出力
2	PV节点	有功出力、电压幅值	电压相角、无功出力
3	PQ节点	有功负荷、无功负荷	电压幅值、相角
4	PQ节点	有功负荷、无功负荷	电压幅值、相角
5	PQ节点	有功负荷、无功负荷	电压幅值、相角

平衡节点是全系统的“锚”。因为潮流方程本质是功率平衡方程，它只约束各节点功率注入和流出的关系，如果不固定某一点的相角，整个系统的相角就会整体漂移，雅可比矩阵奇异，根本解不出来。工程上常把容量最大、带调频能力的电厂设为平衡节点，我这个算例里也是如此。

1.2 为什么选Simulink而不是其他工具

做潮流仿真可选的工具其实很多，比如PSASP、BPA、PSS/E这些专业电力系统分析软件，还有纯手写MATLAB脚本实现牛顿-拉夫逊迭代。那为什么我优先推荐Simulink？

答案在于“仿真链路”。潮流计算在Simulink里不是孤立的，它最大的价值是给后续的电磁暂态仿真提供初始运行点。用Simulink搭建2机5节点系统，你先通过潮流计算得到各节点电压和功率初值，然后直接在这个模型里扔一个三相短路故障，切换到时域仿真，就能看暂态稳定的机间功角摇摆曲线。这套流程在PSASP/BPA里是“两个独立软件、两套数据”，而在Simulink里是“一套模型、连续工作”。

另外，Simulink的图形化建模方式对新手极其友好。你可以像画电气一次接线图一样，把三相电压源、变压器、线路、负荷这些模块拖到画布上连线，双击设置参数就能跑。这对理解网络拓扑和功率流向很有帮助，不像手写代码那样需要先建立节点导纳矩阵的抽象思维。

1.3 你需要的软件环境

做这个模型，我建议用MATLAB R2020a及以上的版本，主要原因是新版本的Simscape Electrical（也就是原来的SimPowerSystems）模块库对潮流求解支持更完善。如果你用的是更老的版本，比如R2016a甚至R2014a，也能搭，但部分模块名称和位置对不上，需要自己摸索。

注意一个关键点：Simulink自带的模块库里只有数学运算和信号处理这些基础模块，电力系统仿真必须依赖Simscape Electrical工具箱。你在MATLAB命令窗口输入 powerlib，如果能弹出模块库窗口，说明工具箱已经装好了；如果提示没有这个命令，那需要重新安装MATLAB并勾选Simscape Electrical组件。

需要模型API调用？免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。

2. 潮流计算原理与Simulink建模选型

2.1 潮流计算到底在算什么

说句实话，很多人在做这个模型的时候，把“搭模型”和“算潮流”搞成了两件事，这是误区。实际上，Simulink里的潮流计算是嵌入在仿真初始化过程中的，核心还是在解一组非线性方程。

每个节点的功率平衡方程为：

code复制Pi = Vi Σ Vj (Gij cosθij + Bij sinθij)
Qi = Vi Σ Vj (Gij sinθij - Bij cosθij)

其中Pi、Qi是节点注入功率（发电机出力减负荷），Vi、Vj是节点电压幅值，Gij、Bij是节点导纳矩阵的实部和虚部，θij是两节点电压相角差。这组方程里，Gij和Bij由网络拓扑和线路参数决定，Pi、Qi由节点类型决定已知量还是未知量。整个潮流计算就是在迭代解这个方程组。

工程上最常用的是牛顿-拉夫逊法。它的原理不复杂：先猜一组电压初值，代入方程算功率不平衡量，再求解修正方程得到电压修正量，反复迭代直到不平衡量小于收敛精度。在Simulink的powergui模块里选择Load Flow工具，它内部用的就是类似的方法。

在Simulink的Simscape Electrical建模中，潮流计算的产物不是一堆数字报表，而是各节点电压和相角的“初值”。它会自动把这些初值写入各个模块的内部状态，这样你切换到暂态仿真模式时，系统是从稳态运行点开始变化，而不是从零开始突然跳变。

2.2 节点类型与Simulink模块的对应关系

前面提到节点类型有平衡节点、PV节点和PQ节点之分，在Simulink里这三类节点对应着不同的模块参数设置方式。

平衡节点用Three-Phase Source模块实现，设置时要指定电压幅值（Voltage magnitude）和相角（Phase angle），频率选50Hz。这里要注意，模块内部还有一个Internal impedance参数，一般设成很小的值代表理想电压源。三相电压源的连接方式也要选对，我习惯选Yg（中性点接地），这样和后续变压器建模能统一。

PV节点同样用Three-Phase Source，但设置方式完全不同。在powergui的Load Flow工具中，你需要把该节点所在的电源指定为PV节点，并设定有功功率P和电压幅值V，此时模块的相角是未知的，由潮流计算自动生成。

PQ节点则不直接连接电源，而是接Three-Phase Series RLC Load模块作为负荷。这个模块直接输入有功功率P、无功功率Q和额定电压Vn，潮流计算会解出该节点的实际电压幅值和相角。需要留意的是，负荷模块的额定电压要和系统基准电压一致，否则潮流结果会受影响。

2.3 输电线路与变压器的建模选型

输电线路在5节点系统里一般用PI型等值电路模块Three-Phase PI Section Line。这个模块需要输入线路的正序电阻R、正序电感L、正序电容C和线路长度。对220kV架空线，每公里的典型参数大概是R=0.04欧、L=1.26mH、C=0.009uF。潮流仿真关心的是工频稳态特性，PI型等值模型在这个场景下精度足够，而且参数直观、调试方便。

变压器用Three-Phase Transformer (Two Windings)模块。这个模块设置参数时有很多坑，常见的问题是搞不清一次侧和二次侧的连接方式。我的原则是：高压侧用Yg连接，低压侧用Delta连接，也就是常见的Yd11接法。这样可以给系统提供接地点，同时阻断零序电流通路，和实际电网的习惯一致。

变压器阻抗参数一般用标幺值表示。比如一台额定容量300MVA、变比230/115kV的变压器，短路电压百分数Uk%=10%，那么在Simulink中短路阻抗Z=0.1pu，X/R比取10，则R=0.00995pu，X=0.0995pu。这些数值在变压器的对话框里直接填标幺值就能用。

2.4 标幺值系统与基准值设定

这是个很多人忽略但极其重要的环节。Simulink的Simscape Electrical模块里，很多参数既可以填国际单位也可以填标幺值，而潮流计算的结果又都是以标幺值形式存在的。如果基准值设得混乱，模型跑得出来结果也是错的。

我建议全系统统一采用以下基准值：

基准容量：100MVA
基准电压：230kV（对应220kV等级）
频率：50Hz

在设置负荷模块的额定电压时填的是kV实际值，在设置变压器时填的是标幺值阻抗，在潮流计算工具里看结果时读的是pu标幺值电压。这套混用习惯看起来不统一，其实是Simscape Electrical的工作方式：物理连接用国际单位，电气参数可以用pu值，潮流解算在内部会统一转换。

一个具体的设置实例：3号节点接了一个100MW、40Mvar的负荷，在Three-Phase Series RLC Load里填P=100e6、Q=40e6、Vn=230e3，注意Simulink里负荷有功无功单位是W和var，不是MW和Mvar，很多人在这里漏掉数量级，导致潮流结果偏得离谱。

3. 实操过程：从零搭建一套可运行的模型

3.1 系统拓扑与典型参数准备

动手搭模型之前，先把网络拓扑图和参数表准备好，这会省下大量排查时间。我下面给出一个典型的2机5节点算例参数，你可以直接照抄。

系统基准值：100MVA，230kV，50Hz。

线路参数（PI型等值，单位长度值）：

线路L12：长度120km
线路L13：长度80km
线路L23：长度100km
线路L24：长度150km
线路L25：长度100km
线路L45：长度120km

220kV线路每公里参数：

正序电阻：0.04Ω/km
正序电感：1.26mH/km
正序电容：0.009μF/km

节点负荷：

节点3：P=100MW，Q=40Mvar
节点4：P=80MW，Q=30Mvar
节点5：P=120MW，Q=50Mvar

发电机参数：

节点1（平衡节点）：额定电压230kV，频率50Hz，电压幅值设1.05pu，相角设0度
节点2（PV节点）：有功出力P=150MW，电压幅值设1.02pu

变压器参数：两台升压变压器，额定容量250MVA，变比230/115kV，短路电压Uk%=10%，X/R=10，连接方式Yd11。注意这里230kV对应发电机出口的115kV升压到220kV输电网络，如果你不想这么复杂，也可以把发电机直接接入230kV母线，省掉变压器，但那就少学了变压器建模这一环。

3.2 Simulink模型画布搭建步骤

打开MATLAB，在命令窗口输入 simulink，新建一个空白模型。依次从Simscape Electrical库中拖入以下模块：

powergui模块（必须放，否则模型无法运行）
两个Three-Phase Source（节点1和节点2的发电机）
两个Three-Phase Transformer (Two Windings)（升压变压器）
五段Three-Phase PI Section Line（对应5条线路）
三个Three-Phase Series RLC Load（3、4、5节点的负荷）
若干Three-Phase V-I Measurement（用于测量节点电压电流）

把模块排列好，按照节点编号连线。连线时有一个容易踩的坑：Simscape Electrical的电气线是粗线，和Simulink普通信号线的细线是两回事，两者不能直接相连。要测量电压和功率，必须通过Three-Phase V-I Measurement模块，把电气量转换成Simulink信号量。

连接完成后，双击powergui模块，在Configuration选项卡里确保仿真类型选的是Phasor simulation，频率设为50Hz。注意，这里Phasor simulation就是相量模式，是潮流计算的基础。如果你选Continuous模式，运行的是电磁暂态仿真，出来的不是稳态潮流结果。

3.3 发电机与负荷参数配置细节

双击节点1的Three-Phase Source，设置如下：

Phase-to-phase voltage (Vrms)：230e3
Phase angle of phase A：0
Frequency：50
Internal connection：Yg
在powergui的Load Flow里指定节点1为平衡节点，电压幅值输入1.05，相角0

节点2的Three-Phase Source类似，但Phase angle先随便填0，后续在Load Flow工具里会重新计算。在Load Flow工具里指定节点2为PV节点，输入P=150e6，V=1.02。

三个负荷模块分别设置：

节点3：P=100e6，Q=40e6，Vn=230e3
节点4：P=80e6，Q=30e6，Vn=230e3
节点5：P=120e6，Q=50e6，Vn=230e3

关于负荷模块的容性/感性设置，在Three-Phase Series RLC Load里如果Q是正值，表示感性负荷（吸收无功）；如果Q是负值，表示容性负荷（发出无功）。传统电力系统负荷大多是感性的，所以填正值没问题。

3.4 使用powergui的Load Flow工具求解潮流

这是整个建模流程里最核心的一步。双击powergui模块，在左侧工具栏中找到Load Flow选项，点击打开Load Flow工具窗口。

在Load Flow工具中，你会看到已经识别出来的所有Bus和机器。确认每个节点的类型设置正确后，点击Compute按钮，系统开始迭代求解潮流。如果参数都正确，几秒钟内就会收敛，显示各节点电压幅值（pu）和相角（度），以及各条支路的功率分布。

求解完成后，一定要点击Apply to model按钮，把潮流结果写入模型状态。这一步很多人会漏掉。如果没有Apply，你后续切到时域仿真模式时，系统还是从初始电压0开始，会出现剧烈的暂态振荡，这会直接导致后续工作出错。

我当时第一次做这个模型时，就是没点Apply to model，然后切到连续仿真模式看短路故障波形，结果电压波形一直在振荡，还以为是模型搭错了，排查了很久才发现是初值没写入。

3.5 查看仿真结果

潮流计算完成后，可以通过两种方式查看结果：

第一种，直接在Load Flow工具的表格里看数值结果。工具会列出每个节点的电压幅值、相角、有功功率和无功功率，以及每条线路的传输功率。这种查看方式适合快速排查数据异常。

第二种，连接Three-Phase V-I Measurement模块，把电压电流信号接到Scope或者MATLAB的To Workspace模块，在MATLAB工作区里分析结果。我用得比较多的是To Workspace方式，因为后续做N-1静态安全分析时，需要对多组故障场景批量处理结果，工作区的方式更适合脚本化操作。

下面是一段简单的MATLAB脚本，用于从工作区读取电压数据并绘制电压分布图：

matlab复制% 假设Scope数据已经存在变量ScopeData中
% 提取各节点A相电压幅值
t = ScopeData.time;
V_a1 = ScopeData.signals(1).values(:, 1);
V_a2 = ScopeData.signals(2).values(:, 1);
V_a3 = ScopeData.signals(3).values(:, 1);

figure;
plot(t, [V_a1 V_a2 V_a3], 'LineWidth', 1.5);
xlabel('Time (s)');
ylabel('Phase A Voltage (V)');
legend('Bus 1', 'Bus 2', 'Bus 3');
grid on;
title('Node Voltage after Load Flow Initialization');

注意，这段脚本假设你已经把To Workspace模块的变量名设成了ScopeData，并且信号排列顺序和连接顺序一致。实际使用时需要根据你自己的信号排列调整索引。

4. 常见问题与调试经验实录

4.1 仿真不收敛，问题出在哪

这是做2机5节点潮流最常碰到的问题，没有之一。Load Flow工具提示“Failed to converge”时，不要慌，按下面顺序排查。

首选检查负荷总量是否过大。潮流迭代不收敛的物理含义是，给定的发电机出力或平衡节点容量不足以满足负荷需求，或者网络阻抗实在太大、送不过去。你可以试着把三个节点的负荷等比缩小10%再算一次，如果收敛了，说明系统参数过于紧凑，需要调整发电机出力或网络参数。

其次检查线路参数的数量级是否正确。PI型线路模块中R、L、C的单位分别是欧姆、亨利、法拉，如果你把0.04Ω/km直接填成了0.04kΩ/km，相当于线路阻抗大了1000倍，功率根本送不出去。这个错误在参数从Excel表格复制到Simulink对话框时特别容易发生。

最后检查发电机的无功出力上限。在Load Flow工具中，PV节点需要设置无功出力范围（最小Q和最大Q）。如果PV节点电压目标值设定过高，导致需要吸收大量无功才能压低电压，但无功下限又不允许负太多，迭代就会发散。把PV节点的电压从1.02降到1.00，往往马上收敛。

4.2 潮流结果与手算对不上

很多时候Simulink能收敛，但结果和你用牛顿-拉夫逊手算或用其他软件算出来的数值对不上，差了千分之几甚至百分之几。这种情况大多是基准值不统一造成的。

比如你在手算时设基准容量为100MVA，但负荷模块里填的是100e6W，那么正好一致。如果负荷填成了100e6MW，等于要求系统输送10^12W的功率，Simulink照样能“收敛”出一个符合功率平衡的结果，只是节点电压会掉到非常低的水平，看起来好像也没报错。这个问题我见过太多次了。

还有一个容易被忽略的点：变压器的阻抗标幺值是基于变压器自身额定容量计算的，而不是系统基准容量。如果一台250MVA变压器的短路阻抗标幺值是0.1pu，换算到100MVA系统基准下应该是0.1*100/250=0.04pu。Simulink的变压器模块输入标幺值时，填写的是以自身额定容量为基准的值，这一点和潮流计算程序的要求不太一样，需要仔细看模块说明。

4.3 切换到时域仿真后电压波形发散

如果你在潮流计算结束后，把powergui从Phasor simulation切换回Continuous模式，运行仿真时发现电压波形不是平滑的正弦波，而是剧烈振荡或发散，那大概率是没点Apply to model。

另一种可能是，你虽然在Load Flow工具里点了Apply，但之后又改动了某个模块的参数，导致系统状态和潮流结果不再匹配。这时候只要重新运行一次Load Flow并再点一次Apply即可。

还有一种情况是求解器配置问题。在Simulink模型窗口的Simulation -> Model Configuration Parameters里，把Solver设成ode23tb，相对误差设成1e-4，最大步长不要超过1e-3秒。我实测下来，这个配置对电力系统电磁暂态仿真的稳定性和速度都比较友好。

4.4 常见问题速查表

现象	可能原因	解决方案
潮流不收敛	负荷过大、线路阻抗数量级错误、PV节点无功越限	检查负荷量级、核对R/L/C单位、放宽PV节点无功范围
电压普遍偏低	基准电压不一致、负荷功率数量级偏大	统一在230kV系统填Vn=230e3、检查P/Q是否漏乘1e6
潮流结果和手算不一致	变压器阻抗基准容量未折算	阻抗值乘以（系统基准/变压器容量）
切到时域波形发散	没有Apply to model	在Load Flow工具中点Apply to model
模块库找不到Three-Phase Source	Simscape Electrical未安装	安装Simscape Electrical工具箱
连线后全是红虚线	Y型连接和Delta连接不匹配	检查所有模块接地方式，确保网络中有接地点

4.5 调试效率提升技巧

调试Simulink潮流模型时，不要每次都从头点GUI操作。你可以用MATLAB脚本自动化整个流程：打开模型、设置参数、运行Load Flow、读取结果。这样批量测试不同负荷场景时效率能提升好几倍。

下面是一个自动化的示例脚本：

matlab复制% 打开模型
open_system('two_machine_5_bus');

% 修改负荷节点3的有功功率
set_param('two_machine_5_bus/Load3', 'ActivePower', '80e6');

% 运行Load Flow
power_loadflow('two_machine_5_bus', 'solved');

% 读取结果
busData = psbpowerflow('two_machine_5_bus');
disp(busData.Bus);

注意这个脚本里模型名和模块名需要和你实际模型完全一致，如果有中文或空格需要特殊处理。psbpowerflow函数在不同MATLAB版本里返回值结构可能略有差异，建议用 doc psbpowerflow 查看本机版本的具体说明。

5. 模型验证与后续扩展方向

5.1 怎么判断模型结果是可信的

建模完成后必须做验证，不能看起来差不多就算完成了。我惯用的验证方法有三个：

第一，功率平衡校验。把潮流结果中各发电机出力加起来，减去各负荷消耗，再减去各线路上的损耗，应当为零（误差在0.01%以内）。在Simulink里可以用Display模块直接看测量到的发电功率和负荷功率之和。

第二，和手算牛顿-拉夫逊对比。取出节点导纳矩阵和潮流结果，用MATLAB写一个简单的牛顿-拉夫逊迭代程序，对比两套结果的节点电压。误差在1e-6量级说明Simulink的实现正确。

第三，和商业软件对比。如果你手头有PSASP或BPA，可以用同一个算例分别计算，对比线路功率和节点电压。这个方法最权威，但不是所有人都有商业软件授权。

第二和第三种方法我在给研究生讲这个模型时都会推荐，尤其是手写一个20行牛顿-拉夫逊迭代脚本，对理解潮流计算本质非常有帮助，顺便还能发现自己在Simulink里设置的网络参数有没有误解。

5.2 静态安全分析与N-1校验

这个模型最直接的应用就是做N-1静态安全分析。所谓N-1，是指系统中任意一条线路或一台变压器退出运行后，系统仍能保证所有负荷供电且各线路不过载、各节点电压不越限。

在Simulink里做N-1分析非常简单：依次断路一条线路（用一个Three-Phase Breaker模块串联在线路中，仿真时断开），重新运行潮流，观察剩余线路的功率和节点电压。你可以把任意一条线路断开，记录结果，形成一张N-1校验表。

我实际测试过，在2机5节点系统里断开支路L24，节点4可能因为失去主供电路径而电压大幅下降，甚至低于0.9pu。这时候可以检查发电机的无功出力和节点PV电压设置，或者考虑增加线路输送能力。这个分析过程完整地展示了静态安全分析的研究思路。

5.3 暂态稳定仿真初值生成

前面反复提到潮流结果是暂态仿真的基础，这里给一个具体的联动流程：先按本文3.4节完成潮流计算并Apply to model，然后做一个“1.0s时线路L13中点发生三相短路，1.2s切除故障”的场景。

具体操作是在线路L13上串联一个Three-Phase Fault模块，在Fault模块里设置Transition times为[1.0 1.2]，Simulink会在1.0s触发短路，1.2s切除。把powergui切到Continuous模式，运行仿真10s，观察两台发电机的功角差。如果从潮流初始化开始，功角差应该在一个稳定的初值附近，短路瞬间突然变化，故障切除后或恢复稳定或持续增大（失稳）。

这套“先潮流初始化，再暂态分析”的流程，正是电力系统稳定分析的标准操作。这也是我为什么坚持用Simulink做2机5节点系统而不是纯手写程序的原因——同一个模型可以贯穿稳态和暂态两条分析主线。

5.4 模型扩展：从2机5节点到多机系统

5节点模型是骨架，往里添肌肉就能变成几十节点的系统。扩展方向主要有三个：增加发电机台数（比如变成3机9节点，这是IEEE标准算例）、增加负荷类型（把恒功率负荷改成恒阻抗或恒电流模型，或者加感应电动机动态负荷）、增加无功补偿设备（在负荷节点装并联电容器或SVC）。

从代码和数据层面，扩展系统最主要的成本在节点导纳矩阵的更新，好在Simulink会通过连线自动确定拓扑，不需要手写节点导纳矩阵。你只需要拖入新的模块、连好线、设置参数，重新运行Load Flow就能得到新系统的潮流结果。但要注意的是，系统规模越大，潮流的初值越敏感，有时候需要对PV节点电压和无功限幅做更细致的调整。

我个人在实际操作中最大的体会是，别嫌5节点系统“太简单”。把2机5节点系统的建模流程吃透，扩展到IEEE 39节点甚至更大的系统只是工作量问题，不是原理问题。那些我在5节点模型里踩过的坑——基准值不统一、变压器阻抗折算搞错、忘记Apply to model——换成大系统照样会遇到，而且大系统里更难排查。所以，认认真真把这个小模型做一遍，比急着上手复杂系统划算得多。