主题:跪求二维快速傅里叶变换程序
方正二两酒
[专家分:0] 发布于 2012-12-02 21:41:00
由于我的ivf怎么也调用不了imsl, 现在只能自己编写一个二维快速傅里叶变换的fortran程序,由于对快速傅里叶变换的算法不怎么懂,现在急着要出结果,求大家行行好给我一个二维快速傅里叶变换的fortran 程序!
最后更新于:2012-12-02 21:53:00
回复列表 (共2个回复)
沙发
方正二两酒 [专家分:0] 发布于 2012-12-02 21:51:00
这个是我自己按照那个numerical recipe 上的算法编的,完全不对,不知道为什么。附件是我编的程序和算法。 请高手帮我看看吧,真的很急啊!
板凳
rockymt [专家分:110] 发布于 2012-12-04 10:29:00
伙计,试试FFTPACK5.1
单精度http://people.sc.fsu.edu/~jburkardt/f_src/fftpack5.1/fftpack5.1.html
双精度http://people.sc.fsu.edu/~jburkardt/f_src/fftpack5.1d/fftpack5.1d.html
下载fftpack5.1.f90或者fftpack5.1d.f90,网页上还有例子
估计你要用到的函数是CFFT2I CFFT2F CFFT2B, 这三个函数处理复数
要么RFFT2I RFFT2F RFFT2B,处理实数
自个好好看看说明和例子
我来回复