频谱图与响应问题

Arick0420

　　[num dof]=size(disp);
　　for i=1:dof;
　　f(:,i)=disp(:,i);
　　dt=0.01;
　　w=pi/dt;
　　Nt=num;
　　Nw=Nt/2;
　　dw=w/Nw;
　　%ff=fft(f)*dt;
　　T=Nt*dt;
　　sf(:,i)=1/(2*pi*T)*(abs(fft(f(:,i)))*dt).^2;
　　%sf1=sum(sf,2);
　　%cc=ifft(sqrt(sf*2*pi*T)/dt);
　　figure(i)
　　plot([0:dw:(Nw-1)*dw],sf(1:Nw,i))
　　title('傅氏频谱图')
　　%xlim([dw Nw*dw])
　　xlabel('Frequency (rad/sec)')
　　ylabel('H(w) (m/N)')
　　%ylabel('|Fx(w)|')
　　%xlim([0*dw 10000*dw])
　　%legend('1-DOf','2-DOF','3-DOF')
　　xlim([0 400])
　　end
　　请问以上这些要怎么解释?? 例如 sf(:,i)=1/(2*pi*T)*(abs(fft(f(:,i)))*dt).^2
　　这是甚么意思?为什么有2*pi*T? 还有 dt 跟 dw 有甚么不同跟影响 ? 拜托请帮我! 还有我用这两组档案下去跑响应与频谱图
　　用newmark法的dt设定有限制吗? 我发现Q1：newmark的dt只有带0.01才是比较像正确的响应,但频谱图只能对应到w1。 Q2：两边 dt 同时带0.0001时，响应图是错的，但可以对应到完整的频谱图。 Q3：
　　要得到一组正确的响应图与可以对应的频谱图该改哪些东西呢? 我
　　知道我问题很多很杂
　　请版友帮帮忙
　　谢谢