gcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function gcall(a)

global Bus DAE

if ~a.n, return, end

V1 = DAE.V(a.bus1);
V2 = DAE.V(a.bus2);
t1 = DAE.a(a.bus1);
t2 = DAE.a(a.bus2);
ss = sin(t1-t2);
cc = cos(t1-t2);

% update B
B = btcsc(a);

P1 = V1.*V2.*ss.*B;

DAE.gp = DAE.gp ...
         + sparse(a.bus1,1,P1,Bus.n,1) ...
         - sparse(a.bus2,1,P1,Bus.n,1);

DAE.gq = DAE.gq ...
         + sparse(a.bus1,1,V1.*(V1-V2.*cc).*B,Bus.n,1) ...
         + sparse(a.bus2,1,V2.*(V2-V1.*cc).*B,Bus.n,1);