num=4;
den=[1 1 4];
G=tf(num, den)
step(G)
css=dcgain(G)
[y,t]=step(G);
[ymax, k]=max(y);
tp=t(k);
Mp=(ymax-css)/css*100;
n=1;while y(n)<=css;
n=n+1;end;tr=t(n);
i=length(t);
while(y(i)>0.98*css)&(y(i)<1.02*css)
i=i-1;end;ts=t(i);
disp(['稳态值:css=',num2str(css)])
disp(['峰值时间:tp=',num2str(tp)])
disp(['上升时间:tr=',num2str(tr)])
disp(['超调量:Mp=',num2str(Mp),'%'])
disp(['稳态时间:ts=',num2str(ts)])