gcall.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 11 行

function gcall(a)

global Bus DAE

if ~a.n, return, end

V = DAE.V(a.bus);
B = bsvc(a);

DAE.gq = DAE.gq + sparse(a.bus,1,-B.*V.*V,Bus.n,1);