计算结果不同

代人发帖

为什么MFFFT(x,y,n,0)的计算结果与标准FFT结果显著不同？

n:=16; //16点信号变换
VARIABLE : x[n]=(15, 32, 9, 222, 118, 151, 5, 7, 56, 233, 56, 121, 235, 89, 98, 111);
VARIABLE : y[n*2]=0;  //以复数方式返回结果，因此需要2倍容量
If IsLastbar Then
Begin
        MFFFT(x,y,n,0);
        for i = 1 to n*2 do
        msgout(1,numtostr(y,3));
End


以下为软件输出



这是标准输出

技术009

不知道你的计算过程是怎样的，但是通过和使用py模块执行的计算结果对照：

(1558+0j)
(-83.16474885677513+182.1210324381293j)
(-423.4213562373095-138.7766952966369j)
(-146.6785066419351+326.1556062575687j)
(256-44j)
(-0.3753172272272991+362.1657075253354j)
(-140.5786437626905-214.7766952966369j)
(66.21857272593752-249.868866294104j)
(-374+0j)
(66.21857272593752+249.868866294104j)
(-140.5786437626905+214.7766952966369j)
(-0.3753172272272991-362.1657075253354j)
(256+44j)
(-146.6785066419351-326.1556062575687j)
(-423.4213562373095+138.7766952966369j)
(-83.16474885677513-182.1210324381293j)

至少是保持了一致性的。


相关py代码：

[Python] 复制代码

import numpy as np


data_list = [15, 32, 9, 222, 118, 151, 5, 7, 56, 233, 56, 121, 235, 89, 98, 111]  

# 对列表做快速傅里叶变换
fft_result = np.fft.fft(data_list)


for i in fft_result:
    print(i)