gcall.m

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

function gcall(a)

if ~a.n, return, end

global DAE

type = a.con(:,4);
a1 = find(type == 1);
a2 = find(type == 2);
a3 = find(type == 3);
a4 = find(type == 4);
a5 = find(type == 5);

Vs = a.u.*DAE.x(a.Vs);

DAE.g = DAE.g + sparse(a.svc,1,Vs(a1),DAE.m,1);
DAE.g = DAE.g + sparse(a.tcsc,1,Vs(a2),DAE.m,1);
DAE.g = DAE.g + sparse(a.statcom,1,Vs(a3).*a.kr,DAE.m,1);
DAE.g = DAE.g + sparse(a.sssc,1,Vs(a4),DAE.m,1);
DAE.g = DAE.g + sparse(a.upfc,1,Vs([a5;a5;a5]).*a.z,DAE.m,1);