function y=Ex0516(z1)
% EX0516    主函數(shù)調(diào)用子函數(shù)，根據(jù)阻尼系數(shù)繪制二階系統(tǒng)時域曲線
t=0:0.1:20;
if (z1>=0)&(z1<1)
    y=plotxy1(z1,t);
elseif z1==1
    y=plotxy2(z1,t);
else
    y=plotxy3(z1,t);
end

function y1=plotxy1(zeta,x)
%畫欠阻尼二階系統(tǒng)時域曲線
y1=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));
plot(x,y1)

function y2=plotxy2(zeta,x)
%畫臨界阻尼二階系統(tǒng)時域曲線
y2=1-exp(-x).*(1+x);
plot(x,y2)

function y3=plotxy3(zeta,x)
%畫過阻尼二階系統(tǒng)時域曲線
y3=1-1/(2*sqrt(zeta^2-1))*(exp(-((zeta-sqrt(zeta^2-1))*x))./(zeta-sqrt(zeta^2-1))...
    -exp(-((zeta+sqrt(zeta^2-1))*x))./(zeta+sqrt(zeta^2-1)));
plot(x,y3)