C++单片机如何使用CMSIS-DSP计算1024点实数FFT并寻找频域基波和谐波-C++-PHP中文网

arm_rfft_fast_init_f32必须在调用arm_rfft_fast_f32前初始化实例，否则因未初始化指针导致乱码或硬故障；需静态分配结构体、传入2的幂点数（如1024），并检查返回状态。

c++单片机如何使用cmsis-dsp计算1024点实数fft并寻找频域基波和谐波

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++14

C++14 对 C++11 的修正与增强版本，适合旧系统维护和较老工具链兼容。

下载

有效频点只有前 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])，别漏了平方和开方