11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
# Copyright (c) P.J.Gawthrop, 1996.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.10 1998/07/27 20:27:25 peterg
## Now get parameters for simpar file.
##
## Revision 1.9 1998/02/26 15:12:26 peterg
## Removed parameter stuff
##
# Revision 1.8 1996/08/30 14:54:36 peter
|
>
>
>
|
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 1998/08/11 13:32:42 peterg
## Lowercase mttLAST etc
##
## Revision 1.10 1998/07/27 20:27:25 peterg
## Now get parameters for simpar file.
##
## Revision 1.9 1998/02/26 15:12:26 peterg
## Removed parameter stuff
##
# Revision 1.8 1996/08/30 14:54:36 peter
|
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
|
if exist('mttwsteps')==0
mttwsteps = 100;
end;
W = logspace(mttwmin,mttwmax,mttwsteps)';
if exist('u0')==0
u0 = zeros(nu,1);
u0(1) = 1;
end;
[n,m]=size(W);
if m>n
W=W';
end;
[n,m]=size(u0);
if m>n
u0=u0';
end;
[A,B,C,D,E] = $1_dm
fr = dm2fr(A,B,C,D,E,W,u0)
lw = log10(W);
lmfr = log10(abs(fr));
pfr = angle(fr)*180/pi;
% Complex frequency response
write_matrix([W fr], '$1_fr');
|
|
<
|
>
>
|
|
|
|
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
|
if exist('mttwsteps')==0
mttwsteps = 100;
end;
W = logspace(mttwmin,mttwmax,mttwsteps)';
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
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');
|