fcall.m

电力系统分析计算程序

function fcall(p)

global DAE

if ~p.n, return, end

q1 = DAE.x(p.q1);
q1_max = p.con(:,8);
q1_min = p.con(:,9);

DAE.f(p.q1) = p.u.*p.con(:,6).*(p.con(:,5)-DAE.y(p.vbus));

% non-windoup limits
idx = find(q1 >= q1_max & DAE.f(p.q1) > 0);
if idx, DAE.f(p.q1(idx)) = 0; end
DAE.x(p.q1) = min(q1,q1_max);
idx = find(q1 <= q1_min & DAE.f(p.q1) < 0);
if idx, DAE.f(p.q1(idx)) = 0; end
DAE.x(p.q1) = max(q1,q1_min);