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