B = 0.01;
V = 380;
X = 15;
Inom= 20;
w_nom = 3600*pi/30;

%% A vazio

Eo = V;
Io = 0;
ko = Eo/w_nom;
tetao = 0;
phio = acos(0);


%% velocidade nominal;
%|V| = |E|

teta_nom = asin((B*X*w_nom^2)/(3*V^2));
I_wn = (380 - (ko*w_nom)*exp(j*teta_nom))/(j*X);
phi_nom = acos (V*sin(teta_nom)/(X*I_wn));

%%
%encontrar w_max -> Ia <= 1.5Inom
% T = Bw = 3Va*Ea*sin(teta)/(X*w)
w = [ 370: 0.01 : 900];
 for i=1:1:length(w)
    teta(i) = asin((B*X*w(i)^2)/(3*(ko*w(i))^2));
    I(i) = (380 - (ko*w(i))*exp(j*teta(i)))/(j*X);
 end
 indice = find (abs(I)>1.5*Inom, 1); %primeiro valor de velocidade maior que a maxima
 wmax = w(indice-1);
 teta_wmax = teta(indice-1);
 I_wmax = I(indice-1);
 phi_wmax = acos(ko*wmax*sin(teta_wmax)/(X*I_wmax));



%%
%try to plot

ph1o = V*exp(i*0);
ph2o = Eo*exp(i*tetao);
ph3o = Io*exp(i*phio);
ph4o = i*X*Io;

p1 = [0; ph1o];
vecs = cumsum([ph2o; ph4o]);
vecs = [p1+vecs ; vecs];
p3 = [0; ph3o];

figure;
plot(real(p1), imag(p1), '->');
%xlim([-30 30]);
%ylim([-30 30]);
xlabel('real part');
ylabel('imaginary part');
grid on;
hold on;
plot(real(vecs), imag(vecs), '-+');
plot(real(p3), imag(p3), '-*');
%plot(real(p4), imag(p4), '--');

%%
ph1nom = V*exp(i*0);                 %Va
ph2nom = ko*w_nom*exp(i*teta_nom);      %Ea
ph3nom = I_wn*exp(i*phi_nom);       %Ia
ph4nom = i*X*I_wn;                 %jX*Ia

p1 = [0; ph1nom];
vecs = cumsum([ph2nom; ph4nom]);
vecs = [0 ; vecs];
p3 = [0; ph3nom];
%p4 = [ph2nom; ph4nom];

figure;
plot(real(p1), imag(p1), '->');
%xlim([-30 30]);
%ylim([-30 30]);
xlabel('real part');
ylabel('imaginary part');
grid on;
hold on;
plot(real(vecs), imag(vecs), '-+');
plot(real(p3), imag(p3), '-*');
%plot(real(p4), imag(p4), '--');
%
%