fm_windup.m

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

function fm_windup(idx,xmax,xmin)
% FM_WINDUP setup non-windup limiter for TDs
%
% FM_WINDUP(IDX,XMAX,XMIN)
%       IDX  = index of state variables
%       XMAX = state variable upper limit
%       XMIN = state variable lower limit
%
%Author:    Federico Milano
%Date:      08-Mar-2006
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
% Copyright (C) 2006 Federico Milano

global DAE

x = DAE.x(idx);
u = find((x >= xmax | x <= xmin) & DAE.f(idx) == 0);

if ~isempty(u)
  k = idx(u);
  DAE.tn(k) = 0;
  DAE.Ac(k,:) = 0;
  DAE.Ac(:,k) = 0;
  DAE.Ac = DAE.Ac - sparse(k,k,1,DAE.m+DAE.n,DAE.m+DAE.n);
end