mfcc源代码

function c=melcepst（sfs）
%MELCEPST Calculate the mel cepstrum of a signal C=（SFSWNCPNINCFLFH）
%
%
% Simple use: c=melcepst（sfs）	% calculate mel cepstrum with 12 coefs 256 sample frames
%				  c=melcepst（sfs‘e0dD‘） % include log energy 0th cepstral coef delta and delta-delta coefs
%
% Inputs:
%     s	 speech signal
%     fs  sample rate in Hz （default 11025）
%     nc  number of cepstral coefficients excluding 0‘th coefficient （default 12）
%     n   length of frame （default power of 2 <30 ms））
%     p   number of filters in filterbank （default floor（3*log（fs）） ）
%     inc frame increment （default n/2）
%     fl  low end of the lowest filter as a fraction of fs （default = 0）
%     fh  high end of highest filter as a fraction of fs （default = 0.5）
%
%		w   any sensible combination of the following:
%
%				‘R‘  rectangular window in time domain
%				‘N‘	Hanning window in time domain
%				‘M‘	Hamming window in time domain （default）
%
%		      ‘t‘  triangular shaped filters in mel domain （default）
%		      ‘n‘  hanning shaped filters in mel domain
%		      ‘m‘  hamming shaped filters in mel domain
%
%				‘p‘	filters act in the power domain
%				‘a‘	filters act in the absolut