mfcc.m

来自「线性预测这一术语是维纳于1947年首次提出的」· M 代码 · 共 28 行

bank=melbankm(32,256,8000,0,0.5,'m');
bank=full(bank);
bank=bank/max(bank(:));
for k=1:16
     n=0:31;
     dctcoef(k,:)=cos((2*n+1)*k*pi/(2*32));
end
w=1+8*cos(pi*[1:16]./16);
w=w/max(w);
m=zeros((s2-s1+1),16);
j=1;
for i=s1:s2
     y=xx(i,:);
     s=y'.*hamming(256);
     t=abs(fft(s));
     t=t.^2;
     c1=dctcoef*log(bank*t(1:129));
     c2=c1.*w';
     m(j,:)=c2';
     j=j+1;
end
dtm=zeros(size(m));
for i=3:size(m,1)-2
      dtm(i,:)=-2*m(i-2)-m(i-1,:)+m(i+1,:)+2*m(i+2,:);
end
dtm=dtm/3;
ccc=[m dtm];
ccc=ccc(3:size(m,1)-2,:);