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.
呼び出し方法
wsave = NumRu::DCL.ezffti(n)
azero, a, b = NumRu::DCL.ezfftf(r,wsave)
r = NumRu::DCL.ezfftb(azero,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を続けて呼ぶと, もとの 値が返される.