企业官网建设流程全解析

简介：本资源是一套基于MATPOWER工具箱的电力系统潮流计算实践代码包，面向电气工程专业本科生、研究生及电力系统仿真初学者，用于掌握多机多节点系统的建模、导纳矩阵构建与稳态潮流分析全流程。压缩包共6个文件，含3个核心MATLAB脚本（用于构建节点导纳矩阵、数据校验及主案例定义）、1个MATPOWER 7.1安装压缩包、1份详细使用说明文档（.docx）和1个文本说明文件（.txt），总大小32.09MB。已有4151人学习下载，反映出其在教学实验与课程设计中的广泛适用性。用户可直接运行Case28_08.m调用runpf完成8机28节点系统的潮流计算，代码中完整实现了发电机出力约束设置、含变压器支路建模、节点无功补偿配置及三类基础数据矩阵（节点、发电机、支路）的规范构造，具备良好的可读性与复现性，是理解MATPOWER案例文件结构与电力系统建模逻辑的实用参考。 Matpower自带案例库里有9节点、14节点、30节点、39节点这些熟面孔，但提到“8机28节点”，老手第一反应通常是：这不是标准型号。这意味着你得自己搭一个可运行的电力系统模型。这不算难，但里面琐碎的坑非常多，尤其是数据结构不熟的时候，一个列填错就能让潮流迭代半天还不收敛。我最近在一个项目里正好用Matpower完整搭建了一套8机28节点系统，从母线数据到发电机组、支路参数、最优潮流，再到结果可视化，整个过程记录下来，给要做论文仿真、毕业设计或者工程测试算例的朋友做个参考。

1. 为什么是8机28节点：自建测试系统的现实需求

1.1 Matpower到底是什么

Matpower是一个基于MATLAB的开源电力系统分析工具箱，核心功能包括潮流计算（power flow）、最优潮流（OPF）、连续潮流（CPF）、状态估计等。它最大的优势是数据接口极其简洁：把整个电力系统定义成一个结构体mpc，包含母线、发电机、支路、成本函数四张核心矩阵，然后调用一句runpf(mpc)就能完成交流潮流计算，runopf(mpc)完成最优潮流。对研究者来说，不需要从零写牛顿-拉夫逊迭代，也不需要自己处理稀疏矩阵的编号问题，省下来的时间可以全部花在算法和策略上。

1.2 28节点为什么值得“手搓”

IEEE标准测试系统里确实没有28节点，实际项目中遇到8机28节点，一般有三种来源：一是从30节点系统删改而来，二是来自某篇论文的算例数据，三是仿真者根据实际电网抽象出来的简化模型。无论哪种来源，核心都指向一个需求：需要一个规模适中、机组数量足够多、能体现多机协调特性的测试系统。

8台发电机意味着系统可以在潮流层面讨论机组出力分配、经济调度，在动态层面可以研究机组间的低频振荡、暂态稳定，甚至故障后的功率转移。28个节点又保证了拓扑结构有一定的复杂度——可以设置不同的电压等级、多条联络线、多个负荷中心，而不像9节点或14节点那样“一眼看穿”。

1.3 这类自建系统最常用的场景

我的实际体验是，这类系统最集中的需求集中在三个方向：

• 算法验证：比如粒子群优化求解经济调度、神经网络做负荷预测后的潮流校验、强化学习做电压控制，都需要一个“比小系统有说服力，又比大系统跑得快”的算例。
• 教学实验：本科生或者研究生做课程设计，需要一个能手动改写参数、观察不同运行方式影响的中等规模系统。
• 工程预研：在真实电网数据不方便公开时，用一个结构类似的抽象系统做方案验证。

所以，8机28节点不是凭空想出来的数字，它恰好是“教学与研究之间”的一个平衡点。

2. bus、gen、branch三张表：Matpower数据结构的底层逻辑

所有Matpower程序设计都绕不开三张表：mpc.bus、mpc.gen、mpc.branch。这三张矩阵就是整个电力系统的“骨架”，搞不清楚它们的列含义，后面所有工作都会像在迷宫里打转。

2.1 全局参数和数据结构

先看一个典型的Matpower数据文件开头：

function mpc = case28_8m mpc.version = '2'; mpc.baseMVA = 100;

baseMVA是系统功率基准值，通常取100 MVA。潮流计算中电压是标幺值，功率是有名值（MW/Mvar），支路阻抗是标幺值。很多人第一次接触时容易把Pd、Qd写成标幺值，这是非常常见的错误——Matpower要求母线负荷直接填有名值MW/Mvar，阻抗才是标幺值。

2.2 bus矩阵：电力系统的“节点档案”

mpc.bus一共13列，每一行代表一个节点。列含义如下：

写文件时最容易犯的错误是type选错。功率平衡节点只能有一个，通常选在装机容量最大的母线上，负责平衡全网功率差额；PV节点接发电机，给定有功出力和电压幅值；PQ节点是普通负荷节点。如果平衡节点缺失，潮流计算会直接报“slack bus not found”之类的错误。

2.3 gen矩阵：发电机的“运行档案”

mpc.gen矩阵每一行对应一台发电机，核心前10列分别是母线号、有功出力Pg、无功出力Qg、无功上限Qmax、无功下限Qmin、电压设定值Vg、机组容量基准mBase、运行状态、有功上限Pmax、有功下限Pmin。

这里有一个非常关键的细节：mpc.bus里的Vm和mpc.gen里的Vg都要给一个合理的电压初值。PV节点的电压设定值，在bus矩阵第8列和gen矩阵第6列都有体现，我个人的习惯是两处保持一致，避免排查问题时被自己搞晕。实际项目里很多人只改了bus表里的电压，忘了gen表里的Vg，最后结果和预期对不上。

2.4 branch矩阵：电力系统的“连接关系”

mpc.branch每一行代表一条输电线路或者变压器支路，前13列的含义是：

纯线路的ratio填0，有载调压变压器则填变比值。计算支路负载率时，Matpower调用的就是rateA这一列，如果全填0，运行结果里支路潮流依然能算出来，但负载率显示为0，这一细节在做安全分析时特别容易踩。

2.5 gencost：最优潮流的经济参数

mpc.gencost是成本函数矩阵，runopf求解最优潮流的依据就在它里面。一般用多项式模型，指定机组有功出力的二次成本函数，形如c2*P^2 + c1*P + c0。由于不同版本列数定义略有差异，我建议直接用loadcase('case30')打开自带案例，仿照它的gencost格式去改，永远比对着文档手搓要稳。

3. 从零到可运行：8机28节点系统的搭建全流程

在搞清楚三张表之后，剩下的工作就是“填格子”。但怎么填、填多少、怎么验证，是有讲究的。

3.1 最稳的起点：从case30改造

我强烈不建议从空白文件开始盲写28节点。实际项目里最快的路径是：以case30为模板，改造成8机28节点。

理由很简单：case30的参数是经过验证的，潮流能收敛，负荷水平、阻抗数量级、电压范围都是合理的。你只需要做三件事：删掉两个节点、增加两台发电机、重新整理编号。

在MATLAB里手动改case30.m比较麻烦，更推荐用代码的方式处理：

% 加载自带30节点系统 mpc = loadcase('case30'); % 删除第29号和第30号节点（以及连接到它们的支路） bus_del = [29; 30]; branch_del = []; for i = 1:size(mpc.branch, 1) if ismember(mpc.branch(i, 1), bus_del) || ismember(mpc.branch(i, 2), bus_del) branch_del = [branch_del; i]; end end mpc.branch(branch_del, :) = []; mpc.bus(bus_del, :) = []; % 删除与这两个节点关联的发电机组（如果有） gen_del = []; for i = 1:size(mpc.gen, 1) if ismember(mpc.gen(i, 1), bus_del) gen_del = [gen_del; i]; end end mpc.gen(gen_del, :) = []; % 修改节点编号：删除中间编号后会产生跳号 % 使用 ext2int 重新编号 mpc = ext2int(mpc);

ext2int这个函数是Matpower内部的编号统一工具，可以把不连续的外部节点编号映射成1到N的连续内部编号。很多人在删除节点后不重新编号，结果bus矩阵第1列出现了29、30这种空洞，虽然Matpower多数情况下能处理，但后续可视化、结果导出到Excel时非常容易出错。

3.2 增加两台发电机

删除节点后，原case30系统是6台发电机，离8机还差两台。加机组的时候要考虑节点类型的变化：

% 在第5号和第12号母线上分别加一台发电机 newgen1 = [5, 40, 0, 60, -60, 1.0, 100, 1, 120, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; newgen2 = [12, 50, 0, 80, -40, 1.0, 100, 1, 150, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; mpc.gen = [mpc.gen; newgen1; newgen2];

注意：mpc.gen一个矩阵有21列（版本不同列数可能略有差异），补齐时必须把后面不用的列都填0。同时要把mpc.bus中对应节点的type从1改为2，表示它变成了PV节点：

mpc.bus(5, 2) = 2; mpc.bus(12, 2) = 2;

这一步很多人会忘，结果新增的发电机在潮流计算里根本没被当成PV节点，系统里的机组数量看着是8台，实际参与电压调节的只有原来的6台。

3.3 从零开始写28节点数据文件

如果手里已经有文献里的系统参数，不愿意从case30改，也可以直接写数据文件。我的建议是先把单线图拓扑画清楚，标好节点号、支路阻抗、负荷大小、发电机接入点，然后按顺序填表。写之前先规划好3个问题：

• 系统的基准电压用什么等级？全部用同一个电压等级最省事（比如230kV），否则变压器支路的变比参数要仔细算。
• 平衡节点放在哪？选一台容量大、有调节能力的机组所在节点，通常放节点1。
• 系统总负荷量级是多少？自建系统建议控制在300~500MW之间，太大会导致机组出力超过Pmax，太小则电压普遍偏高。

实际编写时，数据文件的格式就是把我前面讲的三张表依次填好，保存成case28_8m.m即可。写完之后用loadcase加载，再调用check_mpc做基础校验，这个函数会检查矩阵尺寸、节点编号、支路端点是否越界等常见问题。不同版本check_mpc调用方式略有差异，但载入数据时如果数据有致命错误，命令窗口通常会直接报错或者给出warning。

3.4 跑通第一个潮流计算

数据文件准备好了，运行潮流只需要几行代码：

mpc = loadcase('case28_8m'); opt = mpoption('verbose', 2); res = runpf(mpc, opt);

verbose设置为2，可以在命令窗口看到牛顿法每次迭代的最大功率失配量。正常情况下3~5次迭代就收敛了，res.success返回1。如果verbose=0，则只显示收敛结果，不显示迭代过程。

跑通之后先别急着高兴，先做一次基础自检：查看发电机总出力、负荷总量、网损、各节点电压幅值，确认数据没有“表面收敛但结果离谱”的情况。

4. 潮流结果怎么读：从res结构体到可视化

runpf返回的结果是一个结构体res，里面包含了完整的潮流计算结果。很多新手跑完只看一个res.success就结束了，这太浪费了。

4.1 结果结构体的关键字段

res里的核心字段和mpc是同构的，也包含bus、gen、branch三张表，但各表内容已经被刷新成了潮流计算后的运行状态：

• res.bus(:, 8)：各节点电压幅值（pu）
• res.bus(:, 9)：各节点电压相角（度）
• res.gen(:, 2)：各发电机实际有功出力（MW）
• res.gen(:, 3)：各发电机实际无功出力（Mvar）
• res.branch(:, 14)和res.branch(:, 15)：支路首末端有功潮流（MW）
• res.branch(:, 16)和res.branch(:, 17)：支路首末端无功潮流（Mvar）

注意不同版本Matpower里，branch矩阵在结果中的实际潮流列位置可能有差异。最稳妥的方式是在命令窗口输入res.branch(1, :)，看一下前几列，确认哪些列是有功、无功，再看结果。

4.2 怎么判断结果是否合理

一个可用的潮流结果必须满足几个基本条件：

• 电压幅值基本在0.95~1.05 pu之间，如果出现0.8或1.2这种数值，说明系统无功分布有问题。
• 相角是逐步变化的，相邻节点相角差通常在5度以内，出现几十度的相角差说明支路参数异常。
• 发电机无功出力在Qmin和Qmax范围内。如果某台机组无功出力打到上限，说明该节点附近无功不足，系统可能处于临界运行状态。
• 总有功出力减总负荷等于网损，网损一般在总负荷的2%~5%之间。如果网损达到负荷的30%以上，支路阻抗的数据几乎可以肯定填错了。

判断网损在Matpower里可以直接算：

P_gen_total = sum(res.gen(:, 2)); P_load_total = sum(res.bus(:, 3)); P_loss = P_gen_total - P_load_total; fprintf('系统总负荷: %.2f MW\n', P_load_total); fprintf('发电机总出力: %.2f MW\n', P_gen_total); fprintf('网损: %.2f MW (%.2f%%)\n', P_loss, P_loss / P_load_total * 100);

4.3 可视化与报告导出

Matpower本身不自带特别花哨的画图工具，但用MATLAB基础绘图功能就足够完成大部分展示需求。我常用的三个可视化手段：

• 节点电压幅值柱状图，快速判断全网电压水平是否正常。

figure; bar(res.bus(:, 1), res.bus(:, 8)); xlabel('节点编号'); ylabel('电压幅值(pu)'); title('28节点系统各母线电压幅值'); ylim([0.9 1.1]);

• 支路负载率排序图，找到系统的重载线路。

loading = res.branch(:, 14) ./ res.branch(:, 6); figure; bar(loading * 100); xlabel('支路编号'); ylabel('负载率(%)'); title('各支路负载率');

（注意：如果rateA填了0，这里的负载率全是Inf或NaN，所以我前面才强调rateA一定要填。）

• 相角分布曲线，能看到以平衡节点为参考的相角变化趋势。

在做课程报告或者论文的时候，我通常还会把res.bus和res.branch导出为Excel或者CSV。

writetable(array2table(res.bus), 'result_bus.xlsx'); writetable(array2table(res.branch), 'result_branch.xlsx');

这样后续在Word或者LaTeX里引用数据很方便，不用每次重新跑一遍再手动抄写。

5. 排错实战：搭建与仿真中最容易踩的坑

这部分是我最想写的，因为8机28节点系统本身不算复杂，真正的门槛全在那些“看似不起眼的错误”上。

5.1 潮流不收敛：先查数据，再调迭代

好几次项目里，我搭建的自定义系统第一次跑runpf都会碰到不收敛。Matpower的提示一般是“Newton's method power flow did not converge in 10 iterations”。碰到这个提示，很多人的第一反应是增加迭代次数，这完全是本末倒置。

我的排查顺序是固定的：

1. 检查bus表里有没有type=4的隔离节点，如果负荷节点被错误标为隔离节点，系统会出现功率缺口。
2. 检查每个节点是否至少连了一条支路，孤立节点会让雅可比矩阵奇异。
3. 检查支路x是否为0或者极小值，如果某条支路电抗填了0.0001，相当于一条短路支路，矩阵条件数急剧恶化。
4. 检查负荷总量和发电机总出力是否匹配，如果负荷远大于发电机出力上限，潮流很难收敛到可行域。
5. 检查是否有节点电压初值设置得过低（比如0.8pu），导致迭代起点太差。

有一次我排查了很久，最后发现是某条支路我复制粘贴时把x=0.06写成了x=0.006，数量级差了一个数量级，潮流就怎么都算不拢。所以数据文件的校验一定要做。

5.2 PV节点无功越限：一个最容易忽略的细节

Matpower的牛顿法在迭代过程中会考虑发电机的无功上下限。当某台发电机需要发出的无功超过Qmax时，它会自动把该节点从PQ节点处理，给定无功在Qmax处，电压由潮流自然决定。

这个机制本身是合理的，但它带来的坑是：你的8机系统里某台发电机可能名义上是PV节点，实际上因为Qmax设置太小，在结果里已经变成了无功封顶的状态，电压远远低于初始设定值。很多人做最优潮流或者静态安全分析时看不到这个现象，是因为默认情况下Matpower输出里只在“Generator Reactive Power Limits”部分显示Limits信息，不仔细看很容易忽略。

我的经验是：新增发电机时，Qmax/Qmin不要随意填，要参考原有case30里的数值量级。比如100MVA基准下，Qmax给60~80Mvar、Qmin给-60~-40Mvar是比较常见的。如果填了±9999这种极端值，潮流是很容收敛，但结果没有任何工程意义。

5.3 节点编号与行号的混淆

bus矩阵的第1列是节点编号，而矩阵的行号是顺序索引。删除节点后如果不重新编号，节点编号可能变成[1,2,3,5,6,7,...]这种不连续序列，当你在后续程序里用res.bus(find(res.bus(:,1)==5), 8)这种方式取数据时容易漏项。

更麻烦的是，gen矩阵和branch矩阵里引用节点编号时，必须和bus表第1列保持一致。如果节点编号不连续，只要有任何一个地方写错，就会引用到错误的节点。所以我强烈建议在删改节点后立刻调用ext2int重新编号，或者在bus表里把编号重新写成连续的1到N。

5.4 支路阻抗的单位和基准问题

branch表里的r、x、b都是标幺值，不是欧姆，也不是西门子。很多从实际工程数据换算过来的人，容易把线路的欧姆值直接填进表里。标幺值和有名值的换算公式是：

Zpu = Z_ohm * baseMVA / (baseKV^2)

一个230kV、100MVA基准下的系统，如果一条线路阻抗是10欧姆，对应的标幺值是10 * 100 / (230^2) ≈ 0.0189 pu。如果直接把10填进去，潮流结果会显示电压崩溃，因为系统的电抗被放大了500多倍。

5.5 建议用自带工具做一次完整检查

Matpower有check_mpc这样的数据校验函数，虽然不同版本调用方法不一致，但它能检查矩阵尺寸、支路端点引用、重复节点编号等明显问题。我自己一般在写完数据文件后，会先跑一下这个函数，确认没有结构性问题后再跑潮流。另外一个笨办法是：把系统节点数、支路数、发电机数打印出来，和设计图对照一遍：

fprintf('节点数: %d\n', size(mpc.bus, 1)); fprintf('支路数: %d\n', size(mpc.branch, 1)); fprintf('发电机数: %d\n', size(mpc.gen, 1));

如果数量都对不上，那大概率是自己写数据文件时少了行或者多了行，先在这里排查，比在潮流结果里排查效率高得多。

6. 从潮流到更多玩法：OPF与多机系统扩展

8机28节点系统的价值不止于潮流计算。把runpf换成runopf，就能做经济调度意义上的最优潮流。

6.1 最优潮流的快速上手

调用方式非常简单：

res_opf = runopf(mpc);

但要做好OPF，mpc.gencost必须填写完整，否则Matpower会按默认成本计算，结果可能与你预期的不符。gencost里定义的是发电机有功出力的成本曲线，单位一般换算成$/h。用二次多项式表示：

% 三台发电机为例：c2, c1, c0 分别为二次项、一次项、常数项 mpc.gencost = [ 2 0 0 3 0.005 20 50; 2 0 0 3 0.010 15 40; 2 0 0 3 0.002 25 60; ];

这里的MODEL=2表示多项式成本模型，NCOST=3表示有三个系数，从最高次数往常数项排列。不同版本对gencost列数的定义有细微差别，我的习惯是打开自带case30.m，把它的gencost直接复制过来，改成自己的机组参数，这样永远不会因为列数不一致报错。

6.2 静态安全分析和N-1校验

当系统搭好并跑通潮流后，我通常还会做一轮N-1静态安全校验：遍历每一条支路，让它退出运行，重新跑一次潮流，看是否出现电压越限或支路过载。

实现方式很简单：

n_branch = size(mpc.branch, 1); for i = 1:n_branch mpc_i = mpc; mpc_i.branch(i, 11) = 0; % 把第i条支路状态设为停运 res_i = runpf(mpc_i); if res_i.success == 0 fprintf('支路%d退出运行后潮流不收敛\n', i); end end

这个循环虽然朴素，但在8机28节点这种小系统上跑得很快，是做“系统薄弱环节分析”时非常实用的工具。

6.3 多机系统动态仿真的延伸

如果下一步要做暂态稳定或者小干扰稳定分析，8机28节点这套数据同样可以作为基础。常见做法是把构建好的mpc对象导出成PSS/E的raw文件，或者直接导入到PSAT、PST等动态仿真工具里。因为在潮流层面已经验证过系统是收敛、可运行的，导入到动态仿真时只需要关心动态参数（惯性时间常数、励磁系统模型、调速器参数）的填写，不需要再回头调潮流初始值。

6.4 代码包的组织习惯

最后说说代码组织。我习惯把“数据文件”和“计算脚本”分开，目录结构大致是这样的：

case28_8m.m % 系统数据文件，只定义mpc结构 main_powerflow.m % 潮流计算主程序 main_opf.m % 最优潮流主程序 static_security.m % N-1校验脚本 plot_results.m % 可视化脚本

数据文件里只有mpc的定义，不放任何计算逻辑。主程序通过loadcase加载数据，再调用Matpower的函数。这样做的好处是换一个系统数据文件，主程序一行都不用改，直接切换loadcase的参数就行。

在实际使用中，我最深的体会是：Matpower项目里最花时间的永远不是runpf这一行命令，而是把数据文件整理清楚。8机28节点这种自建系统，关键就是先把bus、gen、branch这三张表的结构吃透，后面的程序调用都是锦上添花。希望这篇分享能帮你少走点弯路。

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