gcall.m

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

function gcall(p)

global DAE

if ~p.n, return, end

DAE.g = DAE.g ...
    - sparse(p.pm,1,p.u.*pmech(p),DAE.m,1) ...
    + sparse(p.wref,1,p.u.*p.con(:,3)-DAE.y(p.wref),DAE.m,1);