home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21fb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / is_detec.m < prev    next >
Text File  |  1999-12-24  |  2KB  |  63 lines
  1. ## Copyright (C) 1993, 1994, 1995 Auburn University.  All Rights Reserved
  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, 59 Temple Place, Suite 330, Boston, MA 02111 USA.
  18.  
  19. ## -*- texinfo -*-
  20. ## @deftypefn {Function File } { [@var{retval}, @var{U}] =} is_detec (@var{a}, @var{c}@{, @var{tol}@})
  21. ## @deftypefnx {Function File } { [@var{retval}, @var{U}] =} is_detec (@var{sys}@{, @var{tol}@})
  22. ## Test for detactability (observability of unstable modes) of (@var{a},@var{c}).  
  23. ## 
  24. ##  Returns 1 if the system @var{a} or the pair (@var{a},@var{c})is 
  25. ##  detectable, 0 if not.
  26. ## 
  27. ## @strong{See} @code{is_stabi} for detailed description of arguments and
  28. ## computational method.
  29. ## 
  30. ##  Default: tol = 10*norm(a,'fro')*eps 
  31. ## 
  32. ## @end deftypefn
  33.  
  34. ## See also: size, rows, columns, length, is_mat, is_scal, is_vec.
  35.  
  36. function [retval, U] = is_detec (a, c, tol)
  37.  
  38.   ## Written by A. S. Hodel (scotte@eng.auburn.edu) August 1993.
  39.   ## Updated by John Ingram (ingraje@eng.auburn.edu) July 1996.
  40.  
  41.   if( nargin < 1) 
  42.     usage("[retval,U] = is_detec(a , c {, tol})");
  43.   elseif(is_struct(a))
  44.     ## system form
  45.     if(nargin == 2)
  46.       tol = c;
  47.     elseif(nargin > 2)
  48.       usage("[retval,U] = is_detec(sys {, tol})");
  49.     endif
  50.     [a,b,c] = sys2ss(a);
  51.   elseif(nargin > 3)
  52.     usage("[retval,U] = is_detec(a , c {, tol})");
  53.   endif
  54.   if(exist("tol"))
  55.     [retval,U] = is_stabi (a', c', tol);
  56.   else
  57.     [retval,U] = is_stabi (a', c');
  58.   endif
  59.   
  60.  
  61. endfunction
  62.  
  63.