gcall.m

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

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

ids = DAE.y(a.ids);
V = DAE.y(a.vbus);
t = DAE.y(a.bus);
st = sin(t);
ct = cos(t);

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

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

DAE.g = DAE.g ...
        - sparse(a.bus,1,ps,DAE.m,1) ...
        - sparse(a.vbus,1,V.*idc./cos(t)+tan(t).*ps,DAE.m,1) ...
        + sparse(a.ids,1,vqs.*ids-vds.*iqs,DAE.m,1);