Program in Scilab:
t=0:0.2:10;
zeta=[0 0.2 0.4 0.6 0.8 1];
cv=[1 2 3 4 5 6];
clf;
for n=1:6
num=poly([1 0 0],'s','coeff');
den=poly([1 2*zeta(n) 1],'s','coeff');
g=syslin('c',num/den);
gs=csim('step',t,g);
plot2d(t,gs,style=cv(n))
end;
xgrid;
xtitle('response to unit-step with values: wn=1
,zeta=0,0.2,0.4,0.6,0.8,1','t(seg)','Output');
legends(['zeta=0','zeta=0.2','zeta=0.4','zeta=0.6'
,'zeta=0.8','zeta=1'],[1,2,3,4,5,6],opt=4);
|
|