gcall.m

来自「电力系统分析计算程序」· M 代码 · 共 20 行

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 = a.con(:,15).*exp(i*DAE.x(a.alpha));

Ss = Vf.*conj((Vf-Vt.*m).*y);
Sr = Vt.*conj((Vt.*m-Vf).*y.*conj(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);