setup.m

电力系统分析计算程序

function a = setup(a)

global Bus

if isempty(a.con), return, end

a.n = length(a.con(:,1));
a.bus = getint(Bus,a.con(:,1));
if length(a.con(1,:)) < a.ncol
  a.con(:,a.ncol) = ones(a.n,1);
end 
a.Pr = zeros(a.n,1);
a.u = a.con(:,a.ncol);
a.store = a.con;