语音信号处理第4讲剖析

4.1概述4.2语音分帧4.3语音信号的时域分析4.4语音信号的频域分析第四章语音信号分析第一页，共三十四页。语音信号分析语音信号处理的前提和基础，只有分析出可表示语音信号特征的参数，才有可能利用这些参数进行高效的语音通信、语音合成和语音识别等处理。贯穿于语音分析全过程的是“短时分析技术”语音信号从整体来看其特征及表征其本质特征的参数均是随时间而变化的，所以它是一个非平衡态过程，不能用处理平衡信号的数字信号处理技术对其进行分析处理。

但是在一个短时间范围内（一般认为在10-30ms的短时间内），其特性基本保持不变即相对稳定，因而可以将其看作是一个准稳态过程，即语音信号具有短时平稳性。

4.1概述第二页，共三十四页。根据所分析出的参数的性质的不同可以分为：

时域分析、频域分析、倒频域分析,线性预测分析等;分析方法的不同：

模型分析方法和非模型分析方法不论是分析怎么样的参数以及彩什么分析方法，在按帧进行语音分析，提取语音参数之前，有一些经常使用的、共同的短时分析技术必须预先进行，如语音信号的数字化、语音信号的端点检测、预加重、加窗和分帧等，这些也是不可忽视的语音信号分析的关键技术。

4.1概述第三页，共三十四页。

4.2语音分帧进行过预加重数字滤波处理后，接下来就要进行加窗分帧处理。一般每秒的帧数约为33-100帧，视实际情况而定。分帧虽然可以采用连续分段的方法，但一般要采用如图3-1所示的交叠分段的方法，这是为了使帧与帧之间平滑过渡，保持其连续性。前一帧和后一帧的交叠部分称为帧移。帧移与帧长的比值一般取为0-1/2。分帧是用可移动的有限长度窗口进行加权的方法来实现的，这就是用窗函数ω(n)来乘s(n)，从而形成加窗语音信号sω(n)=s(n)ω(n)。第四页，共三十四页。语音分帧帧长与帧移：

取的点数为帧长，乘以采样周期得帧长时间

两帧间隔为帧移两帧一定要有重叠，否则有参数突变在语音信号数字处理中常用的窗函数是矩形窗和汉明窗等，它们的表达式如下（其中N为帧长）：矩形窗：汉明窗：第五页，共三十四页。第六页，共三十四页。第七页，共三十四页。语音分帧1.窗口的形状虽然，不同的短时分析方法(时域、频域、倒频域分析)以及求取不同的语音特征参数可能对窗函数的要求不尽一样，一般来讲，一个好的窗函数的标准是：在时域因为是语音波形乘以窗函数，所以要减小时间窗两端的坡度，使窗口边缘两端不引起急剧变化而平滑过渡到零，这样可以使截取出的语音波形缓慢降为零，减小语音帧的截断效应；在频域要有较宽的3dB带宽以及较小的边带最大值。第八页，共三十四页。语音分帧2.窗口的长度采样周期Ts=1/fs，窗口长度N和频率分辨率Δf之间存在下列关系：Δf=1/NTs可见，采样周期一定时，Δf随窗口宽度N的增加而减小，即频率分辨率相应得到提高，但同时时间分辨率降低；如果窗口取短，频率分辨率下降，而时间分辨率提高，因而二者是矛盾的。应该根据不同的需要选择合适的窗口长度。第九页，共三十四页。

语音分帧有时窗口长度的选择，更重要的是要考虑语音信号的基音周期。通常认为在一个语音帧内应包含1～7个基音周期。然而不同人的基音周期变化很大，从女性和儿童的2ms到老年男子的14ms(即基音频率的变化范围为500～70Hz)，所以N的选择比较困难。通常在10kHz取样频率下，N折中选择为100～200点为宜(即10～20ms持续时间)。这样，经过上面介绍的处理过程，语音信号就已经被分割成一帧一帧的加过窗函数的短时信号，然后再把每一个短时语音帧看成平稳的随机信号，利用数字信号处理技术来提取语音特征参数。在进行处理时，按帧从数据区中取出数据，处理完成后再取下一帧，等等，最后得到由每一帧参数组成的语音特征参数的时间序列。第十页，共三十四页。语音信号的时域分析就是分析和提取语音信号的时域参数。进行语音分析时，最先接触到并且也是最直观的是它的时域波形。语音信号本身就是时域信号，时域分析是最早使用，也是应用最广泛的一种分析方法，这种方法直接利用语音信号的时域波形。时域分析通常用于最基本的参数分析及应用，如语音的分割、预处理、大分类等。特点是：①表示语音信号比较直观、物理意义明确。②实现起来比较简单、运算量少。③可以得到语音的一些重要的参数。④只使用示波器等通用设备，使用较为简单等。4.3语音信号的时域分析第十一页，共三十四页。

第十二页，共三十四页。

短时平均能量及含义短时平均能量：每区一个n，得到一个第十三页，共三十四页。

短时能量及短时平均幅度分析En是一个度量语音信号幅度值变化的函数，但它有一个缺陷，即它对高电平非常敏感（因为它计算时用的是信号的平方)。短时平均幅度函数Mn，它定义为：Mn也是一帧语音信号能量大小的表征，它与En的区别在于计算时小取样值和大取样值不会因取平方而造成较大差异。第十四页，共三十四页。（1）区分清/浊音：

大，对应浊音，小，对应清音。（2）在SNR高的情况下，能进行有声/无声判决

无声时，背景噪声的小

有声时，显著增大。判决时可设置一个门限（3）大致能定出浊音变为清音的时刻，或反之。存在的问题：短时能量函数对信号电平值过于敏感，在实际应用中（如定点设备）很容易溢出。此时可以用平均幅度函数代替，单对清/浊音，有/无声

的幅度差不如短时能量明显。

短时平均幅度函数和能量函数的作用第十五页，共三十四页。

短时过零率（ZCR）分析短时过零率：一帧语音中语音信号波形穿过横轴(零电平)的次数。对于连续语音信号，过零即意味着时域波形通过时间轴；对于离散信号，如果相邻的取样值改变符号则称为过零。过零率就是样本改变符号的次数。第十六页，共三十四页。

短时过零率分析定义语音信号xn(m)的短时过零率Zn为：式中，sgn[]是符号函数，即第十七页，共三十四页。

短时过零率分析(1)区分清/浊音：

清音平均过零率高，集中在高频端

浊音平均过零率低，集中在低频端(2)从背景噪声中找出是否有语音。以及语音的起点。清音——随机噪声，高频浊音——周期信号，低频第十八页，共三十四页。

短时过零率分析在实际应用中，短时平均过零率容易受到A/D转换的直流偏移、50Hz交流电源的干扰以及噪声的影响。减少这些干扰可以有两种方法:(1)采用带通滤波器消除信号中的直流和50Hz低频分量(2)用过门限率来修改过零率，减少随机噪声的影响过门限率反应了穿过正负门限的次数，如果存在随机噪声，只要信号没有超过[-T,T]的范围，就没有过零率的产生第十九页，共三十四页。

短时相关分析相关分析是一种常用的时域波形分析方法，并有自相关和互相关之分。这里主要讨论自相关函数。自相关函数具有一些性质，如它是偶函数；假设序列具有周期性，则其自相关函数也是同周期的周期函数等。可以应用于语音信号的时域分析中。在语音信号分析中，采用的是短时自相关函数。第二十页，共三十四页。

1．短时自相关函数定义语音信号xn(m)的短时自相关函数Rn(k)的计算式如下：这里K是最大的延迟点数。短时自相关函数具有以下性质：(1)如果xn(m)是周期的(设周期为N)，则自相关函数是同周期的周期函数，即Rn(k)=Rn(k+Np)。(2)Rn(k)是偶函数，即Rn(k)=Rn(-k)。(3)当k=0时，自相关函数具有最大值，即Rn(0)≥|Rn(k)|，并且Rn(0)等于确定性信号序列的能量或随机性序列的平均功率。第二十一页，共三十四页。

1．短时自相关函数第二十二页，共三十四页。短时自相关函数存在的问题：随着k的变化，参加运算的项减少。极限k=N-1时无运算！改进方法：两个不同长度的窗2.修正的短时自相关函数矩形窗1：矩形窗2：第二十三页，共三十四页。

3．相关函数的作用（1）区分清/浊音

浊音语音的自相关函数具有一定的周期性

清音语音的自相关函数不具有周期性，类似噪声。（2）估计浊音语音信号的周期，即估计基音周期第二十四页，共三十四页。

短时平均幅度差函数短时自相关函数是语音信号时域分析的重要参量。但是，计算自相关函数的运算量很大，其原因是乘法运算所需要的时间较长。为了避免乘法，一个简单的方法就是利用差值。为此常常采用另一种与自相关函数有类似作用的参量，即短时平均幅度差函数(AMDF)。平均幅度差函数能够代替自相关函数进行语音分析的原理：如果信号是完全的周期信号(设周期为Np)，则相距为周期的整数倍的样点上的幅值是相等的，差值为零。第二十五页，共三十四页。

短时平均幅度差函数短时平均幅度差函数的定义：对于周期性的语音信号，也呈周期性，与相反的是：在谷点时，对应的是峰值平均幅度差函数和自相关函数有密切的关系可由下式表达：第二十六页，共三十四页。短时平均幅度差函数第二十七页，共三十四页。语音信号的频域分析就是分析语音信号的频域特征。从广义上讲，语音信号的频域分析包括语音信号的频谱、功率谱、倒频谱、频谱包络分析等，而常用的频域分析方法有带通滤波器组法、傅里叶变换法、线性预测法等几种。本章介绍的是语音信号的傅里叶分析法。因为语音波是一个非平稳过程，因此适用于周期、瞬变或平稳随机信号的标准傅里叶变换不能用来直接表示语音信号，而应该用短时傅里叶变换对语音信号的频谱进行分析，相应的频谱称为“短时谱”。4.4语音信号的频域分析第二十八页，共三十四页。

利用语音的短时傅里叶变换求语音的短时谱对第n帧语音信号，短时傅里叶变换，其定义如下：可知函数是时间n的离散函数，又是角频率的连续函数，因此，可以从两个角度解释短时傅里叶变换的含义：（1）标准傅里叶变换（2）滤波器角度是窗口函数序列，不同的窗口函数序列，将得到不同的傅里叶变换结果第二十九页，共三十四页。讨论（1）当n取固定值时，记

，标准傅里叶变换为：（2）ω固定不变，记ω=L，相当于一个滤波器当n取不同值时，窗沿着序列滑动标准的FT第三十页，共三十四页。离散的短时傅里叶变换令则在语音信号数字处理中，都是采用短时傅里叶变换也叫短时谱特点：（1）时变性：即是角频率的函数也是时间n的函数（2）周期性：是关于角频率的周期函数，周期为短时功率谱

功率谱:函数的功率谱=该函数自相关+傅里叶变换第三十一页，共三十四页。根据功率谱定义，可以写出短时功率谱与短时傅里叶变换之间的关系：或者：式中*表示复共轭运算。并且功率谱Sn(ejω)是短时自