ts=0 te=10 dt=0.01

sys=tf([1 6 4],[1 7 10])

t=ts:dt:te

f=4*sin(t)

y=lsim(sys,f,t)

plot(t,y) xlabel( t(sec) ) ylabel( y(t) ) grid on

三

ts=0 te=5 dt=0.01

sys=tf([1 6 4],[1 7 10])

t=ts:dt:te

y1=impulse(sys,t)

y2=step(sys,t)

plot(t,y1) xlabel( t(sec) ) y