fm_out.m

电力系统分析计算程序

function  fm_out(flag,t,k)
% FM_OUT define output variable vector during time domain simulations
%
% FM_OUT(FLAG,T,K)
%   T = actual time
%   K = time step number
%   FLAG = 0 ->  Output Structure Initialization
%	   1 ->  Memory Allocation
%          2 ->  k-th Step Assignment
%          3 ->  Vector Redimensioning
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
% Copyright (C) 2002-2008 Federico Milano

global Settings Bus DAE Varout Varname Fig

persistent compute_pijs
nb = Bus.n;

switch flag

 case 0

  if DAE.n > 5000 | nb > 5000
    Settings.chunk = 10;
  end
  if length(Varname.idx) > Settings.maxvar
    fm_disp('No more than 1500 variables will be stored for plotting.')
    Varname.idx = Varname.idx(1:Settings.maxvar);
  end

  Varout.t = zeros(Settings.chunk,1);
  Varout.idx = Varname.idx;
  idx0 = length(Varout.idx);
  Varout.vars = zeros(Settings.chunk,idx0);
  idx1 = DAE.n+DAE.m+2*nb;
  if ~isempty(find(Varout.idx > idx1))
    compute_pijs = idx1;
  else
    compute_pijs = 0;
  end
  % network visualization
  fm_threed('init')

 case 1

  Varout.t = [Varout.t; zeros(Settings.chunk,1)];
  idx0 = length(Varout.idx);
  Varout.vars = [Varout.vars; zeros(Settings.chunk,idx0)];

 case 2

  Varout.t(k) = t;

  fm_call('series')

  vars = [];

  if DAE.n
    vars = DAE.x;
  end

  if DAE.m
    vars = [vars; DAE.y; DAE.g(1:2*nb)];
  end

  if compute_pijs

    [ps,qs,pr,qr,bfr,bto] = fm_flows;

    ss = abs(ps + i*qs);
    sr = abs(pr + i*qr);
    is = ss./DAE.y(bfr);
    ir = sr./DAE.y(bto);

    vars = [vars; ps; pr; qs; qr; is; ir; ss; sr];

  end

  if Settings.vs
    idx0 = DAE.n+DAE.m+2*nb+8*Settings.nseries;
    Varout.vars(k,[1:length(Varout.idx)]) = ...
        vars(Varout.idx(find(Varout.idx <= idx0))).';
    fm_vs(k,idx0)
  else
    Varout.vars(k,:) = vars(Varout.idx).';
  end
  fm_threed('newframe')

 case 3

  Varout.t = Varout.t(1:k);
  Varout.vars = Varout.vars(1:k,:);
  fm_threed('finish')

end