MFCC特征提取(可用程序).docx

《MFCC特征提取(可用程序).docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、/******************************************************MFCC特征提取程序*读取一个音频文件(.wav)，将根据帧长分割后的每帧2阶MFCC*系数写在输出文件中，以","为间隔*******************************************************/#include#include#include#include#include#include#include

2、lex>#include//MFCC的特征选取，最后提取了13维左右最后结果保存在文件中usingnamespacestd;#defineSP_EMPHASIS_FACTOR0.97f/*预加重的系数*/typedefstruct_TWavHeader{intrId;//标志符（RIFF）intrLen;//数据大小,包括数据头的大小和音频文件的大小intwId;//格式类型（"WAVE"）intfId;//"fmt"intfLen;//Sizeof(WAVEFORMATEX)shortwFormatTag;//编码格式，包括WAVE_FORM

3、AT_PCM，WAVEFORMAT_ADPCM等shortnChannels;//声道数，单声道为1，双声道为2intnSamplesPerSec;//采样频率intnAvgBytesPerSec;//每秒的数据量shortnBlockAlign;//块对齐shortwBitsPerSample;//WAVE文件的采样大小intdId;//"data"intwSampleLength;//音频数据的大小}TWavHeader;constintFS=16;/*修改帧长*/constlongFrmLen=1024;//可修改帧长constunsignedlongF

4、FTLen=512;//参与FFT运算的512个数据constdoublePI=3.;constintFiltNum=40;//滤波器组数，一共40组constintPCEP=13;//最后得到的关于的13个MFCC的系数doubleHamming[FrmLen];inttemp_1;//计算次数的vectorMFCCcoefficient;staticdoublelast=0;//一窗数据最后一个点的值，此点用于预加重voidpreemphasis(double*buf,double*result,shortFrmLen);//预加重void

5、InitHamming();voidHammingWindow(double*result,double*data);voidcompute_fft(double*buffer,vector>&vecList);voidFFT(constunsignedlong&ulN,vector>&vecList);//FFT的实际程序voidInitFilt(double*FiltCoe1,double*FiltCoe2,int*Num);//初始化滤波器voidCFilt(double*spdata,do

6、uble*FiltCoe1,double*FiltCoe2,int*Num,double*En,vector>&vecList);//计算每个滤波器组内的总能量voidMFCC(double*En,double*Cep);//计算MFCC的13个系数intmain(){//TWavHeaderwaveheader;FILE*sourcefile,*MFCCFile;shortbuffer[FrmLen];doubledBuff[FrmLen];doubleresult[FrmLen];//预加重结果doubledata[FrmL

7、en];//加窗后得到的数据doubleFiltCoe1[FFTLen/2+1];//左系数doubleFiltCoe2[FFTLen/2+1];//右系数intNum[FFTLen/2+1];//一般而言，每个点会包含在相邻的两个滤波器中，这里是与该点相关的第二个滤波器doubleEn[FiltNum+1];//频带能量doubleCep[PCEP];//MFCC结果//对定义的变量赋初值temp_1=0;inti=0;for(i=0;i

8、(i=0;i