RCEPS

Purpose:

Calculates the real cepstrum.

Syntax:

RCEPS(series, n)

series	-	An input series or array.
n	-	Optional. An integer, the number of samples to use. If n > length(series), the series is zero padded. Defaults to length(series).

Returns:

A real series or array.

Example:

W1: gtri(100, 1, 1/100)^3

W2: W1-delay(W1, 60)/2

W3: rceps(W1)

W4: rceps(W1, 512)

W2 adds a synthesized echo at 60 seconds to the data of W1.

W3 displays a small peak at t == 60 indicating the presence of the echo. W4 performs the same calculation with the data padded to 512 samples.

Remarks:

The complex cepstrum of a series is essentially ifft(log(fft(s))). However, the complex log calculation requires the evaluation of the continuous phase component. RCEPS ignores the phase component and calculates: