1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
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
|
-
+
+
+
+
|
#! /bin/sh
######################################
##### Model Transformation Tools #####
######################################
# Bourne shell script: ode_r2m
# Reduce ODE to simulab ODE
# Reduce ODE to matlab ODE
# P.J.Gawthrop 14 June 1991, 12 Jan 1994, April 1994, Jan 95.
# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
# Revision 1.6 1997/01/05 19:34:35 peterg
# Don't write globals which are already assigned to a number.
#
## Revision 1.5 1996/09/13 19:41:39 peter
## *** empty log message ***
##
## Revision 1.4 1996/09/12 18:33:32 peter
## Put back under rcs
##
## Revision 1.3 1996/08/30 11:04:28 peter
|
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
|
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
-
+
-
+
-
-
+
+
-
-
+
+
|
write "% Read in the arguments";
write "$1_args";
write "% Set up the State variables";
FOR i := 1:MTTNx DO
BEGIN
write "MTTx", i, " = x(", i, ");";
write "mttx", i, " = x(", i, ");";
END;
write "% Set up the Input variables";
IF MTTNu>0 THEN
FOR i := 1:MTTNu DO
BEGIN
write "MTTu", i, " = u(", i, ");";
write "mttu", i, " = u(", i, ");";
END;
END;
% Firstly do the dx = f(x,t) function.
OUT "$1_ode.m";
write "function MTTdX = $1_ode(x,t);";
write "% dX = $1_ode(x,t);";
write "function mttdx = $1_ode(x,t);";
write "% mttdx = $1_ode(x,t);";
write "%ODE in Simulab form for system $1;;";
write "%File $1_ode.m;;";
write "%Generated by MTT;;";
common();
%Fortran switches - one line expressions
OFF echo;
ON fort$
cardno!* := 1$
fortwidth!* := 10000$
OFF period$
MTTdx := MTTdx;
SHUT "$1_ode.m";
OFF fort;
% Now do the y = g(x,t) function.
OUT "$1_odeo.m";
write "function MTTy = $1_odeo(x,t);";
write "% dX = $1_odeo(x,t);";
write "function mtty = $1_odeo(x,t);";
write "% mtty = $1_odeo(x,t);";
write "%ODE in Simulab form for system $1;;";
write "%File $1_odeo.m;;";
write "%Generated by MTT;;";
common();
%Fortran switches - one line expressions
|