XML
Laplace.T Root-locus Transient response Frecuency response
       
 
 
 

10 Program 5.18 Ogata 4edition (Response to mechanical vibratory system)



    Let us draw the solution of Example, response to mechanical vibratory system by Scilab:

     

    Program in Scilab:
    num=poly([100 10 0],'s','coeff');
    
    den=poly([100 10 1],'s','coeff');
    
    t1=0:0.001:0.537;
    
    t2=0.538:0.001:1.5;
    
    x1=2.452*(t1^2);
    
    x2=0.707*ones(t2);
    
    y=syslin('c',num/den);
    
    x=[x1 x2];
    
    t=[t1 t2];
    
    g=csim(x,t,y);
    
    clf;
    
    plot2d(t,-x,2);
    
    plot2d(t,-g,5);
    
    xgrid;
    
    xtitle('Response to mechanical
     vibratory system','t(seg)'
    ,'Input negative X  and output
     negative Y');legends(['X','Y']
    ,[2,5],opt=4);
    

    Image Programa5_18