clear;clc;f=50;%系统额定频率Tj=8.47;%归算后的发电机惯性时间常数PT=1;%正常运行时发电机向无穷大系统传输的有功功率P2M=0.48;%故障存在时发电机的最大功率P3M=1.38;%故障切除后发电机的最大功率%下面是利用改进欧拉法进行逐段计算%需要注意故障切除前后电磁功率有跃变h=0.05;%设置步长0.05sDuration=2;%设置计算时段长度2sDelta(1)=33.92;%初始功角Omega(1)=1;%初始转速t(1)=0;Delta_h=