fcall.m

电力系统分析计算程序

function fcall(p)

global DAE

if ~p.n, return, end

V1 = DAE.y(p.v1);
vp = DAE.x(p.vp);
vq = DAE.x(p.vq);
iq = DAE.x(p.iq);

Kr = p.con(:,7);
Tr = p.con(:,8);

vp_max = p.con(:,9);
vp_min = p.con(:,10);
vq_max = p.con(:,11);
vq_min = p.con(:,12);
iq_max = p.con(:,13);
iq_min = p.con(:,14);

u = ~(vp >= vp_max & DAE.f(p.vp) > 0) & ...
    ~(vp <= vp_min & DAE.f(p.vp) < 0) & p.u;
DAE.f(p.vp) = u.*(DAE.y(p.vp0)-vp)./Tr; 

u = ~(vq >= vq_max & DAE.f(p.vq) > 0) & ...
    ~(vq <= vq_min & DAE.f(p.vq) < 0) & p.u;
DAE.f(p.vq) = u.*(DAE.y(p.vq0)-vq)./Tr;

u = ~(iq >= iq_max & DAE.f(p.iq) > 0) & ...
    ~(iq <= iq_min & DAE.f(p.iq) < 0) & p.u;
DAE.f(p.iq) = u.*(Kr.*(DAE.y(p.vref)-V1)-iq)./Tr;

DAE.x(p.vp) = min(vp_max,DAE.x(p.vp));
DAE.x(p.vp) = max(vp_min,DAE.x(p.vp));
DAE.x(p.vq) = min(vq_max,DAE.x(p.vq));
DAE.x(p.vq) = max(vq_min,DAE.x(p.vq));      
DAE.x(p.iq) = min(iq_max,DAE.x(p.iq));
DAE.x(p.iq) = max(iq_min,DAE.x(p.iq));