matLab Gerador sincrono

Simular um gerador Sincrono e obter os diagramas Fasoriais, classificar quanto a excitação
Entradas: Vt, Xs, S e fp.
Simular no matlab. Mostrar gráfico.

vt = input(' Digite o valor de Vt: '); xs = input(' Digite o valor de Xs: '); s = input(' Digite o valor de S: '); fp = input(' Digite o valor do f.p. (negativo para capacitivo, positivo caso contrário): '); p = s * abs(fp) % potência ativa; ang = acos(fp); % Angulo do fp em radianos. angd = acosd(fp) % Angulo do fp em graus. I = s / 3*vt; jxsi = xs * I; temp = vt+ xs * I * cosd(90 - angd); ef = sqrt((temp)^2 + (xs*I*sind(90-angd)^2))

if ef == vt texto = 'Excitação Normal'; else if ef > vt texto = 'Sobrexcitado'; else texto = 'Subexcitado'; end; end; ang = 3 * pi/2 - ang; % inverte o angulo para o gráfico. % transformando grandezas escalares em vetores. vtg =0:vt/100:vt; xsg = 0:0.01:xs; sg = 0:s/50:s; pg = -p/10:p/50:p; ig = 0:0.01:I; temp = size(ig); tam = temp(1,2); % quantidade de elementos do vetor ig. angg = linspace(ang,ang,tam); efg = 0:0.05:ef; temp = size(efg); tam = temp(1,2) dirEfg = 0+p/tam:p/tam:p;

jxsg = vt:0.08:jxsi+vt; temp = size(jxsg); tam = temp(1,2); djxs = 0+p/tam:p/tam:p; hold on plot(p,-pg,'black'); % pot ativa. plot(0,-vtg,'blue'); % tensão V polar(angg,ig,'k'); % corrente I. plot(dirEfg,-efg,'red'); % Ef plot(djxs,-jxsg,'m'); % jxsi. plot(sg,0,'red'); % Linha horizontal à direita plot(-sg,0,'red'); % Linha horizontal à esquerda hold off

text(0,-vt/2,'V');
%text(-2*vt,-2*i,'I');
text(p/2,-ef/2,'Ef'); text(p/2-10,-vt-jxsi/2,'jXsI'); text(-s,0,texto); text(p,-p/2,'P const'); legend('azul (V)','preto ( I )','vermelho (Ef)','rosa (jXsi)',0);

Programa exemplo: vt = 10; xs = 12.5; s = 100; fp = 0.60;