独孤一1984 发表于 2012-7-26 16:25

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\,请问为什么横轴频率的坐标没变,可从图形可以明显看出一幅图属于下一幅图的局部,那为什么横轴坐标不变,我的程序哪里出了问题?

rogen 发表于 2012-7-28 13:41

从你的程序来看,你的分辨率就没有变化,y=fft(x,N)才能通过改变N的值,实现分辨率的变化!而在你的程序中你的N没有设置,所以出现这种情况。你的程序其实是有问题的,当你的N取大于20000的时候,你的F与Ayy的数据点就不对应了,程序就错了。现在你的程序能运行,主要是你的F数据点数刚好小于Ayy的数据点数。主要问题就在你的fft不合适!

独孤一1984 发表于 2012-7-28 14:57

回复 2 # rogen 的帖子

昨天我用MATLAB运行,仔细看每个变量的值,也发现了这个问题。谢谢你!
页: [1]
查看完整版本: fft横坐标与分辨率的关系