11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
+
+
+
|
# Copyright (c) P.J.Gawthrop, 1996.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.11 1997/08/29 07:56:54 peterg
## Minor updates
##
# Revision 1.10 1997/01/07 09:16:03 peterg
# Added step_factor parameter - gives that number of integration steps
# per sample.
#
## Revision 1.9 1997/01/06 21:36:44 peterg
## Fixed bug mtt_error --> mtt_error.txt
## Replaced lsode by Euler integration.
|
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
81
82
83
84
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
|
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
81
82
83
84
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
|
+
+
+
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
-
-
+
+
+
+
+
+
+
+
-
-
+
-
-
-
+
+
+
+
-
+
-
+
-
+
|
##
## Revision 1.1 1996/08/15 11:56:38 peter
## Initial revision
##
###############################################################
echo Creating $1_odes.m
echo Creating $1_odeso.m
rm -f ode2odes_m.log
rm -f mtt_error.txt
if [ "$2" = "" ];
then
PARAMS='T=[0:0.1:10]; x0=zeros(nx,1);'
echo Using default parameter $PARAMS
echo $PARAMS>$1_args.m
else
PARAMS=$2;
fi
#if [ "$2" = "" ];
#then
# PARAMS='T=[0:0.1:10]; x0=zeros(nx,1);'
# echo Using default parameter $PARAMS
# echo $PARAMS>$1_args.m
#else
# PARAMS=$2;
#fi
# PARAMS="$PARAMS ;"
$MATRIX << EOF > ode2odes_m.log 2>mtt_error.txt
step_factor = 1; %Set up default step factor
$1_numpar;
$1_numpar;
%Read in state
x = $1_state
%Read in simulation parameters
$1_simpar
T = [0:DT:LAST]
[nx,ny,nu,nz,nyz] = $1_def;
t=0; %Just in case it appears in the parameter list.
$PARAMS
T=[0:0.1:10]
T=[0:1:100]
end;
if exist('x0')==0
x0 = zeros(nx,1);
end;
% if exist('x0')==0
% % x0 = zeros(nx,1);
% x0 = x;
% end;
[n,m]=size(T);
if m>n
T=T';
end;
if nx>0
% x = lsode('$1_ode', x0, T);
%Euler integration
x = x0;
% x = x0;
X=[];
dt = (T(2)-T(1))/step_factor;
dt = (T(2)-T(1))/STEPFACTOR;
for t=T'
X = [X x];
ts = t;
for i=1:step_factor
for i=1:STEPFACTOR
dx = $1_ode(x,ts);
ts = ts + dt;
x = x + dx*dt;
end;
end;
X = X';
|