gycall.m

这是一个很适合研究和学习用的电力系统仿真软件

function Gycall(a)

global Bus DAE

if ~a.n, return, end

V = DAE.V(a.bus);

DAE.J22 = DAE.J22 + sparse(a.bus,a.bus,-2*a.Be.*V,Bus.n,Bus.n);