|
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 |
|
7016
|
6 ## under the terms of the GNU General Public License as published by |
|
|
7 ## the Free Software Foundation; either version 3 of the License, or (at |
|
|
8 ## your option) any later version. |
|
3430
|
9 ## |
|
7016
|
10 ## Octave is distributed in the hope that it will be useful, but |
|
|
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
|
|
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
|
13 ## General Public License for more details. |
|
3430
|
14 ## |
|
|
15 ## You should have received a copy of the GNU General Public License |
|
7016
|
16 ## along with Octave; see the file COPYING. If not, see |
|
|
17 ## <http://www.gnu.org/licenses/>. |
|
3430
|
18 |
|
|
19 ## -*- texinfo -*- |
|
5016
|
20 ## @deftypefn {Function File} {@var{sys} =} sysappend (@var{syst}, @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 |
|
5016
|
25 ## @item syst |
|
3430
|
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 ## |
|
5016
|
48 ## @strong{Outputs} |
|
|
49 ## @table @var |
|
|
50 ## @item sys |
|
3430
|
51 ## @example |
|
|
52 ## @group |
|
5016
|
53 ## sys.b := [syst.b , b] |
|
|
54 ## sys.c := [syst.c ] |
|
3430
|
55 ## [ c ] |
|
5016
|
56 ## sys.d := [syst.d | D12 ] |
|
|
57 ## [ D21 | D22 ] |
|
3430
|
58 ## @end group |
|
|
59 ## @end example |
|
3502
|
60 ## where @math{D12}, @math{D21}, and @math{D22} are the appropriate dimensioned |
|
3430
|
61 ## blocks of the input parameter @var{d}. |
|
|
62 ## @itemize @bullet |
|
3502
|
63 ## @item The leading block @math{D11} of @var{d} is ignored. |
|
3430
|
64 ## @item If @var{inname} and @var{outname} are not given as arguments, |
|
|
65 ## the new inputs and outputs are be assigned default names. |
|
|
66 ## @item @var{yd} is a binary vector of length rows(c) that indicates |
|
|
67 ## continuous/sampled outputs. Default value for @var{yd} is: |
|
5016
|
68 ## @itemize @minus |
|
|
69 ## @item @var{sys} is continuous or mixed |
|
3430
|
70 ## @var{yd} = @code{zeros(1,rows(c))} |
|
|
71 ## |
|
5016
|
72 ## @item @var{sys} is discrete |
|
3430
|
73 ## @var{yd} = @code{ones(1,rows(c))} |
|
|
74 ## @end itemize |
|
5016
|
75 ## @end itemize |
|
|
76 ## @end table |
|
3430
|
77 ## @end deftypefn |
|
|
78 |
|
|
79 ## Author: John Ingram <ingraje@eng.auburn.edu> |
|
|
80 ## Created: August 1996 |
|
|
81 |
|
|
82 function retsys = sysappend (sys, b, c, d, outname, inname, yd) |
|
|
83 |
|
5568
|
84 ## check input arguments |
|
|
85 if ( (nargin < 2) | (nargin > 7) | (!isstruct(sys))) |
|
6046
|
86 print_usage (); |
|
5568
|
87 elseif(!isstruct(sys)) |
|
|
88 error("sys must be a system data structure"); |
|
|
89 endif |
|
3430
|
90 |
|
5568
|
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 |
|
5568
|
95 ## default c |
|
|
96 if(nargin < 3) c = []; endif |
|
3430
|
97 |
|
5568
|
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 |
|
5568
|
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 |
|
5443
|
109 if(!ischar(inname)) |
|
4460
|
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 |
|
5568
|
115 else |
|
4460
|
116 inname = __sysdefioname__(Bm,"u",(Am+1)); |
|
5568
|
117 endif |
|
4771
|
118 |
|
5568
|
119 if(Am) |
|
|
120 Ain = __sysconcat__(Ain,inname); |
|
|
121 else |
|
|
122 Ain = inname; |
|
|
123 endif |
|
3430
|
124 |
|
5568
|
125 ## default b matrix |
|
|
126 if(isempty(b)) b = zeros(Ann+Anz,(Bm-Am)); |
|
|
127 elseif(rows(b) != Ann+Anz | columns(b) != (Bm-Am)) |
|
4460
|
128 error(sprintf("b(%dx%d); should be (%dx%d)", rows(b), columns(b), ... |
|
|
129 (Ann+Anz), (Bm-Am))); |
|
3430
|
130 endif |
|
|
131 |
|
5568
|
132 ## append new b matrix |
|
|
133 Ab = [Ab,b]; |
|
|
134 endif |
|
3430
|
135 |
|
5568
|
136 ## Append new output(s) if any |
|
|
137 Bp = max(rows(d),rows(c)+Ap); |
|
|
138 if(Bp != Ap) |
|
|
139 |
|
|
140 ## construct new signal names, output classification |
|
|
141 if(nargin >= 5) # new names were passed |
|
5443
|
142 if(!ischar(outname)) |
|
4460
|
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 |
|
5568
|
148 else |
|
4460
|
149 outname = __sysdefioname__(Bp,"y",(Ap+1)); |
|
5568
|
150 endif |
|
|
151 if(Ap) Aout = __sysconcat__(Aout,outname); |
|
|
152 else Aout = outname; endif |
|
3430
|
153 |
|
5568
|
154 ## construct new yd entries |
|
|
155 if(nargin == 7) |
|
4460
|
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 |
|
5568
|
162 else |
|
4460
|
163 ## default yd values |
|
|
164 yd = ones(1,Bp)*( (Ats > 0) & (Ann == 0) & isempty(find(Ayd == 0)) ) ; |
|
5568
|
165 endif |
|
|
166 Ayd = [vec(Ayd);vec(yd)]; |
|
3430
|
167 |
|
5568
|
168 ## default c matrix |
|
|
169 if(isempty(c)) c = zeros((Bp-Ap),Ann+Anz); |
|
|
170 elseif(columns(c) != Ann+Anz | rows(c) != (Bp-Ap)) |
|
4460
|
171 error(sprintf("c(%dx%d); should be (%dx%d)", rows(c), columns(c), ... |
|
|
172 (Bp-Ap), (Ann+Anz) )); |
|
3430
|
173 endif |
|
|
174 |
|
5568
|
175 ## append new c matrix |
|
|
176 Ac = [Ac;c]; |
|
|
177 endif |
|
3430
|
178 |
|
5568
|
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 |
|
5568
|
185 ## Splice in original D matrix |
|
|
186 if(Am & Ap) d(1:Ap, 1:Am) = Ad; endif |
|
|
187 Ad = d; |
|
3430
|
188 |
|
5568
|
189 ## construct return system |
|
|
190 retsys = ss(Aa,Ab,Ac,Ad,Ats,Ann,Anz,Ast,Ain,Aout,find(Ayd == 1)); |
|
3430
|
191 |
|
|
192 endfunction |
