mfcc程序

clear all;
close;
clc;
[x,fs]=audioread('bj.wav');
%x1=zeros(60001,1);
%x1(1:60001,1)=x(30000:90000,1);
x1=x(:,1);
l=length(x1);
%x1=x1(1:60000)-0.97*x1(2:60001);%预加重
%x1=x1(1:l-1)-0.95*x1(2:l);
xx=double(x1);
x1=filter((1-0.9375),1,xx);
x1max=max(x1);
x1=x1/x1max;
sound(x1,fs);
figure(1);
plot(x1);
x2=x1(:,1);
inc=512;%针移
win=hamming(1024);%汉明窗
frameout=enframeself(x1,win,inc);
[h,w]=size(frameout);
%短时傅里叶变换
for i=1:h
    frameout(i,:)=fft(frameout(i,:));
end
%画语谱图
figure(2);
tt=(0:inc:(l-w))/fs;
ff=(0:(w/2))*fs/w;
yuputu=frameout(:,1:1+inc)';
imagesc(tt/1000,ff/1000,20*log10(abs(yuputu)));
%笛卡尔轴模式，原点在左下角，y轴是竖直的，
%由底至顶标数，x轴是水平的ÿ