home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / h2norm.m < prev    next >
Text File  |  1999-03-05  |  2KB  |  59 lines
  1. # Copyright (C) 1996,1998 A. Scottedward Hodel 
  2. #
  3. # This file is part of Octave. 
  4. #
  5. # Octave is free software; you can redistribute it and/or modify it 
  6. # under the terms of the GNU General Public License as published by the 
  7. # Free Software Foundation; either version 2, or (at your option) any 
  8. # later version. 
  9. # 
  10. # Octave is distributed in the hope that it will be useful, but WITHOUT 
  11. # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 
  12. # FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License 
  13. # for more details.
  14. # 
  15. # You should have received a copy of the GNU General Public License 
  16. # along with Octave; see the file COPYING.  If not, write to the Free 
  17. # Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
  18.  
  19. function h2gain = h2norm(sys)
  20.   # Usage: h2gain = h2norm(sys)
  21.   #
  22.   # Computes the H2 norm system data structure (continuous time only)
  23.   # sys = system data structure [see ss2sys()]
  24.   # returns h2gain = Inf if system is unstable
  25.   #
  26.   # Reference:
  27.   # Doyle, Glover, Khargonekar, Francis, "State Space Solutions to Standard
  28.   # H2 and Hinf Control Problems", IEEE TAC August 1989
  29.   #
  30.  
  31.   # A. S. Hodel Aug 1995
  32.   # updated for system data structure by John Ingram November 1996
  33.  
  34.   if((nargin != 1))
  35.     usage("h2gain = h2norm(sys)");
  36.   elseif(!is_struct(sys))
  37.     error("Sys must be in system data structure");
  38.   end
  39.   dflg = is_digital(sys);
  40.  
  41.   if(!is_stable(sys))
  42.     warning("h2norm: unstable input system; returning Inf");
  43.     h2gain = Inf;
  44.   else
  45.     # compute gain
  46.     [a,b,c,d] = sys2ss(sys);
  47.     if(dflg)
  48.       M = dlyap(a,b*b');
  49.     else
  50.       M = lyap (a,b*b');
  51.     endif
  52.     if( min(real(eig(M))) < 0)
  53.       error("h2norm: grammian not >= 0 (lightly damped modes?)")
  54.     endif
  55.  
  56.     h2gain = sqrt(trace(d'*d + c*M*c'));
  57.   endif
  58. endfunction
  59.