clc
clear
close all
%開環傳遞函數描述
num=[20];
den=[1 8 36 40 0];
%求閉環傳遞函數
[numc,denc]=cloop(num,den);
%繪制閉環系統的階躍響應曲線
t=0:0.1:10;
y=step(numc,denc,t);
[y1,x,t1]=step(numc,denc);
%對于傳遞函數調用，狀態變量x返回為空矩陣
plot(t,y,'r:',t1,y1)
title('the step responce')
xlabel('time-sec')
%求穩態值
disp('系統穩態值dc為：')
dc=dcgain(numc,denc)