简介：PSAT（电力系统分析工具箱）2.1.11版本是MATLAB环境下最新发布的电力系统分析工具箱，面向电力系统科研与工程人员，适用于潮流计算、连续潮流、最优潮流、小扰动稳定性分析、时域仿真以及自定义模型搭建等应用场景。该压缩包共包含1713个文件，以1432个m格式源代码为核心，辅以121个bmp、45个jpg图像资源、33个mdl仿真模型及少量演示视频，整体大小约19.23MB，目录结构清晰，便于按模块检索学习。目前已有1753人浏览下载，为PSAT时隔三年的重要版本更新。该版本支持MATLAB R2019a及之前版本，新增功能丰富，源代码完全公开，用户可按研究兴趣修改扩展；同时附带多种电力软件格式转换接口，可配合MATLAB控制、信号处理等工具箱，将控制科学新思想与电力系统分析有机结合，适合开展上千节点规模的系统仿真与算法验证。解压后除psat主程序目录外，还包含辅助文件夹，便于跨平台查看与整理。

1. 电力系统仿真工具箱 PSAT：潮流计算与暂态稳定分析一次装齐

做电力系统分析的人对 PSAT 应该不陌生，但很多初学者卡在第一步——版本匹配和安装路径上。这份 psat-2.1.11-mat.zip 是目前较新的 MATLAB 版本打包，解压即用，涵盖了潮流计算、连续潮流、最优潮流、暂态稳定和小信号分析等完整功能链。我拆完这份资源后最直接的感受是：它把 toolbox 里散落的函数和 GUI 统一封装好了，省去了从旧版本手动迁移数据格式的麻烦。适合正在做配电网潮流、微电网稳定性分析的研究生，也适合需要快速搭建 IEEE 标准节点测试模型的从业者。

2. 安装与环境配置：PSAT 在 MATLAB 里的正确打开方式

2.1 路径添加与初始化顺序

PSAT 的运行机制依赖 MATLAB 的搜索路径。解压后第一件事不是双击 .m 文件，而是先把整个 psat 目录加入路径。我习惯用命令行的方式操作，比 setpath 对话框更可控。

% 将 PSAT 主目录及其子目录全部加入搜索路径 addpath(genpath('C:\Users\yourname\psat-2.1.11')); % 保存路径配置，避免每次重开 MATLAB 都要重新添加 savepath; % 清空工作区中的旧变量和缓存函数，这一步容易漏 clear all; close all; clc;

addpath的genpath参数会递归处理所有子目录，包括@class目录下的类定义文件。PSAT 2.1.11 引入了面向对象的元件模型，如果只用addpath不加genpath，运行d_psat时会报类未定义的错误。clear all在这里不是走过场——如果之前加载过低版本的 PSAT 函数，MATLAB 的缓存机制会优先命中旧函数，导致新版本不生效。

验证安装是否成功，直接输入：

d_psat

如果弹出 PSAT 的主界面窗口，说明路径和类文件都正常。这一步卡住的话，多半是 MATLAB 版本低于 R2016a，类定义语法不兼容。

2.2 数据文件格式与模型导入

PSAT 的模型数据有两种接入方式：一是通过 GUI 里的File -> Load读取.m格式的数据文件；二是用loadfile函数直接在命令行加载。数据文件结构是结构体数组，每个字段对应一类元件。

% 在命令行中直接加载 psat 格式的母线数据文件 [mpc, success] = loadfile('case9.m'); % 如果返回的 success 为 0，多半是文件路径不对或格式不兼容 if success == 0 error('数据文件加载失败，请检查 case9.m 内容'); end % 查看母线数据表前两行，确认 bus 类型列是否有效 disp(mpc.Bus(1:2, :));

mpc.Bus的第二列是母线类型标识：1 表示 PQ 节点、2 表示 PV 节点、3 表示平衡节点。很多数据文件的这部分是从旧版 PSAT 或直接手写转换来的，数值偏移一位就很常见。我一般加载后先检查size(mpc.Bus, 2) >= 13，因为新版格式的母线表至少有 13 列，少了说明版本转换不完整。loadfile返回的success标志位能帮忙快速定位问题，比用exist判断文件是否存在更可靠。

2.3 仿真参数初始化与求解器选择

PSAT 用psat结构体管理全局仿真参数，不用像别的工具箱那样逐个 set 参数。常用的是通过命令行设置容差和迭代上限。

% 设置潮流计算的收敛容差为 1e-8，迭代上限 50 次 psat = setparam(psat, 'tolerance', 1e-8); psat = setparam(psat, 'maxiter', 50); % 选择 Newton-Raphson 求解器，比默认的高斯-赛德尔更适合大系统 psat = setparam(psat, 'Alg', 1); % 执行潮流计算，nf 为计算结果结构体 nf = runpf(psat);

Alg参数有 1、2、3 三个取值，分别对应 Newton-Raphson、快速分解法和 Gauss-Seidel。IEEE 30 节点以下规模用哪个差别不大，超过 300 节点时 Newton-Raphson 收敛速度快一倍不止，但内存占用也高一些。setparam第二个参数直接用点号访问结构体里嵌着的字段，写错字段名会静默失败，不报错但也不生效，所以设置完建议disp(psat.PowerFlow.MaxIter)确认一下。

3. 潮流计算实战：从搭建母线模型到结果提取

3.1 用代码构建单线图模型

没有现成数据文件时，手动搭建模型是绕不开的。PSAT 支持通过addbus、addline等函数逐条添加元件，比直接改 .m 数据文件更直观且不易出错。

% 初始化一个新的数据模型 psat = initpsat; % 添加三组母线：一组平衡节点，两组 PQ 节点 psat = addbus(psat, struct('Bus', 1, 'Vmin', 0.9, 'Vmax', 1.1, ... 'type', 3, 'P', 0, 'Q', 0)); psat = addbus(psat, struct('Bus', 2, 'Vmin', 0.9, 'Vmax', 1.1, ... 'type', 1, 'P', 0.5, 'Q', 0.2)); psat = addbus(psat, struct('Bus', 3, 'Vmin', 0.9, 'Vmax', 1.1, ... 'type', 1, 'P', 0.3, 'Q', 0.1)); % 添加两条输电线路，电阻电抗以标幺值表示 psat = addline(psat, struct('bus1', 1, 'bus2', 2, 'r', 0.01, 'x', 0.10)); psat = addline(psat, struct('bus1', 2, 'bus2', 3, 'r', 0.02, 'x', 0.15));

initpsat不只是清空变量，它会把所有元件计数器归零，并重置默认参数，避免上一次运行残留的模型混入新模型。addbus中的type字段就是前面说的母线类型，3 为平衡节点表示电压幅值和相角已知，1 为 PQ 节点表示有功无功已知。P、Q字段的单位是标幺值，基准功率在psat.BaseMVA里统一设置。如果写了P = 500而基准是 100 MVA，相当于注入 5 倍的功率，潮流计算直接不收敛，这个单位换算错位的坑后面专门说。

3.2 潮流计算运行与结果结构体解析

模型搭好后运行runpf，结果会返回一个结构体，关键是读懂里面的字段含义。很多人看结果只看有没有收敛标志，忽略了详细数据。

% 运行潮流计算并解析结果 results = runpf(psat); % 检查收敛状态，0 表示未收敛 if results.convergence == 0 error('潮流计算未收敛，请调整初值或检查母线参数'); end % 提取节点电压幅值和相角 V_mag = results.Bus.V; % 电压幅值，标幺值向量 V_ang = results.Bus.theta; % 相角，弧度 % 计算某条线路的传输功率 line_power = results.Line.Pf; % 输出母线 2 的电压幅值 disp(['母线 2 电压幅值: ', num2str(V_mag(2))]);

results.Bus是结构体，里面的V和theta按母线编号 1、2、3 的顺序对应，不是按物理位置排列的。results.Line.Pf是线路首端有功功率，负值表示功率方向与定义方向相反。这里有个容易混淆的点：results.Bus.V在潮流计算里是电压幅值，但在暂态稳定计算里会被替换成包含时间和动态响应的数组，类型完全不同。提前判断isfield(results, 'Bus')再访问字段，能避免后期在两种模式间来回切换时踩坑。

3.3 潮流结果的可视化验证

数值算完了，最好把电压分布画出来做直观验证。PSAT 没内置绘图函数，但结果结构体和 MATLAB 原生的绘图接口无缝对接。

% 画节点电压幅值柱状图 figure; bar(V_mag); xlabel('母线编号'); ylabel('电压幅值 (pu)'); ylim([0.9, 1.1]); % 设置坐标范围，方便异常值一眼看出 grid on; % 标出电压越限的母线 idx_over = find(V_mag > 1.05 | V_mag < 0.95); hold on; plot(idx_over, V_mag(idx_over), 'rx', 'MarkerSize', 10); % 输出越限提示 if ~isempty(idx_over) disp(['警告：母线 ', num2str(idx_over'), ' 电压越限']); end

ylim([0.9, 1.1])这个范围设置是有讲究的，很多变电站的实际允许电压范围就是 0.95 到 1.05，画图时直接按正常运行范围显示，越限点用红叉标记后一眼就能发现。idx_over的查找逻辑里用了|而不是||，因为需要同时比较两个边界条件，写||会导致只判断一个条件。

4. 暂态稳定与最优潮流：PSAT 的高级分析模块实战

4.1 暂态稳定仿真的时间步长与故障设置

暂态稳定模块是 PSAT 的招牌功能。设置短路故障时主要关注两个参数：故障时刻和清除时刻。PSAT 通过addfault函数定义，时间单位是秒。

% 在母线 2 上设置三相短路故障 psat = addfault(psat, struct('bus', 2, 'type', 1, ... 't_start', 0.1, 't_end', 0.2, 'R', 0.001)); % 设置仿真总时长与步长 psat = setparam(psat, 'tspan', [0, 1.0]); psat = setparam(psat, 'dt', 0.005); % 运行暂态稳定仿真 td_out = runtd(psat);

addfault的type字段：1 表示三相短路、2 表示单相接地、3 表示两相相间。R是故障电阻，三相短路时理想值为 0，但设成 0 会导致导纳矩阵奇异，实际算例一般取 0.0001 到 0.001 之间。dt步长直接决定仿真精度，0.005 秒适合 1 秒内的故障切除过程模拟，如果观察对象是低频振荡（0.5-2 Hz），步长可以放宽到 0.01 秒，兼顾速度。t_end - t_start要大于至少一个步长，否则故障瞬间不会触发。

4.2 求解器在暂态模式下的切换逻辑

PSAT 在暂态稳定模式下会自动切换求解器，不再使用潮流计算时的 Newton-Raphson，而是采用隐式梯形积分法。这个切换是自动的，但有一个前提——潮流计算必须先收敛。

% 先做潮流计算获取初始运行点 results_pf = runpf(psat); % 确认收敛后再跑暂态分析 if results_pf.convergence == 1 % 在这之后调用 runtd 才会使用正确的初始值 td_out = runtd(psat); else error('潮流未收敛，暂态初始值不准确'); end

这个顺序是关键。很多人直接跳过潮流计算调用runtd，PSAT 会用默认电压初值（通常都是 1.0 pu）来计算，如果实际运行点偏离 1.0 pu 较远，暂态仿真的响应曲线会带着一个虚假的初始振荡，等于是从错误起点开始跑。弧光接地等不对称故障还要依赖潮流计算后的负序和零序网络参数，没有收敛结果这些参数不确定，故障引起的非周期分量会被错误估计。

4.3 最优潮流：目标函数选择与约束配置

OPF 模块适合处理经济调度和有功无功协调优化问题。PSAT 支持多种目标函数，通过参数切换。

% 设置目标函数为发电成本最小化 psat = setparam(psat, 'OPF', 'objective', 'generation_cost'); % 设置各发电机成本系数（二次成本函数 a + bP + cP^2） psat = addgen(psat, struct('bus', 1, 'Pg', 0.5, ... 'Pmax', 1.0, 'Pmin', 0.1, 'cost', [0.02, 2.1, 0])); % 运行最优潮流计算 opf_out = runopf(psat);

cost数组的三位分别对应二次项、一次项和常数项系数。Pmax和Pmin是有功出力上下限，单位是标幺值。OPF 的约束系统默认包含母线电压上下限、线路热极限、发电机无功极限，这些都是从潮流计算模型中继承的。跑 OPF 失败最常见的情况是可行域为空——发电机出力下限设得比负荷需求还高，或者线路热极限设得太小，导致没有任何运行点能满足所有约束。这种情况 PSAT 不会给出明确的错误信息，只会在opf_out.convergence返回 0，排查时会花不少时间。

4.4 连续潮流与 PV 曲线绘制

连续潮流（CPF）用于分析电压稳定性，核心输出是 PV 曲线。PSAT 通过runpf(psat, 'cpf')触发，也可以用独立函数。

% 设置逐步增加负荷方向 psat = setparam(psat, 'CPF', 'base', 1.0); psat = setparam(psat, 'CPF', 'step', 0.05); % 运行连续潮流计算，最多迭代 20 步 cpf_out = runcpf(psat, 20); % 提取节点电压对负荷增长的响应 lambda = cpf_out.Lambda; % 负荷增长因子 V_at_node = cpf_out.V(:, 3); % 母线 3 电压随负荷变化 % 绘制 PV 曲线 figure; plot(lambda, V_at_node, 'b-', 'LineWidth', 1.5); xlabel('负荷增长因子 \lambda'); ylabel('母线 3 电压幅值 (pu)'); grid on;

连续潮流里最关键的就是步长选择。step设成 0.05 适合 20 步左右逼近鞍结分岔点，如果步长过大（比如 0.2），在接近分岔点时会越过转折点导致发散。PSAT 的 CPF 模块内置了步长自适应逻辑，但初始步长如果超出极限范围，自适应来不及调整就会失败。拐点附近电压变化剧烈，绘制的 PV 曲线有尖锐的弯曲时，适当减小步长后曲线会变得更平滑，可以清楚看到鼻尖点对应的电压稳定极限。

5. 避坑避坑：PSAT 安装与仿真的高频翻车现场

5.1 现象：MATLAB 提示函数名冲突

命令行执行d_psat，系统报错说d_psat与 MATLAB 自带的某个函数库重名，或者直接报 undefined function。检查后发现which d_psat定位到的路径不是刚解压的新目录。

原因有两点可能：一是旧版 PSAT 安装目录没删除干净，还在搜索路径的前端；二是 MATLAB 的matlab.toolbox路径优先级比用户目录高，新加的路径排在后面没生效。

解决方法是执行pathtool，手动把 PSAT 目录向上移动到路径列表的顶部，然后savepath。如果同时存在多个版本的 PSAT 目录，只保留一个有效的，其余全部从路径中移除，避免which命中错误版本。

5.2 现象：潮流计算迭代不收敛

runpf返回的convergence字段一直是 0，但数据文件明明是标准的 IEEE 节点模型。这个问题在 2.1.11 版本里仍然常见。

排查路径有三条：先看容差是不是设太严了，tolerance设成 1e-10 时对 300 节点以上的系统很苛刻，放宽到 1e-8 一般就能收敛；再检查 P、Q 的单位，某些手动整理的数据文件里 P 的单位可能是 MW 而不是标幺值，直接写进模型等于放大了基准值，试试用psat.BaseMVA核实一下基准功率；最后检查平衡节点的电压初值，V字段设成 0 时会影响到雅可比矩阵的初始生成，改成 1.0 之后再运行。

5.3 现象：连续潮流在分岔点卡死

连续潮流在逼近鞍结分岔点时迭代不收敛，无论怎么减小步长都无解，可能一直循环到maxiter上限。

根本原因是接近分岔点时雅可比矩阵接近奇异，数值计算本身条件数很差，单纯减小步长不能根治。解决方式是用预测-校正机制的参数调整：把预测步设为切线预测而不是线性预测，在 PSAT 的 CPF 参数中找到pred参数，从默认值改为 2（对应切线预测法）。这种预测法在分岔点附近追踪能力明显增强，但计算量会大一些，每步需要多解一次线性方程组。

5.4 现象：暂态稳定仿真的波形有明显畸变

故障切除后的电压波形出现高频振荡，看起来像数值发散，但是又存在明显的周期。关掉 PSAT 的变步长控制，强制使用定步长后情况有改变。

这类问题通常是dt与系统最小时间常数不匹配。PSAT 暂态模型里的发电机励磁绕组时间常数可能只有 0.05 秒，如果dt用了 0.02 秒，采样间隔接近时间常数的二分之一，就会造成一定程度上的数值混叠。解决方式是检查和系统内最快的动态环节时间常数对比，确保dt不超过该时间常数的四分之一。同步电机的励磁模型时间常数通常在 1 到 3 秒，但调速器和 PSS 的时间常数可能小到 0.01 秒，这类环节存在时要相应缩小步长。

5.5 现象：OPF 结果明显不合理

最优潮流算完了，但某台发电机的出力到了边界值、同时另一台发电机接近零出力，与凭借物理直觉的判断不符，甚至出现功率倒送。

排查重点放在约束是否激活上。PSAT 的 OPF 默认会激活线路热约束，如果线路容量设置过小，优化器可能牺牲发电机经济性来迁就线路容量。在psat.Line.Con里检查线路rate字段，确认各线路的热极限值和实际功率方向。另一种可能是电压约束过于严格，母线电压下限太接近 1.0，导致无功在局部区域强制平衡，有功经济调度因此被扭曲。建议把电压约束放宽到 0.95 到 1.05 后再看结果是否回归正常。

6. 进阶技巧：用脚本批量分析与自定义 PSAT 元件模型

PSAT 最被低估的不是计算能力，而是它面向二次开发的接口。熟练之后不需要每次都在 GUI 里手动点击，完全可以写一套脚本，自动化处理批量的仿真任务。

% 批量扫描不同负荷水平下的潮流结果 load_levels = [0.5, 0.8, 1.0, 1.2, 1.4]; % 负荷系数 voltage_results = zeros(length(load_levels), 3); % 每个负荷系数下重新运行潮流，记录低压母线电压 for k = 1:length(load_levels) % 用 psat 自带的负荷系数调整函数修改负荷 psat = setload(psat, 'scaling', load_levels(k)); % 运行潮流并暂存结果 r_k = runpf(psat); % 取母线电压幅值，暂存到矩阵 voltage_results(k, :) = r_k.Bus.V(:)'; end % 统一输出对比表格 disp(array2table(voltage_results, ... 'VariableNames', {'母线1', '母线2', '母线3'}, ... 'RowNames', string(load_levels)));

setload是 PSAT 里常被忽略的功能接口，它不改母线的基准结构，只调整所有 PQ 节点的负荷功率，相当于给整个系统乘一个系数。批量算电压时比逐母线修改 P、Q 高效很多。跑完一个场景后setload的 scaling 会自动重置为 1.0，需要回到当前场景关键校验，防止上一次的负荷放大系数污染下一个算例。

自定义元件需要改 PSAT 的库文件时，建议复制原文件再改，不要直接动@class目录下的原始类。PSAT 元件的核心逻辑集中在ybus_update和power_residual两个方法里，前者补充该元件对导纳矩阵的贡献，后者补充对节点注入功率不平衡量的影响。具体说，要在现有 PV 节点模型上增加一个本地电压控制器的话，先在@PV目录下加一个方法additional_control，然后修改power_residual在计算增量时调用新增方法。接上自定义负载模型时，利用 admittance matrix 和 Jacobian 的稀疏结构化修改能保持系统整体的稀疏性，避免大规模矩阵变稠。

我的习惯是用脚本化工作流跑任何结构大于 30 节点的模型：一、加载数据；二、setparam设好收敛容差和迭代上限；三、跑基准潮流并把 CV 值存档；四、批量扫描运行点并自动标注越限母线。从那以后我每次跑暂态稳定前都强制走一遍潮流初值校验，不然初始振荡分不清是真实响应还是默认初值带来的误差。希望帮到你。

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