%系統校正后性能指標及階躍響應曲線
global y t;
k0=1000;n1=1;d1=conv([1 0],[1,2]);
s1=tf(k0*n1,d1);
n2=[0.04916 1];d2=[0.008438 1]；
s2=tf(n2,d2);
sope=s1*s2;
sys=feedback(sope,1);
step(sys)   %繪制階躍響應曲線
[y,t]=step(sys);   %求出階躍響應的函數值及其對應時間
[sigma,tp,ts]=ste(y,t)  %調用函數ste()