fcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function fcall(a)

global DAE Pod

if ~a.n, return, end

V = DAE.V(a.bus);
ist = DAE.x(a.ist);
Kr = a.con(:,5);
Tr = a.con(:,6);
ist_max = a.con(:,7);
ist_min = a.con(:,8);

% Updating reference voltages
Vref = a.Vref;
if a.pod
  Vref(a.pod) = Vref(a.pod) + DAE.x(Pod.Vs(Pod.statcom));
end

u = ~(ist >= ist_max & DAE.f(a.ist) > 0) & ...
    ~(ist <= ist_min & DAE.f(a.ist) < 0);
DAE.f(a.ist) = u.*(Kr.*(Vref-V)-ist)./Tr;  
DAE.x(a.ist) = min(ist_max,DAE.x(a.ist));
DAE.x(a.ist) = max(ist_min,DAE.x(a.ist));