戻る|進む

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

2.3.3 sinti/sint

1.

機能

奇の周期データのSINE変換をおこなう. sintiは初期化をおこなう; sintはSINE変換をおこなう.

2.

定義

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

$$X_{k} = 2\sum_{i=1}^{N}x_{i}\sin \frac{\pi ik}{N+1},\mbox{\hspace{1em}}( k = 1, \ldots, N ).$$

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

3.

呼び出し方法

wsave = NumRu::DCL.sinti(n)
NumRu::DCL.sint(x,wsave)

4.

パラメーターの説明

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

5.

備考

(a)

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