home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / parallel.m < prev    next >
Text File  |  1999-03-05  |  2KB  |  65 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 sysp = parallel(Asys,Bsys)
  20. # function sysp = parallel(Asys,Bsys)
  21. # Forms the parallel connection of two systems.
  22. #
  23. #              ____________________
  24. #              |      ________    |
  25. #     u  ----->|----> | Asys |--->|----> y1
  26. #         |    |      --------    |
  27. #         |    |      ________    |
  28. #         |--->|----> | Bsys |--->|----> y2
  29. #              |      --------    |
  30. #              --------------------
  31. #                   Ksys
  32.  
  33. # Written by David Clem August 15, 1994
  34. # completely rewritten Oct 1996 a s hodel
  35. # SYS_INTERNAL accesses members of system structure
  36.  
  37.   if(nargin != 2)
  38.     usage("sysp = parallel(Asys,Bsys)");
  39.   endif
  40.   if(! is_struct(Asys) )
  41.     error("1st input argument is not a system data structure")
  42.   elseif (! is_struct(Bsys) )
  43.     error("2nd input argument is not a system data structure")
  44.   endif
  45.   [Ann,Anz,mA] = sysdimensions(Asys);
  46.   [Bnn,Bnz,mB] = sysdimensions(Bsys);
  47.   if(mA != mB)
  48.     error(["Asys has ",num2str(mA)," inputs, Bsys has ",num2str(mB)," inputs"]);
  49.   endif
  50.  
  51.   # save signal names
  52.   Ain = sysgetsignals(Asys,"in");
  53.  
  54.   # change signal names to avoid warning messages from sysgroup
  55.   Asys = syssetsignals(Asys,"in",sysdefioname(length(Ain),"Ain_u"));
  56.   Bsys = syssetsignals(Bsys,"in",sysdefioname(length(Ain),"Bin_u"));
  57.  
  58.   sysp = sysgroup(Asys,Bsys);
  59.   sysD = ss2sys([],[],[],[eye(mA);eye(mA)]);
  60.   
  61.   sysp = sysmult(sysp,sysD);
  62.   sysp = syssetsignals(sysp,"in",Ain);
  63.   
  64. endfunction
  65.