huitu.m

matlab实用教程

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')