home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / sys2ss.m < prev    next >
Text File  |  1999-04-29  |  3KB  |  76 lines
  1. # Copyright (C) 1996, 1998 A. Scottedward Hodel 
  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 [a,b,c,d,tsam,n,nz,stname,inname,outname,yd] = sys2ss(sys)
  20.   # function [a,b,c,d(,tsam,n,nz,stname,inname,outname,yd)] = sys2ss(sys)
  21.   # Convert from system data structure to state space form
  22.   # inputs:
  23.   #    sys: system data structure for the state space system
  24.   #
  25.   #      x' = Ax + Bu
  26.   #      y  = Cx + Du
  27.   #
  28.   #  or a similar discrete-time system.  
  29.   #
  30.   # outputs:
  31.   #    a,b,c,d: state space matrices for sys
  32.   #    tsam: sampling time of sys (0 if continuous)
  33.   #    n, nz: number of continuous, discrete states (discrete states come
  34.   #          last in state vector x)
  35.   #    stname, inname, outname: signal names (strings);  names of states,
  36.   #          inputs, and outputs, respectively
  37.   #    yd: binary vector; yd(ii) is nonzero if output y is discrete.
  38.   # 
  39.   # A warning message is printed if the system is a mixed 
  40.   # continuous/discrete system.
  41.  
  42.   # Written by David Clem August 19, 1994
  43.   # Updates by John Ingram July 14, 1996
  44.  
  45.   if(nargin != 1)
  46.     usage("[a,b,c,d,tsam,n,nz,stname,inname,outname,yd] = sys2ss(sys)")
  47.   endif
  48.  
  49.   if (nargout > 11)
  50.     warning(["sys2ss: ",num2str(nargout)," out arguments exceeds max=11"])
  51.     usage("[a,b,c,d,tsam,n,nz,stname,inname,outname,yd] = sys2ss(sys)")
  52.   endif
  53.  
  54.   if( ! is_struct(sys) )
  55.     error("input argument must be a system data structure");
  56.   endif
  57.  
  58.   sys = sysupdat(sys,"ss");        # make sure state space data is there
  59.   [n,nz,m,p] = sysdimen(sys);
  60.   [stname,inname,outname,yd] = sysgetsg(sys);
  61.   tsam = sysgetts(sys);
  62.  
  63.   cont = sum(yd == 0) + n;
  64.   dig = sum(yd != 0) + nz + tsam;
  65.   if(cont*dig)
  66.     warning("sys2ss: input system is mixed continuous/discrete");
  67.   endif
  68.  
  69.   a = sys.a;
  70.   b = sys.b;
  71.   c = sys.c;
  72.   d = sys.d;
  73.  
  74. endfunction
  75.  
  76.