#include?
#include?
#include?
#include?
#include?“l(fā)pc.h“

#define?frameLEN		240		//幀長
#define?frameSHIFT		80		//幀移
#define?LPCFACTOR		10		//lpc階數(shù)
#define?LPCCFACTOR		10		//lpcc階數(shù)
#define?REALMAX			100000000

//////////////////////////////////////////////////////////////////
///		函數(shù)功能：計(jì)算lpc系數(shù)		
///		參數(shù)含義：	fData----一幀數(shù)據(jù)，輸入?yún)?shù)
///					fResult---返回結(jié)果，輸出參數(shù)
//////////////////////////////////////////////////////////////////
void?lpc（float*?fData?float*?fResult）
{
	int?ijk;
	float?r[LPCFACTOR+1];
	float?l[LPCFACTOR];
	float?a[LPCFACTOR][LPCFACTOR];
	float?e[LPCFACTOR];
	float?c[LPCFACTOR];
	float?d[LPCFACTOR];
		
	for（i=0;?i	{
		float?sum=0.0;
		for（j=0;?jameLEN-i;?j++）
		{
			sum?+=?（*（fData+j））?*?（*（fData+j+i））;
		}
		*（r+i）?=?sum;
	}


	l[0]?=?r[1]/r[0];
	a[0][0]?=?l[0];
	e[0]?=?（1-squre（l[0]））*r[0];
	
	for（i=1;?i	{
		float?sum=0.0;
		for（j=0;?j		{
			sum?+=?（a[i-1][j]）*（r[i-j]）;
		}
		c[i]?=?sum;
		l[i]?=?（r[i+1]-c[i]）/e[i-1];
		a[i][i]?=?l[i];
		for（j=0;?j		{
			a[i][j]?=?a[i-1][j]-l[i]*a[i-1][i-j-1];
		}
		e[i]?=?（1-squre（l[i]））*e[i-1];
	}
	for（k=0;?k	{
		d[k]?=?a[LPCFACTOR-1][k];
		d[k]?=?（-1）*d[k];
	}
	memcpy（fResult?d?sizeof（float）*LPCFACTOR）;
	return;
}


/////////////////////////////////////////////////////////////////
//		函數(shù)功能：由lpc得到lpcc
//		參數(shù)含義：	lpcdata---lpc系數(shù)所在的內(nèi)存首地址輸入?yún)?shù)
//					lpccdata---返回的lpcc系數(shù)，輸出參數(shù)
/////////////////////////////////////////////////////////////////
void?lpctolpcc（float*?lpcdata?float*?lpccdata）
{
	int?mk;
	float?sum=0.0;
	float?tempData[LPCCFACTOR];
	for（m=0;?m	{
		sum?=?*（lpcdata+m）;
		for（k=0;?k		{
			sum?+=?（*（lpcdata+k））*tempData[m-k-1]*（1.0-（float）（k+1）/（m+1））;
		}
		tempData[m]=sum;
	}
	for（m=LPCFACTOR;?m	{
		sum=0.0;
		for（k=0;?k		{
			sum?+=?（*（lpcdata+k））*tempData[m-k-1]*（1.0-（float）（k+1）/（m+1））;
		}
		tempData[m]=sum;
	}
	memcpy（lpccdata?tempData?sizeof（float）*LPCCFACTOR）;
	return;
}

////////////////////////////////////////////////////////////////////////////////////
//		函數(shù)功能：對(duì)模板參數(shù)和樣本參數(shù)進(jìn)行dtw比較
//		參數(shù)含義：	mobandata---模板參數(shù)數(shù)組，輸入?yún)?shù)，大小為mbfrmcnt*LPCCFACTOR
//					yangbendata---樣本參數(shù)數(shù)組，輸入?yún)?shù)，大小為ybfrmcnt*LPCCFACTOR
//					mbfrmcnt---模板的幀數(shù)，輸入?yún)?shù)
//					ybfrmcnt---樣本幀數(shù)，輸入?yún)?shù)
//					cmpresult---比較結(jié)果，輸出參數(shù)，大小為mbfrmcnt*ybfrmcnt
////////////////////////////////////////////////////////////////////////////////////
float?dtw（float*?mobandata?float*?yangbendata?int?mbfrmcnt?int?ybfrmcnt）
{
	int?ijk;
	float?sumd1d2d3retvalue;
	float?paccd[50][50];
	//生成幀匹配距離矩陣
	float?pd[50][50];
//	float*?pd?=?（float*）malloc（sizeof（float）*mbfrmcnt*ybfrmcnt）;
	for（i=0;?i	{
		for（j=0;?j		{
			sum=0.0;
			for（k=0;?k			{
				sum?+=?squre（（*（mobandata+i*LPCCFACTOR+k））-（*（yangbendata+j*LPCCFACTOR+k）））;
			}
			pd[i][j]=sum;
//			*（pd+i*ybfrmcnt+j）?=?sum;
		}
	}

	//
//	paccd?=?（float*）malloc（sizeof（float）*mbfrmcnt*ybfrm

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????4041??2006-11-04?13:37??DSP實(shí)現(xiàn)簡單語言系統(tǒng)\用DSP實(shí)現(xiàn)的一個(gè)簡單的語音識(shí)別系統(tǒng)\lpc.c

?????目錄??????????0??2007-05-09?11:57??DSP實(shí)現(xiàn)簡單語言系統(tǒng)\用DSP實(shí)現(xiàn)的一個(gè)簡單的語音識(shí)別系統(tǒng)

?????目錄??????????0??2010-04-02?15:49??DSP實(shí)現(xiàn)簡單語言系統(tǒng)

-----------?---------??----------?-----??----

?????????????????4041????????????????????3