home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / starp.m < prev    next >
Text File  |  1999-04-29  |  3KB  |  120 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] = starp(P, K, ny, nu);
  20. #
  21. # [sys] = starp(P, K, ny, nu)
  22. #
  23. # Redheffer star product or upper/lower LFT, respectively.
  24. #
  25. #
  26. #               +-------+
  27. #     --------->|       |---------> 
  28. #               |   P   |
  29. #          +--->|       |---+  ny
  30. #          |    +-------+   |
  31. #          +-------------------+
  32. #                           |  |
  33. #          +----------------+  |
  34. #          |                   |
  35. #          |    +-------+      |
  36. #          +--->|       |------+ nu 
  37. #               |   K   |
  38. #     --------->|       |--------->
  39. #               +-------+
  40. #
  41. # If ny and nu "consume" all inputs and outputs of K then the result
  42. # is a lower fractional transformation. If ny and nu "consume" all
  43. # inputs and outputs of P then the result is an upper fractional
  44. # transformation.
  45. #
  46. # ny and/or nu may be negative (= negative feedback)
  47.  
  48. # Written by Kai Mueller May 1998
  49.  
  50.   if((nargin != 2) && (nargin != 4))
  51.     usage("[sys] = starp(P, K, ny, nu)");
  52.   endif
  53.   if (!is_struct(P))
  54.     error("---> P must be in system data structure");
  55.   endif
  56.   if (!is_struct(K))
  57.     error("---> K must be in system data structure");
  58.   endif
  59.  
  60.   P = sysupdat(P, "ss");
  61.   [n, nz, mp, pp] = sysdimen(P);
  62.   np = n + nz;
  63.   K = sysupdat(K, "ss");
  64.   [n, nz, mk, pk] = sysdimen(K);
  65.   nk = n + nz;
  66.   ny_sign = 1;
  67.   nu_sign = 1;
  68.   if (nargin == 2)
  69.     # perform a LFT of P and K (upper or lower)
  70.     ny = min([pp, mk]);
  71.     nu = min([pk, mp]);
  72.   else
  73.     if (ny < 0)
  74.       ny = -ny;
  75.       ny_sign = -1;
  76.     endif
  77.     if (nu < 0)
  78.       nu = -nu;
  79.       nu_sign = -1;
  80.     endif
  81.   endif
  82.   if (ny > pp)
  83.     error("---> P has not enough outputs.");
  84.   endif
  85.   if (nu > mp)
  86.     error("---> P has not enough inputs.");
  87.   endif
  88.   if (ny > mk)
  89.     error("---> K has not enough inputs.");
  90.   endif
  91.   if (nu > pk)
  92.     error("---> K has not enough outputs.");
  93.   endif
  94.   nwp  = mp - nu;
  95.   nzp  = pp - ny;
  96.   nwk  = mk - ny;
  97.   nzk  = pk - nu;
  98.   if ((nwp + nwk) < 1)
  99.     error("---> no inputs left for star product.");
  100.   endif
  101.   if ((nzp + nzk) < 1)
  102.     error("---> no outputs left for star product.");
  103.   endif
  104.  
  105.   # checks done, form sys
  106.   if (nzp)  Olst = [1:nzp];  endif
  107.   if (nzk)  Olst = [Olst, pp+nu+1:pp+pk];  endif
  108.   if (nwp)  Ilst = [1:nwp];  endif
  109.   if (nwk)  Ilst = [Ilst, mp+ny+1:mp+mk];  endif
  110.   Clst = zeros(ny+nu,2);
  111.   for ii = 1:nu
  112.     Clst(ii,:) = [nwp+ii, nu_sign*(pp+ii)];
  113.   endfor
  114.   for ii = 1:ny
  115.     Clst(nu+ii,:) = [mp+ii, ny_sign*(nzp+ii)];
  116.   endfor
  117.   sys = buildssc(Clst,[],Olst,Ilst,P,K);
  118.  
  119. endfunction
  120.