www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第3章/li3_22.m

    t=[0:0.1:18]';
for x=0.2:0.2:0.8
    b=sqrt([1-x^2]);
    z=atan(b/x);
    y1=-t*x;
    y2=t*b+z;
    y=1-exp(y1).*sin(y2)/b;
    plot(t,y),hold on
end
xlabel('t(秒)');ylabel('y');
title('二阶系统阶跃响应');
text(3.3,0.9,'{\xi}=0.8')
text(4.3,1.4,'{\xi}=0.2')