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.13 1998/09/03 09:12:32 peterg
## Removed debugging lines
##
## Revision 1.12 1998/09/02 12:02:14 peterg
## Added mttinput paramater
##
## Revision 1.11 1998/08/11 13:32:42 peterg
## Lowercase mttLAST etc
##
## Revision 1.10 1998/07/27 20:27:25 peterg
|
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
123
124
125
|
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
123
124
125
126
127
128
|
-
+
-
+
-
-
+
+
-
-
+
+
-
-
+
+
-
+
-
+
|
#else
# PARAMS=$2;
#fi
$MATRIX << EOF > dm2fr_m.log 2>mtt_error.txt
# Data
$1_numpar;
mttpar=$1_numpar;
# Simulation parameters
$1_simpar;
simpar=$1_simpar;
[nx,ny,nu,nz,nyz] = $1_def
%Defaults
if exist('mttwmin')==0
mttwmin = -1;
if exist('simpar.wmin')==0
simpar.wmin = -1;
end;
if exist('mttwmax')==0
mttwmax = 2;
if exist('simpar.wmax')==0
simpar.wmax = 2;
end;
if exist('mttwsteps')==0
mttwsteps = 100;
if exist('simpar.wsteps')==0
simpar.wsteps = 100;
end;
W = logspace(mttwmin,mttwmax,mttwsteps)';
W = logspace(simpar.wmin,simpar.wmax,simpar.wsteps)';
if exist('mttinput')==0
mttinput = 1;
end;
mttu0 = zeros(nu,1);
mttu0(mttinput) = 1;
[n,m]=size(W);
if m>n
W=W';
end;
[n,m]=size(mttu0);
if m>n
mttu0=mttu0';
end;
[A,B,C,D,E] = $1_dm;
[A,B,C,D,E] = $1_dm(mttpar);
fr = dm2fr(A,B,C,D,E,W,mttu0);
lw = log10(W);
lmfr = log10(abs(fr));
pfr = angle(fr)*180/pi;
% Complex frequency response
write_matrix([W fr], '$1_fr');
|