gcall.m

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

function gcall(a)

global DAE

if ~a.n, return, end

V = DAE.y(a.vbus);
Dw = domega(a);

DAE.g = DAE.g + sparse(a.bus,1,a.u.*a.con(:,2).*(V.^a.con(:,3)).* ...
                         (1+Dw).^a.con(:,4),DAE.m,1);
DAE.g = DAE.g + sparse(a.vbus,1,a.u.*a.con(:,5).*(V.^a.con(:,6)).* ...
                         (1+Dw).^a.con(:,7),DAE.m,1);