home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-bin / lib / octave / 1.1.1 / m / control / is_controllable.m < prev    next >
Encoding:
Text File  |  1996-10-12  |  2.0 KB  |  74 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_controllable (a, b, tol)
  20.  
  21. # Usage: is_controllable (a, b {,tol})
  22. #
  23. # Returns 1 if the pair (a, b) is controllable, or 0 if not.
  24. #
  25. # See also: size, rows, columns, length, is_matrix, is_scalar, is_vector
  26. #
  27. # This should really use the method below, but I'm being lazy for now:
  28. #
  29. # Controllability is determined by applying Arnoldi iteration with
  30. # complete re-orthogonalization to obtain an orthogonal basis of the
  31. # Krylov subspace.
  32. #
  33. # (FIX ME... The Krylov subspace approach is not done yet!)
  34. #                      n-1
  35. #   span ([b,a*b,...,a^   b]).
  36. #
  37. # tol is a roundoff paramter, set to 2*eps if omitted.
  38.  
  39. # Written by A. S. Hodel (scotte@eng.auburn.edu) August, 1993.
  40.  
  41.   if (nargin == 2 || nargin == 3)
  42.  
  43.     n = is_square (a);
  44.     [nr, nc] = size (b);
  45.  
  46.     if (n == 0 || n != nr || nc == 0)
  47.       retval = 0;
  48.     else
  49.  
  50.       m = b;
  51.       tmp = b;
  52.       for ii = 1:(n-1)
  53.         tmp = a * tmp;
  54.         m = [m, tmp];
  55.       endfor
  56.  
  57. # If n is of any significant size, m will be low rank, so be careful!
  58.  
  59.       if (nargin == 3)
  60.         if (is_scalar (tol))
  61.           retval = (rank (m, tol) == n);
  62.         else
  63.           error ("is_controllable: tol must be a scalar");
  64.         endif
  65.       else
  66.         retval = (rank (m) == n);
  67.       endif
  68.     endif
  69.   else
  70.     usage ("is_controllable (a, b)");
  71.   endif
  72.  
  73. endfunction
  74.