remove.m

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

function a = remove(a,idx)

if ~a.n, return, end
if isempty(idx), return, end

for i = 1:length(idx)
  jdx = find(a.ty1 == idx(i));
  if jdx 
    a.ty1(jdx) = [];
  end
  jdx = find(a.ty2 == idx(i));
  if jdx 
    a.ty2(jdx) = [];
    a.x2(jdx) = [];
  end  
  jdx = find(a.pod == idx(i));
  if jdx, a.pod(jdx) = []; end
end

a.con(idx,:) = [];
a.bus1(idx) = [];
a.bus2(idx) = [];
a.line(idx) = [];
a.n = a.n - length(idx);
a.x1(idx,:) = [];
a.B(idx,:) = [];
a.Cp(idx,:) = [];
a.x0(idx,:) = [];
a.Pref(idx,:) = [];
a.y(idx,:) = [];