MATPOWER二机五节点建模与Simulink联合仿真实战解析

简介：MATPOWER是电力系统潮流计算与优化分析的常用开源工具箱，该压缩包专注二机五节点与五机二节点两类典型教学模型，适合电力系统专业学生、科研人员以及MATPOWER初学者快速上手。包体共2个文件，包含1个MATLAB脚本和1个Simulink模型，整体仅26KB，轻量精炼，便于下载后直接运行和修改学习。脚本可调用MATPOWER完成潮流计算，输出节点电压、功率分布等关键结果；Simulink模型则用于构建系统动态仿真，帮助理解多发电机间的同步与交互特性。压缩包虽小，但覆盖了从静态潮流到动态仿真的基础流程，可配合教学演示或课程设计使用。已有579人学习下载，是入门电力系统分析、验证MATPOWER用法的实用素材。

1. 为什么二机五节点是MATPOWER入门绕不开的最小算例

最近翻到一个case2_5.rar，解压后只有case2_5.m和cljs.slx两个文件。初看像普通作业包，但对做过机电暂态仿真的人而言，二机五节点（Two Machine Five Bus）恰好是判断“会不会用 MATPOWER 和 Simulink 配合”的试金石。它比 IEEE 9 节点更容易手算验证，又比单机无穷大系统多一台发电机，能清晰看到两机之间的功角摇摆、功率分配和电压相互牵制。本文从这两个文件出发，先拆解 MATPOWER 中二机五节点的建模方法，再讲runpf潮流计算和cljs.slx的联合仿真衔接，最后给出扩展成五机二节点时的重构思路。文中所有代码都可以直接落到 MATLAB 环境里跑，重点在参数边界和排错习惯，不写空泛概念。

2. 从 case2_5.m 看 MATPOWER 数据结构与二机五节点参数边界

2.1 两机五节点模型的拓扑结构与物理含义

二机五节点系统通常由两个发电机节点和三个负荷节点组成，发电机节点可以是平衡节点和 PV 节点，负荷节点是 PQ 节点。拓扑上，发电机节点通过较长线路连接到负荷区域，两个发电机节点之间可能直接相连，也可能通过中间负荷母线间接耦合。这种结构模拟的是“两个电厂向同一片负荷区供电”的场景。系统中每一台发电机的有功、无功出力会直接影响另一端母线的电压幅值和相角。相比三节点模型，二机五节点多了无功分配的维度，因此很适合验证 PV 节点无功上下限对系统电压稳定性的作用。

在 MATPOWER 里，模型不是 Simulink 里的图形，而是由mpc.bus、mpc.gen、mpc.branch、mpc.gencost四个核心矩阵构成的 case 文件。分析case2_5.m时，必须同时理解这四个矩阵之间的索引关系，否则后续改负荷或扩节点很容易因为行号错位导致潮流不收敛。

2.2 一个可直接运行的 case2_5.m 骨架

实际的case2_5.m可能有更多细节，但核心结构一定是下面这样的：

function mpc = case2_5 % 二机五节点 MATPOWER 算例 % 节点1: 平衡节点 (Slack) % 节点2: PV 节点 % 节点3~5: PQ 负荷节点 mpc.version = '2'; mpc.baseMVA = 100; % 基准功率 100 MVA mpc.bus = [ 1 3 0 0 0 0 1 1 0 230 1 1.1 0.9; 2 2 0 0 0 0 1 1 0 230 1 1.1 0.9; 3 1 125 50 0 0 1 1 0 230 1 1.1 0.9; 4 1 90 30 0 0 1 1 0 230 1 1.1 0.9; 5 1 60 20 0 0 1 1 0 230 1 1.1 0.9; ]; mpc.gen = [ 1 150 0 100 -100 1 1 1 0 0 0; 2 100 0 100 -100 1 1 1 0 0 0; ]; mpc.branch = [ 1 3 0.0100 0.0600 0.0600 170 0 0 0 0 1 -360 360; 2 4 0.0150 0.0800 0.0500 170 0 0 0 0 1 -360 360; 3 4 0.0200 0.1000 0.0400 130 0 0 0 0 1 -360 360; 3 5 0.0250 0.1200 0.0300 130 0 0 0 0 1 -360 360; 4 5 0.0300 0.1500 0.0200 130 0 0 0 0 1 -360 360; ]; mpc.gencost = [ 2 1500 0 3 0.0300 0.2000 0.0100; 2 1500 0 3 0.0200 0.1500 0.0200; ];

mpc.bus每一行固定是 13 列，其中第 2 列是节点类型，3表示平衡节点，2表示 PV 节点，1表示 PQ 节点。第 3、4 列是负荷有功和无功，这里所有负荷都放在 3、4、5 号节点，1、2 号节点纯发电，这样两机的相互作用更直观。mpc.gen中第 1 列是发电机所在母线的编号，第 2、3 列是初始有功和无功出力，第 4、5 列是 Qmax / Qmin，第 6、7 列是电压幅值和机组状态。mpc.branch中第 1、2 列是首末端母线号，第 3~5 列是线路电阻、电抗、对地电纳，单位全部是标幺值。这里的baseMVA = 100，意味着分支参数是基于 100 MVA 折算的，扩容或改造算例时不能直接套用，否则电压降落会失真。

2.3 关键列向量的含义和参数边界

下面这张表总结了case2_5.m中每个字段在后续runpf中的实际作用。

修改负荷或发电时，gen的 Qmax / Qmin 必须和bus类型联动。比如把 2 号节点改成容量更大的发电机，如果 Qmax 还是 100 Mvar，高负荷水平下无功可能顶到上限，MATPOWER 会把它从 PV 节点变成 PQ 节点，电压就会脱离设定值。很多人遇到“为什么发电机电压设了 1.0，潮流结果却是 0.96”的疑问，原因就在这里。

2.4 常见误用：把发电机节点全设为平衡节点

刚接触 MATPOWER 的人常常把两个发电机节点都设成类型3，以为这样更“真实”。实际上一旦存在两个平衡节点，系统功率基准就失去了参照，潮流计算会出现一个自由度冗余，导致结果无法收敛或某个平衡节点的有功出力被算法强行置零。正确做法是只保留一个平衡节点，另一个发电机节点用 PV 节点，给定有功出力和电压幅值，让 MATPOWER 通过调节无功和注入电流来平衡系统。case2_5.m里节点 1 作为 slack，节点 2 作为 PV，就是为了让两台发电机在迭代过程中自然形成功角差。

另一种常见误用是把负荷直接写成负的发电机功率，而不放在bus的 PD / QD 列。这样虽然也能算，但负荷随电压变化的特性没了，后续接 Simulink 动态仿真时，负荷模块和潮流初值对不上，经常出现初始化失败。因此建议始终用标准的 PD / QD 定义负荷。

3. 用 runpf 迭代二机五节点潮流：结果校验与灵敏度调整

3.1 第一次运行：先看收敛状态而不是只读电压

拿到case2_5.m后，第一个动作不是直接看结果，而是确认潮流迭代是否收敛。MATPOWER 的runpf支持迭代过程打印和输出的精细控制。打开 MATLAB，将case2_5.m放到当前路径，执行以下脚本：

mpc = case2_5; opt = mpoption('verbose', 2, 'out.all', 0); res = runpf(mpc, opt); fprintf('success=%d, iterations=%d\n', res.success, res.iterations);

mpoption('verbose', 2)会让runpf每迭代一次就打印一次最大功率失配量，方便观察是电压迭代震荡还是无功越限。out.all设置为 0 表示不保存输出文件，避免每个算例留一堆垃圾文件。res.success为 1 代表牛顿法收敛，res.iterations是实际迭代次数。二机五节点这种规模，正常情况 3~6 次就能收敛；如果超过 15 次，优先检查mpc.gen中是否有 Qmax / Qmin 设置过紧导致节点类型反复切换。

收敛后，res.bus里存着完整的节点状态。下面这个脚本可以很清楚地展示每个节点的稳态电压：

V = res.bus(:, 8); Va = res.bus(:, 9); for i = 1:length(V) fprintf('bus %d: V=%.4f p.u., angle=%.4f deg\n', ... i, V(i), Va(i)); end

res.bus的第 8 列是电压幅值标幺值，第 9 列是相角度。注意节点 1 是平衡节点，相角固定为 0，节点 2 的相角会随负荷变化。运行后如果发现 V(3)、V(4) 明显低于 0.95，说明无功支撑不足，应增大节点 2 的电压设定值或降低线路电抗。这里不是看一个电压值合格就行，要同时检查节点 2 和节点 3、4 之间的相角差，因为二机系统最关心的就是两机间的功角差。

3.2 从结果矩阵里提取功率流并做功率平衡校验

潮流算完以后，除了电压，还要核对线路潮流和发电出力，否则无法确定结果可用。res.branch中第 14、15、16、17 列分别存储线路首端有功、首端无功、末端有功和末端无功。下面这段代码提取所有支路的有功功率，并计算总网损：

Pflow = res.branch(:, 14); % 首端有功，正方向为 from -> to Ploss_total = sum(res.branch(:, 14) + res.branch(:, 16)); fprintf('total line loss = %.2f MW\n', Ploss_total * mpc.baseMVA);

这里把标幺值乘以baseMVA得到兆瓦数。res.branch(:, 14) + res.branch(:, 16)首端有功加末端有功就是该支路损耗，取所有支路损耗之和即可得到整个网络的网损。用这个值可以验证潮流结果的合理性：外面最热的天气下，负荷 300 MW 左右时，三条配电线路网损如果超过 10 MW，通常说明线路电阻参数偏大或无功流动太多。runpf返回的res.gen(:, 2)是发电机最终有功，两台机组出力之和应该等于所有负荷有功加上网损。

做功率平衡校验还有一个隐藏好处：能提前发现case2_5.m的节点负荷是否把厂用电算重了。如果总发电和总负荷相差超过 2%，多半是mpc.bus里某个节点的 PD 列和mpc.gen的厂用电重复计算，这种情况在网上下载的算例里很常见。

3.3 负荷水平扫描：观察电压稳定边界

二机五节点最适合做负荷爬坡实验。把 3、4、5 号节点的负荷统一乘以一个系数，反复调用runpf，可以得到系统的最低电压曲线。这是一种不需要搭建动态模型就能提前发现电压稳定边界的方法：

loads = [0.5, 0.8, 1.0, 1.2, 1.5]; V_min = zeros(1, length(loads)); for k = 1:length(loads) mpc = case2_5; mpc.bus(3:5, 3) = mpc.bus(3:5, 3) * loads(k); mpc.bus(3:5, 4) = mpc.bus(3:5, 4) * loads(k); res = runpf(mpc, mpoption('verbose', 0)); if res.success V_min(k) = min(res.bus(:, 8)); else V_min(k) = NaN; end end [loads; V_min]

这里直接修改mpc.bus的 PD 和 QD 列。注意同时缩放有功和无功，负荷功率因数保持不变。runpf失败时返回的res.success为 0，V_min设为 NaN 方便后续绘图。实际运行结果一般是：负荷倍率在 0.5 到 1.0 之间时最低电压缓慢下降；超过 1.2 后电压会突然加速下跌，这就是所谓的“电压崩溃前兆”。二机系统中，节点 4 和 5 的电压通常最先掉，因为它们离两台发电机都远，又承担了较多负荷。

典型结果如下表：

这个表对应一个“正常阻抗、无功上限足够”的二机五节点算例。如果你拿到的case2_5.m在 1.2 倍负荷时就发散，优先检查节点 2 的 Qmax 是否小于 200 Mvar。

3.4 不收敛时的检查顺序

潮流不收敛时，按下面的顺序排查，不要直接改线路参数。第一，检查是否只有一个平衡节点。第二，检查 PV 节点的无功限幅是否被顶住，runpf的 verbose 输出会出现“Qmax limit reached for generator at bus 2”之类的信息，这时要把 Qmax 调大，或者把该节点改成 PQ 节点观察无功需求量。第三，检查线路电抗 x 是否太大，x 超过 0.3 p.u. 时二机五节点系统的潮流极难收敛，因为线路两端相角差过大。第四，检查负荷 PD / QD 是否超过发电机总容量，超了以后算法无论如何都无法满足功率平衡。这些检查步骤的目的不是改到收敛为止，而是要找到系统真正需要调整的参数边界。

4. cljs.slx 与 MATPOWER 的联合仿真：从潮流状态到动态响应

4.1 为什么动态仿真必须有稳态潮流初值

cljs.slx如果是 Simulink 模型，多半是二机五节点的机电暂态仿真模型。Simulink 里的同步发电机模块需要初始转速、初始转子角、初始励磁电压和初始机械功率，这些初值不是随便填的，而是必须由潮流计算给出。许多人在 Simulink 里手搭两机模型，仿真一开始功角就乱跳，问题出在初值不是稳态解。MATPOWER 算出的稳态电压、相角和发电机出力，正是 Simulink 中同步电机模块 initialization 所需的参数。

常见做法是：在 Simulink 模型里使用powergui的 Machine Initialization 工具，把潮流结果填进去。但为了批量实验，更推荐用 MATLAB 脚本直接生成并写入初值。下面这段脚本展示了从res结构体提取初值并写回 Simulink 的思路：

res = runpf('case2_5', mpoption('verbose', 0)); Vmag = res.bus(:, 8); Vang = res.bus(:, 9) * pi / 180; Pgen = res.gen(:, 2) * res.baseMVA; set_param('cljs/Generator1', 'Pm', num2str(Pgen(1))); set_param('cljs/Generator1', 'Vf', num2str(1.0)); % 初始励磁需要根据内部状态计算

上面代码中res.baseMVA自动取到 case 文件里的基准功率，Pgen单位转为 MW。set_param的第二个参数Pm是同步电机的机械功率输入，具体名称取决于模块库版本，实际使用时先用get_param('cljs/Generator1', 'DialogParameters')查看可配置项。更稳妥的方式是在powergui的 Machine Initialization 对话框里点击 “Compute from Initial State”，它会读取潮流数据并自动计算励磁和机械功率。

4.2 cljs.slx 中模块参数与潮流结果的双向同步

cljs.slx里如果包含两个同步电机模块、五条母线和对应的负荷模块，那么 MATPOWER 的潮流结果必须映射到三个地方：电压源模块的幅值和相角、同步电机的初始转子和励磁、负荷模块的初始吸收功率。下面给出一个封装成函数的方法，方便每次修改 case 后重新同步：

function init_cljs_sim() % 重新计算潮流并初始化 Simulink 模型 cljs.slx mpc = case2_5; res = runpf(mpc, mpoption('verbose', 0)); if ~res.success error('潮流不收敛，无法初始化 Simulink 模型'); end assignin('base', 'res_case2_5', res); set_param('cljs', 'SimulationCommand', 'update'); end

assignin('base', ...)把潮流结果写到 MATLAB 基准工作区，Simulink 模型里的Goto/From或 mask 参数可以直接引用工作区变量。set_param('cljs', 'SimulationCommand', 'update')是告诉 Simulink 重新检查模型参数，这样如果在模型里写了res_case2_5.bus(3, 8)之类的表达式，它会立刻刷新。实际工程中，我会把init_cljs_sim()放在模型打开前的InitFcn回调里，每修改一次case2_5.m，运行一次函数，Simulink 初值就自动更新。

4.3 跑动态过程：典型短路故障观察功角摇摆

Simulink 模型准备好后，可以在任意一条线路上加三相短路故障，观察两机之间的转子角差变化。对二机五节点系统来说，最有价值的输出是两台发电机的转子角差曲线。执行下面的脚本可以跑一次 5 秒仿真并画出功角动态：

simOut = sim('cljs', 'StopTime', '5'); t = simOut.tout; delta1 = simOut.logsout.get('Generator1.delta').Values.Data; delta2 = simOut.logsout.get('Generator2.delta').Values.Data; figure; plot(t, (delta1 - delta2) * 180 / pi, 'LineWidth', 1.2); grid on; xlabel('Time (s)'); ylabel('Rotor angle difference (deg)');

Generator1.delta是日志记录信号名，实际名称由cljs.slx内部定义，可以在信号线上点右键选择 Export Logged Data 查看。曲线如果在前 1 秒内出现超过 180° 的摆开，说明系统阻尼不足或者故障时间过长。二机五节点结构简单，一般故障持续 0.1~0.2 秒后切除，功角应该逐渐收敛。如果曲线持续振荡不衰减，检查同步发电机的阻尼绕组参数，机械阻尼系数 D 是否太小。

4.4 联合仿真的参数同步陷阱

二机五节点联合仿真最容易出问题的点是基准值匹配。MATPOWER 里所有结果都是标幺值，而 Simulink 同步电机模块的基准容量可能默认是某个具体值，两者的baseMVA不一致，初值填进去就会造成瞬态偏移。确认case2_5.m里的baseMVA是否与同步电机模块的 nominal power 一致。另一个坑是负荷模型的初始功率：Simulink 中的三相并联 RLC 负荷模块需要输入初始有功和无功，如果直接键入 125 MW 和 50 Mvar，但潮流结果经过变压器变比折算后实际消耗与设定不符，动态仿真仍会有一个持续的过渡过程。最好通过powerlib中的 Three-Phase Dynamic Load 模块，把电压和有功、无功按照潮流结果自动计算恒定阻抗，避免初值偏差。

联合仿真本身不是一次“能不能够用”的问题，而是一个不断缩小稳态误差的过程。检验方法是先不设故障，只跑 0.1 秒，观察功角和电压是否水平不动。如果一启动就有直线上升，说明初值里有一项参数没有同步好，回到 4.2 节重新检查 init 函数。

5. 从二机五节点扩展到五机二节点：重构思路与排错技巧

五机二节点算例平时不常见，但它在“多机并联到同一条母线”的研究里非常有用。扩展思路很直接：让多个发电机挂在同一节点，把 bus 数量降下来，gen 数量增上去。MATPOWER 允许mpc.gen中出现多台发电机使用同一个 bus 编号，但要注意保持潮流可解。下面这段示例把五台发电机分配到两个母线：

mpc.bus = [ 1 3 0 0 0 0 1 1 0 230 1 1.1 0.9; 2 2 300 100 0 0 1 1 0 230 1 1.1 0.9; ]; mpc.gen = [ 1 120 0 100 -100 1 1 1 0 0 0; 1 100 0 80 -80 1 1 1 0 0 0; 1 90 0 60 -60 1 1 1 0 0 0; 2 120 0 100 -100 1 1 1 0 0 0; 2 110 0 80 -80 1 1 1 0 0 0; ]; mpc.gencost = [ 2 1500 0 3 0.0200 0.1500 0.0100; 2 1500 0 3 0.0250 0.1800 0.0150; 2 1500 0 3 0.0300 0.2000 0.0200; 2 1500 0 3 0.0180 0.1200 0.0120; 2 1500 0 3 0.0240 0.1600 0.0140; ];

这里节点 1 是平衡节点，节点 2 是 PV 节点。同一个母线下的多台发电机，它们之间的有功分配由gencost决定。如果不关心经济调度，可以把所有发电机的成本曲线设成一样，否则 MATPOWER 会在同一母线内的机组间按成本优化分配功率。这个扩展的常见坑是：mpc.bus中的 PD / QD 仍然写在节点 2 上，但节点 2 本身是 PV 节点，负荷不能超过五台机的无功支撑能力。若潮流失败，把节点 2 的 Qmax 累加值调大是第一步。

排错技巧：无论二机五节点还是五机二节点，遇到不收敛先做“单机验证”。把非平衡节点的所有发电机逐个改成 PQ 节点，即把gen对应行的母线类型改为1，同时把gen的无功上下限设成该机实际发出的无功值。如果这样收敛，说明原来 PV 节点电压设定值太高，导致无功越限；如果这样也不收敛，问题在节点负荷或线路阻抗本身。这个方法可以快速定位“是发电机的问题还是网络的问题”，比反复调 Qmax 效率高得多。

最后一个实用验证：把五机二节点的节点 2 电压幅值从 1.0 逐步降到 0.97，同时观察潮流迭代次数变化。如果电压降低反而让迭代次数下降，说明之前接近节点无功极限，扩展后需要重新设计无功补偿或调整发电机无功上限。一套参数调顺以后，这套 case 文件和 Simulink 模型还可以继续用来做小干扰稳定分析和短路电流校验，不局限于单纯算潮流。

本文还有配套的精品资源，点击获取