二阶系统单位阶跃响应曲线

clear all;

t=0:0.01:12;

for s=0:0.1:0.9 %欠阻尼二阶系统单位阶跃响应h=1-exp(-s*t).*sin(t*(1-s*s)^0.5+acos(s));

plot(t,h,'g');%绿线表示欠阻尼响应曲线

axis([0,12,0,2]);

hold on;

end

h=1-exp(-t).*(1+t);%临界阻尼二阶系统单位阶跃响应plot(t,h,'b');%蓝线表示临界阻尼响应曲线

axis([0,12,0,2]);

title('二阶系统单位阶跃响应');

hold on;

for s=1.1:0.1:2 %过阻尼二阶系统单位阶跃响应T1=s-(s*s-1)^0.5;

T2=s+(s*s-1)^0.5;

h=1+exp(-t*T1)/(T1/T2-1)+exp(-t*T2)/(T2/T1-1);

plot(t,h,'r');%红线表示过阻尼响应曲线

axis([0,12,0,2]);

hold on;

end

hold off;