[Scilab-users] Model ID of First order delay and Dead Time Process

[Serge Steer]

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 at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
>