clc
clear
close
%開環傳遞函數描述
numo=20;
deno=[1 8 36 40 0];
%求閉環傳遞函數
[numc,denc]=cloop(numo,deno,-1);
%繪制閉環系統的脈沖激勵響應曲線
t=1:0.1:10;
[y,x]=impulse(numc,denc,t);
plot(t,y)
title('the impulse responce')
xlabel('time-sec')