为什么这组数据做不出傅立叶变换....
一共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
数据直流量过大,所以fft后最大为直流量分量
而其他频率振动分量很小,所以在图上看不出来 谢谢,我来取中间的频率段画画看
Ref: http://forum.vibunion.com/thread-121006-1-1.html ChaChing 发表于 2012-11-3 21:08 static/image/common/back.gif
Ref: http://forum.chinavib.com/thread-121006-1-1.html
谢谢...
页:
[1]