1条回答 默认 最新
- joel_1993 2021-10-30 09:53关注
你好代码供参考
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('欠阻尼','临界阻尼','过阻尼')
有帮助望采纳
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 1无用