MATLAB MAQUINAS ELECTRICAS

f=60;

Vta=1;

Xd=0.78;

Xq=0.63;

Xbarra=0.09;

%a. Corriente de armadura vs excitación

Eafmin=Vbus-XD=0.04;

Eafmax=Vbus+XD=1.96;

for Eaf=0.04:0.01:1.96

    Ia=1/(sqrt(Xd*Xd+Xq*Xq))*(Eaf-Vta);

    plot(Eaf,Ia)

end

xlabel('Excitación [p.u.]')

ylabel('Corriente de armadura [A]')

title('Corriente de armadura vs Excitación')

                        axis([0 2 0 1])

   b)

f=60;

Vta=1;

Xd=0.78;

Xq=0.63;

Xbarra=0.09;

%a. Corriente de armadura vs excitación

Eafmin=Vbus-XD=0.04;

Eafmax=Vbus+XD=1.96;

for Eaf=0.04:0.01:1.96

    Ia=1/(sqrt(Xd*Xd+Xq*Xq))*(Eaf-Vta);

    plot(Eaf,Ia)

end

xlabel('Excitación [p.u.]')

ylabel('Corriente de armadura [A]')

title('Corriente de armadura vs Excitación')

axis([0 2 0 1])

%b. Gráfica adicional corriente de campo

hold on

P=0.25;

Ia=P/(3*(Xd*Xd+Xq*Xq));

plot(Ia,If)

c)

f=60;

Vta=1;

Xd=0.78;

Xq=0.63;

Xbarra=0.09;

%a. Corriente de armadura vs excitación

Eafmin=Vbus-XD=0.04;

Eafmax=Vbus+XD=1.96;

for Eaf=0.04:0.01:1.96

    Ia=1/(sqrt(Xd*Xd+Xq*Xq))*(Eaf-Vta);

    plot(Eaf,Ia)

end

xlabel('Excitación [p.u.]')

ylabel('Corriente de armadura [A]')

title('Corriente de armadura vs Excitación')

axis([0 2 0 1])

%c. Gráfica adicional

hold on

P1=0.25;

P2=0.50;

P3=0.75;

Ia1=P1/(3*(Xd*Xd+Xq*Xq));

Ia2=P2/(3*(Xd*Xd+Xq*Xq));

Ia3=P3/(3*(Xd*Xd+Xq*Xq));

plot(Ia1,If)

plot(Ia2,If)

plot(Ia3,If)