← Back to branch summary

~ubuntu-branches/ubuntu/karmic/scilab/karmic

« back to all changes in this revision

Viewing changes to man/control/time_id.cat

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2002-03-21 16:57:43 UTC
  • Revision ID: james.westby@ubuntu.com-20020321165743-e9mv12c1tb1plztg
Tags: upstream-2.6
Import upstream version 2.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
man/control/time_id.cat
 
1
time_id           Scilab Group           Scilab Function            time_id
 
2
NAME
 
3
   time_id - SISO least square identification
 
4
  
 
5
CALLING SEQUENCE
 
6
 [H [,err]]=time_id(n,u,y)
 
7
PARAMETERS
 
8
 n          : order of transfer
 
9
            
 
10
 u          : one of the following
 
11
            
 
12
           u1           : a vector of inputs to the system
 
13
                        
 
14
           "impuls"     : if y is an impulse response
 
15
                        
 
16
           "step"       : if y is a step response.
 
17
                        
 
18
 y          : vector of response. 
 
19
            
 
20
 H          : rational function with degree n denominator and  degree n-1
 
21
            numerator if y(1)==0 or rational function with degree n
 
22
            denominator and  numerator if y(1)<>0.
 
23
            
 
24
 err        : ||y - impuls(H,npt)||^2, where impuls(H,npt) are the npt
 
25
            first coefficients of impulse response of H
 
26
            
 
27
DESCRIPTION
 
28
   Identification of discrete time response. If y is strictly proper
 
29
  (y(1)=0) then time_id computes the least square solution of the linear
 
30
  equation:  Den*y-Num*u=0 with the constraint  coeff(Den,n):=1. if y(1)~=0
 
31
  then the algorithm first computes the proper part solution and then add 
 
32
  y(1) to the solution
 
33
  
 
34
EXAMPLE
 
35
 z=poly(0,'z');
 
36
 h=(1-2*z)/(z^2-0.5*z+5)
 
37
 rep=[0;ldiv(h('num'),h('den'),20)]; //impulse response
 
38
 H=time_id(2,'impuls',rep)
 
39
 //  Same example with flts and u
 
40
 u=zeros(1,20);u(1)=1;
 
41
 rep=flts(u,tf2ss(h));        //impulse response
 
42
 H=time_id(2,u,rep)
 
43
 //  step response
 
44
 u=ones(1,20);
 
45
 rep=flts(u,tf2ss(h));     //step response.
 
46
 H=time_id(2,'step',rep)
 
47
 H=time_id(3,u,rep)    //with u as input and too high order required
 
48
AUTHOR
 
49
   Serge Steer INRIA
 
50
  
 
51
SEE ALSO
 
52
   imrep2ss, arl2, armax, frep2tf
 
53