gycall.m

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

function Gycall(a)

global DAE

if ~a.n, return, end

m = DAE.x(a.m);
V = DAE.y(a.vbus);
a1 = a.con(:,11);
a2 = a.con(:,12);

DAE.Gy = DAE.Gy + sparse(a.bus,a.vbus, ...
          a.u.*a.con(:,9).*a1.*(V.^(a1-1))./(m.^a1),DAE.m,DAE.m);

DAE.Gy = DAE.Gy + sparse(a.vbus,a.vbus, ...
          a.u.*a.con(:,10).*a2.*(V.^(a2-1))./(m.^a2),DAE.m,DAE.m);