clear;
vel = 80*10^3 / 3600; %m/s
r = 0.3302; %13 polegadas em metros
w = vel / r;

Ptotal = 10000;
Pmot = Ptotal/4; %potencia de cada motor

Va = 380;
X = 15;
Inom = 20;
w_nom = 3600*pi/30;
ko = Va/w_nom;

Ea = ko*w;

%P = 3*Va*Ea*sen(delta)/X

delta = asin(Pmot*X / (3*Va*Ea));
I = (380 - (ko*w)*exp(j*delta))/(j*X);
fp =((ko*w)*sin(delta)/(X*I));
phi = acos(fp);

%%
figure (1)

Ea = ko*w*exp(-j*delta);
Vx = Va-Ea;
Ia = Vx/(j*X);

Vap = quiver(0,0,real(Va),imag(Va),'b','AutoScale','off')
hold on
Eap = quiver(0,0,real(Ea),imag(Ea),'g','AutoScale','off')
Vxp = quiver(real(Ea),imag(Ea),real(Vx),imag(Vx),'k', 'AutoScale','off')
Iap = quiver(0,0,real(Ia*X),imag(Ia*X),'r','AutoScale','off')
hold off