remove.m

电力系统分析计算程序

function a = remove(a,k)

global DAE

if ~a.n, return, end

if isempty(k), return, end
if ~isempty(DAE.Gk)
  DAE.Gk(a.bus(k)) = 0;
end
a.con(k,:) = [];
a.bus(k) = [];
a.vbus(k) = [];
a.u(k) = [];
a.pq(k) = [];
a.qg(k) = [];
a.n = a.n - length(k);
a.qmax(k) = [];
a.qmin(k) = [];