gcall.m

来自「电力系统分析计算程序」· M 代码 · 共 14 行

function gcall(a)

global DAE

if ~a.n, return, end

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

DAE.g = DAE.g ...
        - sparse(a.vbus,1,DAE.y(a.q),DAE.m,1) ...
        + sparse(a.q,1,B.*V.*V-DAE.y(a.q),DAE.m,1) ...
        + sparse(a.vref,1,a.con(:,8)-DAE.y(a.vref),DAE.m,1);