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