◆ Rfft1b()

1次元実フーリエ逆変換

目的: 本ルーチンは実数配列中の周期数列の1次元フーリエ変換を計算する. この変換はフーリエ逆変換あるいはフーリエ合成と呼ばれ, 数列をスペクトル空間から物理空間に変換する.
Nが偶数の場合:
R(K) = R(0) + (-1)^k R(N-1) + ΣR(2j-1)Cos(2πjk/N) + ΣR(2j)Sin(2πjk/N) (Σは j = 1 〜 N/2-1) (k = 0 〜 N-1)

Nが奇数の場合:
R(K) = R(0) + ΣR(2j-1)Cos(2πjk/N) + ΣR(2j)Sin(2πjk/N) (Σは j = 1 〜 (N-1)/2) (k = 0 〜 N-1)

この変換は正規化されており, Rfft1fに続くRfft1bの呼び出し(あるいはその逆)により, アルゴリズム上の制約, 丸め誤差などを除き, 元の配列を復元する.

引数

[in]	N	入力データ列の長さ. (N >= 1) (Nが小さな素数の積で表されると効率が良い)
[in,out]	R()	配列 R(LR - 1) (LR >= Inc*(N - 1) + 1) [in] 入力データ列. [out] フーリエ逆変換されたデータ列.
[in]	Wsave()	配列 Wsave(LWsave - 1) (LWsave >= N + ln(N)/ln(2) + 4) 作業データ. 入力データ列の長さNごとに, Rfft1fあるいはRfft1bを最初に呼び出す前にRfft1iにより初期化しておかなければならない.
[out]	Info	= 0: 正常終了. = -1: パラメータ N の誤り. (N < 1) = -2: パラメータ R() の誤り. (配列R()の大きさが不足) = -3: パラメータ Wsave() の誤り. (配列Wsave()の大きさが不足) = -5: パラメータ Inc の誤り. (Inc < 1)
[in]	Inc	(省略可) データ列の連続する要素の配列R()内でのインデックスの間隔. (Inc >= 1) (省略時 = 1)