Differences From Artifact [72a8466d78]:

To Artifact [c4b4dcedcc]:


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



13
14
15
16
17
18
19
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.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
% %% Now has u0 input vector
% %%
% %% Revision 1.2  1996/08/15 10:24:28  peter
50
51
52
53
54
55
56
57

58
59
60
61
62
63
64
53
54
55
56
57
58
59

60
61
62
63
64
65
66
67







-
+







end;


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






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