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
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.000
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
Unified Diff 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

#! /bin/sh

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

# Bourne shell script: dm2fr_m2dat

# Transformation: descriptor matrix to frequency response

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

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$



## 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





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




 echo $PARAMS





$MATRIX << EOF  | awk '$1 !~ /[A-Za-cf-z]/ && NF>0 {print} '  >$1_fr.dat




  $PARAMS;





  [A,B,C,D,E] = $1_dm($1_numpar);
  fr = dm2fr(A,B,C,D,E,w);














  rfr = real(fr);
  ifr = imag(fr);






  [Ny,Nu]=size(D);

  NN = Ny*Nu;

  xy = [rfr(:,1), ifr(:,1)];
  for i=2:NN
    xy = [xy rfr(:,i) ifr(:,i)];
  end;
xy


EOF








|










>
>
>








|
>
>
>
>



|





>
>
>
|

>
>
>
>
|
>
>
>

<
>
>
>
>


|
>
>
>

>
>
>
>
>
>
>
>
>
>
|
|
>
>
>
>
>

<
>
|
|
|
|
|

<
>




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_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.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); 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;
  $PARAMS

  %Defaults
  if exist('w')==0
    w = logspace(-1,2,100);
  end;

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


  [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,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
  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');


% Nichols style
  re = lw;
  im = pfr;
  nic = [im(:,1) re(:,1)];
  for i = 2:ny
    nic = [nic im(:,i) re(:,i)]
  end;

  write_matrix(nic, '$1_nifr');

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