gycall.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 14 行

function Gycall(a)

global DAE

if ~a.n, return, end

V = a.u.*DAE.y(a.vbus);

DAE.Gy = DAE.Gy ...
         - sparse(a.vbus,a.q,1,DAE.m,DAE.m) ...
         + sparse(a.q,a.vbus,2*a.Be.*V,DAE.m,DAE.m) ...
         - sparse(a.q,a.q,1,DAE.m,DAE.m) ...
         - sparse(a.vref,a.vref,1,DAE.m,DAE.m);