实验1_实验一语音信号MFCC特征提取_MFCC_

clear; tic [signal,FS] = audioread('D:\1.wav');%读取音频信号 x=signal';%将sianal的转置矩阵赋给x y=filter([1,-0.97],1,x);%预加重 N=400;%每帧400采样点（25ms*16kHZ=400）,帧移10ms*16kHZ=160 N_FFT=512;%样本数 i=1; start=1; framesample=zeros(ceil(length(y)/N),N_FFT);%成帧 while start+400-1<=length(y) framesample(i,1:N)=y(start:start+400-1) ; start=start+160; i=i+1; end framesample(i,1:length(y)-start+1)=y(start:end); H=hamming(N_FFT);%加窗 AW_sample=framesample.*H'; z=fft(AW_sample');%傅里叶变换 zz=abs(z); plot(zz); AFFT_sample=(z.^2)./N_FFT;%功率谱 f_max=8000;%采样频率*2为最大频率 mel_max=2595*log10(1+f_max/700);%求解mel的最大值，最小值为0 for i=1:40 %40个滤波器 mel(i)=mel_max/41*i; fmel(i)=round(700*(10^(mel(i)/2595)-1));%依次算出每个点的值 end fmel(41)=f_max; f_interval=f_max/(N_FFT/2); fbank=[1 round(fmel./f_interval)]; for i=1:40 Hm(i,:)=mytriangle(fbank(i),fbank(i+1),fbank(i+2),N_FFT/2);%通过三角滤波器 end fbankm=Hm*AFFT_sample(1:N_FFT/2,:); Fm_sample=10.*log10(fbankm); DCT=dct(Fm_sample);%倒谱系数 MFCC=DCT(1:13,:); toc