◆ rfft1b()

1次元実フーリエ逆変換

目的: rfft1bは実数配列中の周期数列の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の呼び出し(あるいはその逆)により, アルゴリズム上の制約, 丸め誤差などを除き, 元の配列を復元する.

戻り値: info (int)
= 0: 正常終了
= -1: 入力パラメータ n の誤り (n < 1)
= -2: 入力パラメータ r の誤り
= -3: 入力パラメータ wsave の誤り
= -4: 入力パラメータ inc の誤り (inc < 1)

引数

[in]	n	入力データ列の長さ. (n >= 1) nが小さな素数の積で表されると効率が良い.
[in,out]	r	Numpy ndarray (1次元配列, float, 長さ inc*(n - 1) + 1) [in] 入力データ列. [out] フーリエ逆変換されたデータ列.
[in]	wsave	Numpy ndarray (1次元配列, float, 長さ n + ln(n)/ln(2) + 4) 作業データ. 入力データ列の長さnごとに, rfft1fあるいはrfft1bを最初に呼び出す前にrfft1iにより初期化しておかなければならない.
[in]	inc	(省略可) データ列の連続する要素の配列r内でのインデックスの間隔. (inc >= 1) (省略時 = 1)