声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1282|回复: 3

[FFT] 二维FFT疑问

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

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

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

x
我对纵波在铝梁上的传播作了2维FFT分析。有点疑问需要大家帮忙解决。
图1是我输入的激振波形,频率为400.4kHz.
input waveform.png
利用有限元计算后,得到位移-时间的二维分布图 sig.png

采样周期8e-8s,空间采样间隔8e-4m.2维FFT分析后得到contour图,频率-波数关系
contour.png

根据该图可以估算出纵波的相速度约为5000m/s.但是为何频率集中在0-200kHz之间。我觉得频率应该分布在400kHz周围,波数应该在80左右。
以下是我的程序
sig=csvread('FFT2D_data.csv',0,0);
ts=8e-8;es=8e-4;
t=linspace(1,250,250);space=linspace(1,120,120);
t=ts*t;space=es*space;
mesh(space,t,sig)
xlabel('Space(m)');ylabel('Time(s)');
max_R=512;max_C=120;
sig=data1;
ts=8e-8;fs=1/ts;
t_nfft=2048;
freq=[fs*(0:t_nfft-1)/t_nfft]';
Nf=round(length(freq)/2);
freq=freq(1:Nf);

s_nfft=1024;
es=8e-4;wl=1/es;
wl=[wl*(0:s_nfft-1)/s_nfft];
Nwl=round(length(wl)/2);wl=wl(1:Nwl);
F=abs(fft2(sig,t_nfft,s_nfft));
mesh(F);contour(F);
F=F(1:length(freq),1:length(wl));
contour(wl,freq,F);
xlabel('wavenumber(1/\lambda(m))');ylabel('Freq(Hz)');

数据文件 sig.mat (120.05 KB, 下载次数: 8)

[ 本帖最后由 xinglong-liu 于 2007-11-27 21:38 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-11-27 22:32 | 显示全部楼层

看不出来!

'FFT2D_data.csv'是倒进去的吧!
 楼主| 发表于 2007-11-28 09:25 | 显示全部楼层
是的。FFT2D_data.csv是我从ansys分析结果中导入的。结果数据就是附件中的sig.mat.
 楼主| 发表于 2007-11-28 22:28 | 显示全部楼层
无人问津?期待中。。。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 23:17 , Processed in 0.073676 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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