DCL:MATH2:FFTLIB : 高速フーリエ変換:メソッド(元サブルーチン)の説明
周期実数値データのフーリエ変換をおこなう. rfftiは初期化をおこなう; rfftfはフーリエ順変換をおこなう; rfftbはフーリエ逆変換をおこなう.
Nが偶数のとき N' = N/2-1, Nが奇数のとき N'' = (N-1)/2とおく.
順変換は次のように定義される.
ただしNが偶数のとき,
逆変換は次のように定義される.
Nが偶数のとき,
Nが奇数のとき,
wsave = NumRu::DCL.rffti(n)
NumRu::DCL.rfftf(r,wsave)
NumRu::DCL.rfftb(r,wsave)
n (I) 処理するデータの長さ. wsave (R) 作業用配列. 長さは少なくとも 2 n+15以上でなければならない. r (R) 処理する実数型配列. 入力パラメータでもあり出力パラメータでもある. (上記定義参照).