以文本方式查看主题

-  金字塔客服中心 - 专业程序化交易软件提供商  (http://222.73.7.161/bbs/index.asp)
--  公式模型编写问题提交  (http://222.73.7.161/bbs/list.asp?boardid=4)
----  傅里叶变换  (http://222.73.7.161/bbs/dispbbs.asp?boardid=4&id=181821)

--  作者:gerrypan
--  发布时间:2020/8/20 16:35:19
--  傅里叶变换
T:=16; //16点信号变换
VARIABLE : x[T]=(15, 32, 9, 222, 118, 151, 5, 7, 56, 233, 56, 121, 235, 89, 98, 111);
VARIABLE : y[T*2]=0;  //以复数方式返回结果,因此需要2倍容量
If IsLastbar Then
 Begin
        MFFFT(x,y,T,0);
          w:2*3.14/T ;
         ft1:=0;
    
        for j=1 to T do begin
        for i = 0 to T - 1 do begin
     
ft1:=ft1 + y[i * 2]*cos(i*w * j)/T + y[i * 2 + 1]*sin(i*w*j)/T;
end;
msgout(1,NUMTOSTR(ft1 / T,5));
end;
End

结果 :
2020/08/20 16:30:57  5.85846
2020/08/20 16:30:57  10.70716
2020/08/20 16:30:57  21.31674
2020/08/20 16:30:57  24.41963
2020/08/20 16:30:57  32.06571
2020/08/20 16:30:57  37.75488
2020/08/20 16:30:57  41.76019
2020/08/20 16:30:57  48.59118
2020/08/20 16:30:57  52.76111
2020/08/20 16:30:57  58.53273
2020/08/20 16:30:57  66.14874
2020/08/20 16:30:57  69.21599
2020/08/20 16:30:57  79.82165
2020/08/20 16:30:57  84.64967
2020/08/20 16:30:57  90.49294
2020/08/20 16:30:57  97.42217

不太对啊。

--  作者:FireScript
--  发布时间:2020/8/20 17:09:18
--  
 抱歉,这种算法我们也不是很了解的。