bruce001. 发表于 2012-11-2 23:36

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

一共60000个点,想做FFT,时间步长是(2*pi/800)/100,所以采样频率我选的100/(2*pi/800),可是画出来的图上没有频率
clc
clear
fin=fopen('plate.txt');
tt=fscanf(fin, '%f\r\n',);
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')
=max(abs(ky));%将abs(ky)最大值点的横坐标、纵坐标分别赋给value和index
text(f(index),value,sprintf('maxpiont=(%f,%f)',f(index),value),'fontsize',8);%显示最大值点坐标
grid on






犟牛 发表于 2012-11-3 06:52

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

bruce001. 发表于 2012-11-3 13:18

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

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

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

bruce001. 发表于 2012-11-5 10:45

ChaChing 发表于 2012-11-3 21:08 static/image/common/back.gif
Ref: http://forum.chinavib.com/thread-121006-1-1.html

谢谢...
页: [1]
查看完整版本: 为什么这组数据做不出傅立叶变换....