fcall.m

用于电力系统的一个很好的分析软件

function fcall(p)

global DAE

if ~p.n, return, end

Vs = DAE.x(p.Vs);
Qgr = p.con(:,7);
Vs_max = p.con(:,8);
Vs_min = p.con(:,9);

DAE.f(p.Vs) = (Qgr.*DAE.y(p.cac)-DAE.y(p.q)).*p.dVsdQ.*p.u;

% non-windup limits
idx = find(Vs >= Vs_max & DAE.f(p.Vs) > 0);
if idx, DAE.f(p.Vs(idx)) = 0; end
DAE.x(p.Vs) = min(Vs,Vs_max);
idx = find(Vs <= Vs_min & DAE.f(p.Vs) < 0);
if idx, DAE.f(p.Vs(idx)) = 0; end
DAE.x(p.Vs) = max(Vs,Vs_min);