如何在STM32F407G中使用实部FFT（RFFT）函数？

我明白你在STM32F407上用CMSIS-DSP库处理实数据FFT时遇到的困扰——复数据的CFFT能正常跑，但实数据的RFFT总是出问题。其实核心是实FFT的调用流程和复FFT有差异，下面我给你拆解正确的实现步骤，附完整代码示例。

关键知识点 #

CMSIS-DSP针对实数据提供了快速实FFT接口（arm_rfft_fast_xxx系列），它利用实数据FFT的共轭对称性，比普通复FFT更高效，完全不需要手动处理虚部输入。

完整实现步骤与代码 #

1. 必要头文件与宏定义 #

先确保包含CMSIS-DSP相关头文件，采样长度必须是2的幂（比如32、64、128等）：

#include "stm32f407xx.h"
#include "arm_math.h"
#include "arm_const_structs.h"
#include "core_cm4.h"

#define TEST_LENGTH_SAMPLES 32  // 必须是2的幂，这里和你的定义保持一致

2. 定义输入输出数组 #

准备实数据输入数组，以及存储FFT复数输出、幅度谱的数组：

float32_t ffttestrealip[TEST_LENGTH_SAMPLES];  // 你的实数据输入
static float32_t fft_output[TEST_LENGTH_SAMPLES];  // 存储复数频谱输出
static float32_t fft_magnitude[TEST_LENGTH_SAMPLES/2 + 1];  // 存储有效幅度谱

3. 初始化RFFT实例 #

这是最容易漏掉的关键步骤，必须先初始化快速实FFT的实例：

arm_rfft_fast_instance_f32 rfft_instance;
arm_status status;

// 初始化RFFT实例，传入实例指针和采样长度
status = arm_rfft_fast_init_f32(&rfft_instance, TEST_LENGTH_SAMPLES);
if (status != ARM_MATH_SUCCESS) {
    // 初始化失败，可添加错误处理逻辑
    while(1);
}

4. 执行实数据FFT #

调用arm_rfft_fast_f32处理实数据，最后一个参数0代表正向FFT（时域转频域）：

// 正向实FFT：输入实数据，输出复数格式的频谱
arm_rfft_fast_f32(&rfft_instance, ffttestrealip, fft_output, 0);

5. 计算幅度谱（可选但常用） #

输出的复数频谱需要转换为幅度值才方便分析，由于实数据FFT的对称性，只有前N/2 + 1个点是有效的：

// 计算幅度谱，只取有效点
arm_cmplx_mag_f32(fft_output, fft_magnitude, TEST_LENGTH_SAMPLES/2 + 1);

常见坑点提醒 #

• 采样长度必须是2的幂：CMSIS-DSP的快速实FFT仅支持2^n长度的输入，否则初始化或FFT执行会失败。
• 输出数组长度要求：arm_rfft_fast_f32的输出数组长度需和输入一致，因为它会把实数据转换为复数格式（每两个float对应一个复数的实部和虚部）存储。
• 正向/反向FFT区分：函数最后一个参数，0是正向FFT（时域转频域），1是反向FFT（频域转时域）。

你可以给ffttestrealip填充已知的实信号（比如正弦波），然后查看fft_magnitude是否出现预期的峰值，以此验证FFT是否正常工作。

内容的提问来源于stack exchange，提问作者user3147192