test4_2_21.m

来自「两点间的头镇试验源码」· M 代码 · 共 19 行

% 实验4第二题的第2小题，绘制的kesai范围改动一下，kesai=0.2:0.3:2 
% 这是第一种方法，不求用不线形来标识，手标

clear
t=linspace(0,18,1000);
figure
for xi=0.2:0.3:2.0
    beta=sqrt(1-xi^2);
    theta=atan(sqrt(1-xi^2)/xi);
    y=1-1./beta.*exp(-xi.*t).*sin(beta*t+theta);
    plot(t,y),hold on
    pause        
end
gtext('{\xi}=0.2');
gtext('{\xi}=2');
xlabel('t'),ylabel('y')
title('二阶系统响应曲线')
hold off