声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 994|回复: 0

[FFT] 频谱图与响应问题

[复制链接]
发表于 2008-11-4 01:33 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-9-13 13:31 编辑

  [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:
  要得到一组正确的响应图与可以对应的频谱图该改哪些东西呢? 我
  知道我问题很多很杂
  请版友帮帮忙
  谢谢

FFT.m

708 Bytes, 下载次数: 8

BDisp3_free4_10000.m

1.62 KB, 下载次数: 8

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 19:43 , Processed in 0.058011 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表