simulink实现标准IEEE33配电网系统50HZ将各节点数据统计起来输出到工作区再matlab中跑出某时刻节点电压分布适合用于观察某时刻节点电压变化情况。 #特别是当用于接入双馈风机时用powergui无法进行潮流计算通过此方法能过很好的解决此问题。 有参考文献。最近在搞双馈风机接入配电网的仿真发现Powergui遇到这种动态元件直接歇菜——潮流算不动啊折腾半天搞了个野路子用Simulink实时导出节点数据再用Matlab脚本直接画电压分布图。实测能动态观察某时刻电压变化特别适合风机接入的场景分享给同样被Powergui坑过的兄弟。第一步搭个IEEE33配电网模型先扒拉出Simulink的Power System Components照着IEEE33节点参数画拓扑。重点注意三点母线电压设为12.66kV/50Hz别用60Hz会挨导师骂负载数据按文献[1]的P/Q值挂载在每个节点接一个Voltage Measurement模块输出勾选Per Unit归一化方便后续处理% 节点负载设置示例节点2 Load2 100e3; % 单位kW QLoad2 60e3; % 单位kVar set_param(IEEE33_model/Load2, ActivePower, num2str(Load2)); set_param(IEEE33_model/Load2, ReactivePower, num2str(QLoad2));第二步把数据怼进工作区Powergui不能算那就暴力导出所有数据在每个Voltage Measurement后挂To Workspace模块变量名填NodeVoltage别用默认的simout会覆盖保存格式选Structure With Time采样时间设0.01秒50Hz对应20ms周期别太小否则卡死跑完仿真后工作区里会有一坨结构体。这时候写个脚本提取某时刻的数据% 提取t5秒时的节点电压 t_target 5; [~, idx] min(abs(NodeVoltage.time - t_target)); VoltageData zeros(33,1); for i 1:33 VoltageData(i) NodeVoltage.signals(i).values(idx); end第三步画个炫酷的电压分布图数据到手直接上可视化figure(Color,[1 1 1]); bar(VoltageData); xlabel(节点编号); ylabel(电压 (p.u.)); title([t,num2str(t_target),s时节点电压分布]); grid on; ylim([0.9 1.05]); % 电压越限一目了然如果看到某个节点电压狂跌比如节点18大概率是风机接入后引起的局部电压塌陷。这时候可以对比有无风机的情况一眼看出影响范围。为什么不用Powergui官方给的Powergui做潮流需要静态模型但双馈风机这种旋转电机变流器的组合自带动态特性一仿真就报错Algebraic Loop。而我们的方法本质上是时间步进数据抓取避开传统潮流计算直接硬核记录瞬时值——虽然土但管用啊参考文献simulink实现标准IEEE33配电网系统50HZ将各节点数据统计起来输出到工作区再matlab中跑出某时刻节点电压分布适合用于观察某时刻节点电压变化情况。 #特别是当用于接入双馈风机时用powergui无法进行潮流计算通过此方法能过很好的解决此问题。 有参考文献。[1] 张某某, 含分布式电源的配电网潮流分析, 《电力系统自动化》, 2020.[2] Simulink User Guide, MathWorks, 2023.