Fortran Coder 程序员聚集地

快速傅里叶变换FFT [我来说两句]

评论列表（评论 8）以下网友评论只代表网友个人观点，不代表本站观点。

2021-12-05 19:03:51 njuxyc(njuxyc)

程序没问题，和mathematica的计算结果一致

2019-10-12 09:35:16 sysslm(sysslm)

Spineless 于 2015-12-14 11:11:32发布

这个代码有误，用序列（1,2,3,4）计算，算出来的结果是[（10,0），（-2,-2），（-2,0），（-2,2）]，但是正确答案应该是[（10,0），（-2,2），（-2,0），（-2,-2）]

第四行FFT_Forward和第五行FFT_Inverse的取值调换一下就对了

回复支持（1）

2019-07-21 13:45:02 aoeo2jam(aoeo2jam)

乘除法次数太多，居然还有指数运算，计算速度不行，如果作为核心程序大量调用，效率是不行的。

回复支持（5）

2019-01-07 09:25:23 hq1162377655(hq1162377655)

你好，请问有维度是奇数的FFT代码？谢谢！

回复支持（7）

2017-09-20 13:49:45 168371286(168371286)

糖盒love玲珑于 2016-03-16 10:24:49发布

这段代码 没有错误!!! 变换的时候,数学公式上有个i(虚数单位), 其实也可以用-i,只要正向的变换和返回的时候一致就可以了. 要是想和Matlab中计算的一致, 你自己修改 x(itmp) = ctmp * exp( forback * cmplx( 0.0_DP , e*(j-1),kind=DP ) ) 为 x(itmp) = ctmp * exp( forback * cmplx( 0.0_DP , -e*(j-1),kind=DP ) ) 即可!

请问，二维数组（矩阵）离散傅里叶变换怎么实现？

回复支持（6）

2016-03-16 10:24:49 糖盒(糖盒love玲珑)

这段代码没有错误!!! 变换的时候,数学公式上有个i(虚数单位), 其实也可以用-i,只要正向的变换和返回的时候一致就可以了. 要是想和Matlab中计算的一致, 你自己修改 x(itmp) = ctmp * exp( forback * cmplx( 0.0_DP , e*(j-1),kind=DP ) ) 为 x(itmp) = ctmp * exp( forback * cmplx( 0.0_DP , -e*(j-1),kind=DP ) ) 即可!

回复支持（3）

2015-12-17 08:35:09 vvt(vvt)

Spineless 于 2015-12-14 11:11:32发布

这个代码有误，用序列（1,2,3,4）计算，算出来的结果是[（10,0），（-2,-2），（-2,0），（-2,2）]，但是正确答案应该是[（10,0），（-2,2），（-2,0），（-2,-2）]

经验证，此代码与matlab计算结果一致。为：（10,0），（-2,-2），（-2,0），（-2,2）

回复支持（4）

2015-12-14 11:11:32 Spineless(Spineless)

这个代码有误，用序列（1,2,3,4）计算，算出来的结果是[（10,0），（-2,-2），（-2,0），（-2,2）]，但是正确答案应该是[（10,0），（-2,2），（-2,0），（-2,-2）]