Matlab - SHE2_PowerInSalientPolesGenerators

clear all
clc

Xd = 1.2;
Xq = 0.6 * Xd;
U = 30000
E = 32000 * sqrt(2)/2 * (1+i);
theta = 0: pi/200: pi;

Pe = U*E/Xd * sin(theta);
Pr = 0.5*U*U*(1/Xq - 1/Xd)*sin(2*theta);
P = Pe + Pr;

plot(theta, Pe);
hold on
plot(theta, Pr);
hold on
plt(theta, P);
legend({'Pe','Pr', 'P'},'Location','southwest');