home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / is_abcd.m < prev    next >
Text File  |  1999-03-05  |  3KB  |  97 lines
  1. # Copyright (C) 1997 Kai P. Mueller
  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_abcd(a, b, c, d)
  20.   # ------------------------------------------------------
  21.   # retval = is_abcd(a [, b, c, d])
  22.   # Returns retval = 1 if the dimensions of a, b, c, d
  23.   # are compatible, otherwise retval = 0.
  24.   # The matrices b, c, or d may be omitted.
  25.   # ------------------------------------------------------
  26.   # 
  27.   # see also: abcddim
  28.  
  29.   # Written by Kai P. Mueller November 4, 1997
  30.   # based on is_controllable.m of Scottedward Hodel
  31.   # modified by
  32.  
  33.   retval = 0;
  34.   switch (nargin)
  35.     case (1)
  36.       # A only
  37.       [na, ma] = size(a);
  38.       if (na != ma)
  39.         disp("Matrix A ist not square.")
  40.       endif
  41.     case (2)
  42.       # A, B only
  43.       [na, ma] = size(a);  [nb, mb] = size(b);
  44.       if (na != ma)
  45.         disp("Matrix A ist not square.")
  46.     return;
  47.       endif
  48.       if (na != nb)
  49.         disp("A and B column dimension different.")
  50.         return;
  51.       endif
  52.     case (3)
  53.       # A, B, C only
  54.       [na, ma] = size(a);  [nb, mb] = size(b);  [nc, mc] = size(c);
  55.       if (na != ma)
  56.         disp("Matrix A ist not square.")
  57.     return;
  58.       endif
  59.       if (na != nb)
  60.         disp("A and B column dimensions not compatible.")
  61.     return;
  62.       endif
  63.       if (ma != mc)
  64.         disp("A and C row dimensions not compatible.")
  65.     return;
  66.       endif
  67.     case (4)
  68.       # all matrices A, B, C, D
  69.       [na, ma] = size(a);  [nb, mb] = size(b);
  70.       [nc, mc] = size(c);  [nd, md] = size(d);
  71.       if (na != ma)
  72.         disp("Matrix A ist not square.")
  73.     return;
  74.       endif
  75.       if (na != nb)
  76.         disp("A and B column dimensions not compatible.")
  77.     return;
  78.       endif
  79.       if (ma != mc)
  80.         disp("A and C row dimensions not compatible.")
  81.     return;
  82.       endif
  83.       if (mb != md)
  84.         disp("B and D row dimensions not compatible.")
  85.     return;
  86.       endif
  87.       if (nc != nd)
  88.         disp("C and D column dimensions not compatible.")
  89.     return;
  90.       endif
  91.     otherwise
  92.       usage("retval = is_abcd(a [, b, c, d])")
  93.   endswitch
  94.   # all tests passed, signal ok.
  95.   retval = 1;
  96. endfunction
  97.