gycall.m

用于电力系统的一个很好的分析软件

function Gycall(a)

global DAE

if ~a.n, return, end

id = a.u.*DAE.x(a.Id);
iq = a.u.*DAE.x(a.Iq);
V = DAE.y(a.vbus);
theta = DAE.y(a.bus);
delta = DAE.x(a.delta);
cdt = cos(delta-theta);
sdt = sin(delta-theta);

DAE.Gy = DAE.Gy + sparse(a.bus,a.bus,V.*cdt.*id-V.*sdt.*iq,DAE.m,DAE.m);
DAE.Gy = DAE.Gy - sparse(a.bus,a.vbus,sdt.*id+cdt.*iq,DAE.m,DAE.m);
DAE.Gy = DAE.Gy - sparse(a.vbus,a.bus,V.*sdt.*id+V.*cdt.*iq,DAE.m,DAE.m);
DAE.Gy = DAE.Gy + sparse(a.vbus,a.vbus,sdt.*iq-cdt.*id,DAE.m,DAE.m);