You made a mistake because the response you compute with
rep=flts(u,tf2ss(h)); k=find(rep<>0,1)
is not exactly a step response because the step does not occur at the
first sample
If you use
H=time_id(10,u,rep)
to identify the response you vill got a batter result
An other an probably better solution is to extract the delay before
starting the identification:
k=find(rep<>0,1) //here the threshold has to be improved in case of
noisy signal
H=time_id(1,"step",rep(k:$))
Serge
Le 11/03/2016 18:44, noguchi a écrit :
Dear Mr.Serge,
Thank you for your advise. I tried to use time_ID.
Please look at the following Scinote:
z=poly(0,'z');
h=(0.07/(z-0.934))*(1/z^10)
u=zeros(1,100);
for i=6:1:100
u(1,i)=1;
end
t=1:1:100;
rep=flts(u,tf2ss(h));
plot(t,rep,t,u)
H=time_id(10,'step',rep)
rep=flts(u,tf2ss(H));
plot(t,rep,'.r')
There is a big difference between the original model equation and Identified
model equation.
The responsed curves are also different each other.
Are there any good command so that I can obtain parameters for first order
time delay constant and its gain and dead time?
Thank you foryour help.
Y. Noguchi
--
View this message in context:
http://mailinglists.scilab.org/Model-ID-of-First-order-delay-and-Dead-Time-Process-tp4033659p4033692.html
Sent from the Scilab users - Mailing Lists Archives mailing list archive at
Nabble.com.
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
