电力系统暂态上机暂态稳定.docxVIP

clear; clc; f=50;Tj=8.47; f=50; Tj=8.47; PT=1; P2M=0.48; P3M=1.38; %系统额定频率 %归算后的发电机惯性时间常数 %正常运行时发电机向无穷大系统传输的有功功率 %故障存在时发电机的最大功率 %故障切除后发电机的最大功率 %下面是利用改进欧拉法进行逐段计算 %需要注意故障切除前后电磁功率有跃变 h=0.01; %设置步长 0.01s Duration=20 %设置计算时段长度 20s CutTime=input( 输入故障切除时间 \n); %故障发生时的功角变化过程 Delta(1)=33.92; %初始功角 Omega(1)=1; % 初始转速 t(1)=0; for i=1:round(CutTime/h) d(i)=(Omega(i)-1)*360*f; o(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj;% 时段初变化率 Delta(i+1)=Delta(i)+d(i)*h; Omega(i+1)=Omega(i)+o(i)*h;% 时段末估计植 dd(i)=(Omega(i+1)-1)*360*f; oo(i)=(PT-P2M*sin(Delta(i+1)*pi/180))/Tj;% 时段末变化率 ddd(i)=(d(i)+dd(i))/2; ooo(i)=(o(i)+oo(i))