|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我对纵波在铝梁上的传播作了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)
[ 本帖最后由 xinglong-liu 于 2007-11-27 21:38 编辑 ] |
|