glcall.m

电力系统分析计算程序

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