gycall.m

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

function Gycall(p)

global DAE Settings

if ~p.n, return, end

DAE.J11(p.bus,:) = 0;
DAE.J12(p.bus,:) = 0;
DAE.J11(:,p.bus) = 0;
DAE.J21(:,p.bus) = 0;
DAE.J11(p.bus,p.bus) = DAE.J11(p.bus,p.bus) + speye(p.n);

if Settings.pv2pq
  idx = p.bus(find(~p.dq));
else
  idx = p.bus;
end

if isempty(idx),return, end

n = length(idx);
DAE.J21(idx,:) = 0;
DAE.J22(idx,:) = 0;
DAE.J12(:,idx) = 0;
DAE.J22(:,idx) = 0;
DAE.J22(idx,idx) = DAE.J22(idx,idx) + speye(n);