setup.m

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

function a = setup(a)

global Bus Settings

if isempty(a.con), return, end

a.n = length(a.con(:,1));
a.bus = uint32(Bus.int(round(a.con(:,1))));
a.gen = uint8(zeros(a.n,1)); 
[b,h,k] = unique(a.bus);
if length(k) > length(h)
  fm_disp(['Error: More than one PQ load ', ...
           'connected to the same bus.'],2)
  Settings.ok = 0;
  return
end
if length(a.con(1,:)) == 5
  a.con = [a.con, 1.2*ones(a.n,1), ...
           0.8*ones(a.n,1), zeros(a.n,1)];
elseif length(a.con(1,:)) == 7
  a.con = [a.con, zeros(a.n,1)];
end
for i = 1:a.n
  if a.con(i,6) == 0, a.con(i,6) = 1.2; end
  if a.con(i,7) == 0, a.con(i,7) = 0.8; end
end
a.P0 = a.con(:,4);
a.Q0 = a.con(:,5);
a.vmax = uint8(ones(a.n,1));
a.vmin = uint8(ones(a.n,1));
a.store = a.con;