发明名称 用于语种识别的约束异方差线性鉴别分析方法
摘要 用于语种识别的约束异方差线性鉴别分析方法涉及一种高维特征矢量的降维和去相关处理方法,其特征在于,所述方法首先由语音信号提取MFCC特征,选取连续M帧MFCC特征并列放置得到倒谱矩阵,然后将其按行展开,形成超矢量,再分块计算超矢量的均值和协方差,并用迭代法分块计算变换矩阵,接着用变换矩阵对超矢量分块进行变换,最后将各块进行降维和拼接处理,可得到新的特征矢量。该方法计算量小,得到的特征矢量去除了各维之间的相关性,鉴别性高且维数低。可以用于语种识别。
申请公布号 CN101546555B 申请公布日期 2011.05.11
申请号 CN200910081899.3 申请日期 2009.04.14
申请人 清华大学 发明人 张卫强;刘加
分类号 G10L15/02(2006.01)I 主分类号 G10L15/02(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 朱琨
主权项 1.用于语种识别的约束异方差线性鉴别分析方法,其特征在于,所述方法是在数字集成电路芯片中依次按以下步骤实现的:步骤(1):把一帧语音信号输入到N维特征提取电路,得到N维Mel频率倒谱系数(MFCC)的基本特征矢量,N=10,步骤(2):用步骤(1)所述的方法提取出连续M帧语音信号的Mel频率倒谱系数(MFCC)的N维基本特征矢量,M=21,并输入到一个倒谱矩阵形成电路,得到一个倒谱矩阵C:<img file="FSB00000365438000011.GIF" wi="880" he="310" />步骤(3):把步骤(2)得到的倒谱矩阵C输入到一个串并转换电路,按行展开,得到一个超矢量特征c:<img file="FSB00000365438000012.GIF" wi="357" he="536" />步骤(4):把步骤(3)得到的超矢量特征c输入到一个分路器,分成N块:<img file="FSB00000365438000013.GIF" wi="357" he="369" />其中c<sup>(n)</sup>=[c<sub>n1</sub> c<sub>n2</sub>…c<sub>nM</sub>]<sup>T</sup>代表超矢量特征中的第n块,T代表转置, 步骤(5):对每个语种的语音数据均按照步骤(1)-(4)进行操作,假设一共有L个语种,第<i>l</i>个语种可得到超矢量特征I<sub><i>l</i></sub>帧,所有语种共得到超矢量特征I帧,即<img file="FSB00000365438000021.GIF" wi="237" he="79" />步骤(6):按以下步骤用一个逻辑电路分块计算各语种的超矢量特征的统计量和所有语种总的超矢量特征的统计量,<i>n</i>=1,2,...,<i>N</i>:步骤(6.1):计算第<i>l</i>个语种的各帧中第n块超矢量特征c<sup>(n)</sup>的均值<img file="FSB00000365438000022.GIF" wi="114" he="61" /><img file="FSB00000365438000023.GIF" wi="350" he="155" /><i>i</i>=1,2,...,<i>I</i>,计算第<i>l</i>个语种的第n块超矢量特征c<sup>(n)</sup>的协方差矩阵<img file="FSB00000365438000024.GIF" wi="119" he="57" /><img file="FSB00000365438000025.GIF" wi="770" he="158" /><i>i</i>=1,2,...,<i>I</i>,其中∑<sub><i>g</i>(<i>i</i>)=<i>l</i></sub>表示对所有帧中属于第<i>l</i>个语种的特征矢量求和;步骤(6.2):计算所有语种第n块超矢量特征c<sup>(n)</sup>的总均值m<sup>(n)</sup>:<img file="FSB00000365438000026.GIF" wi="389" he="153" />计算所有语种第n块超矢量特征c<sup>(n)</sup>的总协方差矩阵T<sup>(n)</sup>:<img file="FSB00000365438000027.GIF" wi="837" he="136" />步骤(7):用一个变换矩阵运算电路按以下步骤分块计算超矢量特征的变换矩阵,<i>n</i>=1,2,...,<i>N</i>:步骤(7.1):初始化第n块变换矩阵为单位阵:A<sup>(n)</sup>=I,步骤(7.2):把上述第n块变换矩阵划分为两部分:<img file="FSB00000365438000028.GIF" wi="460" he="166" />其中,<img file="FSB00000365438000029.GIF" wi="92" he="66" />代表第n块变换矩阵A<sup>(n)</sup>的前U<sup>(n)</sup>行,<img file="FSB000003654380000210.GIF" wi="141" he="66" />代表剩下的M-U<sup>(n)</sup> 行,M=21,步骤(7.3):按下式计算第n块变换矩阵的第m行<img file="FSB00000365438000031.GIF" wi="93" he="57" /><img file="FSB00000365438000032.GIF" wi="912" he="171" />其中<img file="FSB00000365438000033.GIF" wi="75" he="74" />是中间变量H<sup>(n)</sup>矩阵的第m行:H<sup>(n)</sup>=det(A<sup>(n)</sup>)(A<sup>(n)</sup>)<sup>-1</sup>,<img file="FSB00000365438000034.GIF" wi="1034" he="293" /><img file="FSB00000365438000035.GIF" wi="471" he="77" />为<img file="FSB00000365438000036.GIF" wi="77" he="57" />的对角线元素:<img file="FSB00000365438000037.GIF" wi="1975" he="283" />其中det代表求矩阵的行列式,diag代表取矩阵的对角线元素,步骤(7.4):当达到20次迭代或两次迭代det(A<sup>(n)</sup>)数值相对变化小于1/1000时,停止迭代,否则转步骤(7.2)继续迭代,步骤(8):用变换矩阵A<sup>(n)</sup>对超矢量特征<img file="FSB00000365438000038.GIF" wi="63" he="57" />分块进行变换,得到新的特征矢量x<sup>(n)</sup>,<i>i</i>=1,2,...,<i>I</i>,<i>n</i>=1,2,...,<i>N</i>:<img file="FSB00000365438000039.GIF" wi="292" he="70" />步骤(9):取出所述特征矢量x<sup>(n)</sup>的前U<sup>(n)</sup>维,得到<img file="FSB000003654380000310.GIF" wi="111" he="62" />n=1,2,....,N,再用一个合路器把各块特征<img file="FSB000003654380000311.GIF" wi="395" he="82" />依次进行拼接,得到新的特征矢量:<img file="FSB000003654380000312.GIF" wi="362" he="300" />
地址 100084 北京市100084-82信箱