戻る|進む

DCL:MATH２:FFTLIB : 高速フーリエ変換:サブルーチンの説明

2.3.4 COSTI/COST

1.

機能

偶の周期データのCOSINE変換をおこなう. COSTIは初期化をおこなう; COSTはCOSINE変換をおこなう.

2.

定義

順変換は次のように定義される.

　$$X_{k} = x_{1} + (-1)^{k-1}x_{N}+ 2\sum_{i=2}^{N}x_{i}\cos \...... (i-1)(k-1)}{N-1},\mbox{\hspace{1em}}( k = 1, \ldots, N ).$$

逆変換は順変換と同じである.

3.

呼び出し方法

CALL COSTI(N,WSAVE)
CALL COST(N,R,WSAVE)

4.

パラメーターの説明

N	(I)	処理するデータの長さ.
WSAVE	(R)	作業用配列. 長さは少なくとも 3 N+15以上でなければならない.
R	(R)	処理する実数型配列. 入力パラメータでもあり出力パラメータでもある (上記定義参照).

5.

備考

(a)

COSTは逆変換でもある. また, この変換では正規化されない. つまりCOSTを2回続けて呼ぶと, もとの 2(N-1)倍の値が返される.