脉冲syms s zetazeta=0.707;num=16;den=1 8*zeta 16;p=roots(den);sys=tf(num, den);t=0:0.01:3;figure(1)impulse(sys,t);gridxlabel(t);ylabel(c(t);title(impulse response);阶跃syms s zetazeta=0.707;num=16;den=1 8*zeta 16;p=roots(den);sys=tf(num, den);t=0:0.01:3;figure(2)step(sys,t);gridxlabel(t);ylabel(c(t);title(step response);斜坡syms s zetazeta=0.707;num=16;den=1 8*zeta 16;p=roots(den);sys=tf(num, den);t=0:0.01:3;figure(3)u=t;lsim(sys,u,t,0);gridxlabel(t);ylabel(c(t);title(ramp response);