Matpower实战：IEEE30潮流计算、数据修改与收敛排查

简介：面向电力系统课程实验、算法验证与科研入门，这份zip压缩包围绕IEEE30标准测试节点系统，提供基于MATPOWER的潮流计算实现，帮助使用者快速理解并完成节点功率、电压幅值相角等核心参数的求解。压缩包共2个文件，打包1个MATLAB脚本文件与1个Simulink模型文件，整体仅56KB，前者用于设置IEEE30系统参数并调用runpf函数执行潮流计算，后者用于以图形化方式查看系统结构。已有346人学习/浏览，适合电力系统相关专业学生、研究人员及MATPOWER初学者作为直接可运行的参考案例。借助该资源可直接获得针对IEEE30节点的完整仿真文件，省去自行搭建与调试的重复过程；结合资源描述中对节点类型、发电机参数设定及潮流结果分析的说明，还可为课程实验、算法对比或更深入电力系统研究提供支撑。

1. 拿到 ieee30.zip 先别急着跑 runpf，先把 Matpower 的数据契约搞清楚

解压 ieee30.zip 后得到的是一个 case30.m 文件，这个文件就是 Matpower 的 IEEE30 节点标准算例。很多人第一次跑runpf('case30')，看到控制台刷出一堆表就以为算完了，但真正的问题是：这些数据从哪来、收敛判据是什么、换一组负荷怎么改。IEEE30 系统有 30 个节点、41 条支路、6 台发电机组，基准功率 100 MVA，是验证潮流算法、做最优潮流和电压稳定性分析时最常用的公开算例。这篇博客按实际动手顺序讲清楚 Matpower 怎么把 IEEE30 变成可计算的 mpc 结构体，以及参数怎么改、结果怎么读、不收敛时看哪里。适合刚接触电力系统仿真计算的研究生，也适合需要做批量工况计算的工程师。

2. Matpower 求解 IEEE30 潮流的原理与 case30.m 数据结构

2.1 潮流方程的牛顿-拉夫逊解法与收敛判据

潮流计算的本质是求解节点功率平衡方程。对每个节点 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)

P_i 和 Q_i 是节点注入功率，G_ij 和 B_ij 是节点导纳矩阵的实部与虚部。Matpower 内部由makeYbus构造导纳矩阵，再用牛顿-拉夫逊法迭代求解。每次迭代解线性方程 J·Δx = -F(x)，J 是雅可比矩阵，Δx 是电压幅值与相角的修正量，F(x) 是功率不平衡量。收敛判据是最大功率不平衡量小于 pf.tol，默认 1e-8 p.u.。对 IEEE30 这种规模，牛顿法通常 3 到 5 次迭代收敛，这是它成为默认算法的原因。

PQ 分解法（pf.alg 设为 'FDXB' 或 'FDBX'）把雅可比矩阵做解耦近似，每次迭代只解两个低维线性方程，单次迭代更快，但迭代次数增多，适合大规模系统。在 IEEE30 上 PQ 分解法的性能优势不明显，更常用的场景是与牛顿法互相验证：两次结果电压差超过 1e-6 p.u.，基本可以断定数据改出了问题。

2.2 mpc 结构体里的 bus、branch、gen 三张表

loadcase('case30')返回的 mpc 结构体包含 mpc.bus、mpc.branch、mpc.gen 三个矩阵，这是 Matpower 所有计算功能的统一数据入口。理解这三张表的列含义，是改数据的前提。

2.2.1 bus 矩阵哪些列决定潮流结果

mpc.bus 是 30×13 的矩阵，每行一个节点。前四列分别是节点编号、节点类型、有功负荷 Pd（MW）、无功负荷 Qd（MVAr）。节点类型里 3 代表平衡节点，IEEE30 中只有一个，位于 bus 1；2 代表 PV 节点，电压幅值给定，IEEE30 里是 2、5、8、11、13 号节点；1 代表 PQ 节点。第 8、9 列是电压幅值初值 Vm 和相角初值 Va，单位分别为 p.u. 和度，第 10 列是基准电压 kV，第 12、13 列是电压运行上下限。做电压越限检查时，改动最多的是第 8、9、12、13 列。

2.2.2 branch 与 gen 矩阵怎么读

mpc.branch 是 41×13 矩阵，第 1、2 列是首末端节点编号，第 3、4 列是电阻 r 和电抗 x（p.u.），第 5 列是充电电纳 b，第 6 到 8 列是长期、短期、紧急容量上限 rateA/rateB/rateC（MVA），第 9、10 列是变压器变比和相移，第 11 列是投运状态。mpc.gen 是 6×21 矩阵，对应 1、2、5、8、11、13 号节点的机组，第 2、3 列是发电机有功、无功出力，第 4、5 列是无功上下限，第 6 列是机端电压设定值。

需要特别指出：发电机节点的 Pg 和 Vg 是输入量，Qg 是计算结果。平衡节点只有 Vm 和 Va 作为输入，Pg 和 Qg 完全由潮流计算决定。改数据时如果动了平衡节点的 Pg，等于同时给了两个约束，雅可比矩阵会出现矛盾，这是新手最容易踩的坑。

2.3 runpf 内部做了什么：从 loadcase 到 int2ext

runpf的执行链路是：先用loadcase读数据，再调ext2int把外部节点编号映射成从 1 连续排列的内部编号——Matpower 求解时要求内部节点连续，否则稀疏矩阵组装会错位——然后makeYbus生成导纳矩阵，进入牛顿迭代，最后用int2ext把结果映射回原始编号。这个环节最常见的坑是手动删节点后没重排编号，branch 里指向已删节点的支路会让 makeYbus 直接报错。改数据后先执行mpc2 = ext2int(mpc);看内部编号有没有洞，是成本最低的检查手段。

3. 在 MATLAB 中用 Matpower 跑通 IEEE30 潮流计算

3.1 最小可运行代码与输出解读

% 加载 IEEE30 算例，返回 mpc 结构体 mpc = loadcase('case30'); % 运行牛顿-拉夫逊潮流，返回 result 结构体 result = runpf(mpc); % 检查收敛标志与迭代次数 fprintf('success = %d, iterations = %d\n', result.success, result.iterations);

逻辑说明：loadcase接受文件名或 mpc 结构体，runpf的返回值 result 里 success 为 1 表示收敛、0 表示失败，iterations 记录实际迭代次数，IEEE30 默认条件下约 4 次。如果不想看控制台输出，用mpopt = mpoption('out.all', 0);后再跑runpf(mpc, mpopt)，批量计算时必须这样做，否则上百个工况会刷屏。参数说明：runpf第二个参数是 mpoption 结构体，不传时使用全局默认选项，mpoption不带参数调用可以列出全部当前选项。

3.2 mpoption 里需要调的 4 个收敛参数

mpopt = mpoption('pf.alg', 'FDXB', 'pf.tol', 1e-10, 'pf.max_it', 50); result = runpf('case30', mpopt);

这里把算法换成快速解耦 FDXB，收敛阈值收紧到 1e-10，迭代上限放宽到 50。FDXB 在 IEEE30 上迭代次数会比 NR 多 1 到 2 次，单次迭代更快。'FDBX' 与 'FDXB' 的差别在 B' 和 B'' 矩阵的构造方式，重负荷下收敛特性有差异，工程上是两个都跑一遍对比电压结果。收敛阈值降到 1e-10 后，迭代次数可能增加 1 次，但对结果数值几乎无影响，做学术计算时按期刊要求设，不做要求保持默认即可。

3.3 修改负荷重算的正确写法

mpc = loadcase('case30'); % bus 3 的有功负荷从 2.4 MW 改成 50 MW idx = find(mpc.bus(:, 1) == 3); mpc.bus(idx, 3) = 50; mpc.bus(idx, 4) = 20; % 无功负荷同步改为 20 Mvar result = runpf(mpc, mpopt);

逻辑说明：用find定位第 3 号节点所在的行索引，再改列。这样写比直接mpc.bus(3,3)=50安全得多，因为 bus 矩阵的行号与节点编号在删除节点后会错位。第 5 行跑完潮流后，result.bus 里对应行的第 8、9 列就是该节点的新电压幅值和相角。常见误区是只改 Pd 不改 Qd，负荷功率因数突变会让无功潮流严重失真，工程做法是先确定目标功率因数，再同时写入第 3、4 列。修改发电出力时同理：mpc.gen(2, 2) = 60;把第 2 台机组的 Pg 设为 60 MW，但注意平衡节点的 Pg 不要手动设，它是自由变量。

4. IEEE30 潮流计算结果怎么读与收敛失败排查

4.1 result 结构体里电压、相角、支路潮流的字段位置

result 结构体里的 bus、branch、gen 与输入 mpc 同构，但数值是计算结果。result.bus 第 8 列是电压幅值（p.u.）、第 9 列是相角（度），result.branch 第 14、15 列是从首端流向末端的 P、Q（MW/MVAr），第 16、17 列是反向潮流，result.gen 第 2、3 列是计算后的 Pg 与 Qg。

读结果按三个顺序：先看 success 是否为 1，再看电压是否越限（IEEE30 默认运行上下限约 1.06 与 0.94 p.u.），最后看支路潮流是否超过 rateA 长期载流容量。

fprintf('Vmax = %.4f p.u.\n', max(result.bus(:, 8))); fprintf('Vmin = %.4f p.u.\n', min(result.bus(:, 8))); % 找出越限支路，第 6 列是 rateA overload = find(result.branch(:, 14) > result.branch(:, 6));

这段代码先做电压分布检查，再列出支路首端有功功率超过长期载流容量的支路编号。IEEE30 案例里出现越限通常意味着某条支路的潮流被刻意加大，例如把一台大机组切除后功率全部转移到邻近线路。

4.2 不收敛时按这 4 步定位问题

排查顺序先说最常见的负荷问题。IEEE30 原始总有功负荷约 189.2 MW，如果把某节点负荷改成几百 MW，发电机出力上限约束下功率严重失衡，牛顿法必然发散。用sum(mpc.bus(:,3))核对总有功，超过 300 MW 时就要检查是否有输入错误。其次是节点类型问题，平衡节点只能有一个，如果手动把 bus 1 类型改成了 PQ，系统失去相角参考，雅可比矩阵奇异，runpf会直接报矩阵奇异错误。

支路参数数量级也要检查。r、x 都是标幺值，IEEE30 里典型电抗在 0.01 到 0.3 之间，如果 x 为 0，导纳矩阵会出现无穷大元素。电压初值的问题往往出现在重负荷工况：平启动mpc.bus(:,8)=1.0在接近极限负荷时不收敛，改成上一次收敛结果做热启动会更稳：

mpc.bus(:, 8) = prev_result.bus(:, 8); mpc.bus(:, 9) = prev_result.bus(:, 9);

最后提一个单位坑：case30 里 Va 的单位是度，但算法内部换算成弧度。往数据文件里填相角时如果直接填弧度值，等于初值偏差 57 倍，热启动时迭代轨迹会完全偏离。

5. 用 Matpower 扩展 IEEE30 场景的三个实用技巧

5.1 批量负荷扫描的循环写法

实际工程很少只算一个工况，通常是全系统按比例升降负荷，观察哪些节点电压先跌破下限。标准做法是：

mpc = loadcase('case30'); mpopt = mpoption('out.all', 0, 'verbose', 0); for k = 0.8:0.1:1.5 mpc_k = mpc; mpc_k.bus(:, 3) = mpc.bus(:, 3) * k; mpc_k.bus(:, 4) = mpc.bus(:, 4) * k; r = runpf(mpc_k, mpopt); fprintf('k=%.1f success=%d Vmin=%.4f Vmax=%.4f\n', ... k, r.success, min(r.bus(:, 8)), max(r.bus(:, 8))); end

循环内每次新建 mpc_k 再缩放负荷，避免污染原始数据，缩放系数从 0.8 逐步加到 1.5，输出里能直接读出电压最低点跌穿 0.94 p.u. 的临界负荷倍数。想更精确就加一个二分法循环，把 k 的步长从 0.1 改成 0.01。

5.2 从潮流到最优潮流与连续潮流

IEEE30 做最优潮流用runopf，目标函数默认是发电成本最小：

mpopt = mpoption('opf.ac.solver', 'MIPS'); r_opf = runopf('case30', mpopt); fprintf('总发电成本: %.2f $/h\n', r_opf.cost);

与 runpf 的区别是 OPF 把 Pg、Vg、变压器变比作为控制变量，在满足潮流方程、发电机出力上下限和线路容量约束的条件下寻优。IEEE30 这个规模 MIPS 求解器通常 1 秒内出结果。连续潮流用runcpf，可以画 PV 曲线找鞍结分岔点，回答"极限负荷是多少"——单一工况的潮流只是静态断面，连续潮流才是电压稳定分析的工具。注意runcpf对负荷缩放的定义和上面循环里不同，需要先读help runcpf确认参数格式。

5.3 用功率平衡校验抓数据改错

最后一个落地技巧是每次都做的校验：全系统有功注入必须平衡，发电机总出力等于总负荷加全网网损。写一句断言：

loss = sum(r.gen(:, 2)) - sum(r.bus(:, 3)) - sum(r.branch(:, 14) - r.branch(:, 16)); assert(abs(loss - sum(r.branch(:, 14) - r.branch(:, 16)) / 1) < 1e-5, '功率不平衡');

实际运行中更简单的等价校验是abs(sum(r.bus(:,3)) + sum(r.branch(:,14)) - sum(r.gen(:,2)))，各项用同一套结果矩阵计算，误差小于 1e-5 就说明数据一致。改完负荷、切掉机组、调整无功补偿后先跑这一段，能一次性过滤掉大多数低级错误，再进入工况分析。这个断言配合上面的批量扫描，基本覆盖了 IEEE30 潮流计算日常使用的全部场景。

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