具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1和图2，本发明公开了一种应用于电力领域中声纹识别的身份认证方法，其不同之处在于，包括以下步骤：

步骤1，采集语音内容形成语音信号，并对所述语音信号进行预处理，包括端点检测处理，采样量化处理，分帧处理，加窗处理和预加重处理。

首先端点检测消除静音信号。然后使用采样量化保证语音文件的精确性，其中采样频率和信道数影响语音文件的大小，由于是基于语音进行身份认证，所以使用高采样频率(即44100Hz)，并且使用单声道采集语音数据。

分帧处理是将N个采样点集合成一个观测单位，称为帧。设分帧后的信号为S(n)，其中n＝0,1,...,N-1,N，设N的值为200，为帧的大小。

加窗处理采用的窗函数为汉明窗，设窗函数为W(n)，加窗后的信号为S'(n)，则：

其中0≤n≤N-1,a＝0.55。

S'(n)＝S(n)×W(n)

预加重处理是使用一阶预加重数字滤波器去实现，设预加重后的信号S₁(n)：

S₁(n)＝S'(n)-μS'(n-1)

其中μ为预加重系数，μ＝0.93。通过预加重处理可以对高频部分进行增强，使信号的频谱变得平坦。

步骤2：对预处理后的语音信号S₁(n)进行快速傅立叶变换，通过信号在时域上的改变分析信号的特性比较困难，因此需要将信号转换到频域上观察其能量分布，能量分布的差异，表明不同的语音特性。

对每一帧音频信号S₁(n)作快速傅立叶变换，

得到其在频谱上的能量分布X(k)，其中k＝1,2,...,N-1,N。

步骤3：对步骤二得到的语音信号频谱做模平方处理，得到语音信号的功率谱|X(k)|²，其中k＝1,2,...,N-1,N。

步骤4：定义一个滤波器组，该滤波器组有M个三角滤波器，其中心频率为f(m)，其中k＝1,2,...,M-1,M，取M＝30。将步骤三得到的功率谱通过滤波器组，对语音信号进行滤波，其中三角滤波器的频率响应为H_m(k)

可以得到M组参数H_m，其中m＝1,2,...,M-1,M。

步骤5：计算每个滤波器组输出的对数能量S(m)；

步骤6：将步骤五的对数能量S(m)求离散余弦变换，求出L阶的梅尔倒谱参数C(n)

其中L取值为16。

步骤7：通过倒谱均值减，得到信道噪声均值，使得带噪语音特征尽可能接近于零，从而消除信道的不利噪声影响。

步骤8：动态特征提取。标准的倒谱参数只反映了语音参数的静态特性，语音的动态特性可以用这些静态特征的差分谱来描述。实验证明：把动、静态特征结合起来才能有效提高系统的识别性能。差分参数的计算可以采用下面的公式

其中,dt表示第t个一阶差分；Ct表示第t个倒谱系数；Q表示倒谱系数的阶数；K表示一阶导数的时间差，取2。将上式中结果再代入就可以得到二阶差分的参数。

步骤9：利用相似度计算完成身份认证，在认证阶段，需要计算验证声纹码特征向量与注册声纹码特征向量之间的相似性。假设验证声纹码和注册声纹码分别用T和R表示，长度分别为N和M。则

T＝t₁,t₂,...,t_N

R＝r₁,r₂,...,r_M

通过计算它们之间的距离D＝[T,R]，来比较声纹之间的相似度，距离越小则相似度越高。设n和m分别是T和R中任意选择的帧号，则d(T(n),R(m))表示这两帧特征矢量之间的距离。则各帧间距离的累积和即为两序列的总体距离。

D(n,m)＝d(T(n),R(m))+min{D(n-1,m),D(n-1,m-1),D(n,m-1)}

求解出来的累加距离D，即为需要计算验证声纹码特征向量与注册声纹码特征向量之间的相似性差值，距离越大，相似度越小。

需要说明的是，在本文中，诸如术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。