为什么这组数据做不出傅立叶变换....

bruce001.

一共60000个点，想做FFT，时间步长是(2*pi/800)/100，所以采样频率我选的100/(2*pi/800)，可是画出来的图上没有频率
clc
clear
fin=fopen('plate.txt');
tt=fscanf(fin, '%f\r\n',[60000,1]);
fclose(fin);
plot(tt)

figure
yy1=tt;
fs=100/(2*pi/800);
N=2048;
ky=fft(yy1,N);
n=0:length(ky)-1;   % 从0到N-1，一共N个数
f=fs*n/length(ky);
subplot(211);
plot(f,abs(ky));
title('幅频特性曲线')
xlabel('频率Hz')
[value,index]=max(abs(ky));%将abs（ky)最大值点的横坐标、纵坐标分别赋给value和index
text(f(index),value,sprintf('maxpiont=(%f,%f)',f(index),value),'fontsize',8);%显示最大值点坐标
grid on



plate1.rar (290.61 KB, 下载次数: 2)

2012-11-2 23:33 上传

点击文件名下载附件

犟牛

数据直流量过大，所以fft后最大为直流量分量
而其他频率振动分量很小，所以在图上看不出来

bruce001.

谢谢，我来取中间的频率段画画看

ChaChing

Ref: http://forum.vibunion.com/thread-121006-1-1.html

bruce001.

ChaChing 发表于 2012-11-3 21:08

                               
登录/注册后可看大图

Ref: http://forum.chinavib.com/thread-121006-1-1.html

谢谢...