gcall.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 26 行

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);