fcall.m

电力系统分析计算程序

function fcall(a)

global DAE

if ~a.n, return, end

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

u = ~(ist >= ist_max & DAE.f(a.ist) > 0) & ...
    ~(ist <= ist_min & DAE.f(a.ist) < 0) & a.u;
DAE.f(a.ist) = u.*(Kr.*(DAE.y(a.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));