Dear Clément,
<br/><br/><br/>Thank you for your suggestion.
<br/><br/>I made following Scilabe code:
<br/><br/>//
<br/>z=poly(0,'z');
<br/>h=(0.065/(z-0.934))*(1/z^10)//    <== 10 Sampling period dead time
<br/>u=zeros(1,100);
<br/>for i=10:1:100
<br/>    u(1,i)=2.0;
<br/>end
<br/>t=1:1:100;
<br/>rep=flts(u,tf2ss(h));
<br/>plot(t,rep,t,u)//  <== We can see the step type process input with amplitude=2 and its process response with 10 sampling period dead time.
<br/>k=find(rep<>0,1) //here the threshold has to be improved in case of noisy signal 
<br/>//H=time_id(1,"step",rep(k:$))
<br/>H=time_id(2,u,rep)
<br/>rep=flts(u,tf2ss(H));     
<br/>plot(t,rep,'.r')// <== We can see the process response by identified model.
<br/>H
<br/><br/><br/>h  =
<br/> 
<br/>       0.065      
<br/>    -----------   
<br/>          10  11  
<br/>  - 0.934z + z
<br/><br/><br/>H  =
<br/> 
<br/>      0.0265880     
<br/>    -------------   
<br/>  - 0.9779092 + z
<br/><br/><br/>h is the discreate transfer function to provide operation data.
<br/>H is the identified transfer function obtained from the opeartion data using time_id.
<br/><br/>I have two issues.
<br/>1. H is not similar to h even the data doesn't include any noise. How I can obtain transfer function nearly same as h?
<br/>2. I would like to have continuous transfer function. How I should convert the discreate transfer function to continuous transfer function.
<br/><br/>May I ask your further advise on these issues?
<br/><br/>Thanks.
<br/><br/><br/><br/>----- 元のメッセージ -----
<br/>差出人: "Clément David-2 [via Scilab / Xcos - Mailing Lists Archives]" <<a href="/user/SendEmail.jtp?type=node&node=4034630&i=0" target="_top" rel="nofollow" link="external">[hidden email]</a>>
<br/>宛先: "Fukashiimo" <<a href="/user/SendEmail.jtp?type=node&node=4034630&i=1" target="_top" rel="nofollow" link="external">[hidden email]</a>>
<br/>送信済み: 2016年9月25日(日曜日) 04:15:09
<br/>件名: Re: System Identification for First order delay and dead time
<br/><br/><br/><br/>Hello, 
<br/><br/>I suggest you to take a look at the `time_id` function [1]. AFAIK this will give you a first idea of the parameters. If you need more, I suggest you to take a look at the "Optimization" topic in the help and implement a custom cost function. 
<br/><br/>@samuel : IMHO `s` is supposed to be the %s Scilab variable (the Laplace variable also named `p`) ; `u` is the input ; `y` is the output. 
<br/><br/>[1]: <a href="https://help.scilab.org/docs/5.5.2/en_US/time_id.html" target="_top" rel="nofollow" link="external">https://help.scilab.org/docs/5.5.2/en_US/time_id.html</a> 
<br/><br/>Thanks, 
<br/><br/>-- 
<br/>Clément 
<br/>_______________________________________________ 
<br/>users mailing list 
<br/>[hidden email] 
<br/><a href="http://lists.scilab.org/mailman/listinfo/users" target="_top" rel="nofollow" link="external">http://lists.scilab.org/mailman/listinfo/users</a> 
<br/><br/><br/><br/><br/><br/>If you reply to this email, your message will be added to the discussion below: <a href="http://mailinglists.scilab.org/System-Identification-for-First-order-delay-and-dead-time-tp4034608p4034612.html" target="_top" rel="nofollow" link="external">http://mailinglists.scilab.org/System-Identification-for-First-order-delay-and-dead-time-tp4034608p4034612.html</a> 
<br/>To unsubscribe from System Identification for First order delay and dead time, click here . 
<br/>NAML
<br/>

        
        
        
<br/><hr align="left" width="300" />
View this message in context: <a href="http://mailinglists.scilab.org/Scilab-System-Identification-for-First-order-delay-and-dead-time-tp4034630.html">Scilab: System Identification for First order delay and dead time</a><br/>
Sent from the <a href="http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html">Scilab users - Mailing Lists Archives mailing list archive</a> at Nabble.com.<br/>