孤立汉语数字语音识别系统

1、孤立汉语数字语音识别系统孤立汉语数字语音识别系统摘要：本文通过提取声音信号的Mel倒谱系数作为特征，利用动态时间规整技 术实现匹配算法，实现了特定人孤立汉语数字语音的识别，并利用Matlab编写 了简单的图形用户界面。关键词：语音识别；MFCC； DTVV一、引言语言是人类所特有的最重要最自然的交流工具，也是人类信息的重要来源之 一。让机器拥有“听懂”人类口述语言的能力，将使得人与计算机之间的沟通 变得更为方便快捷。自从1952年AT&T贝尔实验室的开发出能识别十个英文数 字的Audry系统以来,语音识别技术已经得到了飞速发展，其中IBM、Microsoft、 Apple等公司在语音识

2、别技术实用化上的巨资投入也使得这项技术在日常生活 中得到广泛应用o计算机技术的发展大大促进了数字信号处理技术的开发与应 用，也使得更多的技术应用在语音识别方面。MATLAB是一种功能强大、效率高、交互性好的数值计算和可视化计算 机高级语言，它将数值分析、信号处理和图形显示有机地融合为一体，使得技 术人员可以在较短的时间里对自己的想法进行实验验证。本文就是在MATLAB 基础上进行汉语数字语言识别的图形用户界面开发与程序设计的o二、语音识别系统概述针对不同的任务，语音识别系统也可以采用对应的多种设计方案，但其大致 的结构和模型思想基本上是相同的。语音识别系统本质上是一种模式识别系统, 它包括了特

3、征提取、模式匹配、参考模型库这三个基本单元。典型系统的基本 结构如图1所示。图1典型语音识别系统目前常用的技术有特征参数匹配法、隐马尔科夫法和神经网络法。其中语音 特征参数可以是能量、基音频率、过量率、共振峰值等，目前常用的参数为基 于发声声道特征模型，通过线性预测分析的线性预测倒谱系数LPCC和基于人 耳的听觉机理,反映听觉特性,模拟人耳对声音频率感知的梅尔倒谱系数MFCC 等。MFCC的频率特性反映了人耳的听觉特性，因而在用于代替人耳来分析语 音时，其性能和鲁棒性比较符合实际听音效果。与LPCC相比，不依赖于全极 点语音产生模型的假定，抗噪能力较强。本系统也将使用MFCC参数作为语音 识别

4、的主要特征参数。三、语音信号的分析与处理1信号采集本实验直接用电脑的声卡进行语音信号的数字化采样量化处理。根据奈奎斯 特采样定理，采样频率必须到达信号的两倍带宽以上时，才不会发生频域上的 混叠失真，能够从采样信号恢复出原始信号。理论上说，通常人声的频率范围 上限可达15kHz或更高，但是即使当带宽限制到低于3kHz附近的频率时语音 仍是较容易听懂的。对于数字电话通信系统，标准采样率为8000Hz。量化是将 时间上离散而幅度依然连续的波形幅度值离散化，其决定了声音的动态范围。 若采用8位，则可将声波分为256级。本系统采用Matlab的audiorecorder对象进行声音信号采集。使用默认参数

5、 8kHz采样率，8位量化精度，单声道。2预加重处理将经采样后的数字语音信号s(ii)通过一个高通漉波器：H(z)= 1 - aXz-1 , 0.9 <a<1.0 (本系统中取0.97)。经过预加重后的信号为：sYn)= s(n) - aX s(n-l)o因为发声过程中声带和嘴唇的效应，使得高频共振峰的振幅低于低频共 振峰的振幅，进行预加重的目的就是为了对信号进行平滑，提升高频部分，补 偿声带和嘴唇的效应，平坦信号频谱，保持整个频带中的信噪比，以便于进行 频谱分析或声道参数分析。3端点检测时域波形图2端点检测结果端点检测即从一段信号中检测出语音信号的起始点和结束点，在孤立词识别 中

6、，可以减少非实时系统中的大量计算。常见的端点检测方法有利用时域特征 的音量、过零率，利用频域特征的频谱变异数、频谱的燧等。时域方法计算量 小，实现简单，且足以满足一般应用°短时能量可以利用帧内数据的平方和进行计算。理论上有声语音的能量值较 大，无声语音（如摩擦音）的过零率较高，可以通过计算短时能量大致判断出 语音的端点，然后使用过零率找到语音端点的相对精确位置。实验时发现汉语 语音通过能量特征已能较好地确定端点，此时通过过零率特征并不能进一步精 化结果，因此系统只使用了能量参数。编程时当能量值大于开启阈值时确定为起始点，小于关闭阈值时确定为语音 结束点。阈值大小根据实验情况确定。图2

7、为汉语5发音的端点检测情况。4 MFCC特征提取数字语音信号经过以上的预加重处理和端点检测之后，将进行梅尔倒谱系数 MFCC的计算。计算方法如下：1）分帧语音整体上是非平稳信号，但其特性在20ms的量级内基本保持不变。取 10ms为一帧，为避免相邻帧间的过大变化，帧之间应设置重叠，选取帧移为帧 长的一半。另外为了方便后面的FFT计算，这里的帧长取为2的整幕指数。2）加窗由于进行离散傅里叶变换时将时域信号进行了周期延拓，在两端一般会有不 连续的现象，为了减少吉布斯效应，通常要进行加窗处理，弱化边缘信号。在 语音信号处理中常用海明窗（hamming windows）,在Matlab中有其实现函数。

8、3）傅里叶变换语音的时域信号变化迅速，不易分析，而频域变化相对缓慢，因此对加窗后 的信号进行快速傅里叶变换FFT,再折半舍弃对称部分后平方计算信号能量谱。4）三角带通滤波器图3梅尔频率滤波器组将每帧能量谱中的线性频谱刻度转化成符合人耳听觉特性的梅尔刻度后，将 其通过如图3的一组N个三角形带通滤波器（N 一般取2030个，本系统取为24),计算通过每个频带的能量并取对数。其中梅尔刻度与线性频率的转换关系为：Mel(/) = 2595*log10(l + )5)离散余弦转换将上面得到的N=24个对数能量Ek带入离散余弦转换公式，求出L阶的梅 尔频率倒谱系数，这里L通常取12。离散余弦转换公式如下：

9、N(1 .= Z与 cos m(km = L2.Lfir2 N)通过离散余弦转换，将能量集中在前面几项中，达到在减少判别参数提高运 算速度的同时又不失其准确性的目的。6)差量倒谱参数很多文献中提到利用差量倒谱参数来对语音的动态特性参数进行描述，可以 提高系统的识别能力。本系统实验发现使用差量倒谱系数后，系统的识别性能 并未有明显提高，有些情况下甚至有所下降。综合考虑，本系统仅使用静态 MFCC参数进行识别。5训练与识别在特征匹配方面有动态时间规整(DTW)技术，有效地解决了语音信号的特征 提取和不等长语音匹配问题。此外还有矢量量化(VQ)和隐马尔科夫模型(HMM) 理论，但这些基于统计的识别技

10、术都需要大量的语音材料来提取统计特征，实 现也比较麻烦。相对而言，DTW实现简单，在特定人孤立词的识别方面也能够 有足够的准确度，为本实验所采用的方法。动态时间规整是将时间规整和距离测度结合起来的一种非线性规整技术。算 法的思想就是把未知量均匀地伸长或缩短，直到与参考模式的长度一致在这 一过程中，未知语音的时间轴要不均匀地扭曲或弯折，以使其特征与模型特征 对正。3 21/M 8762313678910 N图4 DTVV算法搜索路径示意图如图%要找时间规整函数M=w(N),使测试矢量的时间轴N非线性地映射到 模板的时间轴M上，并满足：MO = minZdr(i)R(w(i)l卬,=式中，d是第1

11、帧测试矢量T和第w帧模板矢量R(w)之间的距离(一般采 用欧氏距离)。一般DTW采用逆向思路，从最后点开始，逆推到起始点，寻找其 中的最优路径。在求解中除了要约束路径的起点和终点外，还可以约束斜率在 02之间，使路径不至于过度倾斜，这样每一个格点向前向后都只有三种可能落 点，以加快算法的计算速度。真正进行语音识别时，首先需要采集09的汉语语音信号，计算各自的 MFCC参数并存储起来作为参考模板。进行识别时，将待识别信号的参数特征 依次与模板特征进行DTW匹配，并从中选取距离最小的作为识别结果。四、Matlab的GUI实现最后采用Matlab的GUI工具箱搭建一个简单的用户界面，能够用前面介绍 的原理进行汉语数字语音的识别，而且便于操作，方便使用。用户界面如图5所示，开始时用户要根据界面提示依次读出0-9的汉语读音, 每次录音结束后，界面上会用两条竖线显示出端点检测的结果，因为有时候检 测的端点是错误的，此时需要用户判断本次录音是否有效，选择“是，重读数 字”或“否，往下进行”。10个样本数据采集完毕后，就可以录入语音进行识别, 并