gycall.m

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

function Gycall(a)

global Bus DAE

if ~a.n, return, end

bus1 = a.bus1;
bus2 = a.bus2;

V1 = DAE.V(bus1);
V2 = DAE.V(bus2);
a1 = DAE.a(bus1);
a2 = DAE.a(bus2);

vp = DAE.x(a.vp);
vq = DAE.x(a.vq);
iq = DAE.x(a.iq);

ss = sin(a1-a2+a.gamma);
cc = cos(a1-a2+a.gamma);

c1 = a.y.*sqrt(vp.*vp+vq.*vq);
U = max(V1.^2+V2.^2-2.*V1.*V2.*cos(a1-a2),1e-6);

L0 = V1.*V2.*cos(a1-a2);
L1 = (U-V2.^2+L0)./U;
L2 = (U-V1.^2+L0)./U;
L5 = sin(a1-a2)./U;
L3 = -V2.*L5;
L4 =  V1.*L5;  

P1a1 =  c1.*L1.*V2.*cc;
P1a2 =  c1.*(L2-1).*V2.*cc;
P2a1 = -c1.*L1.*V2.*cc;  
P2a2 = -c1.*(L2-1).*V2.*cc;  
P1v1 =  c1.*L3.*V2.*cc;
P1v2 =  c1.*(ss+V2.*L4.*cc);
P2v1 = -c1.*L3.*V2.*cc;
P2v2 = -c1.*(ss+V2.*L4.*cc);  
Q1a1 = -c1.*(L1-1).*V1.*sin(a.gamma);
Q1a2 = -c1.*L2.*V1.*sin(a.gamma);
Q2a1 =  c1.*L1.*V2.*ss;
Q2a2 =  c1.*(L2-1).*V2.*ss;
Q1v1 =  c1.*(cos(a.gamma)-V1.*L3.*sin(a.gamma))-iq; 
Q1v2 = -c1.*L4.*V1.*sin(a.gamma);
Q2v1 =  c1.*L3.*V2.*ss;
Q2v2 = -c1.*(cc-V2.*L4.*ss);

DAE.J11 = DAE.J11 + sparse(bus1,bus1,P1a1,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(bus1,bus2,P1a2,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(bus2,bus1,P2a1,Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(bus2,bus2,P2a2,Bus.n,Bus.n);  
DAE.J12 = DAE.J12 + sparse(bus1,bus1,P1v1,Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(bus1,bus2,P1v2,Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(bus2,bus1,P2v1,Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(bus2,bus2,P2v2,Bus.n,Bus.n);  
DAE.J21 = DAE.J21 + sparse(bus1,bus1,Q1a1,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(bus1,bus2,Q1a2,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(bus2,bus1,Q2a1,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(bus2,bus2,Q2a2,Bus.n,Bus.n);  
DAE.J22 = DAE.J22 + sparse(bus1,bus1,Q1v1,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(bus1,bus2,Q1v2,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(bus2,bus1,Q2v1,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(bus2,bus2,Q2v2,Bus.n,Bus.n);