快速傅里叶变换（FFT）

最新推荐文章于 2025-02-21 08:22:13 发布

JerryDung

最新推荐文章于 2025-02-21 08:22:13 发布

阅读量1.3k

点赞数 1

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/JerryDung/article/details/8884782

本文介绍了一种基于复数的快速傅立叶变换（FFT）算法实现，该算法使用内存池进行临时数组分配，并通过递归的方式实现了复数序列的FFT计算。文中详细展示了核心函数的伪代码，包括复数乘法、加法运算等关键步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

贴过来做个模板。

C为复数类，newArr和at为内存池。

void FFT(int n, C *x)
{
    C W = C(cos(2 * M_PI / n), sin(2 * M_PI / n)), w = C(1), 
        *u = newArr(n >> 1), *v = newArr(n >> 1);
    for (int i = 0; i < n >> 1; ++i)
        u[i] = x[i << 1], v[i] = x[i << 1 | 1];
    if (n > 1) FFT(n >> 1, u), FFT(n >> 1, v);
    for (int i = 0; i < n >> 1; ++i, w *= W)
        x[i] = u[i] + w * v[i], x[i + (n >> 1)] = u[i] - w * v[i];
    at -= n + 5;
}