Mercurial > hg > octave-nkf

comparison scripts/control/base/lsim.m @ 4174:aabb8307021f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2002-11-13 00:45:22 by jwe]
author jwe
date Wed, 13 Nov 2002 00:45:23 +0000
parents 0d411821682c
children bdbee5282954
comparison
equal deleted inserted replaced
4173:2aea727f516f 4174:aabb8307021f
71 71
72 Ts = 0; 72 Ts = 0;
73 t(2)-t(1); 73 t(2)-t(1);
74 u=u'; 74 u=u';
75 n = max(size(t)); 75 n = max(size(t));
76
76 for ii = 1:(n-1) 77 for ii = 1:(n-1)
77 78
78 ## check if step size changed 79 ## check if step size changed
79 if (t(ii+1) - t(ii) != Ts) 80 ## XXX FIXME XXX -- this is probably not the best test, but it is
81 ## better than a test for exact equality.
82 if (abs (t(ii+1) - t(ii) - Ts) > 10 * eps)
80 Ts = t(ii+1) - t(ii); 83 Ts = t(ii+1) - t(ii);
81 ## [F,G] = c2d(a,b,Ts); 84 ## [F,G] = c2d(a,b,Ts);
82 dsys = c2d(sys, Ts); 85 dsys = c2d(sys, Ts);
83 [F,G] = sys2ss(dsys); 86 [F,G] = sys2ss(dsys);
84 endif 87 endif