百科问答小站 logo
百科问答小站 font logo



如何获取FFT序列中每个点的频率值? 第1页

  

user avatar   iosonosole 网友的相关建议: 
      

先来看看一段简单的MATLAB代码

       dt = 0.002; t  = 0:dt:2; f0 = 100; x  = sin(2*pi*f0*t); X  = fft(x,length(x)); plot(abs(X),'linewidth',2); set(gcf,'color','white') axis tight      

很清楚,这是给一个100Hz正弦波做傅立叶变换,得到

很明显得到一个单频信号,不过下标没有对准100Hz。不过没关系。换成下面的代码再跑一次:

       dt = 0.002; t  = 0:dt:2; f0 = 100; N  = length(x); x  = sin(2*pi*f0*t); X  = fft(x,N);  f  = (0:N-1)/N * (1/dt); plot(f,abs(X),'linewidth',2); set(gcf,'color','white') axis tight      


看频谱的时候,只要看前一半就好了。要得到真实频率的下标其实很简单。只要先将频率归一化

       f = (0:N-1)/N;     

然后,再乘以采样率

       f = f * fs;     

这个采样率就是你输入的数据里,相邻两个点之间的间隔的倒数。比如说在你采集的数据里

那么采样率就是

至于为什么要这样做,这是因为在连续傅立叶变换转换到离散傅里叶变换的过程中,频率间隔是这样的




  

相关话题

  Matlab/NumPy/C++Eigen 速度差距为什么很大? 
  如何通过多帧影像进行超分辨率重构? 
  一个函数经过傅立叶变换,再经过傅立叶逆变换得到的函数还是原函数吗? 
  怎么在simulink里面搭建如图所示的模块? 
  一个初中生,做数学作业老是偷偷地用 MATLAB 计算,如何让他知道用 MATLAB 做数学题的危害? 
  控制中拉普拉斯变换是如何推出来的? 
  机器人的视觉伺服控制能在20年内应用于工业界吗? 
  为什么中国程序员工资那么高,连一个 MATLAB 的替代品都开发不出来? 
  想问下专业人士 OpenCv会被深度学习进一步取代吗进一步取代吗? 
  怎样把扭曲的直线图像修直? 

前一个讨论
怎样用非数学语言讲解贝叶斯定理(Bayes's theorem)?
下一个讨论
耳机插头类型不同会影响音质吗?





© 2024-12-26 - tinynew.org. All Rights Reserved.
© 2024-12-26 - tinynew.org. 保留所有权利