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! | l******n
发帖数: 9344 | 2 try fftw
【在 e*r 的大作中提到】
: 【 以下文字转载自 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.
| c****t
发帖数: 19049 | 3 numerical recipe是个软件包？哈哈我孤陋寡闻了。试试IMSL吧或者matlab也行
【在 e*r 的大作中提到】
: 【 以下文字转载自 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.
|
|
