rceps (Signal Processing Toolbox)

実数ケプストラムと最小位相復元

表示

y = rceps(x)
[y,ym] = rceps(x)

詳細

実数ケプストラムは、データ列のフーリエ変換の大きさの自然対数を逆フーリエ変換したものです。

rceps(x) は、実数列xの実数ケプストラムを出力します。実数ケプストラムは、実数値関数です。

[y,ym] = rceps(x) は、入力列の実数ケプストラムyと最小位相復元バージョンymの両方を出力します。

アルゴリズム

rcepsは、参考文献[2]のアルゴリズム7.2をM-ファイルで実現したものです。すなわち、つぎの操作を行います。

y = real(ifft(log(abs(fft(x)))));

ケプストラム領域で適切なウィンドウを適用すると、復元された最小位相信号が作成されます。.

w = [1; 2*ones(n/2-1,1); ones(1 - rem(n,2),1); zeros(n/2-1,1)];
ym = real(ifft(exp(fft(w.*y))));

参考

cceps
複素ケプストラム解析

fft
1次元高速フーリエ変換

hilbert
Hilbert変換

icceps
逆複素ケプストラム

unwrap
位相角の連続性

参考文献

[1] Oppenheim, A.V., and R.W. Schafer, Digital Signal Processing, Englewood Cliffs, NJ, Prentice-Hall, 1975.

[2] Programs for Digital Signal Processing, IEEE Press, New York, 1979.

rc2poly rectpuls