function [mfcc] = compute_mfcc(path)
[wav,Fs,bits] = wavread(path)

sof = zeros(1,23248)
for i=1:23248
    sof(i) = wav(i) - mean(wav)
end

logE = zeros(1,144);
for j=0:143
    s = 0
    for i=1:400
        if(j*160+i <= length(sof)) then
        logE(j+1) = logE(j+1) + (sof(j*160+i))**2
        end
    end
    logE(j+1) = max(-50,log(logE(j+1)));
end


spe = zeros(1,23248);
for i=1:23248
    if i==1 then
        spe(i) = sof(i)
    else
        spe(i) = sof(i) - 0.97*sof(i-1)
    end
end

hamming = window('hm',400)
spec_rep = zeros(144,512);
for i=0:143
    if(i*160+400>length(spe)) then
//        spec_rep(i+1,:) = [conv([spe(i*160+1:length(spe)),zeros(1,400-length(spe)+(i*160))],hamming,'same'),zeros(1,112)]
        spectral_rps(i+1,:) = [[spe(i*160+1:length(spe)),zeros(1,400-length(spe)+(i*160))].*hamming,zeros(1,112)]
    else
//        spec_rep(i+1,:) = [conv(spe(i*160+1:i*160+400),hamming,'same'),zeros(1,112)]
        spectral_rps(i+1,:) = [ (spe(i*160+1:i*160+400) .* hamming) , zeros(1,112)]
    end
    spectral_rps(i+1,:) = abs(fft(spectral_rps(i+1,:)))
end


loadmatfile('C:\Users\freespace\Desktop\MFCC\mel_filters.mat')
energy = zeros(144,23)
for i=1:144
    energy(i,:) =  spectral_rps(i,1:257)*mel_filters;
end


loadmatfile('C:\Users\freespace\Desktop\MFCC\gas_station_mfcc.mat')
mfcc = zeros(14,144);
mfcc(1,:) = logE
for i=0:12
    for f=1:144
       for j=1:23
            mfcc(i+2,f) = mfcc(i+2,f) + max(-50,log(energy(f,j)))*cos(%pi*i/23*(j-1/2))
       end
    end
end

endfunction