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
|
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
|
-
-
+
+
-
+
-
+
-
+
-
+
-
+
-
+
-
+
-
+
-
-
+
+
-
-
+
+
|
LET slin(gain_causality, gain, other_causality, 1, input, causality, 1)
= lin(gain_causality, gain, other_causality, 1, input,
causality, 1);
% Two port version for sC,sI,sR
% Port 1 - standard
FOR ALL gain_causality, gain, sgain, causality, input, sinput, other_causality
LET slin(gain_causality, gain, sgain, other_causality, 1,
FOR ALL COMPONENT, gain_causality, gain, sgain, causality, input, sinput, other_causality
LET slin(COMPONENT,gain_causality, gain, sgain, other_causality, 1,
input, causality, 1,
sinput, causality, 2)
= lin(gain_causality, gain, other_causality, 1,
= lin(COMPONENT,gain_causality, gain, other_causality, 1,
input, causality, 1);
% Port 2 - sensitivity
FOR ALL gain_causality, gain, sgain, causality, input, sinput, other_causality
FOR ALL COMPONENT, gain_causality, gain, sgain, causality, input, sinput, other_causality
SUCH THAT sgain NEQ 0
LET slin(gain_causality, gain, sgain, other_causality, 2,
LET slin(COMPONENT,gain_causality, gain, sgain, other_causality, 2,
input, causality, 1,
sinput, causality, 2)
= lin(gain_causality, gain, other_causality, 1,
= lin(COMPONENT,gain_causality, gain, other_causality, 1,
sinput, causality, 1)
+ sgain*df(lin(gain_causality, gain, other_causality, 1,
+ sgain*df(lin(COMPONENT,gain_causality, gain, other_causality, 1,
input, causality, 1), gain);
%% Don't compute the derivative in this case
FOR ALL gain_causality, gain, sgain, causality, input, sinput, other_causality
FOR ALL COMPONENT, gain_causality, gain, sgain, causality, input, sinput, other_causality
SUCH THAT sgain EQ 0
LET slin(gain_causality, gain, sgain, other_causality, 2,
LET slin(COMPONENT,gain_causality, gain, sgain, other_causality, 2,
input, causality, 1,
sinput, causality, 2)
= lin(gain_causality, gain, other_causality, 1,
= lin(COMPONENT,gain_causality, gain, other_causality, 1,
sinput, causality, 1);
%% This is the version to go with sEMTF
%% It is the CR passed to the AE3 components
%DESCRIPTION four port component - effort i/o
FOR ALL gain, input, junk, m_input, dm_input
LET slin(gain, effort, 2,
FOR ALL COMPONENT, gain, input, junk, m_input, dm_input
LET slin(COMPONENT,gain, effort, 2,
input, effort, 1,
junk, flow, 2,
m_input, effort, 3,
dm_input,effort, 4)
= gain*dm_input*input;
%DESCRIPTION four port component - flow i/o
FOR ALL gain, input, junk, m_input, dm_input
LET slin(gain, flow, 2,
FOR ALL COMPONENT, gain, input, junk, m_input, dm_input
LET slin(COMPONENT,gain, flow, 2,
input, flow, 1,
junk, effort, 2,
m_input, effort, 3,
dm_input,effort, 4)
= gain*dm_input*input;
END;;
|