dfpek.m

市面上电力系统的范例 直得收藏学习的ㄧ本书

% State variable representation of the multimachine system
% during fault. (for use with trstab)
% Copyright (c) 1998 by H. Saadat

function xdot = dfpek(t,x)
global Pm f H E  Y th ngg
Pe=zeros(1, ngg);
for ii = 1:ngg
for jj = 1:ngg
Pe(ii) = Pe(ii) + E(ii)*E(jj)*Y(ii, jj)*cos(th(ii, jj)-x(ii)+x(jj));
end, end
for k=1:ngg
xdot(k)=x(k+ngg);
xdot(k+ngg)=(pi*f)/H(k)*(Pm(k)-Pe(k));
end
xdot=xdot';        % use with MATLAB 5 (remove for MATLAB 4)