简介：loadshed.zip是一份面向电力系统潮流计算与安全分析的MATLAB例程包，适合电力工程专业学生、研究人员及电网调度人员学习切负荷策略的编程实现。资源共18个文件，其中17个为m脚本或函数文件，1个为asv自动保存文件，压缩包整体大小仅18KB，结构紧凑、易读易用。该例程覆盖了电网建模、潮流求解、越限检测与切负荷优化等完整环节：先通过节点与支路参数建立系统模型，再使用牛顿-拉弗森或高斯-塞德尔迭代计算潮流分布；当发现节点电压或线路功率越限时，基于贪心、粒子群等算法智能选择切除负荷点，并反复迭代直至恢复安全运行。包内还包含负荷比例调节、支路故障概率等辅助脚本，便于进一步开展可靠性分析与紧急控制研究。目前已有196人学习/下载，对于想要掌握MATLAB在电力系统紧急控制中具体应用的学习者，具有较高的参考价值。
1. loadshed.zip 不是普通例程包，它是电网减载策略的 MATLAB 沙盘
先纠正一个预期：loadshed.zip 里的“例程”，和 stm32f407、esp32 那类嵌入式例程完全是两码事。解压后看不到寄存器配置，也不存在交叉编译，纯粹是一组.m脚本加一份电网数据文件。它解决的是甩负荷（load shedding）问题：发电机跳机、线路断开造成功率缺额时，在哪些节点、按什么顺序切除多少负荷，才能让系统频率止跌回升，同时尽量少切。下载它的人，多是电力系统方向的研究生、做微电网仿真的工程师，以及毕业设计选了低频减载的本科生。拿到包后别急着点 Run，先搞清楚两件事：电网数据文件在哪改，减载判据是谁在调用。
2. 用 MATLAB 跑通 loadshed 例程：目录结构、路径设置与工具箱检查
2.1 解压后先认四类文件
这类减载例程包解压后的目录结构大同小异，常见做法是主脚本放根目录，算法函数和数据文件混在一起，少数版本用 lib、data 子目录分开。动手运行前先把文件按职责归档一遍，比直接双击主脚本要少踩很多坑：
loadshed/ ├── main_loadshed.m % 主入口：潮流 → 故障 → 减载 → 画图 ├── case14.m % 节点/线路/负荷数据（有的包叫 bus_data.m） ├── pf_newton.m % 牛顿-拉夫逊潮流求解 ├── ufls_stages.m % 低频减载轮次逻辑 └── plot_freq.m % 结果可视化这个清单不是某个固定版本的忠实还原，而是这类例程的惯用骨架，文件命名可能不同，职责基本一致。打开主脚本和数据文件，先核对节点数、基准容量这两个量，后续所有计算都建立在它们上面。
| 文件类型 | 典型文件名 | 运行前要确认的点 |
|---|---|---|
| 主入口脚本 | main_loadshed.m、run_shedding.m | 头部是否写死了 case 数据文件名 |
| 电网数据 | case14.m、bus33.m | 节点数、基准容量 MVA 与主脚本一致 |
| 算法函数 | ufls_stages.m、opt_shed.m | 返回值个数与主脚本调用行匹配 |
| 绘图函数 | plot_freq.m | 是否依赖附加工具箱 |
| Simulink 模型 | shed_model.slx（部分版本带） | 本机 MATLAB 版本能否打开 |
2.2 最小可复现：解压、加路径、跑主脚本
先把 zip 解开，路径里不要出现空格和中文，这点在 Windows 上尤其容易踩。mkdir -p保证目录存在，unzip -d指定解压目标目录，ls -R递归列出全部文件用于核对：
mkdir -p ~/projects/loadshed unzip loadshed.zip -d ~/projects/loadshed cd ~/projects/loadshed ls -R接着打开 MATLAB，把主脚本所在目录连同全部子目录加入搜索路径。很多刚跟着 matlab 下载安装教程装完环境的人，第一次运行例程报 Undefined function or variable，就是漏了这一步：文件在磁盘上，但不在 MATLAB 的 path 里。
cd('C:\Users\你的用户名\projects\loadshed'); addpath(genpath(pwd)); % 当前目录及所有子目录加入搜索路径 savepath; % 写入 pathdef.m，重启后仍生效 main_loadshed; % 运行主例程addpath(genpath(pwd))里pwd返回当前目录，genpath递归找出全部子目录，addpath一次性挂到搜索路径最前面。savepath负责持久化，否则重启 MATLAB 后路径全部丢失；如果系统提示没有权限改写安装目录，就把前两行放进用户目录下的startup.m，每次启动自动执行，这是不污染系统目录的更干净做法。
提示：
savepath报权限错误时不要用管理员权限强行写 MATLAB 根目录，优先改用startup.m方案。
2.3 MATLAB 版本与工具箱：运行前先验三件事
减载例程的能力横跨三个工具箱：最优减载那段通常调intlinprog，属于 matlab 优化工具箱；如果包内带 Simulink 模型，需要 Simulink 本体；绘图里若出现fft或滤波函数，Signal Processing Toolbox 就成了隐藏依赖。命令窗口输入ver，逐项确认这三者在已安装列表里，比运行到一半报工具箱缺失再回头装省事。版本上，R2020a 之后intlinprog的求解器行为很稳定，R2016a 之前的版本对整数变量的预处理有差异，同一份数据可能得出不同的减载节点组合。如果机器上同时装了多个 MATLAB 版本，启动时留意命令行窗口顶部打印的版本号，别用旧版跑新格式的数据文件。
3. 例程核心算法：把减载写成整数规划，用 intlinprog 求解
3.1 为什么是整数规划而不是普通线性规划
最优减载要回答的问题很直接：已知故障后的功率缺额，求每个节点被切除的负荷量，让总切除量等于缺额，同时选中的节点组合代价最小。如果负荷可以无限细分，用linprog就能解；但实际减载靠断路器动作完成，负荷要么整段切除、要么不动，这是天然离散的 0/1 决策，正是intlinprog的适用范围。几十个节点的算例，求解时间通常在毫秒级，嵌进时域仿真循环里完全跑得动。例程包里的 opt_shed.m 或 load_shed_opt.m，去掉注释和画图代码后，核心骨架基本就是下面这十几行。
3.2 一个可直接替换的 opt_loadshed 函数
function [P_shed, status] = opt_loadshed(P_load, delta_P, w) % 输入: % P_load : n 维列向量，各节点有功负荷，p.u. % delta_P : 故障后的系统功率缺额，p.u. % w : n 维列向量，切负荷代价权重，省略时用负荷占比 % 输出: % P_shed : 各节点实际切除量，p.u. % status : 各节点是否被选中(1/0) n = length(P_load); if nargin < 3 w = P_load / sum(P_load); % 权重默认取负荷占比，负荷大切起来代价高 end f = [w; zeros(n, 1)]; % 目标：只最小化切负荷加权和 intcon = (n+1):(2*n); % 后 n 个变量是 0/1 选择标志 A = [eye(n), -diag(P_load)]; % P_shed(i) <= u(i)*P_load(i) b = zeros(n, 1); Aeq = [ones(1, n), zeros(1, n)]; % 总切除量等于功率缺额 beq = delta_P; lb = zeros(2*n, 1); ub = [P_load; ones(n, 1)]; opts = optimoptions('intlinprog', 'Display', 'off', ... 'IntegerTolerance', 1e-5); [x, ~, exitflag] = intlinprog(f, intcon, A, b, Aeq, beq, ... lb, ub, opts); if exitflag ~= 1 warning('intlinprog 返回 exitflag=%d，减载方案不可用', exitflag); end P_shed = x(1:n); status = x(n+1:2*n); end决策变量一共 2n 个：前 n 个是连续量P_shed(i)，后 n 个是 0/1 量status(i)。约束 A 展开后是P_shed(i) - status(i)*P_load(i) <= 0，含义是节点没被选中就不能切负荷，被选中时最多切掉该节点全部负荷。等式约束 Aeq 强制各节点切除量之和等于故障缺额。目标向量 f 前 n 项是代价权重，后 n 项置零——选择标志本身不产生代价，代价全部体现在被切掉的有功上。
参数里容易被略过的是IntegerTolerance，它控制 0/1 变量被判定为整数的误差窗口。默认 1e-5 对工程算例足够；调到 1e-8 会让求解时间成倍增加，而减载方案几乎没有变化。只有论文复现、要求结果严格可比较时，才需要把IntegerTolerance和随机种子一起固定。
3.3 三个必调参数：权重、缺额与单节点上界
| 参数 | 位置 | 典型取值范围 | 影响 |
|---|---|---|---|
w权重向量 | opt_loadshed入参 | 负荷占比 0~1，或人工重要度 | 决定先切哪些节点 |
delta_P功率缺额 | 主脚本故障设置处 | 系统总负荷的 5%~15% | 决定总切除量 |
ub上界 | 求解器边界 | ub(1:n)=P_load | 限制单节点最大切除量 |
3.3.1 权重那行代码的物理含义
w = P_load / sum(P_load)的意思是：负荷越大的节点，切除代价越高。实际调度里医院、炼钢负荷的代价和普通居民负荷完全不同，把 w 换成人工定义的重要度向量[0.1; 0.9; 0.5; ...]再归一化传入即可。只改这一处，最优减载的节点顺序立刻变化，这也是整套例程里最值得做灵敏度分析的旋钮。
4. 把例程改成自己的电网数据：CSV 负荷导入、UFLS 轮次与结果验证
4.1 用 readtable 把写死的负荷数组替换成 CSV 数据
多数 loadshed 例程把负荷写死在 case 文件里，比如 case14.m 的 Pd 列。工程上要做某条馈线的减载策略验证，负荷曲线通常来自能量管理系统导出的 CSV。常见做法是保留 case 文件的网络拓扑，只把负荷数组换成 CSV 内容：
T = readtable('load_profile_2025.csv'); % 两列：Time_s, P_MW t_load = T.Time_s; P_load_mw = T.P_MW; P_load_pu = P_load_mw / 100; % 基准容量 100 MVA，MW 转 p.u. load_system_total = 2.58; % 潮流数据里的系统总有功，p.u. if abs(sum(P_load_pu) - load_system_total) > 0.05 warning('CSV 总负荷与电网数据基准不一致，请检查基准容量'); endreadtable自动识别表头并生成 table 变量，T.Time_s直接取列，省掉了xlsread那类旧接口的文件格式兼容问题。这里最阴的坑是基准容量：case 文件里所有功率都是 p.u. 标幺值，CSV 里通常是 MW 有名值，不除以基准容量就混用，减载结果会整体偏离一个数量级。例程跑出离谱的电压或减载量，先查这一处换算。
4.2 低频减载（UFLS）轮次：阈值、延时与切除比例
另一条并行的技术路线是传统低频减载（UFLS），按频率阈值分级动作，例程里通常对应 ufls_stages.m。默认参数一般按电网稳定导则的典型值配置，拿到例程后第一件事就是和所在地区整定单核对：
| 轮次 | 频率阈值 (Hz) | 动作延时 (s) | 切除比例（占系统总负荷） |
|---|---|---|---|
| 第 1 轮 | 49.0 | 0.2 | 5% |
| 第 2 轮 | 48.8 | 0.3 | 7% |
| 第 3 轮 | 48.5 | 0.3 | 8% |
| 加速轮 | 49.2 且 df/dt < -0.3 | 0.1 | 4% |
function shed_P = ufls_logic(f_meas, df_dt, t_now, P_total) thr = [49.0 48.8 48.5]; % 各轮频率阈值，Hz frac = [0.05 0.07 0.08]; % 各轮切除比例 delay = [0.2 0.3 0.3 ]; % 各轮动作延时，s shed_P = 0; for k = 1:3 if f_meas < thr(k) && t_now >= delay(k) shed_P = shed_P + frac(k) * P_total; end end if f_meas < 49.2 && df_dt < -0.3 && t_now >= 0.1 shed_P = shed_P + 0.04 * P_total; % 加速轮，防频率快速塌落 end endufls_logic的输入是实时频率、频率变化率和当前仿真时间，输出是这一步应切的负荷量，主循环每步调用一次，累计量从系统功率平衡方程里扣除。改例程时最常见的错误是只改阈值不改延时：延时决定了轮次配合，第 1 轮动作后频率仍在跌，第 2 轮才该带着稍长的整定延时跟上。把阈值、比例、延时放进同一个数组而不是散落在主循环里，可读性和可维护性会好得多。加速轮看的是 df/dt 而不是频率绝对值，目的是在频率快速单调下跌时提前动作，代价是可能多切，实际工程里会单独统计加速轮动作次数。
提示：UFLS 整定值必须以当地电网调度部门下发的整定单为准，例程里的 49.0/48.8/48.5 只是教学典型值，直接用于现场口径会出大问题。
4.3 验证减载效果：三条曲线加一个 FFT 指标
跑完例程别只看一个减载总量。我一般把频率、切负荷累积量和最低节点电压画在同一张 figure 里，检查三件事：频率最低点是否高于首轮阈值、稳态频率是否回到额定附近、过切量是否小于预设上限。
t = t_sim; f = f_sim; % 主脚本输出的时域仿真结果 shed_cum = cumsum(P_shed_seq); % 各时刻切负荷累积量, MW figure('Name', '减载效果验证'); subplot(3,1,1); plot(t, f, 'LineWidth', 1.5); hold on; yline(49.0, '--r'); % 低频减载第 1 轮阈值 xlabel('时间 (s)'); ylabel('频率 (Hz)'); subplot(3,1,2); stairs(t, shed_cum, 'LineWidth', 1.2); xlabel('时间 (s)'); ylabel('累积切负荷 (MW)'); subplot(3,1,3); plot(t, V_min, 'LineWidth', 1.5); xlabel('时间 (s)'); ylabel('最低节点电压 (p.u.)');进一步确认减载后不存在振荡失稳风险，可以对频率瞬态段做 FFT，看机电振荡主频是否落在 0.5~2 Hz 之外。把仿真频率按固定步长写成 CSV，再用readtable读回来，Fs = 1/dt; spec = abs(fft(f_trace));就得到主振荡频率。这个流程和常说的 CSV 导入 MATLAB 做 FFT 仿真一致，只是分析对象从信号换成了频率轨迹。
5. 例程跑到一半报错：三个高频报错与多场景加速技巧
5.1 维度不一致：按数据文件到主脚本的顺序对三个数
修改 case 数据后最常见的报错是 Matrix dimensions must agree。比如把 14 节点换成 33 节点，只换了数据文件，主脚本里写死的 nb=14 没同步，潮流迭代到雅可比矩阵那一步就会崩。我一般按数据文件到主脚本的顺序核对节点数、线路数、基准容量三处，命令行whos看 bus 和 branch 矩阵的行数，通常一分钟定位。
5.2 intlinprog 返回负 exitflag 时先松哪条约束
| ExitFlag | 含义 | 先做检查 |
|---|---|---|
| 1 | 找到最优解 | 检查 P_shed 之和是否略大于缺额 |
| 0 | 达到迭代上限 | 增大 MaxTime 或减少整数变量 |
| -2 | 无可行解 | 确认 delta_P 小于最大单节点负荷 |
| -3 | 目标无界 | 检查 lb/ub 维度与 P_load 一致 |
修复顺序是先松 beq，允许 2% 的欠切或过切，再松单节点 ub，仍不行才去掉一半 0/1 变量改用连续模型。这个顺序保证每一步都清楚是哪条约束在卡方案。
5.3 用 parfor 把单场景例程改成灵敏度实验流水线
delta_list = 0.05 : 0.01 : 0.20; % 缺额比例序列 results = cell(length(delta_list), 1); parfor k = 1:length(delta_list) [P_shed, ~] = opt_loadshed(P_load, delta_list(k) * P_total, w); results{k} = [delta_list(k), sum(P_shed .* w)]; end writetable(array2table(cell2mat(results), ... 'VariableNames', {'deltaPU', 'shedCost'}), 'sensitivity.csv');parfor 要求 P_load、P_total、w 是只读广播变量，opt_loadshed 必须在搜索路径上且函数内不能写全局变量。把 w 替换成节点重要度的随机向量循环多次，用同一张表收集输出，就能画出减载策略的风险边界曲线。
本文还有配套的精品资源，点击获取
