STM32 FFT算法实现

然后，打开工程，新建 DSP_LIB 分组，并将 arm_cortexM4lf_math.lib 添加到工程里面，如图：

2， 添加头文件包含路径
添加好.lib 文件后，我们要添加头文件包含路径，将第一步拷贝的 Include 文件夹和 DSP_LIB文件夹，加入头文件包含路径。
3， 添加全局宏定义
最后，为了使用 DSP 库的所有功能，我们还需要添加几个全局宏定义：
1，__FPU_USED
2，__FPU_PRESENT
3，ARM_MATH_CM4
4，__CC_ARM
5，ARM_MATH_MATRIX_CHECK
6，ARM_MATH_ROUNDING
添加方法：点击C/C++选项卡，然后在 Define 里面进行设置，如图

这里，两个宏之间用“,”隔开。并且，上面的全局宏里面，我们没有添加__FPU_USED，因为这个宏定义在 Target 选项卡设置 Code Generation 的时候，选择了：Use FPU（如果没有设置 Use FPU，则必须设置！！），故 MDK 会自动添加这个全局宏，因此不需要我们手动添加了。
至此，STM32F4 的 DSP 库运行环境就搭建完成了。

照一定大小采样频率 F 去采集信号，采集 N 个点，那么通过对这 N 个点进行 FFT 运算，就可以得到这个信号的频谱特性。 这里还涉及到一个采样定理的概念：在进行模拟/数字信号的转换过程中，当采样频率 F 大于信号中最高频率 fmax的 2 倍时(F>2*fmax)，采样之后的数字信号完整地保留了原始信号中的信息，采样定理又称奈奎斯特定理。举个简单的例子：比如我们正常人发声，频率范围一般在 8KHz以内，那么我们要通过采样之后的数据来恢复声音，我们的采样频率必须为 8KHz 的 2 倍以上，也就是必须大于 16KHz 才行。

模拟信号经过 ADC 采样之后，就变成了数字信号，采样得到的数字信号，就可以做 FFT变换了。N 个采样点数据，在经过 FFT 之后，就可以得到 N 个点的 FFT 结果。为了方便进行FFT 运算，通常 N 取 2 的整数次方。

假设采样频率为 F，对一个信号采样，采样点数为 N，那么 FFT 之后结果就是一个 N 点的复数，每一个点就对应着一个频率点（以基波频率为单位递增），这个点的模值（sqrt(实部虚部2)）就是该频点频率值下的幅度特性。具体跟原始信号的幅度有什么关系呢？假设原始信号的峰值为 A，那么 FFT 的结果的每个点（除了第一个点直流分量之外）的模值就是 A 的 N/2倍，而第一个点就是直流分量，它的模值就是直流分量的 N 倍。

这里还有个基波频率，也叫频率分辨率的概念，就是如果我们按照 F 的采样频率去采集一个信号，一共采集 N 个点，那么基波频率（频率分辨率）就是 fk=F/N。 这样,第 n 个点对应信号频率为：F*(n-1)/N；其中 n≥1，当 n=1 时为直流分量。
关于 FFT 就介绍到这。

arm_status arm_cfft_radix4_init_f32( 
arm_cfft_radix4_instance_f32 * S, 
uint16_t fftLen,uint8_t ifftFlag,uint8_t bitReverseFlag)  
void arm_cfft_radix4_f32(const arm_cfft_radix4_instance_f32 * S,float32_t * pSrc) 
void arm_cmplx_mag_f32(float32_t * pSrc,float32_t * pDst,uint32_t numSamples) 

第一个函数 arm_cfft_radix4_init_f32，用于初始化 FFT 运算相关参数，其中：fftLen 用于指定 FFT 长度（16/64/256/1024/4096），本章设置为1024；ifftFlag 用于指定是傅里叶变换(0)还是反傅里叶变换(1)，本次设置为0；bitReverseFlag 用于设置是否按位取反，本章设置为 1；最后，所有这些参数存储在一个 arm_cfft_radix4_instance_f32 结构体指针 S 里面。

第二个函数 arm_cfft_radix4_f32 就是执行基 4 浮点 FFT 运算的，pSrc 传入采集到的输入信号数据（实部+虚部形式），同时 FFT 变换后的数据，也按顺序存放在 pSrc 里面，pSrc 必须大于等于 2 倍 fftLen 长度。另外，S 结构体指针参数是先由 arm_cfft_radix4_init_f32 函数设置好，然后传入该函数的。

第三个函数 arm_cmplx_mag_f32 用于计算复数模值，可以对 FFT 变换后的结果数据，执行取模操作。pSrc 为复数输入数组（大小为 2*numSamples）指针，指向 FFT 变换后的结果；pDst为输出数组（大小为 numSamples）指针，存储取模后的值；numSamples 就是总共有多少个数据需要取模。
通过这三个函数，我们便可以完成 FFT 计算，求出模值后便可以求出幅值、相位，知道采样率和FFT长度就可求出频率。

fft.c 代码

#include "arm_math.h"
#include "common.h"
#include "fft.h"


/*********************************************************************************
************************STM32F407核心开发板******************************
**********************************************************************************
* 文件名称: fft.c                                                             *
* 文件简述：DSP FFT使用                                                           *
* 创建日期：2020.08.21                                                           *
* 版    本：V1.0                                                                 *
* 作    者：近视未看清人心                                                              *
* 说    明：完成FFT信号分析需要用到的函数集，包括求模，幅值、频率、相位
注意：FFT的长度，也就是采样点的个数在fft.h 宏定义中完成* 
**********************************************************************************
*********************************************************************************/	

arm_cfft_radix4_instance_f32 scfft;	//设置FFT参数的结构体变量（基4）

/****************************************************************************
* 名    称: void FFTx4_init(u8 ifftFlag,u8 bitReverseFlag)
* 功