setx0.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 36 行

function a = setx0(a)

global DAE PQ Bus

if ~a.n, return, end

% parameter initialization

% dat:
%  col #1: P0
%  col #2: Q0
%  col #3: V0

a.dat(:,3) = DAE.y(a.vbus);

for i = 1:a.n
  idx = findbus(PQ,a.bus(i));
  if isempty(idx)
    fm_disp(['No PQ load found for initializing Exponential ', ...
             'Recovery Load at bus ',Bus.names{a.bus(i)}])
  else
    PQ = pqactual(PQ,idx);
    a.dat(i,1) = PQ.P0(idx);
    a.dat(i,2) = PQ.Q0(idx);
    if a.u(i), PQ = remove(PQ,idx,'force'); end
  end
end

% state variable initialization
DAE.x(a.xp) = 0;
DAE.x(a.xq) = 0;

% message
fm_disp('Initialization of Exponential Recovery Loads completed.')