ex0516.m

matlab课本程序

function y=Ex0516(z1)
% EX0516    主函数调用子函数，根据阻尼系数绘制二阶系统时域曲线
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)
%画欠阻尼二阶系统时域曲线
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)
%画临界阻尼二阶系统时域曲线
y2=1-exp(-x).*(1+x);
plot(x,y2)

function y3=plotxy3(zeta,x)
%画过阻尼二阶系统时域曲线
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)