gycall.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 24 行

function Gycall(p)

global Bus DAE

if ~p.n, return, end

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)
  h = p.bus(a);
  DAE.J12(h,h)= DAE.J12(h,h) + ...
      diag(2*p.con(a,4).*DAE.V(h)./p.con(a,7)./p.con(a,7));
  DAE.J22(h,h)= DAE.J22(h,h) + ...
      diag(2*p.con(a,5).*DAE.V(h)./p.con(a,7)./p.con(a,7));
end
if ~isempty(b)
  h = p.bus(b);
  DAE.J12(h,h)= DAE.J12(h,h) + ...
      diag(2*p.con(b,4).*DAE.V(h)./p.con(b,6)./p.con(b,6));
  DAE.J22(h,h)= DAE.J22(h,h) + ...
      diag(2*p.con(b,5).*DAE.V(h)./p.con(b,6)./p.con(b,6));
end