gycall.m

电力系统分析计算程序

function Gycall(p)

global DAE Settings

if ~p.n, return, end
if Settings.forcepq, return, end

a = find((DAE.y(p.vbus) < p.con(:,7) & p.con(:,8) & p.u) | p.shunt);
b = find(DAE.y(p.vbus) > p.con(:,6) & p.con(:,8) & p.u);
if ~isempty(a)
  h = p.bus(a);
  k = p.vbus(a);
  DAE.Gy = DAE.Gy + ...
      sparse(h,k,2*p.con(a,4).*DAE.y(k)./p.con(a,7)./p.con(a,7),DAE.m,DAE.m);
  DAE.Gy = DAE.Gy + ...
      sparse(k,k,2*p.con(a,5).*DAE.y(k)./p.con(a,7)./p.con(a,7),DAE.m,DAE.m);
end
if ~isempty(b)
  h = p.bus(b);
  k = p.vbus(b);
  DAE.Gy = DAE.Gy + ...
      sparse(h,k,2*p.con(b,4).*DAE.y(k)./p.con(b,6)./p.con(b,6),DAE.m,DAE.m);
  DAE.Gy = DAE.Gy + ...
      sparse(k,k,2*p.con(b,5).*DAE.y(k)./p.con(b,6)./p.con(b,6),DAE.m,DAE.m);
end