arm_rfft_fast_init_f32必须在调用arm_rfft_fast_f32前初始化实例，否则因未初始化指针导致乱码或硬故障；需静态分配结构体、传入2的幂点数（如1024），并检查返回状态。
CMSIS-DSP 的 arm_rfft_fast_init_f32 必须在调用前初始化实例
单片机上跑 CMSIS-DSP 的实数 FFT，最容易卡住的地方不是算法逻辑，而是忘了初始化 RFFT 实例。CMSIS-DSP 不像 PC 端 FFT 库那样“传入数组就开算”，它要求显式初始化一个 arm_rfft_fast_instance_f32 结构体——否则 arm_rfft_fast_f32 会读取未初始化的内部指针，结果全乱或直接硬故障。
正确做法是静态分配并初始化：
arm_rfft_fast_instance_f32 rfft_inst;
arm_status status = arm_rfft_fast_init_f32(&rfft_inst, 1024);
if (status != ARM_MATH_SUCCESS) {
// 初始化失败，常见原因：未链接 CMSIS-DSP 库或编译时没开浮点支持
}
- 必须用
arm_rfft_fast_init_f32（不是arm_rfft_init_f32），后者是旧接口，不支持自动缩放且输出格式不同 - 1024 是点数，必须是 2 的幂；CMSIS-DSP 不检查非法值，传 1000 会静默越界
- 初始化只需做一次，在
main()开头或模块初始化函数里完成即可
输入数据要补零对齐、输出频谱是复数交错格式
1024 点实数 FFT 输入是 1024 个 float32_t，但 CMSIS-DSP 要求输入缓冲区长度为 1024（不是 512 或 2048），且数据需是实数序列——不能是复数，也不能是未归一化的 ADC 原始值（比如 12-bit ADC 输出 0–4095，得先转成 -1.0 ~ +1.0 范围）。
输出缓冲区也必须是 1024 个 float32_t，但布局是“复数交错”：out[0] 是直流分量实部，out[1] 是直流虚部（恒为 0），out[2] 是 f1 实部，out[3] 是 f1 虚部……直到 out[1022] 和 out[1023] 是 Nyquist 频率分量（实部+虚部，但虚部恒为 0）。
立即学习“C++免费学习笔记（深入）”；
- 有效频点只有前 513 个复数（索引 0 到 512），对应 0 Hz 到 fs/2 Hz；后半部分是共轭冗余，不用处理
- 基波频率 =
fs / 1024，假设采样率fs = 10.24 kHz，则基波是 10 Hz，第 k 个谐波在索引k（k=1 对应 10 Hz，k=2 对应 20 Hz……） - 计算幅值用
sqrtf(out[2*k]*out[2*k] + out[2*k+1]*out[2*k+1])，别漏了平方和开方
FFT 输出幅度受缩放影响，需手动补偿
CMSIS-DSP 的 arm_rfft_fast_f32 默认做 1/N 缩放（N=1024），即输出幅值是理论值的 1/1024。这意味着：一个满幅正弦波（幅值 1.0）经 1024 点 FFT 后，其基波 bin 幅值约为 0.5，而不是理论上的 512 —— 因为缩放已生效。
如果你需要原始能量值（比如做信噪比或与理论模型对比），得反向乘回 1024；如果只找最大值位置（如锁相、频率检测），缩放不影响相对大小，可跳过。
- 缩放行为由初始化时的
arm_rfft_fast_init_f32决定，不可 runtime 修改 - 若想关掉缩放，得改用底层函数
arm_rfft_f32+ 手动管理 twiddle 表，复杂度陡增，一般没必要 - 注意：ADC 采集后的归一化系数（如 ×2/4095）和 FFT 缩放是两层独立因子，别混在一起补偿
在资源受限单片机上查基波和谐波要避开动态内存和浮点除法
STM32F4/F7/H7 类 MCU 虽有 FPU，但 sqrtf 和 atan2f 仍吃周期；而实时系统中频繁调用可能挤占控制环时间。更糟的是，有些老型号（如 F0/F1）连硬件浮点都没有，sqrtf 会链接软件实现，慢到无法接受。
- 替代方案：用幅值平方代替幅值比较（避免开方），只在最终确认目标频点时才算真实幅值
- 基波搜索范围建议限制在 20–200 Hz（对应索引 2–20，假设 fs=10.24 kHz），别扫全 513 点
- 谐波检测建议用“倍频约束”：若基波在索引
k，则二谐波应在2*k±2范围内找峰值，三谐波在3*k±3，避免被噪声误触发 - 慎用
malloc分配 FFT 缓冲区——单片机堆空间小且易碎片化；一律用静态数组或栈分配（确保栈够大，1024×4B = 4 KB）
真正麻烦的从来不是“怎么调用 FFT 函数”，而是如何让幅值稳定、相位可用、且不把中断延迟拖垮。尤其当 ADC 采样和 FFT 计算共享同一个定时器触发时，缓存对齐、DMA 传输边界、以及 __DSB() 内存屏障的插入时机，都得一个个实测验证。
