不同阻尼比情况下，频率响应函数实部、虚部和模值随外激励频率变化曲线

clear
jt=0.04;
k=100;m=1;c=2*sqrt（m*k）*jt;
omegan=sqrt（k/m）;omega=0:0.1:2*omegan;
h=1./（k-m*omega.^2+i*c*omega）;
reh=real（h）;imh=imag（h）;absh=abs（h）;
figure（1）
plot（omegareh‘linewidth‘2）;
text（10real（1./（k-m*（0.95*omegan）.^2+i*c*0.95*omegan））sprintf（‘%2.2f‘jt）‘fontsize‘10）;
title（‘频率响应函数实部响应曲线‘‘fontsize‘10）;
xlabel（‘\omega（\omega_n=10）‘）;ylabel（‘ReH（\omega）‘）;
hold on;grid on;
figure（2）
plot（omegaimh‘linewi