function [res]=obtainmfcc(x,filter_col)
%Apply hamminwindow and followed by fft
y=hamming(length(x))'.*x;
x1=abs(fft(y,1998));
x2=x1(1:1:1000);
res1=log10(sum(abs(fft(y)).^2));
res2=[];
for i=1:1:24
res2=[res2 sum((x2.*filter_col{i}).^2)];
end
res3=log10(res2);
%Suppose the energy values are scaled by some constant number, taking log
%will have the corresponding DC shift. This is further removed using dct as
%follows.
res4=dct(res3);
res5=res4(2:1:13);
res6=diff(res5);
res7=diff(res6);
res=[res1 res5 res6 res7];
