Overview
Comment: | Corrected mistake in step-response calculation. |
---|---|
Downloads: | Tarball | ZIP archive | SQL archive |
Timelines: | family | ancestors | descendants | both | origin/master | trunk |
Files: | files | file ages | folders |
SHA3-256: |
88386835e829bc42cedab089d0068e18 |
User & Date: | gawthrop@users.sourceforge.net on 1996-08-11 10:37:40 |
Other Links: | branch diff | manifest | tags |
Context
1996-08-11
| ||
10:38:30 | Minor debugging changes. check-in: 5495cab3b8 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
10:37:40 | Corrected mistake in step-response calculation. check-in: 88386835e8 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
10:06:50 | Now allows e format and NaN numbers. check-in: fa5c69c7a5 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Modified mttroot/mtt/bin/trans/m/dm2sr.m from [dcd7ecbaa1] to [27a1ed2c88].
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 | function sr = dm2sr(A,B,C,D,E,T); % sr = dm2sr(A,B,C,D,E,T); % Descriptor matrix to impulse response. % NB At the moment - this assumes that E is unity ..... % A,B,C,D,E - descriptor matrices % T vector of time points % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [N,M] = size(T); if M>N T = T'; N = M; end; [Ny,Nu] = size(D); NN=Ny*Nu; [N_y,N_x] = size(C); one = eye(N_x); sr = zeros(N,NN); i = 0; for t = T' i=i+1; | > > > > | | 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 | function sr = dm2sr(A,B,C,D,E,T); % sr = dm2sr(A,B,C,D,E,T); % Descriptor matrix to impulse response. % NB At the moment - this assumes that E is unity ..... % A,B,C,D,E - descriptor matrices % T vector of time points % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ % %% Revision 1.1 1996/08/11 09:42:40 peter % %% Initial revision % %% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [N,M] = size(T); if M>N T = T'; N = M; end; [Ny,Nu] = size(D); NN=Ny*Nu; [N_y,N_x] = size(C); one = eye(N_x); sr = zeros(N,NN); i = 0; for t = T' i=i+1; SR = C*( A\(exp(A*t)-1) )*B + D*ones(size(t)); sr(i,:) = reshape(SR, 1,NN); end; |