fm_dynlf.m

电力系统的psat

function fm_dynlf
% FM_DYNLF define state variable indices for components which are
%          included in power flow analysis
%
% FM_DYNLF
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2005 Federico Milano
%
% This toolbox is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2.0 of the License, or
% (at your option) any later version.
%
% This toolbox is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANDABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
% General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this toolbox; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307,
% USA.

global DAE

DAE.n = 0;

global Mot
if Mot.n
  Mot.slip = zeros(Mot.n,1);
  Mot.e1r = zeros(Mot.n,1);
  Mot.e1m = zeros(Mot.n,1);
  Mot.e2r = zeros(Mot.n,1);
  Mot.e2m = zeros(Mot.n,1);
end
for i = 1:Mot.n
  mot_ord = Mot.con(i,5);
  switch mot_ord
   case 1
    Mot.slip(i) = DAE.n + 1;
    DAE.n = DAE.n+1;
   case 3
    Mot.slip(i) = DAE.n + 1;
    Mot.e1r(i) = DAE.n + 2;
    Mot.e1m(i) = DAE.n + 3;
    DAE.n = DAE.n+3;
   case 5
    Mot.slip(i) = DAE.n + 1;
    Mot.e1r(i) = DAE.n + 2;
    Mot.e1m(i) = DAE.n + 3;
    Mot.e2r(i) = DAE.n + 4;
    Mot.e2m(i) = DAE.n + 5;
    DAE.n = DAE.n+5;
  end
end

global Ltc
if Ltc.n
  for  i = 1:Ltc.n
    Ltc.m(i) = DAE.n + i;
  end
  DAE.n = DAE.n + Ltc.n;
end

global Tap
if Tap.n
  for i = 1:Tap.n
    Tap.m(i) = DAE.n + i;
  end
  DAE.n = DAE.n + Tap.n;
end

global SAE1
if SAE1.n
  for i = 1:SAE1.n
    SAE1.m(i) = DAE.n + i;
  end
  DAE.n = DAE.n + SAE1.n;
end

global SAE2
if SAE2.n
  for i = 1:SAE2.n
    SAE2.m1(i) = DAE.n + 1;
    SAE2.m2(i) = DAE.n + 2;
    DAE.n = DAE.n + 2;
  end
end

global SAE3
if SAE3.n
  for i = 1:SAE3.n
    SAE3.m1(i) = DAE.n + 1;
    SAE3.m2(i) = DAE.n + 2;
    SAE3.m3(i) = DAE.n + 3;
    DAE.n = DAE.n + 3*SAE3.n;
  end
end

global Hvdc
if Hvdc.n
  Hvdc.Id = zeros(Hvdc.n,1);
  Hvdc.xr = zeros(Hvdc.n,1);
  Hvdc.xi = zeros(Hvdc.n,1);
  for i = 1:Hvdc.n
    Hvdc.Id(i) = DAE.n + 1;
    Hvdc.xr(i) = DAE.n + 2;
    Hvdc.xi(i) = DAE.n + 3;
    DAE.n = DAE.n + 3;
  end
end

global Statcom
if Statcom.n
  for i = 1:Statcom.n
    Statcom.Vdc(i) = DAE.n + 1;
    Statcom.alpha(i) = DAE.n + 2;
    Statcom.m(i) = DAE.n + 3;
    Statcom.Vmdc(i) = DAE.n + 4;
    Statcom.Vmac(i) = DAE.n + 5;
    DAE.n = DAE.n + 5;
  end
end

global Sssc
if Sssc.n
  for i = 1:Sssc.n
    Sssc.Vdc(i) = DAE.n + 1;
    Sssc.beta(i) = DAE.n + 2;
    Sssc.m(i) = DAE.n + 3;
    Sssc.Pmac(i) = DAE.n + 4;
    Sssc.Vmdc(i) = DAE.n + 5;
    DAE.n = DAE.n + 5;
  end
end

global Upfc
if Upfc.n
  for i = 1:Upfc.n
    Upfc.Vmdc(i) = DAE.n + 1;
    Upfc.Vmac(i) = DAE.n + 2;
    Upfc.Vdc(i) = DAE.n + 3;
    Upfc.alpha(i) = DAE.n + 4;
    Upfc.m(i) = DAE.n + 5;
    Upfc.x1(i) = DAE.n + 6;
    Upfc.x2(i) = DAE.n + 7;
    Upfc.Imd(i) = DAE.n + 8;
    Upfc.Imq(i) = DAE.n + 9;
    DAE.n = DAE.n + 9;
  end
end

global RLC
if RLC.n
  RLC.vd = zeros(DAE.n,1);
  RLC.vq = zeros(DAE.n,1);
  RLC.id = zeros(DAE.n,1);
  RLC.iq = zeros(DAE.n,1);
  for i = 1:RLC.n
    if RLC.con(i,7) == 0; RLC.con(i,7) = 1e-6; end
    if RLC.con(i,8) == 0; RLC.con(i,8) = 1e-6; end
    if RLC.con(i,9) == 0; RLC.con(i,9) = 1e-6; end
    ord = RLC.con(i,5);
    if ord == 1
      RLC.vd(i) = DAE.n + 1;
      RLC.vq(i) = DAE.n + 2;
      DAE.n = DAE.n + 2;
    elseif ord == 2
      RLC.id(i) = DAE.n + 1;
      RLC.iq(i) = DAE.n + 2;
      DAE.n = DAE.n + 2;
    else
      RLC.vd(i) = DAE.n + 1;
      RLC.vq(i) = DAE.n + 2;
      RLC.id(i) = DAE.n + 3;
      RLC.iq(i) = DAE.n + 4;
      DAE.n = DAE.n + 4;
    end
  end
end

global Phs
if Phs.n
  for i = 1:Phs.n
    Phs.alpha(i) = DAE.n + 1;
    Phs.Pm(i) = DAE.n + 2;
    DAE.n = DAE.n + 2;
  end
end