fft横坐标与分辨率的关系

独孤一1984

1. clc;
   
2. clear;
   
3. x=load('data.txt');
   
4. 
   
5. figure;
   
6. subplot(2,1,1);
   
7. plot(x);
   
8. 
   
9. x=detrend(x)
   
10. 
    
11. fs=20000
    
12. N=5000
    
13. 
    
14. Y=fft(x)
    
15. Ayy = (abs(Y)); %取模
    
16. 
    
17. Ayy=Ayy/(N/2); %换算成实际的幅度
    
18. Ayy(1)=Ayy(1)/2;
    
19. F=([1:N]-1)*fs/N; %换算成实际的频率值
    
20. subplot(212)
    
21. plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
    
22. title('幅度-频率曲线图');
    
23. 
    

复制代码

程序如上，数据文本为data.txt，采样频率设置为20K，采样点数分别取1000,5000,10000,20000做测试，分辨率随之变化：20Hz，4Hz，2Hz，1Hz。如图figure1\figure2\figure3\figure4\，请问为什么横轴频率的坐标没变，可从图形可以明显看出一幅图属于下一幅图的局部，那为什么横轴坐标不变，我的程序哪里出了问题？

rogen

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

独孤一1984

回复 2 # rogen 的帖子

昨天我用MATLAB运行，仔细看每个变量的值，也发现了这个问题。谢谢你！