add.m

来自「电力系统分析计算程序」· M 代码 · 共 32 行

function a = add(a,data,varargin)

global Settings

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

if isempty(data), return, end

% check data size
[nrow,ncol] = size(data);
if ncol < a.ncol
  data = [data, zeros(nrow,a.ncol-ncol)];
  if ncol < a.nu
    data(:,a.nu) = ones(nrow,1);
  end
elseif ncol > a.ncol
  data = data(:,[1:a.ncol]);
end

a.n = a.n + nrow;
a.con = [a.con; data];
a.u = [a.u; data(:,a.nu)];
[a.fr,a.vfr] = getbus(Bus,a.con(:,1));
[a.to,a.vto] = getbus(Bus,a.con(:,2));

Settings.nseries = Settings.nseries + nrow;