home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / lsim.m < prev    next >
Text File  |  1999-03-05  |  3KB  |  96 lines
  1. # Copyright (C) 1996 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 [y,x] = lsim(sys,u,t,x0)
  20. # lsim: Produce output for a linear simulation of a system
  21. #
  22. # lsim(sys,u,t,[x0])
  23. # Produces a plot for the output of the system, sys.
  24. #
  25. # U is an array that contains the system's inputs.  Each column in u 
  26. # corresponds to a different time step.  Each row in u corresponds to a 
  27. # different input.  T is an array that contains the time index of the 
  28. # system.  T should be regularly spaced.  If initial conditions are required
  29. # on the system, the x0 vector should be added to the argument list.
  30. #
  31. # When the lsim function is invoked with output parameters:
  32. # [y,x] = lsim(sys,u,t,[x0])
  33. # a plot is not displayed, however, the data is returned in y = system output
  34. # and x = system states.
  35.  
  36. # Written by David Clem, A. S. Hodel July 1995
  37. # modified by John Ingram for system format August 1996
  38.  
  39.  
  40.   if((nargin < 3)||(nargin > 4))
  41.     usage("[y,x] = lsim(sys,u,t[,x0])");
  42.   endif
  43.  
  44.   if(!is_struct(sys))
  45.     error("sys must be in system data structure");
  46.   endif
  47.  
  48.   sys = sysupdate(sys,"ss");
  49.  
  50.   [ncstates, ndstates, nin, nout] = sysdimensions(sys);
  51.   [a,b,c,d] = sys2ss(sys);
  52.   
  53.   if (nargin == 3)     x0 = zeros(columns(a),1);        endif
  54.  
  55.   if(rows(u) ~= length(t))
  56.     error("lsim: There should be an input value (row) for each time instant");
  57.   endif
  58.   if(columns(u) ~= columns(d))
  59.     error("lsim: U and d should have the same number of inputs");
  60.   endif
  61.   if(columns(x0) > 1)
  62.     error("lsim: Initial condition vector should have only one column");
  63.   endif
  64.   if(rows(x0) > rows(a))
  65.     error("lsim: Initial condition vector is too large");
  66.   endif
  67.  
  68.   Ts = 0;
  69.   t(2)-t(1);
  70.   u=u';
  71.   n = max(size(t));
  72.   for ii = 1:(n-1)
  73.  
  74.     # check if step size changed
  75.     if (t(ii+1) - t(ii) != Ts)
  76.       Ts = t(ii+1) - t(ii);
  77.       # [F,G] = c2d(a,b,Ts);
  78.       dsys = c2d(sys, Ts);
  79.       [F,G] = sys2ss(dsys);
  80.     endif
  81.  
  82.     x(:,ii) = x0;
  83.     x0 = F*x0 + G*u(:,ii);
  84.   endfor
  85.  
  86.   # pick up last point
  87.   x(:,n) = x0;
  88.  
  89.   y = c*x + d*u;
  90.   if(nargout == 0)
  91.    plot(t,y);
  92.    y=[];
  93.    x=[];
  94.   endif
  95. endfunction
  96.