#include?“MFCC.h“
#include?
#include?
#include?
#include?“fftw3.h“
using?namespace?std;


void?MFCC（const?short*?waveData?int?numSamples?int?sampleRate）
{
	float?*preemp?=?new?float[numSamples];
	preemp[0]?=?waveData[0];
	for（int?i?=?1;?i?	{
		preemp[i]?=?waveData[i]?-?0.95*waveData[i-1];
		//printf（“%f?“?preemp[i]）;
	}
	//printf（“\n“）;
	//calculate?the?total?frames?when?overlaps?exist
	int?numframes?=??ceil（（numSamples-frameS_PER_BUFFER）/NOT_OVERLAP）+1;
	//printf（“%d\n“?numframes）;
	float?hammingWindow[frameS_PER_BUFFER];
	float?afterWin[512]?=?{0.0};
	float?energySpectrum[512]?=?{0.0};
	float?mel[NUM_FILTER]?=?{0};
	float?c[13]?=?{0};
	for（int?i?=?0;?iameS_PER_BUFFER;?i++）
	{
		hammingWindow[i]?=?0.54?-?0.46*cos（2*3.14*i/（frameS_PER_BUFFER-1））;
	}
	//handle?all?frames?one?by?one
	for（int?i?=?0;?i?ames;?i++）
	{
		int?j;
		//windowing
		for（j=0;?jameS_PER_BUFFER&&（j+（i-1）*NOT_OVERLAP）		{
			afterWin[j]?=?preemp[j+（i-1）*NOT_OVERLAP]*hammingWindow[j];??
		}
		//Zero?Padding
		for（int?k?=?j-1;?k?		//for（int?a?=?0;?a?		//printf（“\n“）;
		//FFT?+?power
		FFT_Power（afterWin?energySpectrum）;
		//for（int?a?=?0;?a?		//printf（“\n“）;
		//Warping?the?frequency?axis?+?FilterBank
		memset（mel?0?sizeof（float）*NUM_FILTER）;
		computeMel（mel?sampleRate?energySpectrum）;
		//for（int?a?=?0;?a?		//printf（“\n“）;
		//DCT
		memset（c?0?sizeof（float）*13）;
		DCT（mel?c）;
		//for（int?k?=?0;?k?		//printf（“\n“）;
	}
	delete[]?preemp;
}

void?FFT_Power（float?*in?float?*energySpectrum）
{
	int?len?=?512;
	fftwf_complex?*out?=?（fftwf_complex?*）fftwf_malloc（sizeof（fftwf_complex）?*?len）;
	fftwf_plan?p?=?fftwf_plan_dft_r2c_1d（len?in?out?FFTW_ESTIMATE）;?
	fftwf_execute（p）;?
	for?（int?i?=?0;?i?	{??
		energySpectrum[i]?=?out[i][0]*out[i][0]?+?out[i][1]*out[i][1];??
	}
	fftwf_destroy_plan（p）;
	fftwf_free（out）;??
}

void?computeMel（float?*mel?const?int?sampleRate?const?float?*energySpectrum）
{
	int?fmax?=?sampleRate/2;
	float?maxMelFreq?=?1127*log（1+fmax/700）;
	float?melFilters[NUM_FILTER][3];
	float?delta?=?maxMelFreq/（NUM_FILTER+1）;
	float?*m?=?new?float[NUM_FILTER+2];
	float?*h?=?new?float[NUM_FILTER+2];
	float?*f?=?new?float[NUM_FILTER+2];
	for（int?i?=?0;?i?	{
		m[i]?=?i?*?delta;
		h[i]?=?700?*?（exp（m[i]/1127）?-?1）;
		f[i]?=?floor（（256+1）?*?h[i]?/?sampleRate）;??//*********//
	}
	//get?start?peak?end?point?of?every?trigle?filter
	for（int?i?=?0;?i?	{
		melFilters[i][0]?=?f[i];
		melFilters[i][1]?=?f[i+1];
		melFilters[i][2]?=?f[i+2];
	}
	delete[]?m;
	delete[]?h;
	delete[]?f;
	//calculate?the?output?of?every?trigle?filter
	for（int?i?=?0;?i?