view dsp/cepstrum.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

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