gcall.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 21 行

function 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);
ss = sin(a1-a2);
cc = cos(a1-a2);

c1 =  DAE.x(a.vcs).*a.y./ssscden(a);  
P1 = c1.*V1.*V2.*ss;  
Q1 = c1.*(V1.^2-V1.*V2.*cc);
Q2 = c1.*(V2.^2-V1.*V2.*cc);

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,Q1,Bus.n,1) + sparse(a.bus2,1,Q2,Bus.n,1);