2.3.2 EZFFTI/EZFFTF/EZFFTB

1.
機能
RFFTI, RFFTF, RFFTBの簡易型サブルーチン. EZFFTIは初期化をおこなう; EZFFTFはフーリエ順変換をおこなう; EZFFTBはフーリエ逆変換をおこなう.
2.
定義
順変換は次のように定義される. (Nが偶数のとき N'' = N/2-1, Nが奇数のとき N'' = (N-1)/2とおく.)


ただしNが偶数のとき,


逆変換は次のように定義される. (Nが偶数のときN'' = N/2, Nが奇数のとき N'' =(N-1)/2 とおく.)


\begin{displaymath}R_{i} = a_{0} + \sum_{k=1}^{N'}( a_{k} \cos \frac{2\pi (i-1......c{2\pi (i-1)k}{N})\mbox{\hspace{1em}}( i = 1, \ldots, N ). \end{displaymath}

3.
呼び出し方法
CALL EZFFTI(N,WSAVE)
CALL EZFFTF(N,R,A0,A,B,WSAVE)
CALL EZFFTB(N,R,A0,A,B,WSAVE)
4.
パラメーターの説明
N (I) 処理するデータの長さ.
WSAVE (R) 作業用配列. 長さは少なくとも 3 N+15以上でなければならない.
R (R) 処理する実数型配列. EZFFTFにおいては入力パラメータ, EZFFTBにおいては出力パラメータである.
A0 (R) 上記定義におけるA0および a0.
A, B (R) Nが偶数のとき , N/2,N が奇数のとき (N-1)/2 の長さの実数型配列(上記定義参照).
5.
備考
(a)
この変換では正規化される. つまりEZFFTF, EZFFTBを続けて呼ぶと, もとの 値が返される.