findpos.m

来自「电力系统分析计算程序」· M 代码 · 共 24 行

function idx = findpos(a)

global CPF

idx = [];
if ~a.n, return, end

if CPF.onlypqgen
  idx = find(a.u & a.gen);
  if isempty(idx)
    fm_disp('No PQ generator found. Expect meaningless results.')
  end
elseif CPF.onlynegload
  idx = find(a.u.*a.con(:,4) < 0);
  if isempty(idx)
    fm_disp('No negative load found. Expect meaningless results.')
  end
elseif CPF.negload
  idx = find(a.u);
else
  idx = find(a.u.*a.con(:,4) >= 0);
end