首先，對微分方程進行拉氏變換，求其傳遞函數為：G(s)=y(s)/e(s)=(s^2+6s+4)/(s^2+6s+10)
(1)Matlab?編程：
sys=tf([1?6?4],[1?7?10]);
figure(1)?%脈沖響應
impulse(sys)
figure(2)?%階躍響應
step(sys)
（2）Simulink中按下圖連接：