我应该为 vDSP 中的真实 FFT 函数提供什么？

Question

vDSP_fft_zrop 的函数签名（实际上是 vDSP 中的大多数 FFT 函数）接受两个 DSPSplitComplex 结构作为输入和输出。然而，对于真正的 FFT 计算，前向 FFT 只接受实数值的输入，而后向 FFT 只产生实数值的输出。那么我应该如何为全实值填充输入/输出DSPSplitComplex 结构？

为了更清楚地声明，如果我要在下面实现一个更常见的 API：

void forward_fft( std::complex<float>* output, const float* input, unsigned size);

很明显，我应该使用vDSP_ztoc 将vDSP 函数的输出交错到output 内存，但是我应该如何处理input 内存？我应该这样做：

DSPSplitComplex input_for_vdsp;
input_for_vdsp.realp = input;
input_for_vdsp.imagp = nullptr

，或者我应该为imagp指针分配内存？

Answer 1

您可以将实数值表示为使用复数值的实部和虚部的奇偶配置。完整讨论请参阅 Apple 的 Data Packing for Fourier Transforms 文章。