gcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function gcall(p)

global DAE

if ~p.n, return, end

DAE.gp(p.bus) = p.con(:,4) + DAE.gp(p.bus);
DAE.gq(p.bus) = p.con(:,5) + DAE.gq(p.bus);
a = find(DAE.V(p.bus) < p.con(:,7) & p.con(:,8));
b = find(DAE.V(p.bus) > p.con(:,6) & p.con(:,8));
if ~isempty(a)
  k = p.bus(a);
  DAE.gp(k) = p.con(a,4).*DAE.V(k).*DAE.V(k)./p.con(a,7)./p.con(a,7) ...
      + DAE.gp(k) - p.con(a,4);
  DAE.gq(k) = p.con(a,5).*DAE.V(k).*DAE.V(k)./p.con(a,7)./p.con(a,7) ...
      + DAE.gq(k) - p.con(a,5);
end
if ~isempty(b)
  k = p.bus(b);
  DAE.gp(k) = p.con(b,4).*DAE.V(k).*DAE.V(k)./p.con(b,6)./p.con(b,6) + ...
      DAE.gp(k) - p.con(b,4);
  DAE.gq(k) = p.con(b,5).*DAE.V(k).*DAE.V(k)./p.con(b,6)./p.con(b,6) + ...
      DAE.gq(k) - p.con(b,5);
end