% exp4_9_.m
% 計算系統(tǒng)瞬態(tài)性能指標（穩(wěn)態(tài)誤差允許正負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)的峰值及對應(yīng)的下標
[yss,n]=max(y);
% 計算超調(diào)量及峰值時間
percentovershoot=100*(yss-finalvalue)/finalvalue
timetopeak=t(n)
% 計算上升時間
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)
% 計算調(diào)整時間
k=length(t);
while (y(k)>0.98*finalvalue)&(y(k)<1.02*finalvalue)
   k=k-1;
end
settlingtime=t(k)