二维FFT疑问

xinglong-liu

我对纵波在铝梁上的传播作了2维FFT分析。有点疑问需要大家帮忙解决。
图1是我输入的激振波形，频率为400.4kHz.

利用有限元计算后，得到位移-时间的二维分布图

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


根据该图可以估算出纵波的相速度约为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)

2007-11-27 21:33 上传

点击文件名下载附件

[ 本帖最后由 xinglong-liu 于 2007-11-27 21:38 编辑 ]

zsyjacky

'FFT2D_data.csv'是倒进去的吧！

xinglong-liu

是的。FFT2D_data.csv是我从ansys分析结果中导入的。结果数据就是附件中的sig.mat.

xinglong-liu

无人问津？期待中。。。