home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / abcddim.m next >
Text File  |  1999-04-29  |  3KB  |  99 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. # Note: n = 0 (pure gain block) is returned without warning.
  31. #
  32. # Returns n = m = p = -1 if the system is not compatible.
  33. #
  34. # See also: is_abcd
  35.  
  36. # Written by A. S. Hodel (scotte@eng.auburn.edu) August 1993.
  37. # a s hodel: modified to accept pure-gain systems aug 1996
  38.  
  39.   if (nargin != 4)
  40.     error ("abcddim: four arguments required");
  41.   endif
  42.  
  43.   n = m = p = -1;
  44.  
  45.   [a,an,am] = abcddims(a);
  46.   [b,bn,bm] = abcddims(b);
  47.   [c,cn,cm] = abcddims(c);
  48.   [d,dn,dm] = abcddims(d);
  49.  
  50.   if ( (!is_sqr(a)) & (!isempty(a)) )
  51.     warning (["abcddim: a is not square (",num2str(an),"x",num2str(am),")"]);
  52.     return
  53.   endif
  54.  
  55.   if( (bm == 0) & (dm == 0) )
  56.     warning("abcddim: no inputs");
  57.   elseif (bn != am)
  58.     warning (["abcddim: a(",num2str(an),"x",num2str(am), ...
  59.       " and b(",num2str(bn),"x",num2str(bm),") are not compatible"]);
  60.     return
  61.   endif
  62.  
  63.   if( (cn == 0) & (dn == 0 ) )
  64.     warning("abcddim: no outputs");
  65.   elseif (cm != an)
  66.     warning (["abcddim: a(",num2str(an),"x",num2str(am), ...
  67.     " and c(",num2str(cn),"x",num2str(cm),") are not compatible"]);
  68.     return
  69.   endif
  70.  
  71.   have_connections = (bn*cn != 0);
  72.  
  73.   if( (dn == 0) & have_connections)
  74.     warning("abcddim: empty d matrix passed; setting compatibly with b, c");
  75.     [d,dn,dm] = abcddims(zeros(cn,bm));
  76.   endif
  77.  
  78.   if(an > 0)
  79.     [dn, dm] = size(d);
  80.     if ( (cn != dn) & have_connections )
  81.       warning (["abcddim: c(",num2str(cn),"x",num2str(cm), ...
  82.     " and d(",num2str(dn),"x",num2str(dm),") are not compatible"]);
  83.       return
  84.     endif
  85.  
  86.     if ( (bm != dm) & have_connections )
  87.       warning (["abcddim: b(",num2str(bn),"x",num2str(bm), ...
  88.       " and d(",num2str(dn),"x",num2str(dm),") are not compatible"]);
  89.       return
  90.     endif
  91.  
  92.     m = bm;
  93.     p = cn;
  94.   else
  95.     [p,m] = size(d);
  96.   endif
  97.   n = an;
  98. endfunction
  99.