% exp4_9_.m
% 計(jì)算系統(tǒng)瞬態(tài)性能指標(biāo)（穩(wěn)態(tài)誤差允許正負(fù)2％）
% 已知系統(tǒng)G(s)=3/(s+1+3i)(s+1-3i)
clc
clear
% 系統(tǒng)模型建立
num=3;
den=conv([1 1+3i],[1 1-3i]);
% 求系統(tǒng)的單位階躍響應(yīng)
[y,x,t]=step(num,den);
% 求響應(yīng)的穩(wěn)態(tài)值
finalvalue=dcgain(num,den)
% 求響應(yīng)的峰值及對(duì)應(yīng)的下標(biāo)
[yss,n]=max(y);
% 計(jì)算超調(diào)量及峰值時(shí)間
percentovershoot=100*(yss-finalvalue)/finalvalue
timetopeak=t(n)
% 計(jì)算上升時(shí)間
n=1;
while y(n)<0.1*finalvalue
   n=n+1;
end
m=1;
while y(m)<0.9*finalvalue
   m=m+1;
end
risetime=t(m)-t(n)
% 計(jì)算調(diào)整時(shí)間
k=length(t);
while (y(k)>0.98*finalvalue)&(y(k)<1.02*finalvalue)
   k=k-1;
end
settlingtime=t(k)