annotate dsp/cepstrum.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.
author samer
date Sun, 11 Oct 2015 10:20:42 +0100
parents c3b0cd708782
children
rev   line source
samer@32 1 % cepstrum - generalised cepstral coefficients from warped power spectrogram
samer@32 2 %
samer@32 3 % cepstrum ::
samer@32 4 % [[M]->[K]] ~'which coefficients to compute',
samer@32 5 % [[K]] ~'edges for (K-2)-band filterbank',
samer@32 6 % [[N]] ~'centre frequencies of spectrum bins',
samer@32 7 % [[N,T]] ~'input spectra',
samer@32 8 % -> [[M,T]] ~'MFCCs for L time slices'.
samer@32 9 function y = cepstrum(MX,edges,fx,x)
samer@32 10 melW = tri_filterbank(edges,fx);
samer@32 11 dctW = dct(eye(size(melW,1)));
samer@32 12 y = dctW(MX,:)*log(melW*x);
samer@32 13 end