samer@32: % cepstrum - generalised  cepstral coefficients from warped power spectrogram
samer@32: %
samer@32: % cepstrum ::
samer@32: %    [[M]->[K]] ~'which coefficients to compute',
samer@32: %    [[K]]      ~'edges for (K-2)-band filterbank',
samer@32: %    [[N]]      ~'centre frequencies of spectrum bins',
samer@32: %    [[N,T]]    ~'input spectra',
samer@32: % -> [[M,T]]    ~'MFCCs for L time slices'.
samer@32: function y = cepstrum(MX,edges,fx,x)
samer@32: 	melW = tri_filterbank(edges,fx);
samer@32: 	dctW = dct(eye(size(melW,1)));
samer@32: 	y = dctW(MX,:)*log(melW*x);
samer@32: end