www.pudn.com > SpeakerRecognitionsystem.rar > mfcc.m

function r = mfcc(s, fs) 
% MFCC 
% 
% Inputs: s  contains the signal to analize 
%         fs is the sampling rate of the signal 
% 
% Output: r contains the transformed signal 
% 
% 
%%%%%%%%%%%%%%%%%% 
% 
% Author:     Amin Koohi(Sohrevardi) 
%              AminSohrevardi@yahoo.com , AminKoohi@yahoo.com 
 
m = 100; 
n = 256; 
l = length(s); 
 
nbFrame = floor((l - n) / m) + 1; 
 
for i = 1:n 
    for j = 1:nbFrame 
        M(i, j) = s(((j - 1) * m) + i); 
    end 
end 
 
h = hamming(n); 
 
M2 = diag(h) * M; 
 
for i = 1:nbFrame 
    frame(:,i) = fft(M2(:, i)); 
end 
 
t = n / 2; 
tmax = l / fs; 
 
m = melfb(20, n, fs); 
n2 = 1 + floor(n / 2); 
z = m * abs(frame(1:n2, :)).^2; 
 
r = dct(log(z));