【详细易懂】快速傅里叶变换(FFT)生成的频率、幅度具体求解过程

快速傅里叶变换的过程可以用一句话概括:

"通过计算,把信号的配方罗列出来"

即:已知一个时域的随时间变化的信号,该信号可以通过FFT拆分该信号,得到由该信号分解后得到的由不同频率及其幅度组成的频域图,其直观图如下图所示。

1、确定三个参数

这三个参数分别为采样频率Fs、信号最高频率F、采样点数N。

在这里我们假定已知信号最高频率F为600Hz,根据采样定理,我们设定采样频率为1400Hz。

采样定理:采样频率Fs > 信号频率的2倍

采样频率1400Hz即为一秒内有1400个采样点。

最后可得:N=1400 F=600Hz Fs=1400Hz

2、FFT求解生成的幅度
①进行FFT

在上述三个参数的基础上,我们进行FFT,得到 N 个复数。

每个复数都包含着各自特定频率的信息,根据这 N 个复数,我们可以得到原始信号拆分后的各个频率及其对应的幅度值,如下图所示(下图为双边振幅谱)

②求出复数的绝对值

但是直接用复数画出的图不是我们想要的,需要求出全部N个复数的绝对值(即模长),如下图所示,上图中的负数值没有了。

③归一化

最后对纵坐标的振幅需要进行归一化处理才能得到最终的幅度值。

上图中,横坐标为0Hz对应的信号强度(信号振幅),称为直流分量。它对应的振幅求解如下所示:

0Hz对应振幅 = 当前值 / 采样点数N

(即为直流分量对应振幅 = 14000 / 1400 = 10)

注意:直流分量以外的分量所对应的信号振幅求解与上述不太一样,如下所示:

其余频率对应的振幅 = 当前值 /(采样点数N / 2)

(即为:200Hz对应振幅 = 5000 / (1400 / 2) ≈ 7.14)

进行归一化后,我们最终得到如下图所示。

3、FFT求解生成的频率
①频率公式

我们进行FFT后得到的第x个(x从0开始)复数值对应的实际频率为

f(x) = x * (Fs / N)

-f(x) 为实际频率值

-x为第几个复数值

-Fs为采样频率

-N为采样点数

举例:

第0个点的频率 f(0) = 0 * (1400 / 1400) = 0

第1个点的频率 f(0) = 1 * (1400 / 1400) = 1

第2个点的频率 f(0) = 2 * (1400 / 1400) = 2

②删除重复值

由于FFT得到的结果是对称的,因此我们要删除一半的值,只有0 ~ N/2 这一半的频率是有效的,如下图所示。

相关推荐
cxylay12 天前
【FFT】信号处理——快速傅里叶变换【通俗易懂】
fft·频率·快速傅里叶变换·傅里叶变换·时域·频域·幅值
cxylay20 天前
【python版】示波器输出的csv文件(时间与电压数据)如何转换为频率与幅值【方法③】
开发语言·python·示波器·频谱·csv文件·时域·频域
时光飞逝的日子3 个月前
STM32MP135裸机编程:定时器内核时钟频率计算方法
定时器·stm32mp135·频率·裸机·内核时钟
大米粥哥哥4 个月前
Qt 配置Eigen矩阵库 - 并简单测试
开发语言·qt·矩阵·eigen·fft·低通滤波
老歌老听老掉牙6 个月前
概率论基础概念
概率论·事件·频率·随机试验
渣渣小码7 个月前
DFT应用:计算线性卷积
算法·fft·线性卷积·圆周卷积
神也在Copy8 个月前
【深度学习】数据归一化/标准化 Normalization/Standardization
人工智能·深度学习·标准化·归一化
二粒砂9 个月前
从信号处理角度彻底理解FFT
信号处理·fft·傅里叶变换
闻缺陷则喜何志丹10 个月前
【贪心算法】【中位贪心】LeetCode:100123.执行操作使频率分数最大
c++·算法·leetcode·前缀和·贪心算法·频率·中位贪心