<< Signal Processing Supplementary toolbox Signal Processing Supplementary toolbox cheb1ord >>

Signal Processing Supplementary toolbox >> Signal Processing Supplementary toolbox > cceps

cceps

Complex cepstrum

Calling Sequence

Xh = cceps(X)
[Xh, Nd] = cceps(X)
[Xh, Nd, Yh] = cceps(X)
[Xh, Nd, Yh] = cceps(X,N)

Parameters

X

Real input sequence.

N

Cepstrum length.

Xh

Complex cepstrum of real sequence X.

Nd

Number of samples Nd of (circular) delay added to X prior to finding the complex cepstrum.

Yh

Complex cepstrum, which can be used to verify Xh for short sequences that can be rooted and do not have zeros on the unit circle.

Description

Computes the complex cepstrum (Xh) of a real sequence X. The desired length of the cepstrum can also be specified in the second argument N. The values of cepstrum are different for different lengths of the cepstrum specified. However, the values of Yh are independent on the values of the specified cepstrum length. The number of samples delay Nd introduced to the real sequence is required for inversion of complex cepstrum in correct sequence.

Note: Log(0) is undefined. Scilab generates error for log(0). To handle this problem zero values in the input of log function are replaced by very small value %eps. In case of Matlab, log(0) gives -inf.

Examples

// real sequence x
X=[2 5 3 1 0 4 7 6 3 5  8 3 ];
// complex cepstrum 
[Xh,nd]=cceps(X)

Authors


Report an issue
<< Signal Processing Supplementary toolbox Signal Processing Supplementary toolbox cheb1ord >>