home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / syssub.m < prev    next >
Text File  |  1999-04-29  |  3KB  |  97 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 sys = syssub(Gsys,Hsys)
  20. # 
  21. # [sys] = syssub(Gsys,Hsys)
  22. #
  23. #
  24. # returns transfer functin sys = Gsys - Hsys
  25. #
  26. # Method: Gsys and Hsys are connected in parallel
  27. # The vector are connected to both systems; the outputs will be 
  28. # subtracted.  The names given to the system will be the G systems names.
  29. #
  30. #                  ________
  31. #             ----|  Gsys  |---
  32. #        u   |    ----------  +|         
  33. #        -----                (_)----> y
  34. #            |     ________   -|
  35. #             ----|  Hsys  |---
  36. #                  --------
  37.  
  38. # Written by John Ingram July 1996
  39.  
  40.   save_val = implicit_str_to_num_ok;    # save for later
  41.   implicit_str_to_num_ok = 1;
  42.  
  43.   if(nargin != 2)
  44.     usage("syssub:  [sys] = syssub(Gsys,Hsys)");
  45.   endif
  46.  
  47.   # check inputs
  48.   if(!is_struct(Gsys) | !is_struct(Hsys))
  49.     error("Both Gsys and Hsys must be a system data structure");
  50.   endif
  51.  
  52.   # check for compatibility
  53.   [n,nz,mg,pg] = sysdimen(Gsys);
  54.   [n,nz,mh,ph] = sysdimen(Hsys);
  55.   if(mg != mh)
  56.     error(sprintf("Gsys inputs(%d) != Hsys inputs (%d)",mg,mh));
  57.   elseif(pg != ph)
  58.     error(sprintf("Gsys outputs(%d) != Hsys outputs (%d)",pg,ph));
  59.   endif
  60.  
  61.   [Gst, Gin, Gout, Gyd] = sysgetsg(Gsys);
  62.   [Hst, Hin, Hout, Hyd] = sysgetsg(Hsys);
  63.  
  64.   # check for digital to continuous addition
  65.   if (Gyd != Hyd)
  66.     error("can not add a discrete output to a continuous output");
  67.   endif
  68.  
  69.   if( strcmp(sysgetty(Gsys),"tf") | strcmp(sysgetty(Hsys),"tf") )
  70.     # see if adding  transfer functions with identical denominators
  71.     [Gnum,Gden,GT,Gin,Gout] = sys2tf(Gsys);
  72.     [Hnum,Hden,HT,Hin,Hout] = sys2tf(Hsys);
  73.     if( (Hden == Gden) & (HT == GT) )
  74.       sys = tf2sys(Gnum-Hnum,Gden,GT,Gin,Gout);
  75.       return
  76.     endif
  77.     # if not, we go on and do the usual thing...
  78.   endif
  79.  
  80.   # make sure in ss form
  81.   Gsys = sysupdat(Gsys,"ss");
  82.   Hsys = sysupdat(Hsys,"ss");
  83.  
  84.   # change signal names to avoid warning messages from sysgroup
  85.   Gsys = syssetsg(Gsys,"in",sysdefio(length(Gin),"Gin_u"));
  86.   Gsys = syssetsg(Gsys,"out",sysdefio(length(Gout),"Gout_u"));
  87.   Hsys = syssetsg(Hsys,"in",sysdefio(length(Hin),"Hin_u"));
  88.   Hsys = syssetsg(Hsys,"out",sysdefio(length(Hout),"Hout_u"));
  89.   
  90.   sys = sysgroup(Gsys,Hsys);
  91.  
  92.   eyin = eye(mg); eyout = eye(pg);
  93.  
  94.   sys = sysscale(sys,[eyout -eyout],[eyin;eyin],Gout,Gin);
  95.  
  96. endfunction
  97.