gcall.m

基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理

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);