gycall.m

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

function Gycall(a)

global DAE

if ~a.n, return, end

omega_m = DAE.x(a.omega_m);
iqs = DAE.x(a.iqs);
idc = DAE.x(a.idc);

V = DAE.y(a.vbus);
t = DAE.y(a.bus);
ids = DAE.y(a.ids);

rs = a.con(:,6);
xd = a.con(:,7);
xq = a.con(:,8);
psip = a.con(:,9);

c1 = cos(t);
t1 = tan(t);

vds = -rs.*ids+omega_m.*xq.*iqs;
vqs = -rs.*iqs-omega_m.*(xd.*ids-psip);
ps = vds.*ids+vqs.*iqs;

k1 = -2*rs.*ids+omega_m.*(xq-xd).*iqs;

DAE.Gy = DAE.Gy ...
         - sparse(a.bus,a.ids,k1,DAE.m,DAE.m) ...
         - sparse(a.vbus,a.vbus,idc./c1,DAE.m,DAE.m) ...
         - sparse(a.vbus,a.bus,t1.*V.*idc./c1+(1+t1.^2).*ps,DAE.m,DAE.m) ...
         - sparse(a.vbus,a.ids,t1.*k1,DAE.m,DAE.m) ...
         + sparse(a.ids,a.ids,omega_m.*(psip-2*xd.*ids),DAE.m,DAE.m);