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://thunderbox.uwaterloo.ca/~fmilano
%
% 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(k,k) = speye(length(u));
end