gycall.m

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

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.J12(:,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);