gycall.m

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

function Gycall(a)

global DAE

if ~a.n, return, end

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

DAE.Gy = DAE.Gy ... 
  + sparse(a.bus,a.vbus,2*DAE.y(a.G).*V,DAE.m,DAE.m) ...
  - sparse(a.G,a.G,a.u,DAE.m,DAE.m) ...
  + sparse(a.bus,a.G,V.*V,DAE.m,DAE.m);