fft横坐标与分辨率的关系
clc;clear;
x=load('data.txt');
figure;
subplot(2,1,1);
plot(x);
x=detrend(x)
fs=20000
N=5000
Y=fft(x)
Ayy = (abs(Y)); %取模
Ayy=Ayy/(N/2); %换算成实际的幅度
Ayy(1)=Ayy(1)/2;
F=(-1)*fs/N; %换算成实际的频率值
subplot(212)
plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
title('幅度-频率曲线图');
程序如上,数据文本为data.txt,采样频率设置为20K,采样点数分别取1000,5000,10000,20000做测试,分辨率随之变化:20Hz,4Hz,2Hz,1Hz。如图figure1\figure2\figure3\figure4\,请问为什么横轴频率的坐标没变,可从图形可以明显看出一幅图属于下一幅图的局部,那为什么横轴坐标不变,我的程序哪里出了问题? 从你的程序来看,你的分辨率就没有变化,y=fft(x,N)才能通过改变N的值,实现分辨率的变化!而在你的程序中你的N没有设置,所以出现这种情况。你的程序其实是有问题的,当你的N取大于20000的时候,你的F与Ayy的数据点就不对应了,程序就错了。现在你的程序能运行,主要是你的F数据点数刚好小于Ayy的数据点数。主要问题就在你的fft不合适! 回复 2 # rogen 的帖子
昨天我用MATLAB运行,仔细看每个变量的值,也发现了这个问题。谢谢你!
页:
[1]