home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / is_stable.m < prev    next >
Text File  |  1999-03-05  |  2KB  |  64 lines
  1. # Copyright (C) 1993, 1994, 1995 John W. Eaton
  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 retval = is_stable (a, tol, disc)
  20.  
  21. # Usage: retval = is_stable (a {,tol,disc})
  22. # or     retval = is_stable (sys{,tol})
  23. #
  24. # Returns retval = 1 if the matrix a or the system a is stable, or 0 if not.
  25. #
  26. # tol is a roundoff paramter, set to 200*eps if omitted.
  27. # disc != 0: stable if eig(a) in unit circle
  28. #         0: stable if eig(a) in open LHP (default)
  29. #
  30. # See also: size, rows, columns, length, is_matrix, is_scalar, is_vector
  31. #     is_observable, is_stabilizable, is_detectable, krylov, krylovb
  32.  
  33. # Written by A. S. Hodel (scotte@eng.auburn.edu) August, 1993.
  34. # Updated by John Ingram (ingraje@eng.auburn.edu) July, 1996 for systems
  35. # Updated to simpler form by a.s.hodel 1998
  36.  
  37.   if( (nargin < 1) | (nargin > 3) )   usage("is_stable(a {,tol,disc})");
  38.   elseif(is_struct(a))
  39.     # system was passed
  40.     if(nargin < 3)            disc = is_digital(a);
  41.     elseif(disc != is_digital(a))
  42.       warning("is_stable: disc =%d does not match system",disc)
  43.     endif
  44.     sys = sysupdate(a,"ss");
  45.     a = sys2ss(sys);
  46.   else
  47.     if(nargin < 3)        disc = 0;        endif
  48.     if(is_square(a) == 0)
  49.       error("A(%dx%d) must be square",rows(A), columns(A));
  50.     endif
  51.   endif
  52.  
  53.   if(nargin < 2)        tol = 200*eps;
  54.   elseif( !is_scalar(tol) )
  55.     error("is_stable: tol(%dx%d) must be a scalar",rows(tol),columns(tol));
  56.   endif
  57.  
  58.   l = eig(a);
  59.   if(disc)    nbad = sum(abs(l)*(1+tol) > 1);
  60.   else        nbad = sum(real(l)+tol > 0);        endif
  61.   retval = (nbad == 0);   
  62.  
  63. endfunction
  64.