二维FFT疑问
我对纵波在铝梁上的传播作了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=';
Nf=round(length(freq)/2);
freq=freq(1:Nf);
s_nfft=1024;
es=8e-4;wl=1/es;
wl=;
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)');
数据文件
[ 本帖最后由 xinglong-liu 于 2007-11-27 21:38 编辑 ]
看不出来！
'FFT2D_data.csv'是倒进去的吧！ 是的。FFT2D_data.csv是我从ansys分析结果中导入的。结果数据就是附件中的sig.mat. 无人问津？期待中。。。
页:
[1]