[发明专利]用于语种识别的时频二维倒谱特征提取方法

权利要求书

1.用于语种识别的时频二维倒谱特征提取方法，其特征在于所述方法是在数字集成电路芯片中按以下步骤实现的：

步骤(1)：对语音信号进行零均值化和预加重，其中零均值化是指整段语音减去其均值，预加重是对语音进行高通滤波，滤波器传输函数为H(z)＝1-0.975z^-1；

步骤(2)：对语音信号按帧长20ms，帧移10ms进行分帧处理；

步骤(3)：按以下步骤建立一个同时反映语音短时平稳性和语种长时信息的二维时频分布矩阵：

步骤(3.1)：对所述语音信号加汉明窗，得到数据{x(m)，m＝0，1，...，M-1}，M为一帧数据点数；

步骤(3.2)：对加加汉明窗的数据做DFT变换(离散傅立叶变换)，得到：

X(ωk)=Σm=0M-1x(m)e-j2πMmk]]>

其中ω_k代表频率，k代表频率标号；

步骤(3.3)：按Mel频标在频率域中用下式计算每帧内F个三角窗的子带能量e_f，F＝24：

ef=1Uf-Lf+1Σk=LfUf|X(ωk)|2]]>

其中U_f和L_f分别为第f个子带的上下边界，再把F个子带能量组成一个矢量e：

e＝[e₀，e₁，…，e_F-1]^T

其中上标T表示转置；

步骤(3.4)：取步骤(3.3)中T帧矢量并列在一起，形成一个二维时频分布矩阵E(t)，T＝19：

E(t)=[e(t),e(t+1),···,e(t+T-1)]]]>

步骤(4)：对矩阵E(t)进行二维DCT(离散余弦变换)，得到二维倒谱系数：

C(p,q)=γpγqΣτ=0T-1Σf=0F-1ef(t+τ-1)cosπ(2τ+1)p2Tcosπ(2f+1)q2F]]>

其中τ和f为求和变量，γ_p和γ_q为归一化系数：

γp=1/T,p=02/T,p≥1,]]>γq=1/F,q=02/F,q≥1]]>

步骤(5)：选取作为矩阵E(t)主要分量的左上角部分的元素作为特征，用TFC表示，则把左上角部分排列为矢量的重排公式为：

TFC((p+q)2+3p+q2)=C(p,q).]]>