最近在搞电力系统暂态稳定分析，拿经典的3机9节点系统练手。试了两种方法：Matlab编程解微分方程和Simulink时域仿真，发现这俩配合起来用挺有意思

3机9节点系统暂态稳定Matlab编程/simulink仿真 1.Matlab编程计算摇摆曲线，得到3机9节点系统中3台发电机的功角曲线以及转速曲线，通过分析各发电机之间的功角差和转速差来分析系统暂态稳定性。 2.基于Simulink平台，搭建3机9节点系统，通过时域仿真，得到三台机组的功角曲线和转速差曲线，以此判断系统的暂态稳定性。 注: 两种方法可以相互验证!

先撸代码解摇摆方程

直接上Matlab解微分方程的核心代码：

function dydt = swing_eq(t,y) % 参数初始化 H = [23.64; 6.4; 3.01]; D = [0.5; 0.3; 0.2]; Pm = [1.63; 0.85; 0.95]; Pe = calculate_Pe(y); % 这个函数要自己写，涉及网络方程求解 dydt = zeros(6,1); % 状态变量顺序：δ1, ω1, δ2, ω2, δ3, ω3 for i = 1:3 dydt(2*i-1) = y(2*i); dydt(2*i) = (Pm(i) - Pe(i) - D(i)*y(2*i)) / (2*H(i)); end end

这里用经典的二阶摇摆方程，注意每个发电机有两个状态变量：功角δ和转速偏差ω。网络方程求解部分需要处理节点导纳矩阵，这里偷个懒没展开（实际得用牛顿法算潮流）。

跑完仿真后的关键分析点：

figure; subplot(2,1,1); plot(t, delta_diff(:,1)-delta_diff(:,2), 'r', t, delta_diff(:,1)-delta_diff(:,3), 'b'); title('功角差变化'); legend('G1-G2','G1-G3'); subplot(2,1,2); plot(t, omega_diff); title('转速偏差');

看功角差曲线要是发散就直接判定不稳定，如果震荡衰减且最终稳定在某个值，系统还能撑住。转速差要是超过±0.2Hz就得警惕了。

Simulink建模更直观

在Simulink里拖拽元件搭建系统时有个坑要注意：同步电机的初始状态必须和潮流计算结果匹配。建议先用Powergui模块做初始化配置。

!3机9节点Simulink模型截图

（假装这里有模型截图）

故障设置推荐用三相短路模块，设置在母线5，0.1秒触发，0.2秒切除。运行后直接导出Workspace数据：

delta_simulink = simout.Data(:,1:3); omega_simulink = simout.Data(:,4:6); % 计算相对功角 delta_rel = delta_simulink - delta_simulink(:,1); plot(t_sim, delta_rel(:,2), t_sim, delta_rel(:,3));

比较编程和仿真结果时发现，当切除时间超过临界值0.25秒时，两种方法都显示G2的功角差突破120度，转速差像脱缰野马一样冲到了0.35Hz，系统崩得明明白白。

交叉验证小技巧

• 编程法参数调整灵活，但网络方程处理麻烦
• Simulink能直接观察电磁功率波动，适合验证动态过程
• 建议把两种方法的功角曲线叠在一起画，像这样：

hold on; plot(t_program, delta1_prog, 'k--'); plot(t_simulink, delta1_sim, 'r'); % 其他机组同理...

要是两条曲线基本重合，说明模型建得靠谱。上次有个师弟的仿真结果差10度，最后发现是变压器阻抗参数单位搞错了，被导师怼了一下午...

两种方法其实互补：编程法适合研究算法改进，Simulink方便做硬件在环测试。下次试试加AVR和PSS，估计又有新乐子。