setup.m

用于电力系统的一个很好的分析软件

function a = setup(a)

global Syn

if isempty(a.con), return, end

a.n = length(a.con(:,1));

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

% machine indexes and parameters
a.syn = a.con(:,1);
a.con(:,18) = 1./getvar(Syn,a.syn,'M');

% the shaft is inactive if the machine is off-line
a.u = a.u.*Syn.u(a.syn);

% setup inverse of inertias
a.con(:,2) = 1./a.con(:,2);
a.con(:,3) = 1./a.con(:,3);
a.con(:,4) = 1./a.con(:,4);
a.con(:,5) = 1./a.con(:,5);

a.store = a.con;