gycall.m

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

function Gycall(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);

a1 = ss.*a.B;
a2 = cc.*a.B;
a3 = V1.*V2;
a4 = a3.*a1;
a5 = a3.*a2;

% dP1/dy
DAE.J12 = DAE.J12 + sparse(a.bus1,a.bus1, V2.*a1,Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(a.bus1,a.bus2, V1.*a1,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(a.bus1,a.bus1, a5,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(a.bus1,a.bus2,-a5,Bus.n,Bus.n);

% dQ1/dy
DAE.J22 = DAE.J22 + sparse(a.bus1,a.bus1,(2*V1-V2.*cc).*a.B,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(a.bus1,a.bus2,-V1.*a2,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus1,a.bus1, a4,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus1,a.bus2,-a4,Bus.n,Bus.n);

% dP2/dy
DAE.J12 = DAE.J12 + sparse(a.bus2,a.bus1,-V2.*a1,Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(a.bus2,a.bus2,-V1.*a1,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(a.bus2,a.bus1,-a5,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(a.bus2,a.bus2, a5,Bus.n,Bus.n);

% dQ2/dy
DAE.J22 = DAE.J22 + sparse(a.bus2,a.bus2,(2*V2-V1.*cc).*a.B,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(a.bus2,a.bus1,-V2.*a2,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus2,a.bus1, a4,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus2,a.bus2,-a4,Bus.n,Bus.n);