文档介绍:clear all;
t=0::12;
for s=0:: %欠阻尼二阶系统单位阶跃响应
h=l-exp(-s*t).*sin(t*(l-s*s)+acos (s)); plot(t,h,'g');%绿线表clear all;
t=0::12;
for s=0:: %欠阻尼二阶系统单位阶跃响应
h=l-exp(-s*t).*sin(t*(l-s*s)+acos (s)); plot(t,h,'g');%绿线表示欠阻尼响应曲线 axis([0,12,0,2]);
hold on;
end
h=l-exp(-t).*(l+t);%临界阻尼二阶系统单位阶跃响应 plot(t,h,'b');%蓝线表示临界阻尼响应曲线 axis([0,12,0,2]);
title('二阶系统单位阶跃响应');
hold on;
for s=::2 %过阻尼二阶系统单位阶跃响应
Tl=s-(s*s-l);
T2=s+(s*s-l); 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;
二阶系统单位阶跃响应