autorun.m

这是一个很适合研究和学习用的电力系统仿真软件

function check = autorun(msg,type)
% AUTORUN properly launch PSAT routine checking for data
%         files and previous power flow solutions
%
% CHECK = AUTORUN(MSG)
%         MSG   message to be displayed
%         TYPE  0 for static analysis, 1 for dynamic analysis
%         CHECK 1 if everything goes fine, 0 otherwise
%
%Author:    Federico Milano
%Date:      29-Oct-2003
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2006 Federico Milano

global Settings File Snapshot Bus
global DAE LIB SNB OPF CPF clpsat Comp

check = 0;

% check for data file
if isempty(File.data),
  fm_disp(['Set a data file before running ',msg,'.'],2)
  return
end

% check for initial power flow solution
if ~Settings.init
  solvepf
  if ~Settings.init, return, end
end

% check for dynamic components if running a static analysis
if ~type & DAE.n & ~clpsat.init
  dynlf = sum(prod(Comp.prop(:,[3 6 9]),2));
  iscpf = strcmp(msg,'Continuation Power Flow');
  if ~Settings.static & ~dynlf
    Settings.ok = 0;
    uiwait(fm_choice('Dynamic components will be be discarded. Continue?'))
    if Settings.ok
      Settings.static = 1;
      solvepf
      Settings.static = 0; % reset initial condition
    else
      return
    end
  elseif ~Settings.static & ~dynlf & iscpf
    Settings.ok = 0;
    uiwait(fm_choice(['Dynamic components can lead to numerical ' ...
                      'problems, discard?']))
    if Settings.ok
      Settings.static = 1;
      solvepf
      Settings.static = 0; % reset initial condition
    end
  elseif iscpf
    Settings.ok = 0;
    uiwait(fm_choice(['Dynamic components can lead to numerical ' ...
                      'problems, continue?']))
    if ~Settings.ok, return, end
  else
    uiwait(fm_choice(['Dynamic components are not supported for ' ...
                      'static analysis'],2))
    return
  end
end

% check for previous CPF & ATC solutions
if strcmp(msg,'SNB Direct Method')
  one = 1;
else
  one = 0;
end

if CPF.init & ~(one & CPF.init == 1)
  switch CPF.init
   case 1, met = 'CPF';
   case 2, met = 'ATC';
   case 3, met = 'N-1 Cont. An.';
   case 4, met = 'Continuation OPF (PSAT-GAMS)';
  end
  Settings.ok = 0;
  if clpsat.init
    Settings.ok = clpsat.refresh;
  else
    uiwait(fm_choice([met,' has been run last. Do you want to' ...
                      ' restore initial PF solution?']))
  end
  if Settings.ok
    solvepf
    fm_disp(['Initial PF solution will be used as ', ...
	     'base case solution.'])
  else
    fm_disp(['Last ',met,' solution will be used as ', ...
	     'base case solution.'])
  end
  CPF.init = 0;
end

% check for previous time domain simulations
if Settings.init == 2
  Settings.ok = 0;
  if clpsat.init
    Settings.ok = clpsat.refresh;
  else
    uiwait(fm_choice(['TD has been run last. Do you want to' ...
                      ' restore initial PF solution?']))
  end
  if Settings.ok
    solvepf
    fm_disp(['Initial PF solution will be used as ', ...
	     'base case solution.'])
  else
    fm_disp('Last TD point will be used as base case solution.')
  end
  Settings.init = 1;
end

% check for SNB direct method
if SNB.init
  Settings.ok = 0;
  if clpsat.init
    Settings.ok = clpsat.refresh;
  else
    uiwait(fm_choice(['SNB direct method has been run last. Do you want to' ...
                      ' restore initial PF solution?']))
  end
  if Settings.ok
    solvepf
    fm_disp(['Initial PF solution will be used as ', ...
	     'base case solution.'])
  else
    fm_disp('SNB solution will be used as base case solution.')
  end
  SNB.init = 0;
end

% check for LIB direct method
if LIB.init
  Settings.ok = 0;
  if clpsat.init
    Settings.ok = clpsat.refresh;
  else
    uiwait(fm_choice(['LIB direct method has been run last. Do you want to' ...
                      ' restore initial PF solution?']))
  end
  if Settings.ok
    solvepf
    fm_disp('Initial PF solution will be used as base case solution.')
  else
    fm_disp('LIB solution will be used as base case solution.')
  end
  LIB.init = 0;
end

% check for OPF solution
if strcmp(msg,'Optimal Power Flow')
  one = 0;
else
  one = 1;
end

if OPF.init & one
  Settings.ok = 0;
  if clpsat.init
    Settings.ok = clpsat.refresh;
  else
    uiwait(fm_choice(['OPF has been run last. Do you want to' ...
                      ' restore initial PF solution?']))
  end
  if Settings.ok
    solvepf
    fm_disp(['Initial PF solution will be used as ', ...
	     'base case solution.'])
  else
    fm_disp('OPF solution will be used as base case solution.')
  end
  OPF.init = 0;
end

check = 1;

% ---------------------------------------------------
function solvepf

global Settings Varname

fm_disp('Solve base case power flow...')
varname_old = Varname.idx;
Settings.show = 0;
fm_set('lf')
Settings.show = 1;
if ~isempty(varname_old)
  Varname.idx = varname_old;
end