ex0506.m

matlab课本程序

function y=Ex0506(zeta)
% EX0506    使用if结构的二阶系统时域响应
x=0:0.1:20;
if (zeta>0)&(zeta<1)
    y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta));
elseif zeta==1
    y=1-exp(-x).*(1+x);
end
plot(x,y)