home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-bin / lib / octave / 1.1.1 / m / control / abcddim.m next >
Encoding:
Text File  |  1996-10-12  |  1.8 KB  |  69 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 [n, m, p] = abcddim (a, b, c, d)
  20.  
  21. # Usage: [n, m, p] = abcddim (a, b, c, d)
  22. #
  23. # Check for compatibility of the dimensions of the matrices defining
  24. # the linear system (a, b, c, d).
  25. #
  26. # Returns n = number of system states,
  27. #         m = number of system inputs,
  28. #         p = number of system outputs.
  29. #
  30. # Returns n = m = p = -1 if the system is not compatible.
  31.  
  32. # Written by A. S. Hodel (scotte@eng.auburn.edu) August 1993.
  33.  
  34.   if (nargin != 4)
  35.     error ("abcddim: invalid number of arguments.  need four.")
  36.   endif
  37.  
  38.   n = m = p = -1;
  39.  
  40.   [an, am] = size(a);
  41.   if (an != am)
  42.     error ("abcddim: a is not square");
  43.   endif
  44.  
  45.   [bn, bm] = size(b);
  46.   if (bn != am)
  47.     error ("abcddim: a and b are not compatible");
  48.   endif
  49.  
  50.   [cn, cm] = size(c);
  51.   if (cm != an)
  52.     error ("abcddim: a and c are not compatible");
  53.   endif
  54.  
  55.   [dn, dm] = size(d);
  56.   if (cn != dn)
  57.     error ("abcddim: c and d are not compatible");
  58.   endif
  59.  
  60.   if (bm != dm)
  61.     error ("abcddim: b and d are not compatible");
  62.   endif
  63.  
  64.   n = an;
  65.   m = bm;
  66.   p = cn;
  67.  
  68. endfunction
  69.