|
3430
|
1 ## Copyright (C) 1996, 1998 Auburn University. All rights reserved. |
|
|
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, 59 Temple Place, Suite 330, Boston, MA 02111 USA. |
|
|
18 |
|
|
19 ## -*- texinfo -*- |
|
3500
|
20 ## @deftypefn {Function File} {} sysappend (@var{sys}, @var{b}, @var{c}, @var{d}, @var{outname}, @var{inname}, @var{yd}) |
|
3430
|
21 ## appends new inputs and/or outputs to a system |
|
|
22 ## |
|
|
23 ## @strong{Inputs} |
|
|
24 ## @table @var |
|
|
25 ## @item sys |
|
|
26 ## system data structure |
|
|
27 ## |
|
|
28 ## @item b |
|
|
29 ## matrix to be appended to sys "B" matrix (empty if none) |
|
|
30 ## |
|
|
31 ## @item c |
|
|
32 ## matrix to be appended to sys "C" matrix (empty if none) |
|
|
33 ## |
|
|
34 ## @item d |
|
|
35 ## revised sys d matrix (can be passed as [] if the revised d is all zeros) |
|
|
36 ## |
|
|
37 ## @item outname |
|
|
38 ## list of names for new outputs |
|
|
39 ## |
|
|
40 ## @item inname |
|
|
41 ## list of names for new inputs |
|
|
42 ## |
|
|
43 ## @item yd |
|
|
44 ## binary vector; @math{yd(ii)=0} indicates a continuous output; |
|
|
45 ## @math{yd(ii)=1} indicates a discrete output. |
|
|
46 ## @end table |
|
|
47 ## |
|
|
48 ## @strong{Outputs} @var{sys} |
|
|
49 ## @example |
|
|
50 ## @group |
|
|
51 ## sys.b := [sys.b , b] |
|
|
52 ## sys.c := [sys.c ] |
|
|
53 ## [ c ] |
|
|
54 ## sys.d := [sys.d | D12 ] |
|
|
55 ## [D21 | D22 ] |
|
|
56 ## @end group |
|
|
57 ## @end example |
|
3502
|
58 ## where @math{D12}, @math{D21}, and @math{D22} are the appropriate dimensioned |
|
3430
|
59 ## blocks of the input parameter @var{d}. |
|
|
60 ## @itemize @bullet |
|
3502
|
61 ## @item The leading block @math{D11} of @var{d} is ignored. |
|
3430
|
62 ## @item If @var{inname} and @var{outname} are not given as arguments, |
|
|
63 ## the new inputs and outputs are be assigned default names. |
|
|
64 ## @item @var{yd} is a binary vector of length rows(c) that indicates |
|
|
65 ## continuous/sampled outputs. Default value for @var{yd} is: |
|
|
66 ## |
|
|
67 ## @item @var{sys} = continuous or mixed |
|
|
68 ## @var{yd} = @code{zeros(1,rows(c))} |
|
|
69 ## |
|
|
70 ## @item @var{sys} = discrete |
|
|
71 ## @var{yd} = @code{ones(1,rows(c))} |
|
|
72 ## @end itemize |
|
|
73 ## @end deftypefn |
|
|
74 |
|
|
75 ## Author: John Ingram <ingraje@eng.auburn.edu> |
|
|
76 ## Created: August 1996 |
|
|
77 |
|
|
78 function retsys = sysappend (sys, b, c, d, outname, inname, yd) |
|
|
79 |
|
4460
|
80 save_warn_empty_list_elements = warn_empty_list_elements; |
|
|
81 unwind_protect |
|
|
82 warn_empty_list_elements = 0; |
|
3430
|
83 |
|
4460
|
84 ## check input arguments |
|
|
85 if ( (nargin < 2) | (nargin > 7) | (!isstruct(sys))) |
|
|
86 usage("retsys = sysappend(sys,b,c[,d,outname,inname,yd]) "); |
|
|
87 elseif(!isstruct(sys)) |
|
|
88 error("sys must be a system data structure"); |
|
|
89 endif |
|
3430
|
90 |
|
4460
|
91 ## default system type must be state space form |
|
|
92 [Aa,Ab,Ac,Ad,Ats,Ann,Anz,Ast,Ain,Aout,Ayd] = sys2ss(sys); |
|
|
93 [Ann,Anz,Am,Ap] = sysdimensions(sys); |
|
3430
|
94 |
|
4460
|
95 ## default c |
|
|
96 if(nargin < 3) c = []; endif |
|
3430
|
97 |
|
4460
|
98 ## default d |
|
|
99 if(nargin < 4) make_d = 1; |
|
|
100 elseif(isempty(d)) make_d = 1; |
|
|
101 else make_d = 0; endif |
|
|
102 if(make_d) d = zeros(rows(c)+Ap,columns(b) + Am); endif |
|
3430
|
103 |
|
4460
|
104 ## Append new input(s) if any |
|
|
105 Bm = max(columns(d),columns(b)+Am); |
|
|
106 if(Bm != Am) |
|
|
107 ## construct new signal names |
|
|
108 if(nargin >= 6) # new names were passed |
|
|
109 if(!isstr(inname)) |
|
|
110 error("inname must be a string"); |
|
|
111 elseif(rows(inname) != (Bm - Am)) |
|
|
112 error(sprintf("%d new inputs requested; inname(%dx%d)", ... |
|
|
113 (Bm-Am),rows(inname),columns(inname))); |
|
|
114 endif |
|
|
115 else |
|
|
116 inname = __sysdefioname__(Bm,"u",(Am+1)); |
|
|
117 endif |
|
4771
|
118 |
|
|
119 if(Am) |
|
|
120 Ain = __sysconcat__(Ain,inname); |
|
|
121 else |
|
|
122 Ain = inname; |
|
|
123 endif |
|
3430
|
124 |
|
4460
|
125 ## default b matrix |
|
|
126 if(isempty(b)) b = zeros(Ann+Anz,(Bm-Am)); |
|
|
127 elseif(rows(b) != Ann+Anz | columns(b) != (Bm-Am)) |
|
|
128 error(sprintf("b(%dx%d); should be (%dx%d)", rows(b), columns(b), ... |
|
|
129 (Ann+Anz), (Bm-Am))); |
|
3430
|
130 endif |
|
|
131 |
|
4460
|
132 ## append new b matrix |
|
|
133 Ab = [Ab,b]; |
|
3430
|
134 endif |
|
|
135 |
|
4460
|
136 ## Append new output(s) if any |
|
|
137 Bp = max(rows(d),rows(c)+Ap); |
|
|
138 if(Bp != Ap) |
|
3430
|
139 |
|
4460
|
140 ## construct new signal names, output classification |
|
|
141 if(nargin >= 5) # new names were passed |
|
|
142 if(!isstr(outname)) |
|
|
143 error("outname must be a string"); |
|
|
144 elseif(rows(outname) != (Bp - Ap)) |
|
|
145 error(sprintf("%d new outputs requested; outname(%dx%d)", ... |
|
|
146 (Bp-Ap),rows(outname),columns(outname))); |
|
|
147 endif |
|
|
148 else |
|
|
149 outname = __sysdefioname__(Bp,"y",(Ap+1)); |
|
3430
|
150 endif |
|
4771
|
151 if(Ap) Aout = __sysconcat__(Aout,outname); |
|
4460
|
152 else Aout = outname; endif |
|
3430
|
153 |
|
4460
|
154 ## construct new yd entries |
|
|
155 if(nargin == 7) |
|
|
156 if(!isvector(yd)) |
|
|
157 error(sprintf("yd(%dx%d) must be a vector",rows(yd),columns(yd))) |
|
|
158 elseif(rows(c) != length(yd) & rows(d) != length(yd)) |
|
|
159 error(sprintf("length(yd) = %d; c(%dx%d), d(%dx%d); mismatch", ... |
|
|
160 length(yd), rows(c), columns(c),rows(d),columns(d))); |
|
|
161 endif |
|
|
162 else |
|
|
163 ## default yd values |
|
|
164 yd = ones(1,Bp)*( (Ats > 0) & (Ann == 0) & isempty(find(Ayd == 0)) ) ; |
|
3430
|
165 endif |
|
4460
|
166 Ayd = [vec(Ayd);vec(yd)]; |
|
3430
|
167 |
|
4460
|
168 ## default c matrix |
|
|
169 if(isempty(c)) c = zeros((Bp-Ap),Ann+Anz); |
|
|
170 elseif(columns(c) != Ann+Anz | rows(c) != (Bp-Ap)) |
|
|
171 error(sprintf("c(%dx%d); should be (%dx%d)", rows(c), columns(c), ... |
|
|
172 (Bp-Ap), (Ann+Anz) )); |
|
|
173 endif |
|
|
174 |
|
|
175 ## append new c matrix |
|
|
176 Ac = [Ac;c]; |
|
3430
|
177 endif |
|
|
178 |
|
4460
|
179 ## check d matrix |
|
|
180 if(isempty(d)) d = zeros(Bp,Bm); |
|
|
181 elseif(rows(d) != Bp | columns(d) != Bm) |
|
|
182 error(sprintf("d(%dx%d) should be (%dx%d)",rows(d), columns(d), Bp, Bp)); |
|
|
183 endif |
|
3430
|
184 |
|
4460
|
185 ## Splice in original D matrix |
|
|
186 if(Am & Ap) d(1:Ap, 1:Am) = Ad; endif |
|
|
187 Ad = d; |
|
3430
|
188 |
|
4460
|
189 ## construct return system |
|
4771
|
190 retsys = ss(Aa,Ab,Ac,Ad,Ats,Ann,Anz,Ast,Ain,Aout,find(Ayd == 1)); |
|
3430
|
191 |
|
4460
|
192 unwind_protect_cleanup |
|
|
193 warn_empty_list_elements = save_warn_empty_list_elements; |
|
|
194 end_unwind_protect |
|
3430
|
195 |
|
|
196 endfunction |
