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，估计又有新乐子。
