Overview
Comment:Replaced 'C*( (E*j*w - A) \ B*u0 ' by C*( inv(E*j*w - A)*B*u0 ). \
seems to have a bug for complex nos (??)
Downloads: Tarball | ZIP archive | SQL archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: 1ac694cfc981da62b1f17f2d3243034f037895c082e22919ec86de09f4ece300
User & Date: gawthrop@users.sourceforge.net on 1997-06-13 13:05:38
Other Links: branch diff | manifest | tags
Context
1997-06-13
14:08:23
Changed default rep.txt check-in: a735ede176 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
13:05:38
Replaced 'C*( (E*j*w - A) \ B*u0 ' by C*( inv(E*j*w - A)*B*u0 ). \
seems to have a bug for complex nos (??) check-in: 1ac694cfc9 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
12:52:18
Removed extra ; from $PARAMS check-in: 9f6c3bbb2f user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes
Hide Diffs Unified Diffs Ignore Whitespace Patch

Modified mttroot/mtt/bin/trans/m/dm2fr.m from [3b12a6a127] to [e3aede7152]. 

1
2
3
4
5
6
7
8
9
10
11
12



13
14
15

16
17
18
19
20
21
22

function  fr = dm2fr(A,B,C,D,E,W,u0)
% fr = dm2fr(A,B,C,D,E,W,u0)
% Descriptor matrix to frequency response.
% A,B,C,D,E - descriptor matrices
% W vector of frequency points
% u0 input gain vector: u = u0*unit phasor

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$



// Revision 1.6  1996/08/24  14:22:23  peter
// Put in a ; to avoid excessive log output.
//

% %% Revision 1.5  1996/08/16 14:26:37  peter
% %% Check and fix size of u0.
% %%
% %% Revision 1.4  1996/08/15 12:50:51  peter
% %% Put in a conj to undo effect of transpose.
% %%
% %% Revision 1.3  1996/08/15 11:53:44  peter













>
>
>
|
|
<
>








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

function  fr = dm2fr(A,B,C,D,E,W,u0)
% fr = dm2fr(A,B,C,D,E,W,u0)
% Descriptor matrix to frequency response.
% A,B,C,D,E - descriptor matrices
% W vector of frequency points
% u0 input gain vector: u = u0*unit phasor

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.7  1996/11/06  16:40:38  peterg
% %% Explicit definition of j
% %%
% %% Revision 1.6  1996/08/24  14:22:23  peter
% %% Put in a ; to avoid excessive log output.

% %%
% %% Revision 1.5  1996/08/16 14:26:37  peter
% %% Check and fix size of u0.
% %%
% %% Revision 1.4  1996/08/15 12:50:51  peter
% %% Put in a conj to undo effect of transpose.
% %%
% %% Revision 1.3  1996/08/15 11:53:44  peter
︙︙
56
57
58
59
60
61
62
63

64
65
66
67
68
69
70

end;

j = sqrt(-1);
fr = zeros(N,Ny);
i = 0;
for w = W'
  i = i+1;
  FR = C*( (E*j*w - A) \ B*u0 ) + D*u0;

  fr(i,:) = conj(FR');
end;













|
>








59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74

end;

j = sqrt(-1);
fr = zeros(N,Ny);
i = 0;
for w = W'
  i = i+1;
%%  FR = C*( (E*j*w - A) \ B*u0 ) + D*u0;
  FR = C*( inv(E*j*w - A)*B*u0 ) + D*u0;
  fr(i,:) = conj(FR');
end;
︙︙

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