どうするのが一番良いかイマイチ分からないが、 とりあえず Ooura's Mathematical Software Packages にある 「汎用 FFT (高速 フーリエ/コサイン/サイン 変換) パッケージ」 を使ってみる。
これはデータ数が の羃でないとダメという制限はあるが、 性能の高さには定評がある。
oyabun% tar xzf fft.tgz oyabun% cd fft oyabun% ls oyabun% cd sample1 oyabun% make mathpc% ./test4g data length n=? (must be 2^m) 1024 cdft err= 4.24682e-16 rdft err= 4.24682e-16 ddct err= 4.95724e-16 ddst err= 4.98625e-16 dfct err= 3.86626e-16 dfst err= 5.69098e-16 oyabun% cd ../sample2 oyabun% make oyabun% ./pi_fft4g oyabun% less pi.dat |
例えば test4g.c を Java プログラムに書き換えると良い。
I君はその作業をスタートしたようだが、上手く行ったかな?
FFTPACK を Java にするのは無理か? あるじゃん! jfftpack.tgz -- でも上手くコンパイル出来ない。なぜだろう。