生产中应优先使用FFTW等成熟FFT库而非手写；实数输入须用fftw_plan_dft_real_1d并配合fftw_malloc；std::vector不能直接传入复数plan，需校验ABI兼容性；模长和相位计算应使用std::abs与std::arg并注意归一化。
用 std::complex 配合现成 FFT 库最实际
纯手写 Cooley-Tukey FFT 在 C++ 里既容易出错又难优化，除非你真在写教学示例或嵌入式受限环境。生产代码中直接调用成熟库是唯一合理选择。主流选项就三个：FFTW（精度高、C 接口但 C++ 可用）、arm_math.h（ARM Cortex 专用）、xtensor-blas 或 boost::compute（偏高层封装）。其中 FFTW 兼容性最好，x86/ARM/Linux/macOS/Windows 全支持。
fftw_plan_dft_1d 怎么配实数输入和复数输出
多数信号处理场景输入是实数数组（比如 std::vector<double></double>），但 FFTW 的复数变换要求输入也按复数布局。别直接传 double* 给复数 plan——会崩溃或结果全零。
- 用
fftw_plan_dft_real_1d：它专为实数输入设计，自动利用共轭对称性，输出长度为n/2 + 1个fftw_complex（即double[2]） - 输入数组必须用
fftw_malloc分配，不能用new或std::vector::data()—— 否则可能因内存未对齐导致性能暴跌（慢 3–5 倍）或 segfault - 输出的
fftw_complex*第一个元素是 DC 分量，最后一个（当n为偶数时）是 Nyquist 频点，中间是正频率分量；负频率不显式存，靠对称性隐含
为什么 std::vector<:complex>></:complex> 不能直接喂给 fftw_plan_dft_1d
std::complex<double></double> 内存布局虽通常与 fftw_complex 一致（两个连续 double），但 C++ 标准不保证其 ABI 兼容。某些编译器（如 MSVC 旧版）或开启强优化时会破坏字段顺序。
- 安全做法：用
reinterpret_cast<fftw_complex>(vec.data())</fftw_complex>前，先断言sizeof(std::complex<double>) == sizeof(fftw_complex)</double>且alignof(std::complex<double>) == alignof(fftw_complex)</double> - 更稳妥的是绕过
std::complex，直接用fftw_complex*指针操作，或用std::vector<double></double>存实部虚部交错数据（[re0, im0, re1, im1, ...]） - 若用 C++17 以上，可借助
std::span<fftw_complex></fftw_complex>封装，避免裸指针
FFT 结果怎么取模长和相位才不出错
FFT 输出是复数，但直接对 fftw_complex[i][0]（实部）和 [i][1]（虚部）调 std::sqrt 容易溢出或精度丢失，尤其高频噪声大时。
立即学习“C++免费学习笔记（深入）”；
- 用
std::abs(std::complex<double>{re, im})</double>算模长，它内部做了防溢出缩放 - 相位统一用
std::arg(std::complex<double>{re, im})</double>，别手写atan2(im, re)—— 虽等价，但std::arg对re=0 && im=0有明确定义（返回 0），而atan2(0,0)是 domain error - 注意：FFT 输出未归一化，模长要除以
N（点数）才对应真实幅度谱；有些库默认不归一化，有些归一化，得查文档或自己验证
真正麻烦的从来不是调哪个函数，而是内存对齐、数据布局、归一化约定这三处——漏查任意一个，结果看起来“差不多”，实则频谱幅度差十倍、相位跳变、或者某平台跑着跑着就崩了。
