fcall.m

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

function fcall(p)

global DAE Settings

if ~p.n, return, end

m = DAE.x(p.m);
h = p.u.*p.con(:,4);
k = p.u.*p.con(:,5);
mmax = p.con(:,6);
mmin = p.con(:,7);

DAE.f(p.m) = -h.*m + k.*(DAE.y(p.vbus)./m - p.con(:,8));

% non-windup limits
idx = find(m >= mmax & (DAE.f(p.m) > 0 | ~Settings.init));
if idx, DAE.f(p.m(idx)) = 0; end
DAE.x(p.m) = min(DAE.x(p.m),mmax);
idx = find(m <= mmin & (DAE.f(p.m) < 0 | ~Settings.init));
if idx, DAE.f(p.m(idx)) = 0; end
DAE.x(p.m) = max(DAE.x(p.m),mmin);