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

最新推荐文章于 2024-03-19 09:17:03 发布

xueyan0096

最新推荐文章于 2024-03-19 09:17:03 发布

阅读量1.6k

点赞数 1

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

版权

本文详细介绍了Cooley-Tukey FFT算法的递归版本，包括当输入向量长度为2的k次幂时的处理方式和非2的k次幂长度的情况。在向量长度为2^k时，通过参考Wiki提供的代码，阐述了DFT和IDFT的计算，并指出两者仅在指数符号和IDFT额外的1/N系数上的区别。对于非2^k长度的输入，讨论了如何计算卷积，并给出了最终的DFT和IDFT实现。

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

DFT

IDFT

FFT

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

参考wiki FFT有如下代码

/// <summary>
/// Cooley–Tukey FFT algorithm,
/// length of <paramref name="data"/> is power of 2.
/// <paramref name="left"/> is inclusive, <paramref name="right"/> is exclusive.
/// </summary>
private static void CooleyTukey(Complex[] data, int left, int right, int exponentSign)
{
    if (left + 1 < right)
    {
        Separate(data, left, right);
        var half = (right - left) >> 1;
        CooleyTukey(data, left, left + half, exponentSign);   // recurse even items
        CooleyTukey(data, left + half, right, exponentSign);  // recurse odd  items                                               
        for (var k = 0; k < half; k++)// combine results of two half recursions
        {
            var even = data[left + k];   // even
            var odd = data[left + k + half];   // odd                                                      
            var t = 2 * Math.PI * exponentSign * k / (right - left);
            var tmp = new Complex(Math.Cos(t), Math.Sin(t)) * odd; // "twiddle-factor" w*odd