home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts / control / sysappend.m < prev    next >
Encoding:
Text File  |  1999-03-05  |  5.8 KB  |  167 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 retsys = sysappend(sys,b,c,d,outname,inname,yd)
  20.   # 
  21.   # function retsys = sysappend(sys,b[,c,d,outname,inname,yd])
  22.   #
  23.   # This function appends new inputs and/or outputs to a system
  24.   # Inputs:
  25.   #    sys:  system data structure
  26.   #    b: matrix to be appended to sys "B" matrix (empty if none)
  27.   #    c: matrix to be appended to sys "C" matrix (empty if none)
  28.   #     d: revised sys d matrix (can be passed as [] if the revised d is all 
  29.   #        zeros)
  30.   #     outname: names for new outputs
  31.   #     inname: names for new inputs
  32.   #     yd: indicator for which new outputs are continuous/discrete 
  33.   #         (yd(i) = 0 or , respectively)
  34.   # result:
  35.   #   sys.b := [sys.b , b]
  36.   #   sys.c := [sys.c  ]
  37.   #            [ c     ]
  38.   #   sys.d := [sys.d | D12 ]
  39.   #            [D21   | D22 ]
  40.   #         where D12, D21, and D22 are the appropriate dimensioned blocks
  41.   #         of the input parameter d.  The leading block D11 of d is ignored.
  42.   # If inname and outname are not given as arguments, the new inputs and 
  43.   # outputs are be assigned default names.  
  44.   # yd is a vector of length rows(c), and indicates which new outputs are
  45.   # discrete (yd(ii) = 1) and which are continuous (yd(ii) = 0).
  46.   # Default value for yd is:
  47.   #     sys = continuous or mixed: yd = zeros(1,rows(c))
  48.   #     sys = discrete:            yd = ones(1,rows(c))
  49.   
  50.   # written by John Ingram August 1996
  51.   
  52.   sav_implicit_str_to_num_ok = implicit_str_to_num_ok;    # save for later
  53.   sav_empty_list_elements_ok = empty_list_elements_ok;
  54.  
  55.   empty_list_elements_ok = 1;                implicit_str_to_num_ok = 1;
  56.   
  57.   # check input arguments
  58.   if ( (nargin < 2) | (nargin > 7) | (!is_struct(sys)))
  59.     usage("retsys = sysappend(sys,b,c[,d,outname,inname,yd]) ");
  60.   elseif(!is_struct(sys))
  61.     error("sys must be a system data structure");
  62.   endif
  63.   
  64.   # default system type must be state space form
  65.   [Aa,Ab,Ac,Ad,Ats,Ann,Anz,Ast,Ain,Aout,Ayd] = sys2ss(sys);
  66.   [Ann,Anz,Am,Ap] = sysdimensions(sys);
  67.  
  68.   #default c
  69.   if(nargin < 3)      c = [];                                endif
  70.   
  71.   #default d
  72.   if(nargin < 4)     make_d = 1;
  73.   elseif(isempty(d)) make_d = 1;
  74.   else               make_d = 0;                             endif
  75.   if(make_d)         d = zeros(rows(c)+Ap,columns(b) + Am);  endif
  76.  
  77.   #
  78.   # Append new input(s) if any
  79.   Bm = max(columns(d),columns(b)+Am);
  80.   if(Bm != Am)    
  81.     # construct new signal names
  82.     if(nargin >= 6)   # new names were passed
  83.       if(!isstr(inname))
  84.         error("inname must be a string");
  85.       elseif(rows(inname) != (Bm - Am))
  86.         error(sprintf("%d new inputs requested; inname(%dx%d)", ...
  87.       (Bm-Am),rows(inname),columns(inname)));
  88.       endif
  89.     else
  90.       inname = sysdefioname(Bm,"u",(Am+1));
  91.     endif
  92.     if(Am)   Ain = append(Ain,inname);
  93.     else     Ain = inname;        endif
  94.  
  95.     # default b matrix
  96.     if(isempty(b))     b  = zeros(Ann+Anz,(Bm-Am));          
  97.     elseif(rows(b) != Ann+Anz | columns(b) != (Bm-Am))
  98.         error(sprintf("b(%dx%d); should be (%dx%d)", rows(b), columns(b), ...
  99.           (Ann+Anz), (Bm-Am)));
  100.     endif
  101.  
  102.     # append new b matrix
  103.     Ab = [Ab,b];    # empty_list_elements_ok=1 makes this ok
  104.   endif
  105.  
  106.   #
  107.   # Append new output(s) if any
  108.   Bp = max(rows(d),rows(c)+Ap);
  109.   if(Bp != Ap)  
  110.  
  111.     # construct new signal names, output classification
  112.     if(nargin >= 5)  # new names were passed
  113.       if(!isstr(outname))
  114.         error("outname must be a string");
  115.       elseif(rows(outname) != (Bp - Ap))
  116.         error(sprintf("%d new outputs requested; outname(%dx%d)", ...
  117.           (Bp-Ap),rows(outname),columns(outname)));
  118.       endif
  119.     else
  120.       outname = sysdefioname(Bp,"y",(Ap+1));
  121.     endif
  122.     if(Ap)   Aout = append(Aout,outname);
  123.     else     Aout = outname;                endif
  124.  
  125.     # construct new yd entries
  126.     if(nargin == 7)
  127.       if(!is_vector(yd))
  128.         error(sprintf("yd(%dx%d) must be a vector",rows(yd),columns(yd)))
  129.       elseif(rows(c) != length(yd) & rows(d) != length(yd))
  130.         error(sprintf("length(yd) = %d; c(%dx%d), d(%dx%d); mismatch", ...
  131.       length(yd), rows(c), columns(c),rows(d),columns(d)));
  132.       endif
  133.     else
  134.       # default yd values
  135.       yd = ones(1,Bp)*( (Ats > 0) & (Ann == 0)  & isempty(find(Ayd == 0)) ) ;
  136.     endif
  137.     Ayd = [vec(Ayd);vec(yd)];
  138.  
  139.     # default c matrix
  140.     if(isempty(c))      c = zeros((Bp-Ap),Ann+Anz);          
  141.     elseif(columns(c) != Ann+Anz | rows(c) != (Bp-Ap))
  142.         error(sprintf("c(%dx%d); should be (%dx%d)", rows(c), columns(c), ...
  143.           (Bp-Ap), (Ann+Anz) ));
  144.     endif
  145.  
  146.     # append new c matrix
  147.     Ac = [Ac;c];    # empty_list_elements_ok=1 makes this ok
  148.   endif
  149.  
  150.   # check d matrix
  151.   if(isempty(d)) d = zeros(Bp,Bm);
  152.   elseif(rows(d) != Bp | columns(d) != Bm)
  153.     error(sprintf("d(%dx%d) should be (%dx%d)",rows(d), columns(d), Bp, Bp));
  154.   endif
  155.  
  156.   # Splice in original D matrix  
  157.   if(Am & Ap)          d(1:Ap, 1:Am) = Ad;       endif
  158.   Ad = d;
  159.   
  160.   # construct return system
  161.   retsys = ss2sys(Aa,Ab,Ac,Ad,Ats,Ann,Anz,Ast,Ain,Aout,find(Ayd == 1));
  162.   
  163.   implicit_str_to_num_ok = sav_implicit_str_to_num_ok;    # restore value
  164.   empty_list_elements_ok = sav_empty_list_elements_ok;
  165.  
  166. endfunction
  167.