EZFFTI/EZFFTF/EZFFTB

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) 上記定義における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` の長さの実数型配列(上記定義参照).