FFT结果不理想是怎么回事呢
请问大家，FFT变换后，想求得信号频率，但是结果峰值并不明显，第一个峰值是直流量，是数据的问题还是我写的代码有问题？数据第一、二列是时间，第三列是气压值clear;
fid=fopen('2016-05-30-10-24-08.txt');
c=textscan(fid,'%s%s%f');
%链接两个字符串
d=strcat(c{1,1},'/',(c{1,2}));
%将时间字符串转为数字
ymd=datenum(d);
%验证转化的时间是否正确
r=datestr(ymd,'dd-mmm-yyyy HH:MM:SS')
subplot(2,2,1)
plot(ymd,c{1,3});
datetick('x','mm-dd/HH');
xlabel('time');
ylabel('f(t)');
%傅里叶变换
fs=1;
L=length(ymd);
N=length(c{1,3});
n=0:L-1;
t=ymd;
dat=c{1,3}-mean(c{1,3});
tp=;
save tp.mat;
Y=fft(dat,N);
mag=abs(Y);%求幅值
%相位
P=angle(Y);
%频率
f=(0:N-1)*fs/N;
%求峰值，得到幅度最大值及下标
=max(mag(1:N/2)*2/N);
fm=f(Q);
subplot(2,2,2);
plot(f(1:N/2),mag(1:N/2)*2/N);
set(gca, 'Xlim',);
xlabel('frequence/Hz');
ylabel('amplitude');
%title('f=(0:N-1)*fs/Lfs=1Hz 2016-05-25');
title('fs=1Hz 2016-05-25');
subplot(2,2,3)
plot(f(1:N/2),P(1:N/2));
xlabel('frequence/Hz');
ylabel('phase');
%set(gca, 'XTick',[-3.14,0,3.14] ); 哪位大神帮我能看下{:{19}:} 改一下坐标轴上下限的值 试试会不会好一些 sovereign 发表于 2016-11-29 08:20
改一下坐标轴上下限的值 试试会不会好一些
第一个峰值是直流分量，没有意义，我想得到信号的频率，但是看图中振幅非常小，您能帮我看下我的代码是不是没问题吧？ xiujuan507 发表于 2016-11-29 10:23
第一个峰值是直流分量，没有意义，我想得到信号的频率，但是看图中振幅非常小，您能帮我看下我的代码是不 ...
你频率范围也太小了 sovereign 发表于 2016-11-29 13:32
你频率范围也太小了
这个信号频率就是很低啊，10的-3次方以下，我的代码没有错吧，大神 加窗过滤吧 放学别走 发表于 2016-11-29 16:09
加窗过滤吧
过滤后的效果能是什么样子的呢？我做出来的很不光滑
xiujuan507 发表于 2016-11-29 17:11
过滤后的效果能是什么样子的呢？我做出来的很不光滑
难道还需要插值吗？ 敷衍会致命 发表于 2016-11-30 08:31
难道还需要插值吗？
窗函数用来平滑的 貌似代码没什么问题加入窗函数之后再看看结果吧估计差别不大 think2015 发表于 2016-11-30 13:25
貌似代码没什么问题加入窗函数之后再看看结果吧估计差别不大
多谢，第一个波峰是直流分量，没什么意义，其它波峰的幅值较小，请问加窗后直流分量能去掉吗？ xiujuan507 发表于 2016-11-30 18:43
多谢，第一个波峰是直流分量，没什么意义，其它波峰的幅值较小，请问加窗后直流分量能去掉吗？
因为我想得到信号的频率，但是傅里叶变换后，除了第一个直流分量的频率，其它的都不是很明显？不知道这是为什么呢？ xiujuan507 发表于 2016-11-30 18:43
多谢，第一个波峰是直流分量，没什么意义，其它波峰的幅值较小，请问加窗后直流分量能去掉吗？
够呛应该去不掉去直流的话 i1=i-mean（i）试试吧 think2015 发表于 2016-12-1 08:16
够呛应该去不掉去直流的话 i1=i-mean（i）试试吧
这个图就是已经减去平均值了
页:
[1]
2