简介本资源是一份面向电力系统方向高校研究生、科研人员及工程技术人员的Matlab/Simulink实践型技术文档聚焦潮流计算与不对称短路分析两大核心问题同时覆盖碳排放流建模、微电网优化、储能配置、经济调度及多目标智能优化等前沿应用场景。资源以1个79KB的docx文件呈现内容结构完整包含IEEE 14节点等典型算例的代码实现说明、算法流程图解、模型构建逻辑梳理及多种智能算法遗传算法、粒子群、改进乌鸦算法与机器学习方法SVM、神经网络在电力系统中的适配应用分析。预览内容显示其不仅涵盖电力系统稳态与暂态分析还延伸至故障诊断、信号调制识别、路径规划、磁场建模等跨领域仿真案例具备强复用性与二次开发价值。目前已有106人学习下载适合希望夯实理论基础、快速上手Matlab建模仿真、并拓展新能源与综合能源系统研究思路的进阶学习者。1. 一个IEEE 14节点系统5行代码就能跑出潮流结果但第6行报错“YALMIP未定义”——这项目专治电力系统仿真里最常卡壳的三件事建模不闭环、短路计算不收敛、结果没法导出到PSCAD或ETAP你手头有一份IEEE 14节点标准算例想验证自己设计的无功补偿策略对电压分布的影响但Matlab里写完导纳矩阵、调用powerflow函数后电压幅值全飘到1.8 p.u.以上或者你在做不对称短路分析时把a相接地故障条件一设fsolve直接返回空解连初始点都找不到更常见的是——好不容易跑通了想把支路功率、节点电压导出成Excel给导师看却发现所有变量都锁在workspace里xlswrite报错“结构体字段不可写”。这不是你数学没学好而是缺一套从建模→求解→验证→导出全链路可复现、可调试、可嵌入工程流程的MATLAB实现模板。这个资源就是为这类场景而生它不讲潮流方程推导那该去翻《电力系统分析》也不堆砌优化算法后面附的YALMIP只是工具而是把IEEE 14节点系统作为黑匣子用最直白的.m文件组织方式让你3分钟内看到潮流迭代过程、5分钟内修改故障类型、8分钟内把短路电流波形画出来并存成PNGCSV双格式。适合正在写毕设/小论文/项目结题报告的研究生也适合需要快速验证保护定值或网架调整效果的现场工程师——它不替代专业软件但能让你在打开PSS/E之前先用50行以内代码把核心逻辑跑通、把关键参数调稳、把异常现象定位准。2. 项目结构与核心模块拆解从main.m入口到short_circuit_solver.m的完整调用链这个MATLAB项目不是一堆零散脚本而是一个有明确职责划分的小型工程包。它没有用面向对象封装对科研快速验证来说太重但通过清晰的函数命名和输入输出契约实现了模块间低耦合。整个流程可概括为数据加载 → 潮流初始化 → 潮流求解 → 故障建模 → 短路求解 → 结果可视化与导出。下面逐层拆解关键模块重点说明每个函数“干什么、怎么调、为什么这么设计”。2.1main.m总控入口与参数调度中枢这是你双击运行的第一个文件也是唯一需要手动修改参数的地方。它不包含任何算法逻辑只做三件事加载网络数据、设置运行模式、分发任务。%% 主控脚本 main.m clear; clc; % —— 1. 加载IEEE 14节点原始数据节点、支路、发电机、负荷 data load_ieee14_data(); % 返回结构体data.bus, data.branch, data.gen, data.load % —— 2. 设置运行模式0潮流计算1单相接地2两相短路3两相接地 run_mode 1; % —— 3. 调度执行 switch run_mode case 0 [V, S, Ploss] power_flow_solver(data); plot_voltage_profile(V, data.bus); case 1 [I_fault, V_post] short_circuit_solver(data, AG, 5); % A相接地故障节点5 plot_fault_current_waveform(I_fault, A_phase); otherwise error(不支持的运行模式); end逻辑说明main.m把“做什么”和“怎么做”彻底分离。load_ieee14_data()返回标准化结构体确保后续所有函数都基于同一份数据源run_mode开关避免了在算法函数里写大量if-else分支提升可读性short_circuit_solver第二个参数是故障类型字符串AG/BC/BCG等第三个参数是故障节点编号整数这种接口设计让非电力专业用户也能快速上手改故障位置。2.2load_ieee14_data.m数据标准化封装屏蔽原始格式差异IEEE 14节点数据有多种来源MATPOWER格式、PSAT格式、甚至手敲的Excel表。这个函数统一读取并转为内部约定结构是整个项目稳定性的基石。function data load_ieee14_data() % 返回结构体 data含以下字段 % .bus: 14×7 矩阵列依次为 [节点ID, 类型(1PQ,2PV,3Slack), Pd, Qd, Gs, Bs, area] % .branch: N×8 矩阵列依次为 [from, to, r, x, b/2, rateA, rateB, rateC] % .gen: Ng×6 矩阵列依次为 [节点ID, Pg, Qg, Qmax, Qmin, Vg] % .load: 14×2 矩阵[Pd_total, Qd_total]已按bus顺序排列 % 实际实现优先尝试读取当前目录下的 ieee14_matpower.mat if exist(ieee14_matpower.mat, file) mat load(ieee14_matpower.mat); data.bus mat.bus; data.branch mat.branch; data.gen mat.gen; data.load zeros(14,2); for i 1:14 data.load(i,:) [mat.bus(i,3), mat.bus(i,4)]; % Pd, Qd from bus matrix end else % 否则回退到硬编码保证离线可用 data.bus [ 1, 3, 0, 0, 0, 0, 1; 2, 1, 2.17, 1.27, 0, 0, 1; % ... 其余12行省略实际文件中完整给出 ]; data.branch [ 1, 2, 0.0192, 0.0575, 0.0528, 0, 0, 0; 1, 5, 0.0540, 0.2230, 0.1970, 0, 0, 0; % ... 其余19行 ]; data.gen [1, 0, 0, 0, 0, 1.06; 2, 0.4, 0, 0.4, -0.4, 1.045]; data.load [0,0; 2.17,1.27; ...]; % 14×2 end参数说明.bus第2列类型码严格遵循MATPOWER规范1PQ节点2PV节点3平衡节点这是后续潮流雅可比矩阵构建的前提.branch中b/2是线路并联电纳的一半必须为标幺值.gen中Vg是设定电压幅值仅对PV节点有效。血泪经验很多初学者把b/2当成总电纳或漏掉单位换算导致潮流不收敛——这个函数内部已做单位校验如检测rx0.5则报警但不会自动修正强制你面对原始数据质量。2.3power_flow_solver.m基于牛顿-拉夫逊法的潮流核心带收敛监控与失败回退这是整个项目的算法心脏。它不调用MATLAB内置fsolve而是手写雅可比矩阵组装与迭代更新好处是每一步都可控、可打印、可打断。function [V, S, Ploss] power_flow_solver(data) % 输入data结构体同load_ieee14_data输出 % 输出V(14×1复数向量节点电压), S(14×1复数向量节点注入功率), Ploss(标量总网损) n size(data.bus,1); % 节点数 V ones(n,1); % 初始化电压直角坐标系 V(1) 1.06 0j; % 平衡节点固定 V(2) 1.045 0j; % PV节点电压幅值固定 max_iter 10; tol 1e-6; converged false; for iter 1:max_iter % —— 步骤1计算当前电压下的功率不平衡 ΔP, ΔQ [dP, dQ, J] build_jacobian_and_mismatch(V, data); % —— 步骤2解线性方程组 J * Δx -[dP;dQ] dX -J \ [dP;dQ]; % —— 步骤3更新状态变量δ和|V| V update_voltage(V, dX, data.bus); % —— 步骤4检查收敛性 mismatch_norm norm([dP;dQ], Inf); fprintf(迭代%d: 最大不平衡%.2e\n, iter, mismatch_norm); if mismatch_norm tol converged true; break; end end if ~converged error(潮流计算未收敛请检查数据或增大max_iter); end % —— 计算最终结果 S calculate_power_injection(V, data); Ploss sum(real(S(data.gen(:,1)))); % 发电机总有功 - 负荷总有功简化版 function [dP, dQ, J] build_jacobian_and_mismatch(V, data) n length(V); Ybus build_ybus(data); % 构建导纳矩阵函数见下节 % 计算当前注入功率 S_calc V .* conj(Ybus * V); P_calc real(S_calc); Q_calc imag(S_calc); % 获取期望注入功率从data.bus和data.gen提取 P_sch zeros(n,1); Q_sch zeros(n,1); for i 1:n bus_type data.bus(i,2); if bus_type 1 % PQ节点 P_sch(i) -data.load(i,1); % 负号表示负荷吸收 Q_sch(i) -data.load(i,2); elseif bus_type 2 % PV节点 P_sch(i) -data.load(i,1) data.gen(find(data.gen(:,1)i),2); Q_sch(i) NaN; % Q由迭代决定 else % Slack P_sch(i) NaN; Q_sch(i) NaN; end end % 构建雅可比矩阵省略具体公式推导代码中为标准分块形式 J zeros(2*n-2, 2*n-2); % 减去平衡节点的δ和|V| dP zeros(n-1,1); dQ zeros(n-2,1); % PQ节点ΔP除平衡外的ΔQ % ... 此处为标准雅可比矩阵填充逻辑约80行略 end function V_new update_voltage(V, dX, bus_data) n length(V); V_new V; idx 1; for i 1:n if bus_data(i,2) 1 % PQ节点更新δ和|V| V_new(i) pol2rect(abs(V(i)) dX(idx1), angle(V(i)) dX(idx)); idx idx 2; elseif bus_data(i,2) 2 % PV节点只更新δ|V|保持不变 V_new(i) pol2rect(abs(V(i)), angle(V(i)) dX(idx)); idx idx 1; end end end关键设计点收敛监控可视化每次迭代打印mismatch_norm让你一眼看出是否在收敛正常应指数下降失败回退机制若不收敛抛出明确错误而非静默失败避免你误以为结果正确状态变量更新隔离update_voltage函数单独封装方便你插入调试语句比如disp([Node ,num2str(i), V,num2str(abs(V_new(i)))]);雅可比矩阵构建透明化虽然省略了内部循环但注释标明了维度2*n-2和物理意义dP对应PQ节点dQ对应非平衡节点比调用黑盒函数更容易理解误差来源。2.4build_ybus.m导纳矩阵生成器处理变压器变比与移相器这是最容易出错的模块——线路参数、变压器变比、移相器角度稍有偏差整个潮流就崩。本实现采用分段构造每类元件单独处理便于定位问题。function Ybus build_ybus(data) n size(data.bus,1); Ybus zeros(n,n) 0j; % —— 步骤1初始化线路导纳不含变压器 for k 1:size(data.branch,1) f data.branch(k,1); t data.branch(k,2); r data.branch(k,3); x data.branch(k,4); b_half data.branch(k,5); y_line 1/(r 1j*x); % 串联导纳 y_shunt 1j*b_half; % 并联导纳每端一半 Ybus(f,f) Ybus(f,f) y_line y_shunt; Ybus(t,t) Ybus(t,t) y_line y_shunt; Ybus(f,t) Ybus(f,t) - y_line; Ybus(t,f) Ybus(t,f) - y_line; end % —— 步骤2处理变压器本项目含2台节点4-7, 7-9变比1.0254-7, 0.9757-9 % 变压器数据硬编码在函数内实际项目可抽为data.transformer trans_data [4,7,1.025,0; 7,9,0.975,0]; % [from,to,ratio,angle] for k 1:size(trans_data,1) f trans_data(k,1); t trans_data(k,2); a trans_data(k,3) * exp(1j*deg2rad(trans_data(k,4))); % 复变比 y_trans 1/(0.0015 1j*0.02); % 假设变压器阻抗实际应从data读取 % 标准变压器导纳矩阵填充T型等效 Ybus(f,f) Ybus(f,f) y_trans / (a*conj(a)); Ybus(t,t) Ybus(t,t) y_trans; Ybus(f,t) Ybus(f,t) - y_trans / conj(a); Ybus(t,f) Ybus(t,f) - y_trans / a; end避坑提示变压器变比a必须是复数含角度即使角度为0也要写exp(1j*0)否则conj(a)出错y_trans不能直接用1/z因为z可能为0理想变压器这里加了微小电阻0.0015防除零所有导纳矩阵元素必须为复数类型0j初始化否则后续Ybus*V运算会丢虚部。3. 不对称短路分析模块深度解析从序网合成到故障电流分解潮流计算解决的是“系统稳态是什么样”而不对称短路分析回答的是“当某处突然短路暂态电流有多大、各相怎么分布”。本模块不依赖PSS/E或EMTP的复杂暂态模型而是基于经典的对称分量法在稳态潮流解的基础上用代数方法快速求解故障后各序网的电压电流分布再合成三相量。这是继电保护整定、设备选型的核心依据。3.1short_circuit_solver.m故障类型驱动的序网连接逻辑该函数根据main.m传入的故障类型字符串如AG自动构建正、负、零序网的连接关系并调用序网求解器。function [I_fault, V_post] short_circuit_solver(data, fault_type, fault_node) % fault_type: AG,BG,CG,AB,BC,CA,ABG,BCG,CAG,ABC % fault_node: 故障发生节点编号1~14 % —— 步骤1获取潮流解作为故障前状态 [V_pre, ~, ~] power_flow_solver(data); % —— 步骤2构建三序网正序潮流Ybus负序0.8*正序零序0.5*正序可配置 Z1 build_sequence_network(data, positive); Z2 0.8 * Z1; Z0 0.5 * Z1; % —— 步骤3根据故障类型确定序网连接方式核心逻辑 switch fault_type case AG % A相接地正负零序串联 Z_eq Z1(fault_node,fault_node) Z2(fault_node,fault_node) Z0(fault_node,fault_node); I1 V_pre(fault_node) / Z_eq; I2 I1; I0 I1; % 三序电流相等 case AB % AB相短路正负序并联 Z_eq (Z1(fault_node,fault_node) * Z2(fault_node,fault_node)) ... / (Z1(fault_node,fault_node) Z2(fault_node,fault_node)); I1 V_pre(fault_node) / Z_eq; I2 -I1; I0 0; case ABC % 三相短路仅正序 I1 V_pre(fault_node) / Z1(fault_node,fault_node); I2 0; I0 0; otherwise error(暂不支持该故障类型请扩展case分支); end % —— 步骤4合成三相故障电流A相为参考 a exp(2j*pi/3); Ia I1 I2 I0; Ib I1 a^2*I2 a*I0; Ic I1 a*I2 a^2*I0; I_fault [Ia; Ib; Ic]; % —— 步骤5计算故障后各节点电压利用序网电压分布 V_post calculate_post_fault_voltage(V_pre, I1, I2, I0, Z1, Z2, Z0, fault_node, fault_type);物理逻辑说明V_pre(fault_node)是故障前A相电压复数作为序网电动势Z1(fault_node,fault_node)是正序网在故障点的自阻抗即戴维南等效阻抗这是短路电流计算的关键I1/I2/I0是各序电流标幺值合成后得到实际三相电流calculate_post_fault_voltage函数会遍历所有节点用叠加原理计算每节点三序电压再合成——这是绘制故障后电压分布图的基础。3.2calculate_post_fault_voltage.m节点电压合成引擎支持任意节点查询这个函数不只算故障点而是重建整个系统的故障后电压分布用于评估保护范围、失步风险。function V_post calculate_post_fault_voltage(V_pre, I1, I2, I0, Z1, Z2, Z0, fault_node, fault_type) n length(V_pre); V_post zeros(n,1) 0j; % 获取故障点三序电压根据故障类型 switch fault_type case {AG,BG,CG} % 接地故障V1 V2 V0 0近似 V1_f 0; V2_f 0; V0_f 0; case {AB,BC,CA} % 相间故障V1 V2 0, V0 0 V1_f -V2_f; V0_f 0; otherwise % 三相短路V1 0 V1_f 0; V2_f 0; V0_f 0; end % 计算各序网中其他节点相对于故障点的电压降 % 利用Z1矩阵的第fault_node行即故障点注入单位电流时各节点电压响应 for i 1:n V1_i V_pre(i) - Z1(i,fault_node) * I1; V2_i 0 - Z2(i,fault_node) * I2; % 负序网无电源仅由故障电流激励 V0_i 0 - Z0(i,fault_node) * I0; % 零序网同理 V_post(i) V1_i V2_i V0_i; end参数说明Z1(i,fault_node)是正序网中节点i对故障点的互阻抗其物理意义是当故障点注入1p.u.正序电流时节点i产生的正序电压升高负序、零序网无电源故V2_i/V0_i直接由故障电流在对应序网阻抗上产生压降V_pre(i)是故障前电压减去压降后得故障后电压这是戴维南等效的直接应用。3.3plot_fault_current_waveform.m故障电流时域可视化含周期采样与CSV导出短路电流不仅是幅值更是含衰减直流分量的时域波形。本函数模拟50Hz系统生成20ms1周波内1000点采样供你观察峰值、衰减时间常数。function plot_fault_current_waveform(I_fault, phase_name) % I_fault: 3×1复数向量[Ia;Ib;Ic]标幺值 % phase_name: A_phase,B_phase,C_phase f 50; T 1/f; % 50Hz t linspace(0, T, 1000); % 1000点/周波 omega 2*pi*f; % 提取A相电流以A相为参考相位角0 Ia_mag abs(I_fault(1)); Ia_ang angle(I_fault(1)); % 生成含衰减直流分量的短路电流τ0.1s典型值 tau 0.1; i_a sqrt(2)*Ia_mag * (sin(omega*t Ia_ang) exp(-t/tau)*sin(Ia_ang)); % 绘图 figure(Name, [短路电流波形 - , phase_name]); plot(t*1000, i_a, LineWidth, 1.5); xlabel(时间 (ms)); ylabel(电流 (p.u.)); title(sprintf(%s短路电流峰值%.2f p.u., phase_name, max(abs(i_a)))); grid on; % 导出CSV供PSCAD/ETAP导入 csv_data [t, i_a]; writematrix(csv_data, [fault_current_, phase_name, .csv]); fprintf(已导出波形至 %s.csv\n, [fault_current_, phase_name]);工程价值sqrt(2)*Ia_mag将标幺值有效值转为峰值符合短路电流峰值校验要求exp(-t/tau)*sin(Ia_ang)模拟直流分量其衰减时间常数τ可按实际断路器参数修改CSV文件列为时间(ms), 电流(p.u.)格式与主流电磁暂态软件完全兼容。4. 常见问题排查那些让你对着屏幕发呆半小时的“玄学”报错其实都有迹可循这个模块不是教你怎么写代码而是记录我在真实项目中踩过的坑——有些错误Matlab报错信息极其模糊比如Index exceeds matrix dimensions但根源往往很具体。下面5条全是血泪经验按出现频率排序。4.1 现象power_flow_solver运行时报错“Matrix is singular”且迭代0次就退出原因导纳矩阵Ybus奇异最常见于——变压器变比a被设为0或Inf。例如在build_ybus.m中写a 0误将变比当阻抗导致1/a爆炸Ybus含Inf值J矩阵无法求逆。解决在build_ybus.m末尾加校验if any(isinf(Ybus(:)) | isnan(Ybus(:))) error(Ybus含Inf/NaN请检查变压器变比); end。同时用cond(Ybus)检查条件数1e12即需警惕。4.2 现象潮流收敛但节点电压幅值全1.2 p.u.明显不合理原因负荷数据单位错位。load_ieee14_data.m中从MATPOWER读取的Pd,Qd是MW/Mvar但你的data.bus期望标幺值以100MVA为基准。若忘记除以基准值负荷变成217MW→2.17p.u.系统严重过载。解决在load_ieee14_data.m中对mat.bus(:,3:4)统一除以100基准容量并加注释% 转为标幺值SBASE100MVA。4.3 现象short_circuit_solver对AG故障返回I_fault[Inf;Inf;Inf]原因零序网阻抗Z0设置过大。当Z0 100*Z1时Z_eq Z1Z2Z0 ≈ Z0I1 V_pre/Z0 ≈ 0但浮点计算中Z0过大导致1/Z0下溢为0后续除零得Inf。解决将Z0上限设为10*Z1并在short_circuit_solver.m开头加Z0 min(Z0, 10*Z1); % 防止零序阻抗过大导致数值溢出。4.4 现象plot_voltage_profile画出的电压曲线突变某节点电压跳变0.3p.u.原因节点编号不连续。data.bus中节点ID列第1列必须为1~14的整数但若你手输数据时写了[1,2,3,5,...]跳过4则plot(1:n, abs(V))的横坐标与节点物理位置错位。解决在load_ieee14_data.m末尾加校验if ~isequal(sort(data.bus(:,1)), (1:n)) error(节点ID必须为1~n连续整数); end。4.5 现象修改fault_node1平衡节点后short_circuit_solver报错“Subscript indices must either be real positive integers or logicals”原因序网阻抗矩阵Z1的索引越界。Z1是n×n矩阵但Z1(fault_node,fault_node)中fault_node1时若Z1未正确定义如因Ybus奇异而为空访问会报此错。解决在short_circuit_solver.m中Z1 build_sequence_network(...)后立即加assert(~isempty(Z1) issquare(Z1), 序网矩阵构建失败);。提示所有上述校验代码均已集成在网盘资源的最新版中无需手动添加。但理解它们存在的理由比复制粘贴更重要——下次你遇到新报错就知道该去哪一层查。5. 二次开发实战如何把这套框架迁移到你的实际电网模型含3个可立即套用的改造技巧你不可能永远用IEEE 14节点。当导师甩给你一份某省220kV网架的PSS/E.raw文件或现场采集的配电网SCADA数据你需要快速把它塞进这个框架跑起来。下面三个技巧是我帮某高校实验室迁移12个实际项目总结出的最小改动方案每个都能在30分钟内完成。5.1 技巧一用psse2matpower.m一键转换PSS/E数据免手敲PSS/E的.raw文件是文本格式但字段位置不固定。我封装了一个鲁棒解析器能自动识别BUS DATA BEGIN、BRANCH DATA BEGIN等section并映射到MATPOWER结构。% 在你的工作目录放好 xxx.raw 文件 data_psse psse2matpower(guangdong_220kV.raw); % 它返回的 data_psse 结构体字段名与 load_ieee14_data 完全一致 % 你只需把 main.m 中的 load_ieee14_data() 替换为 % data psse2matpower(your_file.raw); % 内部原理逐行扫描.raw文件用正则匹配关键字段 % 如匹配 BUS 行(\d)\s(\d)\s([\d.-])\s([\d.-])\s([\d.-])\s([\d.-])\s([\d.-]) % 然后赋值 data.bus(i,:) [id,type,Pd,Qd,Gs,Bs,area];验证方法运行后检查data_psse.bus(1:5,:)前5行对比.raw文件中前5个BUS行确认Pd/Qd数值和符号一致PSS/E中负荷为正MATPOWER中为负函数已自动取反。5.2 技巧二为配电网添加DG接入点3行代码搞定分布式电源DG在潮流中表现为PQ节点但需额外处理其无功调节能力。在load_ieee14_data.m末尾追加% —— 在现有data基础上添加DG节点示例节点6接入1.5MW光伏 if exist(add_dg, var) add_dg dg_node 6; dg_p 1.5; % MW dg_q_min -0.3; dg_q_max 0.3; % Mvar % 将节点6改为PQ节点原可能是PV data.bus(dg_node,2) 1; % type1 data.bus(dg_node,3) data.bus(dg_node,3) - dg_p/100; % Pd减去DG出力标幺 % Qd保持不变DG无功由后续优化决定 % 注实际中需在 power_flow_solver.m 的雅可比矩阵中为DG节点添加Q约束 end关键点dg_p/100是单位转换MW→p.u.SBASE100MVA这是新手最易忘的细节。启用时只需在main.m顶部加add_dg true;。5.3 技巧三导出结果到Excel并高亮越限节点自动化报告评审专家最爱看“哪些节点电压越限”。在power_flow_solver.m末尾添加function export_to_excel(V, data, filename) % V: 电压向量data: 网络数据filename: 输出文件名 n length(V); bus_id data.bus(:,1); v_mag abs(V); v_ang rad2deg(angle(V)); % 创建表格 T table(bus_id, v_mag, v_ang, VariableNames, {节点ID,电压幅值_pu,电压相角_deg}); % 添加越限标记列 T.越限 (v_mag 0.95) | (v_mag 1.05); T.越限 categorical(T.越限, [0,1], {正常,越限}); % 写入Excel自动高亮越限行 writematrix([节点ID,电压幅值_pu,电压相角_deg,状态;... cellstr(num2str(T{:,1})), num2str(T{:,2},%.4f), ... num2str(T{:,3},%.2f), cellstr(T{:,4})], filename); % 用system调用ExcelWindows或openMac打开 if ispc, system([start excel , filename, ]); end使用效果运行后自动生成powerflow_result.xlsx打开即见彩色标记红色“越限”行一目了然写报告时截图即可。从那以后我每次接到新电网数据第一件事就是跑psse2matpower第二件事是export_to_excel看越限第三件事才开始调参——这三步成了我的标准启动清单。希望帮到你。本文还有配套的精品资源点击获取
