数字语音处理及MATLAB仿真-第三章

数字语音处理及MATLAB仿真张雪英编著1第三章语音信号的短时时域分析3.1概述13.2语音信号的预处理3.3短时平均能量33.4短时平均幅度函数43.5短时平均过零率53.6短时自相关分析623.7基于能量和过零率的语音端点检测73.8基音周期估值8数字语音处理及MATLAB仿真张雪英编著23.1概述语音信号是一种非平稳的时变信号，它携带着各种信息。在语音编码、语音合成、语音识别和语音增强等语音处理中都需要提取语音中包含的各种信息。语音处理的目的：对语音信号进行分析，提取特征参数，用于后续处理；加工语音信号。总之，语音信号分析的目的就在于方便有效的提取并表示语音信号所携带的信息。数字语音处理及MATLAB仿真张雪英编著3根据所分析的参数类型，语音信号分析可以分成时域分析和变换域（频域、倒谱域）分析。其中时域分析方法是最简单、最直观的方法，它直接对语音信号的时域波形进行分析，提取的特征参数主要有语音的短时能量和平均幅度、短时平均过零率、短时自相关函数和短时平均幅度差函数等。数字语音处理及MATLAB仿真张雪英编著43.2语音信号的预处理在对语音信号进行数字处理之前，首先要将模拟语音信号s(t)离散化为s(n).实际中获得数字语音的途径一般有两种，正式的和非正式的。正式的是指大公司或语音研究机构发布的被大家认可的语音数据库，非正式的则是研究者个人用录音软件或硬件电路加麦克风随时随地录制的一些发音或语句。数字语音处理及MATLAB仿真张雪英编著5语音信号的频率范围通常是300~3400Hz，一般情况下取采样率为8kHz即可。本书的数字语音处理对象为语音数据文件，是已经数字化了的语音。有了语音数据文件后，对语音的预处理包括：预加重、加窗分帧等。数字语音处理及MATLAB仿真张雪英编著63.2.1语音信号的预加重处理预加重目的：为了对语音的高频部分进行加重，去除口唇辐射的影响，增加语音的高频分辨率。可通过一阶FIR高通数字滤波器来实现：1()1Hzz设n时刻的语音采样值为x(n)，经过预加重处理后的结果为()()(1)ynxnxn数字语音处理及MATLAB仿真张雪英编著7高通滤波器的幅频特性和相频特性如下数字语音处理及MATLAB仿真张雪英编著8预加重前和预加重后的一段语音信号时域波形数字语音处理及MATLAB仿真张雪英编著9预加重前和预加重后的一段语音信号频谱数字语音处理及MATLAB仿真张雪英编著103.2.2语音信号的加窗处理由于发音器官的惯性运动，可以认为在一小段时间里（一般为10ms~30ms）语音信号近似不变，即语音信号具有短时平稳性。这样，可以把语音信号分为一些短段（称为分析帧）来进行处理。数字语音处理及MATLAB仿真张雪英编著11语音信号的分帧实现方法：采用可移动的有限长度窗口进行加权的方法来实现的。一般每秒的帧数约为33～100帧。分帧一般采用交叠分段的方法，这是为了使帧与帧之间平滑过渡，保持其连续性。前一帧和后一帧的交叠部分称为帧移，帧移与帧长的比值一般取为0~1/2。数字语音处理及MATLAB仿真张雪英编著12图3.3给出了帧移与帧长示意图。数字语音处理及MATLAB仿真张雪英编著13加窗常用的两种方法：矩形窗，窗函数如下：其它0101)(Nnnw汉明(Hamming)窗，窗函数如下0.540.46cos2π/10()0nNnNwn其它，数字语音处理及MATLAB仿真张雪英编著14矩形窗及其频谱如下数字语音处理及MATLAB仿真张雪英编著15汉明窗及其频谱如下思考：两种窗效果有何异同？数字语音处理及MATLAB仿真张雪英编著16加窗方法示意图：数字语音处理及MATLAB仿真张雪英编著17窗长的选择一般选取100～200。原因如下：当窗较宽时，平滑作用大，能量变化不大，故反映不出能量的变化。当窗较窄时，没有平滑作用，反映了能量的快变细节，而看不出包络的变化。数字语音处理及MATLAB仿真张雪英编著18语音信号的分帧处理，实际上就是对各帧进行某种变换或运算。设这种变换或运算用T[]表示，x(n)为输入语音信号，w(n)为窗序列，h(n)是与w(n)有关的滤波器，则各帧经处理后的输出可以表示为：[()]()nmQTxmhnm数字语音处理及MATLAB仿真张雪英编著19几种常见的短时处理方法是：22[()]()()()Txmxmhnwn，nQ1.对应于能量；[()]sgn[()]sgn[(1)]()()Txmxmxmhnwn，2.，nQ对应于平均过零率；[()]()()Txmxmxmkhnwnwnk，3.nQ对应于自相关函数；数字语音处理及MATLAB仿真张雪英编著203.3短时平均能量1．短时平均能量定义定义n时刻某语音信号的短时平均能量En为：)]()([)]()([)1(22nNnmmnmnwmxmnwmxE当窗函数为矩形窗时，有)()1(2nNnmnmxE数字语音处理及MATLAB仿真张雪英编著21若令则短时平均能量可以写成：)()(2nwnh)()()()(22nhnxmnhmxEmnx(n)x2(n)Enh(n)(.)2图3.7语音信号的短时平均能量实现方框图数字语音处理及MATLAB仿真张雪英编著222.En特点：En反映语音信号的幅度或能量随时间缓慢变化的规律。3.窗的长短对于能否由短时能量反映语音信号的幅度变化，起着决定性影响。如果窗选得很长，En不能反映语音信号幅度变化。窗选得太窄，En将不够平滑。通常，当取样频率为10kHz时，选择窗宽度N=100~200是比较合适的。数字语音处理及MATLAB仿真张雪英编著23不同矩形窗长N时的短时能量函数数字语音处理及MATLAB仿真张雪英编著24短时平均能量的主要用途如下：1）可以作为区分清音和浊音的特征参数。2）在信噪比较高的情况下，短时能量还可以作为区分有声和无声的依据。3）可以作为辅助的特征参数用于语音识别中。数字语音处理及MATLAB仿真张雪英编著25MATLAB的具体实现如下：1、用Cooledit读入语音“我到北京去”。2、将读入的语音文件wav保存为txt文件，设置采样率为8kHz，16位，单声道。3、把保存的文件zqq.txt读入Matlab。fid=fopen('zqq.txt','rt');x=fscanf(fid,'%f');fclose(fid);4、对采集到的语音样点值进行分帧。数字语音处理及MATLAB仿真张雪英编著263.4短时平均幅度函数为了克服短时能量函数计算x2(m)的缺点，定义了短时平均幅度函数:mnmnwmxM)(|)(|数字语音处理及MATLAB仿真张雪英编著27x(n)Mnw(n)图3.9短时平均幅度|x(n)||.|Mn与En的比较：1.Mn能较好地反映清音范围内的幅度变化；2.Mn所能反映幅度变化的动态范围比En好；3.Mn反映清音和浊音之间的电平差次于En。数字语音处理及MATLAB仿真张雪英编著28短时平均幅度函数随矩形窗窗长N变化的情况数字语音处理及MATLAB仿真张雪英编著293.5短时平均过零率1.定义在离散时间语音信号情况下，如果相邻的采样具有不同的代数符号就称为发生了过零。单位时间内过零的次数就称为过零率。短时平均过零率的定义为sgnsgn1sgnsgn1*nmZxmxmwnmxnxnwn数字语音处理及MATLAB仿真张雪英编著3010120nNwnN其它1()0sgn()1()0xnxn-xn及在上式中，用1/2N作为幅值，是考虑了对该窗口范围内的过零数取平均的意思。数字语音处理及MATLAB仿真张雪英编著31考虑到w(n-m)的非零值范围为n-m≥0，即m≤n，以及n-m≤N-1,故m≥n-N+1，因此短时平均过零率可以改写为：nNnmnmxmxNZ)1(|1-sgnsgn|21sgnsgn1sgnsgn1*nmZxmxmwnmxnxnwn（定义式）数字语音处理及MATLAB仿真张雪英编著322.实现短时平均过零率Nnx(n)Sgn[x(n)]1一阶差分sgn[xn]-sgn[xn-1]低通滤波h(n)=w(n)-1+1图3.11语音信号的短时平均跨零数Sgn[.]取绝对值│.│数字语音处理及MATLAB仿真张雪英编著33女声“我到北京去”的短时平均过零次数的变化曲线:数字语音处理及MATLAB仿真张雪英编著343.应用清音过零率高，浊音过零率低。局限性：浊音和清音重叠区域只根据短时平均过零率不可能明确地判别清、浊音。清音浊音01020304050607080每10ms内的过零数过零率概率分布数字语音处理及MATLAB仿真张雪英编著35端点检测端点检测目的：从包含语音的一段信号中确定出语音的起点及结束点。有效的端点检测不仅能使处理时间减到最少，而且能抑制无声段的噪声干扰，提高语音处理的质量。数字语音处理及MATLAB仿真张雪英编著363.6短时自相关分析3.6.1短时自相关函数时域离散确定信号的自相关函数定义为：时域离散随机信号的自相关函数定义为：周期为P的周期信号满足：∑∞∞)(mkmxmxkRNNmNkmxmxNkR121limPkRkR数字语音处理及MATLAB仿真张雪英编著37自相关函数具有下述性质：(1)对称性R(k)=R(-k)(2)在k=0处为最大值，即对于所有k来说，|R(k)|≤R(0)(3)对于确定信号，R(0)对应于能量对于随机信号，R(0)对应于平均功率数字语音处理及MATLAB仿真张雪英编著38采用短时分析方法，定义语音信号短时自相关函数为mnmknwkmxmnwmxkR因为nnRkRk所以nnmRkRkxmxmkwnmwnmk3.6.2语音信号的短时自相关函数数字语音处理及MATLAB仿真张雪英编著39定义(3-18)那么短时自相关函数可以写成：knwnwnhkmnhkmxmxkRmkn)()(knxnx)(nhk)(kRn上式表明，序列经过一个冲激响应为的数字滤波器滤波即得到短时自相关函数数字语音处理及MATLAB仿真张雪英编著40Rn(k)z-kx(n)x(n–k)hk(n)mnhkmxmxkRmkn数字语音处理及MATLAB仿真张雪英编著41'mnm)()(mwmw[()()][()()]nmmRkxnmwmxnmkwkmxnmwmxnmkwkm也可采用直接运算的方法，令则可得：kNmnmkwkmnxmwmnxkR10上式可以写成数字语音处理及MATLAB仿真张雪英编著42浊音的短时自相关函数数字语音处理及MATLAB仿真张雪英编著43清音的短时自相关函数数字语音处理及MATLAB仿真张雪英编著44浊音和清音的短时自相关函数有如下几个特点：1）短时自相关函数可以很明显的反映出浊音信号的周期性。2）清音的短时自相关函数没有周期性，也不具有明显突出的峰值，其性质类似于噪声。3）不同的窗对短时自相关函数结果有一定的影响。数字语音处理及MATLAB仿真张雪英编著45图3.16不同矩形窗长时的短时自相关函数数字语音处理及MATLAB仿真张雪英编著463.6.3修正的短时自相关函数修正的短时自相关函数，其定义如下k