b=1;a=[1 5 6];sys=tf(b,a);y1=step(sys);y2=impulse(sys);plot(y1,'linewidth',2)hold onplot(y2,'linewidth',2)legend('单位阶跃响应','单位脉冲响应')set(gcf,'Color','w')title('$\frac{1}{s^2+5s+6}$','interpreter','latex', 'FontSize', 18);