ezffti/ezfftf/ezfftb

DCL:MATH２:FFTLIB : 高速フーリエ変換:メソッド(元サブルーチン)の説明

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) 上記定義におけるA₀および a₀.

a, b (R) nが偶数のとき , n/2,n が奇数のとき (N-1)/2 の長さの実数型配列(上記定義参照).

5.

備考

(a): この変換では正規化される. つまりezfftf, ezfftbを続けて呼ぶと, もとの値が返される.

`n`	`(I)`	処理するデータの長さ.
`wsave`	`(R)`	作業用配列. 長さは少なくとも 3 `n`+15以上でなければならない.
`r`	`(R)`	処理する実数型配列. `ezfftf`においては入力パラメータ, `ezfftb`においては出力パラメータである.
`a0`	`(R)`	上記定義におけるA₀および a₀.
`a`, `b`	`(R)`	`n`が偶数のとき , `n/2,n` が奇数のとき (`N-1)/2` の長さの実数型配列(上記定義参照).