Cepstrum

Following abbreviations are used in the formulas to explain the cepstrum in more detail:

Abbreviation	Explanation
${\displaystyle f(t)}$	Signal, which is a function of time
${\displaystyle C}$	Cepstrum
${\displaystyle {\mathcal {F}}}$	Fourier transform: The abbreviation can stand i.e. for a continuous fourier transform, a discrete fourier transform (DFT) or even a z-transform, as the z-transform is a generalization of the DFT.[3]
${\displaystyle {\mathcal {F}}^{-1}}$	Inverse of the fourier transform
${\displaystyle \left|{\mathcal {F}}\{f(t)\}\right|^{2}}$	Power spectrum
${\displaystyle \log(x)}$	Logarithm of x. The choice of the base b depends on the user. In some articles the base is not specified, others prefer base 10 or e. The choice of the base has not impact on the basic calculation rules, but sometimes base e leads to simplifications (see "Complex Cepstrum").
${\displaystyle \left|x\right|}$	Absolute value, or magnitude of a complex value, which is calculated from real- and imaginary part using the Pythagorean theorem.
${\displaystyle \phi }$	Phase angle of a complex value

The "Cepstrum" was originally defined as "Power Cepstrum" by following relationship:[1][3]

${\displaystyle C_{p}=\left|{\mathcal {F}}^{-1}\left\{\log \left(\left|{\mathcal {F}}\{f(t)\}\right|^{2}\right)\right\}\right|^{2}}$

The power cepstrum has main applications in analysis of sound and vibration signals. It is a complementary tool to spectral analysis.[2]

Sometimes it is also defined as:[2]

${\displaystyle C_{p}=\left|{\mathcal {F}}\left\{\log \left(\left|{\mathcal {F}}\{f(t)\}\right|^{2}\right)\right\}\right|^{2}}$

Due to this formula, the cepstrum is also sometimes called the spectrum of a spectrum. It can be shown that both formulas are consistent with each other as the frequency spectral distribution remains the same, the only difference being a scaling factor [2] which can be applied afterwards. Some articles prefer the second formula.[2][4]

Other notations are possible due to the fact, that the log of the power spectrum is equal to the log of the spectrum, if a scaling factor 2 is applied:[5]

${\displaystyle \log |{\mathcal {F}}|^{2}=2\log |{\mathcal {F}}|}$

and therefore:

${\displaystyle C_{p}=\left|{\mathcal {F}}^{-1}\left\{2\log |{\mathcal {F}}|\right\}\right|^{2}}$, or

${\displaystyle C_{p}=4\cdot \left|{\mathcal {F}}^{-1}\left\{\log |{\mathcal {F}}|\right\}\right|^{2}}$, which provides a relationship to the real cepstrum (see below).

Further, it shall be noted, that the final squaring operation in the formula for the power spectrum ${\displaystyle C_{p}}$ is sometimes called unnecessary[3] and therefore sometimes omitted.[4][2]

The complex cepstrum was defined by Oppenheim in his development of homomorphic system theory.[6][7] The formula is provided also in other literature.[2]

${\displaystyle C_{c}={\mathcal {F}}^{-1}\left\{\log({\mathcal {F}}\{f(t)\})\right\}}$

As ${\displaystyle {\mathcal {F}}}$ is complex the log-term can be also written with ${\displaystyle {\mathcal {F}}}$ as a product of magnitude and phase, and subsequently as a sum. Further simplification is obvious, if log is a natural logarithm with base e:

${\displaystyle \log({\mathcal {F}})=\log({\mathcal {|F|\cdot e^{i\phi }}})}$

${\displaystyle \log _{e}({\mathcal {F}})=\log _{e}({\mathcal {|F|}})+\log _{e}(e^{i\phi })=\log _{e}({\mathcal {|F|}})+i\phi }$

Therefore: The complex cepstrum can be also written as:[8]

${\displaystyle C_{c}={\mathcal {F}}^{-1}\left\{\log _{e}({\mathcal {|F|}})+i\phi \right\}}$

The complex cepstrum retains the information about the phase. Thus it is always possible to return from the quefrency domain to the time domain by the inverse operation:[2][3]

${\displaystyle f(t)={\mathcal {F}}^{-1}\left\{b^{\left({\mathcal {F}}\{C_{c}\}\right)}\right\}}$, where b is the base of the used logarithm

Main application is the modification of the signal in the quefrency domain (liftering) as an analog operation to filtering in the spectral frequency domain.[2][3] An example is the suppression of echo effects by suppression of certain quefrencies.[2]

The real cepstrum is derived from the complex cepstrum by setting the phase of the spectrum to zero.[4] It has a focus on periodic effects in the amplitudes of the spectrum:

${\displaystyle C_{r}={\mathcal {F}}^{-1}\left\{\log({\mathcal {|{\mathcal {F}}\{f(t)\}|}})\right\}}$

Thus, it is also directy related to the power spectrum, as (see above):

${\displaystyle C_{p}=4\cdot C_{r}^{2}}$

Steps in forming cepstrum from time history

A short-time cepstrum analysis was proposed by Schroeder and Noll for application to pitch determination of human speech.[9][10][11]

The phase cepstrum is related to the complex cepstrum as

phase spectrum = (complex cepstrum − time reversal of complex cepstrum)².

The kepstrum, which stands for "Kolmogorov-equation power-series time response", is similar to the cepstrum and has the same relation to it as expected value has to statistical average, i.e. cepstrum is the empirically measured quantity, while kepstrum is the theoretical quantity. It was in use before the cepstrum.[12][13]

The cepstrum can be seen as information about the rate of change in the different spectrum bands. It was originally invented for characterizing the seismic echoes resulting from earthquakes and bomb explosions. It has also been used to determine the fundamental frequency of human speech and to analyze radar signal returns. Cepstrum pitch determination is particularly effective because the effects of the vocal excitation (pitch) and vocal tract (formants) are additive in the logarithm of the power spectrum and thus clearly separate.[11]

The autocepstrum is defined as the cepstrum of the autocorrelation. The autocepstrum is more accurate than the cepstrum in the analysis of data with echoes.

The cepstrum is a representation used in homomorphic signal processing, to convert signals combined by convolution (such as a source and filter) into sums of their cepstra, for linear separation. In particular, the power cepstrum is often used as a feature vector for representing the human voice and musical signals. For these applications, the spectrum is usually first transformed using the mel scale. The result is called the mel-frequency cepstrum or MFC (its coefficients are called mel-frequency cepstral coefficients, or MFCCs). It is used for voice identification, pitch detection and much more. The cepstrum is useful in these applications because the low-frequency periodic excitation from the vocal cords and the formant filtering of the vocal tract, which convolve in the time domain and multiply in the frequency domain, are additive and in different regions in the quefrency domain.

Recently cepstrum based deconvolution was used to remove the effect of the stochastic impulse trains, which originates an sEMG signal, from the power spectrum of sEMG signal itself. In this way, only information on motor unit action potential (MUAP) shape and amplitude were maintained, and then, used to estimate the parameters of a time-domain model of the MUAP itself.[14]

The independent variable of a cepstral graph is called the quefrency.[15] The quefrency is a measure of time, though not in the sense of a signal in the time domain. For example, if the sampling rate of an audio signal is 44100 Hz and there is a large peak in the cepstrum whose quefrency is 100 samples, the peak indicates the presence of a fundamental frequency that is 44100/100 = 441 Hz. This peak occurs in the cepstrum because the harmonics in the spectrum are periodic and the period corresponds to the fundamental frequency, since harmonics are integer multiples of said fundamental frequency.

Note that a pure sine wave can not be used to test the cepstrum for its pitch determination from quefrency as a pure sine wave does not contain any harmonics and does not lead to quefrency peaks. Rather, a test signal containing harmonics should be used (such as the sum of at least two sines where the second sine is some harmonic (multiple) of the first sine, or better, a signal with a square or triangle waveform, as such signals provide many overtones in the spectrum.).

Playing further on the anagram theme, a filter that operates on a cepstrum might be called a lifter. A low-pass lifter is similar to a low-pass filter in the frequency domain. It can be implemented by multiplying by a window in the quefrency domain and then converting back to the frequency domain, resulting in a modified signal, i.e. with signal echo being reduced.

A very important property of the cepstral domain is that the convolution of two signals can be expressed as the addition of their complex cepstra:

${\displaystyle x_{1}*x_{2}\mapsto x'_{1}+x'_{2}.}$

• D. G. Childers, D. P. Skinner, R. C. Kemerait, "The Cepstrum: A Guide to Processing", Proceedings of the IEEE, Vol. 65, No. 10, October 1977, pp. 1428–1443.
• "Speech Signal Analysis"
• "Speech analysis: Cepstral analysis vs. LPC", www.advsolned.com
• "A tutorial on Cepstrum and LPCCs"
• Alan V. Oppenheim and Ronald W. Schafer, "From Frequency to Quefrency: A History of the Cepstrum", IEEE SIGNAL PROCESSING MAGAZINE, Sep. 2004, pp. 95–99