Differences From Artifact [aa5007007b]:

To Artifact [6f8f6f0736]:


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  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
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('u0')==0
  if exist('mttinput')==0
    u0 = zeros(nu,1);
    u0(1) = 1;
    mttinput = 1;
  end;
  mttu0 = zeros(nu,1);
  mttu0(mttinput) = 1;


  [n,m]=size(W);
  if m>n
    W=W';
  end;

  [n,m]=size(u0);
  [n,m]=size(mttu0);
  if m>n
    u0=u0';
    mttu0=mttu0';
  end;

  [A,B,C,D,E] = $1_dm
  fr = dm2fr(A,B,C,D,E,W,u0)
  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');


MTT: Model Transformation Tools
GitHub | SourceHut | Sourceforge | Fossil RSS ]