setup.m

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

function a = setup(a)

global Bus DAE Settings

if isempty(a.con)
  fm_disp('Error: No slack bus found.',2)
  Settings.ok = 0;
  return
end

a.n = length(a.con(:,1));
a.bus = Bus.int(round(a.con(:,1)));
b = unique(a.bus);
if a.n > length(b)
  fm_disp(['Error: More than one slack generator ', ...
           'connected to the same bus.'],2)
  Settings.ok = 0;
  return
end
if length(a.con(1,:)) < 8
  a.con = [a.con(:,1),100*ones(a.n,1), ...
            ones(a.n,1),a.con(:,[2:end])];
end
if length(a.con(1,:)) == 5
  a.con = [a.con, 999*ones(a.n,1), -999*ones(a.n,1), ...
            1.1*ones(a.n,1), 0.9*ones(a.n,1)];
end
if length(a.con(1,:)) == 10
  a.con = [a.con, ones(a.n,1)];
end
DAE.V(a.bus) = a.con(:,4);
if ~sum(DAE.a) & a.n == 1
  DAE.a(:) = a.con(1,5);
else
  DAE.a(a.bus) = a.con(:,5);
end
a.pg = a.con(:,10);
a.dq = uint8(zeros(a.n,1));
a.qg = zeros(a.n,1);
a.qmax = uint8(ones(a.n,1));
a.qmin = uint8(ones(a.n,1));
a.store = a.con;