{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

decay-crippled

# decay-crippled - q-axis current i IT= Terminal current...

This preview shows page 1. Sign up to view the full content.

function xdot=decay(t,x) % Global parameters again global Ra Xdp Xd Xq Td0p Efd H KD Pm ... % Machine parameters ws EB XE RE XTq XTd D; % Network parameters, EB is infinite bus i % called as [t,x]=ode45('decay', [ 0 tend], [wd0 delta0 EQP0]'); % Name state variables for readability wd=x(1); delta=x(2); Eqp=x(3); E % Solve algebraic equations - this you'll need to figure out id= % d-axis current i iq=
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: % q-axis current i IT= % Terminal current referred to inf. bus ET= % Terminal voltage referred to inf. bus (Ra outside) Et= % Terminal voltage magnitude Pe= % Electrical power output P % Now write differential equations xdot(1)=ws/(2*H)*(Pm-Pe-KD*wd/ws); % Acceleration xdot(2)=wd; % Speed deviation xdot(3)=1/Td0p*(Efd-Eqp-(Xd-Xdp)*id); % Eqp derivative xdot=xdot'; x return;...
View Full Document

{[ snackBarMessage ]}