setup.m

电力系统分析计算程序

function a = setup(a,varargin)

global Settings

switch nargin
 case 2
  Bus = varargin{1};
 otherwise
  global Bus
end

if isempty(a.con), return, end

a.n = length(a.con(:,1));
[a.bus1,a.v1] = getbus(Bus,a.con(:,1));
[a.bus2,a.v2] = getbus(Bus,a.con(:,2));

if length(a.con(1,:)) < a.ncol
  a.u = ones(a.n,1);
else
  a.u = a.con(:,a.ncol);
end

a.iters = zeros(a.n,1);
a.delay = zeros(a.n,1);
a.store = a.con;

% initialize delay function
a.iters = 0;
a.delay = 0;

% fix remote control bus number
a.vr = a.v2;
idx = find(a.con(:,16) == 3);
if ~isempty(idx)
  a.vr(idx) = getvint(Bus,a.con(idx,15));
  %a.con(idx,16) = 1;
end

% fix nominal tap ratio
idx = find(a.con(:,6) == 0);
if ~isempty(idx)
  a.dat(idx,6) = 1;
end

Settings.nseries = Settings.nseries + a.n;