简介：PSAT（Power System Analysis Toolbox）2.1.3版是面向电力系统工程师、高校师生及科研人员的开源MATLAB分析工具箱，用于在MATLAB环境中完成发电机、变压器、线路与负荷等元件的建模，并开展稳态潮流、动态仿真与控制策略评估等工作。压缩包共含1616个文件，约1.67MB，以1351个.m源码文件为主体，辅以bmp、jpg、png等界面与示例图片、30个mdl仿真模型，以及ieee2psat、psse2psat、digsilent2psat等一批数据格式转换脚本，便于导入不同来源的电网数据。资源已有546人学习下载，适合需要快速搭建电力系统仿真环境、研究算法实现或进行课程设计的读者参考。解压后可直接集成到MATLAB中调用，借助源码与示例模型理解PSAT内部结构，为二次开发与排错提供便利。
1. 从一份压缩包说起：PSAT 2.1.3 到底能帮你算什么
如果你在电力系统仿真这个圈子里待过一阵，大概率会在某个时刻遇到一个后缀为.zip的压缩包，名字里带着psat-2.1.3-mat这样的字样。它不是一个普通的资料包，而是一套基于 MATLAB 环境运行的开源电力系统分析工具箱。很多人第一次拿到它的时候，面对满屏的.m文件和几个示例数据文件，完全不知道从哪里下手。更常见的场景是：导师丢过来一句“你先用 PSAT 跑个潮流看看”，然后你就开始了漫长的环境配置和报错排查之旅。
PSAT 全称 Power System Analysis Toolbox，它的核心价值在于把潮流计算、连续潮流、最优潮流、小信号稳定性分析、时域仿真这些电力系统分析里的经典任务，封装成了一套可以在 MATLAB 里直接调用的函数和脚本。和那些商业软件相比，它的优势是源码开放、算法透明、方便二次开发；劣势也很明显——文档零散、版本兼容性玄学、报错信息不够友好。2.1.3 这个版本在圈子里流传较广，主要是因为它对 MATLAB 版本的兼容区间相对宽，而且自带了不少经典的测试系统算例。
这篇文章面向的是手里已经拿到这个压缩包、或者正准备去找这个包的人。我会从解压后的目录结构讲起，把环境配置、第一个潮流算例、参数怎么调、报错怎么查这几个环节串起来。目标很明确：让你在本地把 PSAT 跑起来，并且知道每一步在干什么。至于那些需要深度定制的场景，比如自己写模型、改算法，我也会在最后一章给出可操作的切入点。
2. 解压之后先别急着运行：目录结构与运行链路拆解
2.1 压缩包里到底有什么
把psat-2.1.3-mat.zip解压到一个纯英文路径下，你会看到类似这样的结构：
psat-2.1.3-mat/ ├── psat.m % 主入口脚本 ├── init.m % 初始化脚本 ├── pfsolver.m % 潮流求解器 ├── cpfsolver.m % 连续潮流求解器 ├── opfsolver.m % 最优潮流求解器 ├── fm_*/ % 各类分析模块目录 ├── lib/ % 底层函数库 ├── models/ % 元件模型定义 ├── data/ % 测试系统数据 └── docs/ % 零散文档这里有几个关键点需要先建立认知。psat.m是整个工具箱的入口，但它不是一个可以直接双击运行的脚本，而是需要在 MATLAB 命令窗口里调用的。init.m负责把 PSAT 的各个子目录添加到 MATLAB 的搜索路径中，同时设置一些全局变量。如果你跳过init.m直接运行psat.m，大概率会遇到“函数未定义”的报错。
data/目录里存放的是测试算例，常见的有 IEEE 14 节点、30 节点、57 节点等标准系统。这些数据文件通常是.m格式，里面用结构体数组的方式定义了母线、支路、发电机等元件参数。models/目录下是各种电力系统元件的数学模型实现，比如同步发电机、励磁系统、负荷模型等。如果你后续要做自定义建模，这个目录是必须啃的。
2.2 环境配置：MATLAB 版本与路径设置
PSAT 2.1.3 对 MATLAB 版本是有要求的。根据我的踩坑经验，R2014a 到 R2018a 之间的版本兼容性最好。太新的版本（比如 R2020b 之后）可能会因为 MATLAB 语言本身的变更导致部分函数报错，尤其是涉及图形界面和符号计算的部分。如果你用的是较新版本，也不是完全跑不了，但需要做好心理准备，可能要手动修几个文件。
配置步骤如下：
% 第一步：将 PSAT 根目录添加到 MATLAB 路径 % 假设解压路径为 D:\psat-2.1.3-mat addpath(genpath('D:\psat-2.1.3-mat')); % 第二步：运行初始化脚本 % 这一步会设置全局变量、加载默认配置 init; % 第三步：验证路径是否配置成功 % 如果返回非空，说明 PSAT 相关函数已可调用 which psataddpath(genpath(...))的作用是把根目录及其所有子目录一次性加入搜索路径。genpath会递归遍历子文件夹，省去手动逐个添加的麻烦。init脚本执行后会设置一些全局变量，比如global PQ、global PV之类的，这些变量在后续的潮流计算中会被频繁引用。
注意：每次重启 MATLAB 后都需要重新执行
addpath和init，除非你把这两步写进startup.m文件。我一般会在项目目录下单独建一个startup.m，把路径配置和初始化命令放进去，这样每次打开 MATLAB 会自动执行。
2.3 第一个算例：用 IEEE 14 节点系统跑通潮流
环境配好之后，最直接的验证方式就是跑一个标准算例。PSAT 自带的数据目录里有 IEEE 14 节点的数据文件，通常命名为d_014_mdl.m或类似的名字。运行流程如下：
% 清除工作区，避免旧变量干扰 clear all; close all; clc; % 重新初始化 PSAT addpath(genpath('D:\psat-2.1.3-mat')); init; % 加载 IEEE 14 节点测试系统 % 具体文件名以 data 目录下实际存在的为准 run('D:\psat-2.1.3-mat\data\d_014_mdl.m'); % 执行潮流计算 % pfsolver 是 PSAT 的潮流求解器入口 pfsolver; % 查看结果 % 母线电压幅值和相角会存储在 BUS 结构体中 disp(BUS);这里需要解释几个关键点。run命令执行数据文件后，会在工作区里生成一系列结构体变量，比如BUS（母线数据）、LINE（线路数据）、GEN（发电机数据）。pfsolver会读取这些变量，组装雅可比矩阵，然后用牛顿-拉夫逊法迭代求解。求解完成后，BUS结构体里的电压列会被更新为计算结果。
如果你在pfsolver这一步遇到报错，最常见的原因是数据文件里的元件参数和求解器期望的格式不匹配。比如某些版本的 PSAT 要求负荷必须用恒功率模型表示，而数据文件里可能混用了恒阻抗模型。这时候需要检查models/目录下对应的负荷模型定义，确认数据文件里的模型类型字段是否正确。
2.4 结果怎么看：电压、功率、收敛曲线
潮流计算完成后，结果主要分布在几个结构体里。BUS结构体包含每个节点的电压幅值、相角、注入有功和无功。LINE结构体包含每条支路的潮流分布。GEN结构体包含发电机的出力情况。
% 提取母线电压幅值 V_mag = BUS(:, 2); % 第二列通常为电压幅值 % 提取母线电压相角（弧度） V_ang = BUS(:, 3); % 第三列通常为相角 % 绘制电压幅值分布图 figure; bar(V_mag); xlabel('母线编号'); ylabel('电压幅值 (p.u.)'); title('IEEE 14 节点系统电压分布'); grid on;电压幅值是判断潮流结果是否合理的第一指标。对于 IEEE 14 节点系统，大部分母线的电压应该在 0.95 到 1.05 p.u. 之间。如果出现某个母线电压低于 0.8 或高于 1.2，要么是数据文件有问题，要么是系统本身处于重载状态。相角的单位是弧度，如果需要转换成角度，乘以180/pi即可。
收敛曲线是另一个重要信息。pfsolver在迭代过程中会输出每次迭代的误差范数。如果误差在 3 到 5 次迭代内降到 1e-6 以下，说明收敛良好。如果迭代次数超过 10 次还没收敛，或者误差出现震荡，那就要检查系统是否存在孤岛、参数是否合理、或者初始值是否给得太离谱。
3. 从潮流到连续潮流：参数扫描与电压稳定分析
3.1 连续潮流在解决什么问题
普通潮流计算有一个固有局限：当系统负荷逐渐增加时，潮流方程的解会在某个临界点附近消失。这个临界点对应的负荷水平就是系统的最大输电能力，也叫电压稳定极限。普通潮流在接近这个点时，雅可比矩阵会趋于奇异，迭代不收敛。连续潮流通过引入一个连续参数，把负荷增长过程和潮流方程的解曲线一起追踪出来，从而绕过奇异点，画出完整的 PV 曲线。
PSAT 里的连续潮流求解器是cpfsolver。它的输入和普通潮流类似，但需要额外指定一个“增长方向”，也就是哪些负荷在增加、哪些发电机在出力。常见做法是选择一个区域的负荷按比例增长，同时指定平衡机来承担功率缺额。
3.2 配置连续潮流的三个关键参数
连续潮流的配置主要集中在几个参数上，这些参数决定了追踪的步长、方向和终止条件。
% 设置连续潮流参数 % 这些参数通常存储在 CPF 结构体中 CPF.step = 0.01; % 初始步长 CPF.step_min = 1e-5; % 最小步长 CPF.step_max = 0.1; % 最大步长 CPF.lambda_min = 0; % 最小负荷增长因子 CPF.lambda_max = 3; % 最大负荷增长因子 CPF.adapt = 1; % 启用自适应步长 % 指定负荷增长方向 % 假设所有 PQ 节点的有功和无功按相同比例增长 CPF.loaddir = 'all'; % 执行连续潮流 cpfsolver;CPF.step是初始步长，太大会导致曲线追踪不光滑，太小会显著增加计算时间。我一般从 0.01 开始试，如果曲线在临界点附近出现明显折角，就减小到 0.005。CPF.adapt控制是否启用自适应步长，开启后求解器会根据曲率自动调整步长，在平坦区域走大步，在临界点附近走小步。CPF.lambda_max是负荷增长因子的上限，设置得太小可能还没到临界点就停了，设置得太大则会在临界点之后继续追踪无意义的解。
3.3 读懂 PV 曲线：鼻点、裕度与工程含义
连续潮流跑完之后，最直观的结果就是 PV 曲线。横轴是负荷增长因子 lambda，纵轴是某个关键母线的电压幅值。曲线通常呈现一个“鼻形”：随着 lambda 增加，电压逐渐下降，到达鼻点后，电压会继续下降但 lambda 开始减小，这部分对应的是不稳定解分支。
% 提取 PV 曲线数据 % CPF 结构体中通常包含 lambda 和电压的追踪结果 lambda = CPF.lambda; V_critical = CPF.V(:, 1); % 假设关注第一个母线的电压 % 绘制 PV 曲线 figure; plot(lambda, V_critical, 'b-', 'LineWidth', 1.5); xlabel('负荷增长因子 \lambda'); ylabel('关键母线电压 (p.u.)'); title('PV 曲线与电压稳定裕度'); grid on; % 找到鼻点对应的 lambda [lambda_max, idx] = max(lambda); V_at_nose = V_critical(idx); fprintf('最大负荷增长因子: %.4f\n', lambda_max); fprintf('鼻点电压: %.4f p.u.\n', V_at_nose);鼻点对应的 lambda 值就是系统的电压稳定裕度指标。如果 lambda_max 是 1.5，意味着系统还能承受 50% 的负荷增长。这个数字在规划阶段很有参考价值。需要注意的是，PV 曲线的形状和增长方向的选择密切相关。如果只增加某个区域的负荷，鼻点位置会和全网等比例增长时不一样。所以做稳定分析时，增长方向的设定要结合实际场景，不能随便选一个就下结论。
3.4 和普通潮流结果的交叉验证
连续潮流跑完之后，建议用普通潮流在 lambda=0 的位置做一次交叉验证。如果两者得到的电压分布不一致，说明连续潮流的初始点设置有问题。常见原因是连续潮流在启动时对系统状态做了一些假设，比如把某些 PV 节点转换成 PQ 节点，而普通潮流没有做这个转换。
% 在 lambda=0 处用普通潮流验证 % 先将负荷恢复到基准值 % 然后运行 pfsolver pfsolver; % 对比两种方法得到的电压幅值 % 如果差异超过 1e-4，需要检查模型定义 V_diff = abs(BUS(:, 2) - V_critical(1)); max_diff = max(V_diff); fprintf('最大电压差异: %.6f p.u.\n', max_diff);这个验证步骤看起来多余，但在实际项目里能帮你省下大量排查时间。我曾经遇到过一次连续潮流结果异常的情况，追查了半天才发现是数据文件里某个变压器的变比在连续潮流初始化时被错误地重置了。如果当时先做了交叉验证，五分钟就能定位到问题。
4. 避坑指南：PSAT 2.1.3 最常见的五类报错与排查
4.1 报错“Undefined function or variable”的三种可能
现象：运行pfsolver或cpfsolver时，MATLAB 命令窗口报错，提示某个函数或变量未定义。
原因一：路径没配好。init.m没有执行，或者addpath添加的路径不完整。解决方法是重新执行addpath(genpath(...))和init，然后用which命令确认关键函数是否在搜索路径中。
原因二：数据文件没加载。pfsolver依赖工作区里的BUS、LINE等结构体变量。如果忘记run数据文件，求解器找不到输入数据。解决方法是在运行求解器之前，先确认工作区里存在这些变量。
原因三：MATLAB 版本不兼容。某些在旧版本中可用的函数在新版本中被移除或改名。解决方法是查看报错信息里具体是哪个函数，然后在 MATLAB 文档里搜索替代函数，手动修改 PSAT 源码中对应的调用。
4.2 潮流不收敛：从雅可比矩阵到初始值
现象：pfsolver迭代超过 20 次仍未收敛，或者误差范数在某个值附近震荡。
原因：最常见的是系统存在孤岛节点，也就是某个母线没有和任何发电机或平衡节点相连。另一个常见原因是初始电压设置不合理，比如把所有母线电压都设成 1.0，但实际系统里某些母线电压偏离较大。
解决方法：先用graph或sparse检查网络连通性，确认没有孤立节点。然后尝试用平启动（flat start）之外的初始值，比如用上一次潮流的结果作为这一次的初始值。如果还是不行，检查负荷模型是否合理，特别是恒功率负荷在低压下可能导致雅可比矩阵奇异。
% 检查网络连通性 % 假设 LINE 结构体中第一列和第二列是首末节点编号 nbus = size(BUS, 1); adj_matrix = sparse(LINE(:, 1), LINE(:, 2), 1, nbus, nbus); % 如果某个节点的度为 0，说明是孤岛 degree = sum(adj_matrix, 2); isolated = find(degree == 0); if ~isempty(isolated) fprintf('发现孤岛节点: %s\n', num2str(isolated')); end4.3 连续潮流在鼻点附近步长崩溃
现象：连续潮流追踪到接近鼻点时，步长急剧减小，计算时间暴增，甚至卡死。
原因：鼻点附近雅可比矩阵接近奇异，自适应步长策略会不断减小步长以维持收敛。如果CPF.step_min设置得太小，求解器会在一个点上反复迭代。
解决方法：适当增大CPF.step_min，比如从 1e-5 改成 1e-4。同时检查CPF.adapt的参数，有些版本的 PSAT 允许设置步长调整的激进程度。如果还是不行，可以手动在鼻点附近切换成固定小步长，跳过自适应逻辑。
4.4 图形界面打不开或显示异常
现象：运行psat.m后，图形界面窗口空白、按钮无响应，或者报错提示某个 GUI 函数不存在。
原因：PSAT 的图形界面依赖 MATLAB 的 GUIDE 框架，较新版本的 MATLAB 对 GUIDE 的支持有变化。另外，如果系统缺少某些 Java 组件，界面也可能无法正常渲染。
解决方法：优先使用命令行方式调用求解器，绕过图形界面。如果必须用界面，尝试在 MATLAB 启动时加上-nojvm参数关闭 Java 虚拟机，或者换用较旧的 MATLAB 版本。我个人的习惯是全程命令行操作，图形界面只在需要快速查看结果时偶尔用一下。
4.5 数据文件格式不匹配导致的静默错误
现象：潮流计算正常收敛，但结果明显不合理，比如某些母线电压接近 0 或发电机出力远超容量。
原因：数据文件里的元件参数顺序和求解器期望的顺序不一致。PSAT 的不同版本对数据格式的要求有细微差别，比如母线类型编码、负荷模型标识等。这种错误不会导致报错，但会让计算结果完全失真。
解决方法：对照docs/目录下的格式说明（如果有的话），或者直接阅读pfsolver.m里读取数据的代码段，确认每个字段的含义。更稳妥的做法是先用自带的测试算例跑通，然后以测试算例的数据格式为模板，修改成自己的系统数据。
5. 进阶用法：自定义模型与批量算例的自动化
5.1 在 PSAT 里挂一个自定义负荷模型
PSAT 的模型定义集中在models/目录下。每个模型通常对应一个.m文件，里面定义了该元件的微分方程或代数方程。如果你想加一个自定义负荷模型，比如电压相关的 ZIP 负荷，可以按照现有模型的模板来写。
% 自定义 ZIP 负荷模型示例 % 文件名建议放在 models/ 目录下，命名为 zip_load.m function [P, Q] = zip_load(V, P0, Q0, alpha_p, alpha_q) % V: 实际电压幅值 (p.u.) % P0, Q0: 额定电压下的有功和无功 % alpha_p, alpha_q: ZIP 系数 % 恒阻抗部分 P_imp = P0 * (V^2); Q_imp = Q0 * (V^2); % 恒电流部分 P_cur = P0 * V; Q_cur = Q0 * V; % 恒功率部分 P_pow = P0; Q_pow = Q0; % 按系数加权 P = alpha_p(1)*P_imp + alpha_p(2)*P_cur + alpha_p(3)*P_pow; Q = alpha_q(1)*Q_imp + alpha_q(2)*Q_cur + alpha_q(3)*Q_pow; end写完之后，需要在数据文件里把对应负荷的模型类型改成这个新模型的标识，并且在init.m或相关配置里注册这个模型。具体注册方式取决于 PSAT 版本的模型管理机制，常见做法是在models/目录下建一个索引文件，把模型名称和文件名对应起来。
5.2 批量运行多个算例的脚本模板
在实际项目里，经常需要对多个场景批量跑潮流或连续潮流。手动一个个加载数据文件效率太低，用脚本自动化是更合理的做法。
% 批量运行算例的脚本模板 clear all; close all; clc; % 配置 PSAT 路径 addpath(genpath('D:\psat-2.1.3-mat')); init; % 定义要运行的算例列表 case_files = {'d_014_mdl.m', 'd_030_mdl.m', 'd_057_mdl.m'}; results = cell(length(case_files), 1); for i = 1:length(case_files) fprintf('正在运行算例: %s\n', case_files{i}); % 清除旧数据 clear BUS LINE GEN; % 加载新算例 run(fullfile('D:\psat-2.1.3-mat\data', case_files{i})); % 运行潮流 pfsolver; % 保存关键结果 results{i}.V = BUS(:, 2); results{i}.angle = BUS(:, 3); results{i}.case_name = case_files{i}; fprintf('算例 %s 完成，最低电压: %.4f p.u.\n', ... case_files{i}, min(BUS(:, 2))); end % 汇总结果 for i = 1:length(results) fprintf('%s: 最低电压 %.4f, 最高电压 %.4f\n', ... results{i}.case_name, ... min(results{i}.V), ... max(results{i}.V)); end这个模板的核心逻辑是：用循环遍历算例文件列表，每次清除旧变量后加载新数据，运行求解器，然后把关键结果存到 cell 数组里。clear BUS LINE GEN这一步很重要，如果不加，上一次算例的变量会残留，导致求解器读到错误的数据。
5.3 用 PSAT 做小信号稳定性分析的切入点
除了潮流和连续潮流，PSAT 还提供了小信号稳定性分析功能。它的基本原理是在稳态运行点附近对系统方程做线性化，然后计算状态矩阵的特征值。特征值的实部对应阻尼，虚部对应振荡频率。
% 在潮流收敛的基础上做小信号分析 % 先跑一次潮流，得到稳态运行点 pfsolver; % 然后调用小信号分析模块 % 具体函数名以 PSAT 版本为准，常见的是 sssa 或 small_signal sssa; % 查看特征值 % 结果通常存储在 SSSA 结构体中 eigenvalues = SSSA.eigenvalues; figure; plot(real(eigenvalues), imag(eigenvalues), 'bo'); xlabel('实部 (阻尼)'); ylabel('虚部 (频率)'); title('特征值分布'); grid on;特征值分布图里，如果所有特征值的实部都小于 0，说明系统在小扰动下是稳定的。如果有实部大于 0 的特征值，对应的振荡模式就是失稳模式。通过参与因子分析，可以进一步定位是哪些发电机或负荷对失稳模式贡献最大。
5.4 一个我反复使用的检查习惯
每次跑完一个不熟悉的算例，我都会做三件事：第一，检查潮流收敛的迭代次数，超过 8 次就去看数据；第二，检查最低电压母线的位置，如果出现在负荷中心之外的地方，多半是网络参数有问题；第三，把结果和同类型的标准算例做对比，偏差超过 10% 就回头查模型。
这个习惯帮我省下了不少后悔药。有一次跑一个改造过的 30 节点系统，潮流收敛了但最低电压出现在一个偏远节点上，后来发现是一条线路的阻抗参数多写了一个零。如果只看收敛标志，这个问题就漏过去了。
希望帮到你。
本文还有配套的精品资源，点击获取
