e*r 发帖数: 1121 | 1 【 以下文字转载自 Programming 讨论区 】
发信人: ear (猪头洗心革面), 信区: Programming
标 题: numerical recipe里的快速傅立叶变换
发信站: BBS 未名空间站 (Tue Jul 1 22:39:16 2014, 美东)
最近需要编个cross correlation 的程序, 调用了numerical recipe里的four1 (快
速傅立叶变换)。用简单的正弦函数做测试,通过。上真实数据,发现傅立叶变换后返
回很多NAN的value。数据点也不是很变态,大概是0.01 到3之间。网上google有些人说
这个four1毛病很多。我不是很明白。
这里有懂的可以讲讲么?另外可以推荐个比较robust 的快速傅立叶变换子程序吗?非
常感谢!我用的是fortran.
xiexie! |