fcall.m

电力系统分析计算程序

function p = fcall(p)

global DAE Line

if ~p.n, return, end

t2 = p.con(:,3) == 2;
ta = p.con(:,4) == 2;
i2 = find(t2);

x0 = DAE.y(p.x0);
x1 = DAE.x(p.x1);
Pref = DAE.y(p.pref);

[Ps,Qs,Pr,Qr] = flows(Line,'pq',p.line);
[Ps,Qs,Pr,Qr] = flows(p,Ps,Qs,Pr,Qr,'tcsc');

Ki = p.con(:,13);

fx2 = -p.u.*t2.*Ki.*(Pref-Ps-ta.*Pr);
DAE.f(p.x2(i2)) = fx2(i2);

x1_max = p.u.*p.con(:,10);
x1_min = p.u.*p.con(:,11);

u = ~(x1 >= x1_max & DAE.f(p.x1) > 0) & ...
    ~(x1 <= x1_min & DAE.f(p.x1) < 0) & p.u;
DAE.f(p.x1) = u.*(x0-x1)./p.con(:,9);
DAE.x(p.x1) = min(x1_max,x1);
DAE.x(p.x1) = max(x1_min,DAE.x(p.x1));

% update B
p.B = btcsc(p);