Overview
Comment:Now creats a number of versions:
fr complex frequency response
lmfr log magnitude of fr
lpfr phase of fr
nyfr real and imag parts - Nyquist style
nifr Nichols style.
Downloads: Tarball | ZIP archive | SQL archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: d6652c7f6890d59eb4a4e883ffa173139ad248131dfba5b2aa911e7f81d89aee
User & Date: gawthrop@users.sourceforge.net on 1996-08-15 11:52:42
Other Links: branch diff | manifest | tags
Context
1996-08-15
11:53:44
Now has u0 input vector check-in: 30eb07fb39 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
11:52:42
Now creats a number of versions:
fr complex frequency response
lmfr log magnitude of fr
lpfr phase of fr
nyfr real and imag parts - Nyquist style
nifr Nichols style. check-in: d6652c7f68 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
10:24:28
Includes u0 argument. check-in: e40d837faa user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes
Hide Diffs Side-by-Side Diffs Ignore Whitespace Patch

Modified mttroot/mtt/bin/trans/dm2fr_m from [d561b3e13c] to [1f2439c8d4]. 

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





27
28
29
30

31
32
33
34
35



36

37




38




39
40




41
42
43




44










45
46







47
48
49
50
51
52
53






54
55

56
57
58

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
27
28

29
30
31
32
33
34
35
36

37
38
39
40
41
42
43
44
45

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







-
+










+
+
+








-
+
+
+
+
+



-
+





+
+
+
-
+

+
+
+
+
-
+
+
+
+

-
+
+
+
+


-
+
+
+
+

+
+
+
+
+
+
+
+
+
+
-
-
+
+
+
+
+
+
+

-
-
-
-
-
-
+
+
+
+
+
+

-
+




#! /bin/sh

     ###################################### 
     ##### Model Transformation Tools #####
     ######################################

# Bourne shell script: dm2fr_m2dat
# Bourne shell script: dm2fr_m

# Transformation: descriptor matrix to frequency response

# Copyright (c) P.J.Gawthrop, 1996.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.3  1996/08/11 19:08:08  peter
## Parameter passing now enabled.
##
## Revision 1.2  1996/08/11 09:32:12  peter
## Now takes the numpar parameters correctly
##
## Revision 1.1  1996/08/10 14:11:11  peter
## Initial revision
##
###############################################################

echo Creating $1_fr.dat
echo Creating $1_fr.m
echo Creating $1_lmfr.m
echo Creating $1_lpfr.m
echo Creating $1_nyfr.m
echo Creating $1_nifr.m

if [ "$2" = "" ]; 
then
  PARAMS='w=logspace(-1,2,100);'
  PARAMS='w=logspace(-1,2,100); u0=ones(nu,1);'
  echo Using default parameter $PARAMS
else
  PARAMS=$2;
fi

$MATRIX << EOF > dm2fr_m.log

  [nx,ny,nu,nz,nyz] = $1_def;
 echo $PARAMS
  $PARAMS

  %Defaults
  if exist('w')==0
    w = logspace(-1,2,100);
  end;
$MATRIX << EOF  | awk '$1 !~ /[A-Za-cf-z]/ && NF>0 {print} '  >$1_fr.dat

  if exist('u0')==0
    u0 = ones(nu,1);
  end;

  $PARAMS;
  [n,m]=size(w);
  if m>n
    w=w';
  end;

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

% Log magnitude v log frequency
  write_matrix([lw lmfr], '$1_lmfr');

% Angle v log frequency
  write_matrix([lw pfr], '$1_lpfr');

% Nyquist style
  rfr = real(fr);
  ifr = imag(fr);
  re = real(fr);
  im = imag(fr);
  nyq = [re(:,1) im(:,1)];
  for i = 2:ny
    nyq = [nyq re(:,i) im(:,i)]
  end;
  write_matrix(nyq, '$1_nyfr');

  [Ny,Nu]=size(D);
  NN = Ny*Nu;

  xy = [rfr(:,1), ifr(:,1)];
  for i=2:NN
    xy = [xy rfr(:,i) ifr(:,i)];
% Nichols style
  re = lw;
  im = pfr;
  nic = [im(:,1) re(:,1)];
  for i = 2:ny
    nic = [nic im(:,i) re(:,i)]
  end;
xy
  write_matrix(nic, '$1_nifr');

EOF

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