Dibujo de la curva de cálculo de cortocircuito sin carga y diagrama de simulación en la salida del generador

%发电机出口空载短路计算曲线、仿真图绘制
% %***********************************************
N=48;
t1=(0:0.02/N:1.00);
m=size(t1);
fai=0*pi/180;
%空载短路全电流表达式
Ia=(-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai)+...
    4*exp(-6.3*t1).*cos(-fai*pi/180)+0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+fai));
%基频分量
Ia1=-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai);
%倍频分量
Ia2=0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+fai);
%非周期分量
Iap=4*exp(-6.3*t1).*cos(-fai*pi/180);

subplot(4,1,1);
plot(t1,Ia);
grid on;
axis([0 1 -10 10]);
ylabel('Ia(p.u.)');

subplot(4,1,2);
plot(t1,Ia1);
grid on;
axis([0 1 -10 10]);
ylabel('Ia1(p.u.)');

subplot(4,1,3);
plot(t1,Ia2);
grid on;
axis([0 1 -1 1]);
ylabel('Ia2(p.u.)');

subplot(4,1,4);
plot(t1,Iap);
grid on;
axis([0 1 -10 10]);
ylabel('Iap(p.u.)');
xlabel('t/s');