快速傅里叶变换(FFT)算法C++实现代码

最新推荐文章于 2024-08-07 06:47:39 发布

simplecoding

最新推荐文章于 2024-08-07 06:47:39 发布

阅读量1w

点赞数

文章标签： fft c++ 算法存储 math.h include

本文链接：https://blog.csdn.net/simplecoding/article/details/1665324

版权

本文介绍了如何用C++实现快速傅里叶变换(FFT)算法，包括正向和逆向变换。代码中详细展示了Danielson - Lanczos引理的应用，并解释了输入数据和输出数据的存储方式以及频率范围。

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

#include <math.h>
#define DOUBLE_PI 6.283185307179586476925286766559

// 快速傅里叶变换
// data 长度为 (2 * 2^n), data 的偶位为实数部分, data 的奇位为虚数部分
// isInverse表示是否为逆变换
void FFT(double * data, int n, bool isInverse = false)
{
    int mmax, m, j, step, i;
    double temp;
    double theta, sin_htheta, sin_theta, pwr, wr, wi, tempr, tempi;
    n = 2 * (1 << n);
    int nn = n >> 1;
    // 长度为1的傅里叶变换, 位置交换过程
    j = 1;
    for(i = 1; i < n; i += 2)
    {
        if(j > i)
        {
            temp = data[j - 1];
            data[j - 1] = data[i - 1];
            data[i - 1] = temp;
            data[j] = temp;
            data[j] = data[i];
            data[i] = temp;
        }
        // 相反的二进制加法
        m = nn;
        while(m >= 2 && j > m)<