gycall.m

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

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