gcall.m

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

function a = gcall(a)

global Bus DAE

if ~a.n, return, end

V1 = DAE.V(a.bus1);
V2 = DAE.V(a.bus2);
a1 = DAE.a(a.bus1);
a2 = DAE.a(a.bus2);

vp = DAE.x(a.vp);
vq = DAE.x(a.vq);
iq = DAE.x(a.iq);

a.gamma = fgamma(a);
ss = sin(a1-a2+a.gamma);
cc = cos(a1-a2+a.gamma);

c1 = a.y.*sqrt(vp.*vp+vq.*vq);

P1 =  c1.*V2.*ss;
Q1 =  c1.*V1.*cos(a.gamma)-iq.*V1;
Q2 = -c1.*V2.*cc;
  
DAE.gp = DAE.gp + sparse(a.bus1,1,P1,Bus.n,1);
DAE.gp = DAE.gp - sparse(a.bus2,1,P1,Bus.n,1);
DAE.gq = DAE.gq + sparse(a.bus1,1,Q1,Bus.n,1);
DAE.gq = DAE.gq + sparse(a.bus2,1,Q2,Bus.n,1);