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);
|
|