Mercurial > octave-nkf
view scripts/control/hinf/wgt1o.m @ 4771:b8105302cfe8
[project @ 2004-02-16 17:45:50 by jwe]
author | jwe |
---|---|
date | Mon, 16 Feb 2004 17:45:50 +0000 |
parents | 7923abdeb4e5 |
children | bdbee5282954 |
line wrap: on
line source
## Copyright (C) 1998 Kai P. Mueller ## ## This file is part of Octave. ## ## Octave 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, or (at your option) any ## later version. ## ## Octave is distributed in the hope that it will be useful, but WITHOUT ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 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 Octave; see the file COPYING. If not, write to the Free ## Software Foundation, 59 Temple Place, Suite 330, Boston, MA 02111 USA. ## -*- texinfo -*- ## @deftypefn {Function File} {} wgt1o (@var{vl}, @var{vh}, @var{fc}) ## State space description of a first order weighting function. ## ## Weighting function are needed by the H2/H_infinity design procedure. ## These function are part of thye augmented plant P (see hinfdemo ## for an applicattion example). ## ## vl = Gain at low frequencies ## ## vh = Gain at high frequencies ## ## fc = Corner frequency (in Hz, *not* in rad/sec) ## @end deftypefn ## Author: Kai P. Mueller <mueller@ifr.ing.tu-bs.de> ## Created: September 30, 1997 function wsys = wgt1o (vl, vh, fc) if (nargin != 3) usage("wsys = wgt1o(vl, vh, fc)"); endif if(nargout > 1) usage("wsys = wgt1o(vl, vh, fc)"); endif if (vl == vh) a = []; b = []; c = []; else a = [-2*pi*fc]; b = [-2*pi*fc]; c = [vh-vl]; endif d=[vh]; wsys = ss(a,b,c,d); endfunction