【笔记】Cooley–Tukey FFT Algorithm - Iterative Edition

最新推荐文章于 2023-05-22 15:59:40 发布

xueyan0096

最新推荐文章于 2023-05-22 15:59:40 发布

阅读量865

点赞数

本文链接：https://blog.csdn.net/xueyan0096/article/details/84857152

版权

本文深入探讨了Cooley-Tukey快速傅里叶变换(FFT)算法的迭代实现，包括输入向量长度为2的幂次时的情况和非2的幂次时的处理。详细介绍了DFT、IDFT的计算，并阐述了如何处理不同长度的输入向量，特别是通过卷积的关系来说明IDFT的计算方法。同时，提到了IDFT过程中涉及的虚实交换操作。

摘要由CSDN通过智能技术生成

DFT

IDFT

FFT

1 当输入向量长度为2^k时

参考《Introduction to Algorithms》chapter 30.3

每次迭代的元素赋值及辅助方法reorder

/// <summary>
/// Radix-2 Step Helper Method
/// </summary>
/// <param name="data">Sample vector.</param>
/// <param name="exponentSign">Fourier series exponent sign.</param>
/// <param name="levelSize">Level Group Size.</param>
/// <param name="k">Index inside of the level.</param>
private static void Radix2Step(Complex[] data, int exponentSign, int levelSize, int k)
{
    // Twiddle Factor
    var exponent = (exponentSign * k) * Math.PI / levelSize;
    var w = new Complex(Math.Cos(exponent), Math.Sin(exponent));

    var step = levelSize << 1;
    for (var i = k; i < data.Length; i += step)
    {
        var ai = data[i];
        var t = w