remove.m

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

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) = [];
    a.bcv(jdx,:) = [];
  end
  jdx = find(a.ty2 == idx(i));
  if jdx
    a.ty2(jdx) = []; 
    a.alpha(jdx,:) = [];
    a.vm(jdx,:) = [];
  end  
  jdx = find(a.pod == idx(i));
  if jdx, a.pod(jdx) = []; end
  jdx = find(a.cluster == idx(i));
  if jdx, a.cluster(jdx) = []; end  
end

a.con(idx,:) = [];
a.bus(idx) = [];
a.n = a.n - length(idx);
a.Be(idx,:) = [];
a.Vref(idx,:) = [];