%計算系統性能指標
s1=tf(2000,conv([1 0],[1 20]));
s2=tf([0.1871 0.09353],[1 0.09353]);
sope=s1*s2;
sys=feedback(sope,1);
step(sys)   %繪制階躍響應曲線
[y,t]=step(sys);   %求出階躍響應的函數值及其對應時間
[sigma,tp,ts]=ste(y,t)  %調用函數ste()