fft と逆fft によるもとの時系列の作成

もとのデータ

これを fftpして、また逆fft する ソースはここ fft に使ったのは、NETLIB の fftpack を f2c したものから、 (1/n) の因子で悩まなくて良い ezfft の系列。

fft と逆 fft により作られた数学的には元のデータと同じものから、 元のデータを引き、差を見た。データを X 軸に、差を Y軸に書くと

となる。もとのデータに対して、エラーは 10の -5 乗くらいであり、 計算誤差として許容できる程度である。

2000.7.17