gcall.m

用于电力系统的一个很好的分析软件

function gcall(a)

global DAE

if ~a.n, return, end

Vf = a.u.*DAE.y(a.v1).*exp(i*DAE.y(a.bus1));
Vt = a.u.*DAE.y(a.v2).*exp(i*DAE.y(a.bus2));  
y = admittance(a);
m = DAE.x(a.m);

Ss = Vf.*conj((Vf-Vt.*m).*y);
Sr = Vt.*conj((Vt.*m-Vf).*y.*m);

DAE.g = DAE.g + ...
        sparse(a.bus1,1,real(Ss),DAE.m,1) + ...
        sparse(a.bus2,1,real(Sr),DAE.m,1) + ...
        sparse(a.v1,1,imag(Ss),DAE.m,1) + ...
        sparse(a.v2,1,imag(Sr),DAE.m,1);