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
|
function structure = TF_eqn(name,bond_number,bonds,direction,cr,args, ...
structure,eqnfile);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.2 1996/09/10 10:41:37 peter
% %% Now used equation.m to write the equations.
% %%
% %% Revision 1.1 1996/08/19 09:05:04 peter
% %% Initial revision
% %%
if check_bonds(bonds,2,'TF')==0
return
end
% There are 2 ports; extract the information
e_1 = bonds(1,1);
f_1 = bonds(1,2);
e_2 = bonds(2,1);
f_2 = bonds(2,2);
% Effort
outcause = 1;
incause = 1;
if e_1==1
outbond = bond_number(2);
inbond = bond_number(1);
outport = 2;
inport = 1;
else
outbond = bond_number(1);
inbond = bond_number(2);
outport = 1;
inport = 2;
end;
eqn = equation(name,cr,args,outbond,outcause,outport, ...
inbond,incause,inport);
fprintf(eqnfile, '%s',eqn);
% Flow
outcause = -1;
incause = -1;
if f_1==-1
outbond = bond_number(2);
inbond = bond_number(1);
outport = 2;
inport = 1;
else
outbond = bond_number(1);
inbond = bond_number(2);
outport = 1;
inport = 2;
end;
eqn = equation(name,cr,args,outbond,outcause,outport, ...
inbond,incause,inport);
fprintf(eqnfile, '%s',eqn);
|
|
>
<
<
<
<
<
<
>
|
|
>
|
|
|
|
|
|
|
|
|
|
>
>
>
|
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
|
function structure = EMTF_eqn(name,bond_number,bonds,direction,cr,args, ...
structure,eqnfile);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if nargin<8
eqnfile = 'stdout';
end;
% Check that there are exactly two bonds.
if check_bonds(bonds,3,'EMTF')==0
return
end
% There are 3 ports; extract the information
e_1 = bonds(1,1);
f_1 = bonds(1,2);
e_2 = bonds(2,1);
f_2 = bonds(2,2);
e_3 = bonds(3,1);
% Effort
outcause = 1;
incause = [1;1];
if e_1==1
outbond = bond_number(2);
inbond = bond_number([1 3]);
outport = 2;
inport = [1;3];
else
outbond = bond_number(1);
inbond = bond_number([2 3]);
outport = 1;
inport = [2;3];
end;
eqn = equation(name,cr,args,outbond,outcause,outport, ...
inbond,incause,inport);
fprintf(eqnfile, '%s',eqn);
% Flow
outcause = -1;
incause = [-1;1];
if f_1==-1
outbond = bond_number(2);
inbond = bond_number([1 3]);
outport = 2;
inport = [1;3];
else
outbond = bond_number(1);
inbond = bond_number([2 3]);
outport = 1;
inport = [2;3];
end;
eqn = equation(name,cr,args,outbond,outcause,outport, ...
inbond,incause,inport);
fprintf(eqnfile, '%s',eqn);
% Modulation: flow on port 3 is always zero
fprintf(eqnfile, '%s := 0;\n', ...
varname(name,bond_number(3), -1));
|