过程包括了语音信号处理,预加重,分帧加窗,端点检测与求MFCC

const float PI1=3.1415926536;

`typedef struct _AudFrame
{
float fltFrame[FRM_LEN];
float fltSte;	 //这一帧的短时能量
DWORD	 dwZcr;	 //这一帧的过零率
bool	 blVad;//判断这帧是否有效
struct _AudFrame * AudFrmNext;//下一帧的地址,不知道python有没有地址
}AudFrame;`

函数名:ReadWav(CString strFileName)

功能:读取音频

函数名:AudPreEmphasize(void)

功能:对所有采样点进行预处理

解释:
spSound:处理前的采样数据
fpPreSound:处理后的采样数据
fpPreSound[i] = (float)(spSound[i]) - float(spSound[i - 1]) * 0.9375;

函数名:AudEnframe(float *Sound,DWORD FrmLen,DWORD FrmSft,DWORD dwSoundLen)

功能:给每一帧的fltFrame[Frmlen]赋采样点的值,个数是帧长

解释: Sound 输入的采样点数据的起始地址 Frmlen是帧的长度 FrmSft是帧移 dwSoundlen是采样点长度 FrmNum表示有多少帧 由采样点长度,帧长,帧移得到帧的个数FrmNum FrmNum = (dwSoundLen - (FrmLen - FrmSft)) / FrmSft;

函数名:Hamming(DWORD FrmLen)

功能:求汉明窗系数,输入的是每一帧的帧长,要用到PI。这个数组是固定值,只有帧长决定

解释:
FrmLen 帧长,固定值
fltHamm[i] = (float)(0.54 - 0.46cos((2i*PI1) / (FrmLen-1)));

函数名:AudHamming(DWORD FrmLen)

功能:输入的是每一帧的帧长,需要利用到求得的汉明窗系数,具体是每个采样点的值乘以汉明窗系数,再把结果赋予fltFrame[]

解释:
需要先求得汉明窗系数fltHamm[i]; (stpWav->fltFrame)[i] //每一帧采样点数据 for(DWORD i = 0; i < FrmLen; i++) { (stpWav->fltFrame)[i] *= fltHamm[i]; }

函数名:AudSte(fltSound *fpFrmSnd, DWORD FrmLen)

功能:求每一帧的短时能量,即将所有这一帧的所有样点值相加,fpFrmSnd是帧第一个样本值

解释:
fltShortEnergy:每一帧短时能量
fpFrmSnd:每个样本点的值 for(int i = 0; i < FrmLen; i++) { fltShortEnergy += fabs(*fpFrmSnd++); }

函数名:AudZcr(fltSound *fpFrmSnd, DWORD FrmLen,fltSound ZcrThresh)

功能:求解一帧的过零率,fpFrmSnd帧第一个采样点地址,FrmLen帧长,ZcrThresh过零率阀值

解释:

##fpFrmSnd样本点的值 DWORD CVad::AudZcr(fltSound *fpFrmSnd, DWORD FrmLen,fltSound ZcrThresh) { DWORD dwZcrRate = 0;

for(int i = 0; i < FrmLen - 1; i++)
{
	if((fpFrmSnd[i]*fpFrmSnd[i + 1] < 0)&&(fabs(fpFrmSnd[i] - fpFrmSnd[i - 1]) > ZcrThresh))
		dwZcrRate++;
}
return dwZcrRate;
}

函数名: AudNoiseEstimate()

功能:计算双门限阀值

解释: fltSteThresh [2] 短时能量阀值,[0]高 [1]低 dwZcrThresh [2] 过零率阀值, [0]高 [1]低 ZcrThresh = 0; StrThresh = 0.0; ZcrThresh = 所有帧的过零率之和 StrThresh = 所有帧的短时能量之和 NoiseFrmLen 信号帧数 dwZcrThresh[0] = (float)(ZcrThresh) / NoiseFrmLen; dwZcrThresh[1] = (float)(ZcrThresh) / NoiseFrmLen2.5; fltSteThresh[0] = (float)StrThresh / NoiseFrmLen0.7; fltSteThresh[1] = (float)(StrThresh / NoiseFrmLen)*0.5;//*0.95;

函数名: AudVadEstimate(void)

功能:端点检测,需要用到估计阀值的函数,最后得出有效起始点和有效截止点