glcall.m

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

function Glcall(a)

if ~a.n, return, end

global DAE

V = DAE.y(a.vbus);

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