addgen.m

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

function a = addgen(a,b)

if ~b.n, return, end

% set generated powers as negative loads
b.gen = ones(b.n,1);
b.P0 = -b.P0;
b.Q0 = -b.Q0;
b.con(:,4) = b.P0;
b.con(:,5) = b.Q0;

% add PQ generators to PQ loads
a.con = [a.con; b.con];
a.n = a.n + b.n;
a.bus = [a.bus; b.bus];
a.gen = [a.gen; b.gen];
a.P0 = [a.P0; b.P0];
a.Q0 = [a.Q0; b.Q0];
a.vmax = [a.vmax; b.vmax];
a.vmin = [a.vmin; b.vmin];
a.store = [a.store; b.store];

[u,h,k] = unique(a.bus);
if length(k) > length(h)
  fm_disp(['Error: it is not allowed to connect a PQ load and a PQ ', ...
           'generator to the same bus.'],2)
  Settings.ok = 0;
  return
end