ex0521.m

matlab课本程序

function y=Ex0521(z1)
% EX0521    利用函数句柄执行二阶系统时域响应
t=0:0.1:20;
h_plotxy1=str2func('plotxy1')  %创建函数句柄
h_plotxy2=str2func('plotxy2')  %创建函数句柄
h_plotxy3=str2func('plotxy3')  %创建函数句柄
if (z1>=0)&(z1<1)
    y=feval(h_plotxy1,z1,t);  %执行函数
elseif z1==1
    y=feval(h_plotxy2,z1,t);  %执行函数
else
    y=feval(h_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)