[14f77c]: inst / @lti / __lti_group__.m  Maximize  Restore  History

Download this file

81 lines (63 with data), 2.4 kB

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
## Copyright (C) 2009 Lukas F. Reichlin
##
## This file is part of LTI Syncope.
##
## LTI Syncope is free software: you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation, either version 3 of the License, or
## (at your option) any later version.
##
## LTI Syncope is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with LTI Syncope. If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## Block diagonal concatenation of two LTI models.
## This file is part of the Model Abstraction Layer.
## For internal use only.
## Author: Lukas Reichlin <lukas.reichlin@gmail.com>
## Created: September 2009
## Version: 0.1
function retlti = __lti_group__ (lti1, lti2)
retlti = lti ();
retlti.inname = [lti1.inname;
lti2.inname];
retlti.outname = [lti1.outname;
lti2.outname];
if (nfields (lti1.ingroup) || nfields (lti2.ingroup))
m1 = numel (lti1.inname);
lti2.ingroup = structfun (@(x) x + m1, lti2.ingroup, "uniformoutput", false);
retlti.ingroup = __merge_struct__ (lti1.ingroup, lti2.ingroup);
endif
if (nfields (lti1.outgroup) || nfields (lti2.outgroup))
p1 = numel (lti1.outname);
lti2.outgroup = structfun (@(x) x + p1, lti2.outgroup, "uniformoutput", false);
retlti.outgroup = __merge_struct__ (lti1.outgroup, lti2.outgroup);
endif
if (lti1.tsam == lti2.tsam)
retlti.tsam = lti1.tsam;
elseif (lti1.tsam == -2)
retlti.tsam = lti2.tsam;
elseif (lti2.tsam == -2)
retlti.tsam = lti1.tsam;
elseif (lti1.tsam == -1 && lti2.tsam > 0)
retlti.tsam = lti2.tsam;
elseif (lti2.tsam == -1 && lti1.tsam > 0)
retlti.tsam = lti1.tsam;
else
error ("lti_group: systems must have identical sampling times");
endif
endfunction
function ret = __merge_struct__ (a, b)
## FIXME: this is too complicated;
## isn't there a simple function for this task?
fa = fieldnames (a);
fb = fieldnames (b);
ca = struct2cell (a);
cb = struct2cell (b);
ret = cell2struct ([ca; cb], [fa; fb]);
## FIXME: name clashes
endfunction