Signal Toolkit - rceps

Function File: [y, ym] = rceps (x)

Return the cepstrum of the signal x.

If x is a matrix, return the cepstrum of each column.

If called with two output arguments, the minimum phase reconstruction of the signal x is returned in ym.

For example:

 f0 = 70; Fs = 10000;            # 100 Hz fundamental, 10kHz sampling rate
 a = poly (0.985 * exp (1i * pi * [0.1, -0.1, 0.3, -0.3])); # two formants
 s = 0.005 * randn (1024, 1);    # Noise excitation signal
 s(1:Fs/f0:length(s)) = 1;       # Impulse glottal wave
 x = filter (1, a, s);           # Speech signal
 [y, ym] = rceps (x .* hanning (1024));
 

Reference: Programs for Digital Signal Processing, IEEE Press, John Wiley & Sons, New York, 1979.