% 单自由度系统的强迫振动不同阻尼比下的传递率曲线
lambda=linspace(0,3,100);
for zeta=0.2:0.1:1
tr=(1+(2*zeta.*lambda).^2)./((1-lambda.^2).^2+(2*zeta.*lambda).^2);
plot(lambda,tr);
grid on;
hold on
end
title('不同阻尼比下的传递率曲线')
xlabel('频率比')
ylabel('传递率')
gtext({'\zeta=0.125';'\zeta=0.25';'\zeta=0.5';'\zeta=0.707';'\zeta=1'})