function out=MFCC_3(in,cofnumber,fs)
%這個函數(shù)是用來計算摩爾倒譜系數(shù)
%輸入為預(yù)處理后的矩陣 摩爾濾波器個數(shù) 還有采樣頻率
%輸出就是輸入矩陣的摩爾倒譜系數(shù)
in=abs(fft(in)).^2;
[a,b]=size(in);
in=in(1:a/2,:);
inter=fs/a;
index=ones(1,cofnumber);
Mef=2595*log10(1+fs/1400);
step_f=Mef/cofnumber;
for i=1:cofnumber
    k=i*step_f/2595;
    fre=(10.^(k)-1)*700;
    index(i)=floor(fre/inter);
end
template=zeros(a/2,cofnumber-2);
for i=1:cofnumber-2
    mid=index(i+1);
    step1=1/(mid-index(i));
    step2=1/(mid-index(i+2));
    for j=index(i):mid
        template(j,i)=(j-index(i))*step1;
    end
    for j=mid:index(i+2)
        template(j,i)=(j-index(i+2))*step2;
    end
end
out=zeros(cofnumber-2,b);
for i=1:b
    for j=1:cofnumber-2
        out(j,i)=sum(in(:,i).*template(:,j));
        k=out(j,i);
        if isnan(k)
            out(j,i) = 0;
        end
    end
end
out=dct(out);    
t=1:size(out,1);
A=2.5*sin(t*pi*1/(2*size(out,1)));
for i=1:size(out,2)
    out(:,i)=out(:,i).*A';
end