快速傅里叶变换算法因单个减号而出错

Question

所以看完这个关于快速傅立叶变换的视频https://www.youtube.com/watch?v=h7apO7q16V0

我分析了伪代码并在 python 中实现了它，发现它产生的输出与许多 fft 计算器网站的输出不同。我的价值观似乎都在那里，只是奇怪，因为订单不合适，任何人都知道为什么。它是一种不同的算法实现还是什么。

import cmath
import math
def FFT(P):
    n= len(P)

    if n == 1:
        return P

    omega = cmath.exp((2 * cmath.pi * 1j)/n)

    p_even = P[::2]
    p_odd = P[1::2]

    y_even = FFT(p_even)
    y_odd = FFT(p_odd)

    y = [0] * n

    
    for i in range(n//2):
        y[i] = y_even[i] + omega**i*y_odd[i]
        y[i+n//2] = y_even[i] - omega**i*y_odd[i]
    return y

    
poly = [0,1,2,3]
print(FFT([0,1,2,3]))

我测试它的网站是https://tonysader.github.io/FFT_Calculator/？ 我在这个站点输入 0,1,2,3 得到：6, -2+2J, -2, -2+-2J

而我的python程序输出：6，-2-2J，-2，-2+2J

我遵循的伪代码：

Answer 1

我认为您正在运行的程序正在执行逆 FFT。尝试 omega = cmath.exp((-2 * cmath.pi * 1j)/n)。注意减号。