% exp4_9_.m
% 計算系統瞬態性能指標（穩態誤差允許正負2％）
% 已知系統G(s)=3/(s+1+3i)(s+1-3i)
clc
clear
% 系統模型建立
num=3;
den=conv([1 1+3i],[1 1-3i]);
% 求系統的單位階躍響應
[y,x,t]=step(num,den);
% 求響應的穩態值
finalvalue=dcgain(num,den)
% 求響應的峰值及對應的下標
[yss,n]=max(y);
% 計算超調量及峰值時間
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)
% 計算調整時間
k=length(t);
while (y(k)>0.98*finalvalue)&(y(k)<1.02*finalvalue)
   k=k-1;
end
settlingtime=t(k)