电力系统暂态稳定性Matlab编程/ Simulink仿真 单机无穷大系统发生各类（三相短路，单相接地，两相接地，两相相间短路）等短路故障，各类（单相断线，两相断线，三相断线）等断线故障，暂态稳定仿真分析 1.Matlab编程进行数值分析（才用欧拉法Euler，改进欧拉法improve Euler，4阶龙格库塔法Runger-Kutta），计算故障后发电机功角-时间曲线，电机转速-时间曲线，临界切除角和临界切除时间 2. Simulink搭建电力系统暂态仿真模型 通过仿真，观察串联电抗器，并联补偿器，自动重合闸，以及故障切除快慢对暂态稳定性的影响
单机无穷大系统的暂态稳定性分析是电力系统研究中的经典问题。今天咱们直接上代码和模型，用Matlab和Simulink搞点实在的。先来个三相短路场景，假设发电机用经典二阶模型，微分方程长这样：
function dydt = generator_fault(t, y, Pmech, E, V, Xd, H) delta = y(1); omega = y(2); Pe = E*V*sin(delta)/Xd; dydt = [omega; (Pmech - Pe)/(2*H)]; end这里H是惯性常数，Xd是暂态电抗。用欧拉法实现的话，时间步长选0.01秒足够对付大多数情况：
h = 0.01; t = 0:h:5; y = zeros(2, length(t)); y(:,1) = [0.1; 0]; % 初始功角10度，转速0 for k=1:length(t)-1 dy = generator_fault(t(k), y(:,k), 0.8, 1.05, 1.0, 0.3, 5); y(:,k+1) = y(:,k) + dy*h; % 欧拉法核心 end但遇到故障剧烈变化时，欧拉法可能翻车。改用改进欧拉法，在预测-校正两步走：
ypred = y(:,k) + dy*h; dy_corr = generator_fault(t(k)+h, ypred, 0.8, 1.05, 1.0, 0.3, 5); y(:,k+1) = y(:,k) + (dy + dy_corr)*h/2;想要更高精度就上四阶龙格库塔，虽然代码多几行，但误差能压到0.001 rad以内：
k1 = generator_fault(t(k), y(:,k), ...); k2 = generator_fault(t(k)+h/2, y(:,k)+k1*h/2, ...); k3 = generator_fault(t(k)+h/2, y(:,k)+k2*h/2, ...); k4 = generator_fault(t(k)+h, y(:,k)+k3*h, ...); y(:,k+1) = y(:,k) + (k1+2*k2+2*k3+k4)*h/6;跑完仿真后重点看功角摇摆曲线——如果超过120度大概率失稳。临界切除时间可以通过二分法暴力搜索，比如从0.2秒开始逐步试探，直到功角曲线不再发散。
转战Simulink搭建模型时，故障模块的配置是关键。比如设置三相短路从0.5秒开始，0.65秒切除，直接影响到转子摇摆幅度。建议用Powergui模块启用Phasor仿真模式，比默认的Continuous模式快三倍不止。
并联补偿器的作用在仿真中一目了然——加上SVC后母线电压波动减少约40%。这里有个骚操作：在故障期间动态调整电抗器值，能硬生生把临界切除时间从0.6秒拉到0.8秒。
电力系统暂态稳定性Matlab编程/ Simulink仿真 单机无穷大系统发生各类（三相短路，单相接地，两相接地，两相相间短路）等短路故障，各类（单相断线，两相断线，三相断线）等断线故障，暂态稳定仿真分析 1.Matlab编程进行数值分析（才用欧拉法Euler，改进欧拉法improve Euler，4阶龙格库塔法Runger-Kutta），计算故障后发电机功角-时间曲线，电机转速-时间曲线，临界切除角和临界切除时间 2. Simulink搭建电力系统暂态仿真模型 通过仿真，观察串联电抗器，并联补偿器，自动重合闸，以及故障切除快慢对暂态稳定性的影响
自动重合闸的仿真最刺激的是时序设置。实测发现重合闸动作比故障切除晚0.3秒时，系统会出现二次冲击。建议搭配FFT分析工具观察谐波成分，防止保护误动作。
最后说个干货：暂态稳定对计算步长极其敏感。对比发现，当故障切除时间在0.5秒附近时，步长从0.01秒改为0.02秒会导致临界时间判断误差达8%，这时候就得掏出龙格库塔法来救场了。
