fcall.m

电力系统分析计算程序

function fcall(p)

global DAE

if ~p.n, return, end

global DAE

DAE.f(p.v) = p.u.*(DAE.y(p.If) - p.con(:,6))./p.con(:,2);

% anti-windup limit
idx = find(DAE.x(p.v) <= 0 & DAE.f(p.v) < 0);
DAE.f(p.v(idx)) = 0;
idx = find(DAE.x(p.v) >= p.con(:,7) & DAE.f(p.v) > 0);
DAE.f(p.v(idx)) = 0;

DAE.x(p.v) = min(DAE.x(p.v),p.con(:,7));
DAE.x(p.v) = max(p.u.*DAE.x(p.v),0);