收起
你好代码供参考
y1 = @(x,xi) 1 - 1/sqrt(1-xi^2)*exp(-xi*x).*sin(sqrt(1-xi^2)*x+acos(xi));
y2 = @(x,xi) 1 - (1+x).*exp(-x);
y3 = @(x,xi) 1 - 1/(2*sqrt(xi^2-1)).*...
(exp(-(xi-sqrt(xi^2-1))*x)./(xi-sqrt(1-xi^2)-...
exp(-(xi+sqrt(xi^2-1))*x)./(xi+sqrt(1-xi^2))));
x = linspace(0,20,101);
plot(x,y1(x,0.5),'r--');hold on;
plot(x,y2(x,1),'b-.')
plot(x,y3(x,2),'k-')
legend('欠阻尼','临界阻尼','过阻尼')
有帮助望采纳
报告相同问题?