fm_dynidx.m

用于电力系统的一个很好的分析软件

function fm_dynidx
% FM_DYNIDX define indices of state variables for components
%           which are not initialized during the power flow
%           analysis
%
% FM_DYNIDX
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Update:    17-Jul-2007
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
% Copyright (C) 2002-2008 Federico Milano

global DAE
global Syn Exc Tg Oxl Pss Fl Thload Svc Statcom Tcsc Sssc Upfc Mass SSR
global Sofc Cac Cluster Exload Wind Cswt Dfig Busfreq Pmu Pod COI Jimma
global Mixload Ddsg

m_old = DAE.m;

Syn = dynidx(Syn);
Exc = dynidx(Exc);
Tg = dynidx(Tg);
Oxl = dynidx(Oxl);
Pss = dynidx(Pss);
Fl = dynidx(Fl);
Thload = dynidx(Thload);
Svc = dynidx(Svc);
Statcom = dynidx(Statcom);
Tcsc = dynidx(Tcsc);
Sssc = dynidx(Sssc);
Upfc = dynidx(Upfc);
Mass = dynidx(Mass);
SSR = dynidx(SSR);
Sofc = dynidx(Sofc);
Cac = dynidx(Cac);
Cluster = dynidx(Cluster);
Exload = dynidx(Exload);
Wind = dynidx(Wind);
Cswt = dynidx(Cswt);
Dfig = dynidx(Dfig);
Ddsg = dynidx(Ddsg);
Busfreq = dynidx(Busfreq);
Pmu = dynidx(Pmu);
Pod = dynidx(Pod);
COI = dynidx(COI);
Jimma = dynidx(Jimma);
Mixload = dynidx(Mixload);

% rebuild algebraic variables and vectors
k = DAE.m-m_old;
DAE.y = [DAE.y;zeros(k,1)];
DAE.g = [DAE.g;zeros(k,1)];
DAE.Gy = [DAE.Gy,sparse(m_old,k);sparse(k,DAE.m)];