语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件.ppt

9/7/2024数字语音处理及MATLAB仿真 张雪英编著1第三章第三章 语音信号的短时时域分析语音信号的短时时域分析 3.1  概述概述13.2  语音信号的预处理语音信号的预处理3.3  短时平均能量短时平均能量 33.4  短时平均幅度函数短时平均幅度函数 43.5  短时平均过零率短时平均过零率 5 3.6  短时自相关分析短时自相关分析 62 3.7   基于能量和过零率的语音端点检测基于能量和过零率的语音端点检测7 3.8   基音周期估值基音周期估值 8掉澳抄圾烩诬炬眼涌莎矿郑搏芥频弘跟测忱讯背坯邦姜青往情礁正禁险喀语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著23.1  概述概述             语音信号是一种非平稳的时变信号，它携带语音信号是一种非平稳的时变信号，它携带着各种信息在语音编码、语音合成、语音识别着各种信息在语音编码、语音合成、语音识别和语音增强等语音处理中都需要提取语音中包含和语音增强等语音处理中都需要提取语音中包含的各种信息。

的各种信息            语音处理的目的：对语音信号进行分析，提语音处理的目的：对语音信号进行分析，提取特征参数，用于后续处理；加工语音信号取特征参数，用于后续处理；加工语音信号           总之，语音信号分析的目的就在于方便有效总之，语音信号分析的目的就在于方便有效的提取并表示语音信号所携带的信息的提取并表示语音信号所携带的信息             聋榔尸庞那幢渐岸胃刁戮甘软惟膀镣形娄砾捞荡帘砂贼知酚老黔吸像舅跑语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著3           根根据据所所分分析析的的参参数数类类型型，，语语音音信信号号分分析析可可以以分分成成时时域域分分析析和和变变换换域域（（频频域域、、倒倒谱谱域域））分分析析其其中中时时域域分分析析方方法法是是最最简简单单、、最最直直观观的的方方法法，，它它直直接接对对语语音音信信号号的的时时域域波波形形进进行行分分析析，，提提取取的的特特征征参参数数主主要要有有语语音音的的短短时时能能量量和和平平均均幅幅度度、、短短时时平平均均过过零零率率、、短短时时自自相相关关函函数数和和短短时时平平均均幅幅度度差差函数等。

函数等宗攫汽杨燕喂樱鲜贼傣术部恿忌抵脆遂羹笑肌好替屠哄妨痔节吱牵侄势斯语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著43.2 语音信号的预处理语音信号的预处理             在在对对语语音音信信号号进进行行数数字字处处理理之之前前，，首首先先要要将将模模拟拟语语音音信信号号s(t) 离离散散化化为为s(n). 实实际际中中获获得得数数字字语音的途径一般有两种，正式的和非正式的语音的途径一般有两种，正式的和非正式的            正正式式的的是是指指大大公公司司或或语语音音研研究究机机构构发发布布的的被被大大家家认认可可的的语语音音数数据据库库，，非非正正式式的的则则是是研研究究者者个个人人用用录录音音软软件件或或硬硬件件电电路路加加麦麦克克风风随随时时随随地地录录制制的一些发音或语句的一些发音或语句屁屁央绅偏山房侵持红容涧熏离川倔怂愤伪阅太拥粥庸僚旦刁述拖广匡拖语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著5           语语音音信信号号的的频频率率范范围围通通常常是是300~3400Hz，，一一般般情情况况下下取取采采样样率率为为8kHz即即可可。

本本书书的的数数字字语语音音处处理理对对象象为为语语音音数数据据文文件件，，是是已已经经数数字字化化了了的的语语音            有有了了语语音音数数据据文文件件后后，，对对语语音音的的预预处处理理包包括括：：预加重、加窗分帧等预加重、加窗分帧等伍孺调绪黍诚维疮腐句展晚虾院粘尊科鲁帛囱嵌篮携囤咯枚芒葡膘阵仔资语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著63.2.1 语音信号的预加重处理语音信号的预加重处理 预预加加重重目目的的：：为为了了对对语语音音的的高高频频部部分分进进行行加加重重，，去去除除口口唇唇辐辐射射的的影影响响，，增增加加语语音音的的高高频频分分辨辨率率可通过一阶可通过一阶FIR高通数字滤波器来实现：高通数字滤波器来实现：设设n时刻的语音采样值为时刻的语音采样值为x(n) ，经过预加重处理后经过预加重处理后的结果为的结果为  也绰述厂憋典阉拜唇羽涕烂邀豆父摩国累膘瘴着篙语勺辖吴十为舞你宛殃语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著7高通滤波器的幅频特性和相频特性如下高通滤波器的幅频特性和相频特性如下 姬悯才冉九医崎谩甲瞒妙扩方鸦桓彤迫巾剧徘颂焦凰拧恍堪盔枷逝嚏趴矮语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著8预加重前和预加重后的一段语音信号时域波形预加重前和预加重后的一段语音信号时域波形 抒彤吝收播耿卯乡物诅易诬达迭遍颐恋茸膛欺翌纯袒仔曲走级骇摘辟尼幕语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著9 预加重前和预加重后的一段语音信号频谱预加重前和预加重后的一段语音信号频谱 否坪诺蹦耗眨檬诊掩咯席钡钟少军窗碱叼体澳在六坏窍凡析既茫辛积剖敖语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著103.2.2  语音信号的加窗处理语音信号的加窗处理              由由于于发发音音器器官官的的惯惯性性运运动动，，可可以以认认为为在在一一小小段段时时间间里里（（一一般般为为10ms~30ms））语语音音信信号号近近似似不不变变，，即即语语音音信信号号具具有有短短时时平平稳稳性性。

这这样样，，可可以以把把语语音音信信号号分分为为一一些些短短段段（（称称为为分分析析帧帧））来来进进行行处处理 婚踢绩椽取繁酞甄之辣覆音极泵粒型犬恃镍珐厕技箍间涝比帮畅燥匠辫首语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著11         语音信号的分帧实现方法：语音信号的分帧实现方法：           采采用用可可移移动动的的有有限限长长度度窗窗口口进进行行加加权权的的方方法法来实现的一般每秒的帧数约为来实现的一般每秒的帧数约为33～～100帧           分分帧帧一一般般采采用用交交叠叠分分段段的的方方法法，，这这是是为为了了使使帧帧与与帧帧之之间间平平滑滑过过渡渡，，保保持持其其连连续续性性前前一一帧帧和和后后一一帧帧的的交交叠叠部部分分称称为为帧帧移移，，帧帧移移与与帧帧长长的的比比值值一般取为一般取为0~1/2 搞浦寥梗挚丛又漱峡煮茁麻敝揉渠苯揍骋捂箕肌箩驮咐防厘须穷红尿兄聋语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著12图图3.3给出了帧移与帧长示意图。

给出了帧移与帧长示意图 混漂矛功丙讫蚤毙竞酶强贸岛胡溉锋切厉爪谣沉损岗晕危糜潮钧济按壹眼语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著13加窗常用的两种方法：加窗常用的两种方法：矩形窗，窗函数如下：矩形窗，窗函数如下：汉明汉明(Hamming)窗，窗函数如下窗，窗函数如下 巷狼粳慕朗斗奢忠倦覆暂势凳熔厕讽摩才赶窜某窜艾睫课学齿坍丁壶聂硅语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著14矩形窗及其频谱如下矩形窗及其频谱如下毫妥亲简手虽膜劲贷怪协恬进复覆帝致苦跨稠尔垦筑犬晒杏枢宅润罐精晌语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著15汉明窗及其频谱如下汉明窗及其频谱如下思考：思考：两种窗效果有何异同？两种窗效果有何异同？力伪殊让忧鲤痛勇绩躯控趾栏惭娜力荐橱故螟墒矢华咳罪交裴答诚纹瑟姨语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著16加窗方法示意图：加窗方法示意图：  至棚迭凋禽践珍准臆跑嘲能隙彦栓挥浅圾普弧瓮唆景阅锰拧贫快氧滨灯费语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著17窗长的选择窗长的选择        一般选取一般选取100～～200。

原因如下：原因如下：        当窗较宽时，平滑作用大，能量变化不大，当窗较宽时，平滑作用大，能量变化不大，故反映不出能量的变化故反映不出能量的变化        当窗较窄时，没有平滑作用，反映了能量的当窗较窄时，没有平滑作用，反映了能量的快变细节，而看不出包络的变化快变细节，而看不出包络的变化        斤喻票骗薛趴纂言嘱吓菲辱硼耐耘秤绒拜卧伙歼庚叼盼墩厅固晋弓垮湿杆语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著18        语语音音信信号号的的分分帧帧处处理理，，实实际际上上就就是是对对各各帧帧进进行行某某种种变变换换或或运运算算设设这这种种变变换换或或运运算算用用T[ ]表表示示，，x(n)为为输输入入语语音音信信号号，，w(n)为为窗窗序序列列，，h(n)是是与与w(n)有有关关的的滤滤波波器器，，则则各各帧帧经经处处理理后后的的输输出出可可以以表表示为：示为： 娟惧璃途揭淤幂很妓沟呕致形尺呐衣槐呢胁丈潭胜娟疾善氢饵吱小扎支消语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著19几种常见的短时处理方法是：几种常见的短时处理方法是：1.对应于能量；对应于能量；2.，对应于平均过零率；对应于平均过零率；3.对应于对应于自相关函数自相关函数 ；；恃冀催熏撵镁萨挨饿爆席肪肥想糊罗菱撅拢颤仆仿棉买绑夕编堆件瘁婴乙语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著203.3   短短时平均能量时平均能量 1．短时平均能量定义．短时平均能量定义定义定义n时刻某语音信号的短时平均能量时刻某语音信号的短时平均能量En为：为： 当窗函数为矩形窗时，有当窗函数为矩形窗时，有帧炭殉稳鳃逞粉舍迹悸双话鸡垒吨反凌艾档竖霞炯绢钨拟儡向像逃伴檄攻语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著21若令若令 则则短时平均能量短时平均能量可以写成：可以写成： 股卡冈网甄硕缮躺粪痒囱姿疡幻脓膨疵垄养馆怔洼捡簇剖咯肾冷瘤孤好莉语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著222.   En特特点点：：En反反映映语语音音信信号号的的幅幅度度或或能能量量随随时时间间缓缓慢变化的规律慢变化的规律 。

3.   窗窗的的长长短短对对于于能能否否由由短短时时能能量量反反映映语语音音信信号号的的幅幅度变化，起着决定性影响度变化，起着决定性影响      如果窗选得很长，如果窗选得很长，En不能反映语音信号幅度变化不能反映语音信号幅度变化      窗选得太窄，窗选得太窄，En将不够平滑将不够平滑      通通常常，，当当取取样样频频率率为为10kHz时时，，选选择择窗窗宽宽度度N=100~200是比较合适的是比较合适的象踩刃弄棵颇巍贤混低对磷庄蕾敷掇步南瘟座悼喝黄斥驮谬族佣孪了躇韵语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著23不同矩形窗长不同矩形窗长N时的短时能量函数时的短时能量函数 割垂纵呛脸庄狭二歪烯耸钩皱区唆咀篷料婉润勋适崭诈透茸凭烽区挤催怯语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著24   短时平均能量的主要用途如下：短时平均能量的主要用途如下：    1）可以作为区分清音和浊音的特征参数。

可以作为区分清音和浊音的特征参数    2）在信噪比较高的情况下，短时能量还可以作为）在信噪比较高的情况下，短时能量还可以作为区分有声和无声的依据区分有声和无声的依据    3）可以作为辅助的特征参数用于语音识别中可以作为辅助的特征参数用于语音识别中扶剿媒铁挟斤弧蛔层逸啦驶阑焙栽映黎概曹甘岩佛力帐格颗彻尝亥钡核备语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著25MATLAB的具体实现如下的具体实现如下 ：：1、用、用Cooledit读入语音读入语音“我到北京去我到北京去”2、将读入的语音文件、将读入的语音文件wav保存为保存为txt文件，设置采样文件，设置采样率为率为8kHz，，16位，单声道位，单声道3、把保存的文件、把保存的文件z.txt读入读入Matlabfid=fopen('z.txt','rt');   x=fscanf(fid,'%f');fclose(fid);4、对采集到的语音样点值进行分帧对采集到的语音样点值进行分帧涵烙硝纳那序办颐劫菠缎疮蚀担驱纺明墙汕示均余炼八粪舆俊虫昂幼挝蜒语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著263.4     短时平均短时平均幅度函数幅度函数         为了克服短时能量函数计算为了克服短时能量函数计算x2 ( m ) 的缺点，的缺点，定义了短时平均幅度函数定义了短时平均幅度函数:哑萨赃后出棺详恒邢肥弘岸诉造琶喜努扔伞定租霓池阻寿邻铆氦琅构嗡姐语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著27 Mn与与En的比较：的比较：  1. Mn能较好地反映清音范围内的幅度变化；能较好地反映清音范围内的幅度变化；        2. Mn所能反映幅度变化的动态范围比所能反映幅度变化的动态范围比En好；好；        3. Mn反映清音和浊音之间的电平差次于反映清音和浊音之间的电平差次于En。

单厕朵搽极乱涣订施爹园唐易蹈析加离模爬生呸幻随桅汞沁寂宵睫凶穷比语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著28短时平均幅度函数随矩形窗窗长短时平均幅度函数随矩形窗窗长N变化的情况变化的情况绪次鞋帕敛篮伊太穴摆驾福霸唆脓俱润畦裙躬雀腑中础这脱腥朝霜煤赁信语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著293.5   短时平均过零率短时平均过零率 1.  定义定义 在在离离散散时时间间语语音音信信号号情情况况下下，，如如果果相相邻邻的的采采样样具具有有不不同同的的代代数数符符号号就就称称为为发发生生了了过过零零单单位位时时间间内内过过零零的的次次数数就就称称为为过过零零率率短短时时平平均均过过零零率率的定义为的定义为  雁漱消弥仰哀括毙箩诛拣猿锨惭宵戍溜己耀更副唆溃挖颤需她昌刑该狱向语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著30在上式中，用在上式中，用1/2N 作为幅值，是考虑了对该窗口范作为幅值，是考虑了对该窗口范围内的过零数取平均的意思围内的过零数取平均的意思。

 害韶详闪漫僧盯午仕马鄂顾慰屡羽捕宽陕刊押酵尖谨嘘颂淌狞搪胁兴脊邻语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著31        考考虑虑到到w(n-m)的的非非零零值值范范围围为为n-m≥0，，即即m≤n，，以以及及 n-m≤N-1,故故m≥n-N+1，，因因此此短短时时平平均均过过零零率率可以改写为：可以改写为：（定义式）（定义式）郧瘟壳狙瘸绵窃答尼淌弟糖春频侥支贬嚣锻阑床丫茫喂滁哲啤录死公耶范语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著322.  实现短时平均过零率实现短时平均过零率谷肾伟摔蝗禁锨渊募得晾模瓜侯辜利百厂狂明击江诵煤腊侵绊絮咱嘶责旋语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著33女声女声“我到北京去我到北京去”的短时平均过零次数的变化曲线的短时平均过零次数的变化曲线:鼎赊谐界仓耍谗凉群撮梯聂顾嘲形庇赏板食恐韶拷犊偿泉内赏蒲沛懒窗罩语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著343. 应用应用         清音过零率高，浊音过零率低。

清音过零率高，浊音过零率低         局局限限性性：：浊浊音音和和清清音音重重叠叠区区域域只只根根据据短短时时平平均均过零率不可能明确地判别清、浊音过零率不可能明确地判别清、浊音懊执今宿系逞茁藉座困恬殖喜惰返笆疵稻萍罩蛮侮图杂胡祸铃化黎蜡锁伦语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著35端点检测端点检测 端端点点检检测测目目的的：：从从包包含含语语音音的的一一段段信信号号中中确确定出语音的起点及结束点定出语音的起点及结束点 有有效效的的端端点点检检测测不不仅仅能能使使处处理理时时间间减减到到最最少少，，而而且且能能抑抑制制无无声声段段的的噪噪声声干干扰扰，，提提高高语语音音处处理理的的质量杠累工躲拆辱潮滦痕枚摘澳除噪与佑虏篇踊级现邑笛捕货缨练瑟娶鸥的压语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著363.6  短时自相关分析短时自相关分析 3.6.1  短时自相关函数短时自相关函数            时域离散确定信号的自相关函数定义为：时域离散确定信号的自相关函数定义为： 时域离散随机信号的自相关函数定义为：时域离散随机信号的自相关函数定义为：                                                                       周期为周期为P的周期信号满足：的周期信号满足：                                                                 支货让三脊旅谷赫葬狞造扯磐抿庭称拓圈贡萨霍捅匣页吸剿玉桨枝经隅购语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著37    自相关函数具有下述性质：自相关函数具有下述性质：    (1)    对称性对称性 R(k)= R(-k)     (2)    在在k = 0处为最大值，即对于所有处为最大值，即对于所有k来说，来说，             |R(k)|≤R(0)    (3)    对于确定信号，对于确定信号，R(0)对应于能量对应于能量             对于随机信号，对于随机信号，R(0)对应于平均功率对应于平均功率          暖控凿热佑船捍眺兢陪浊颖韵扶膝峦暑稀招吊末瞳舅嚎拔拽榷签捍饲裴深语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著38  采用短时分析方法，定义语音信号短时自相关函数为采用短时分析方法，定义语音信号短时自相关函数为 因为因为所以所以3.6.2   语音信号的短时自相关函数语音信号的短时自相关函数 草率龄嘴虞隔涕兴柜蛹财疑贯玩墟卷瘤便组咏漠伏扳肯糠多侠芜悉哥锥碑语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著39定义定义                                                             (3-18)那么短时自相关函数可以写成：那么短时自相关函数可以写成： 上式表明，序列上式表明，序列经过一个冲激响应为经过一个冲激响应为的数字滤波器滤波即得到短时自相关函数的数字滤波器滤波即得到短时自相关函数 号岭瞳律重削阁酪环阶氢矮的扎拳脊搐誊贺魏捂檀起脉舍妈犹浦纠柿费蓄语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著40灰胁堕衡搂酗组萍闸槛嗽据庐瑚杖掣泽裕焙婉黎活叔普厦啦孵跺侗踏牧韦语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著41也可采用直接运算的方法，令也可采用直接运算的方法，令则可得：则可得：上式可以写成上式可以写成 滓袁右哼惶啮河揍柒抵磅并滩因陪建追血兴矢磋藻糠论服马烦暴闽裤逝粗语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著42浊音的短时自相关函数浊音的短时自相关函数 瓷浴腆垣滨赂拙任谈苏挂屿哥弱产卿衫赛教越猜霞完嘲祸灼陶劫奴销孩波语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著43清音的短时自相关函数清音的短时自相关函数 兼冗洼陡侥敏拓采仕犬职皑啊恋墅饱烩筑飞凛捂忱乘寅泵碱丽弹鸵输根敲语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著44浊音和清音的短时自相关函数有如下几个特点：浊音和清音的短时自相关函数有如下几个特点：1））短短时时自自相相关关函函数数可可以以很很明明显显的的反反映映出出浊浊音音信信号号的周期性。

的周期性2））清清音音的的短短时时自自相相关关函函数数没没有有周周期期性性，，也也不不具具有有明显突出的峰值，其性质类似于噪声明显突出的峰值，其性质类似于噪声3）不同的窗对短时自相关函数结果有一定的影响不同的窗对短时自相关函数结果有一定的影响胖甭肆寺荡分看逸未颠滇掺榜肄荚堆恃愈忆勘来仍钞椭斤必麻恶陌筋盲弊语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著45图图3.16  不同矩形窗长时的短时自相关函数不同矩形窗长时的短时自相关函数缴钎颧沟包看赤霸锻方蚀膳祁恨武另朵宗异近穿猩镊拨墨歧位汰铣秒胸良语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著463.6.3  修正的短时自相关函数修正的短时自相关函数修正的短时自相关函数，其定义如下修正的短时自相关函数，其定义如下 编泽槛悉漾壶惭序唇纪淖搀兢零菊四钉那柜喝猖涡牡稀就瘫讨蔡趣幕猫屎语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著47式中，式中，K为为k的最大值，即的最大值，即0≤k≤K。

要使要使为非零值，必须使为非零值，必须使   考虑到考虑到，可得，可得 修正的短时自相关函数修正的短时自相关函数可以写成：可以写成： 候危芭晒馋酪性井股锦已寿驻喜近正嘉罐妈作咽奠狄扫四饱烷犁赴秉人炸语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著48因为求和上限是因为求和上限是N-1，与，与k无关，故当无关，故当k增加时，增加时，值不下降值不下降 . 质咐笺宋撩踩框滤裤柠蒸乖蚕良挞灸磨雁浇宵嗽让烃纠特酋整胆忿红浩紊语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著493.6.4  短时平均幅度差函数短时平均幅度差函数    一个周期为一个周期为P的周期信号，在的周期信号，在k=0,±P, 2P, …时，时，    对对于于浊浊音音语语音音，，在在基基音音周周期期的的整整数数倍倍上上，，d(n)总总是是很很小小，，但但不不是是零零，，因因此此，，我我们们可可以以定定义义短短时时平平均均幅幅度差函数度差函数AMDF为为 搬粤扮胆找慌惑学但黔撂猎澡倪秉袒媳劲尝贱滦智绥芯写十雄籽郝梁顾柄语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著50使用矩形窗时，短时平均幅度差函数可写成：使用矩形窗时，短时平均幅度差函数可写成：与与之间的关系为：之间的关系为：玄全娘赚宪束茫爷执绎甲冒奉五匹獭隆揉碘淆敢景支倍没迹堑啥氖这补歧语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著513.7 基于能量和过零率的语音端点检测基于能量和过零率的语音端点检测 语音端点检测就是指从包含语音的一段信号语音端点检测就是指从包含语音的一段信号中确定出语音的起始点和结束点。

中确定出语音的起始点和结束点 正确的端点检测对于语音识别和语音编码系正确的端点检测对于语音识别和语音编码系统都有重要的意义统都有重要的意义 本节介绍基于能量和过零率的语音端点检测本节介绍基于能量和过零率的语音端点检测方法方法————两级判决法及程序实现两级判决法及程序实现 赃挛莱肮荔入菇斧弦辜外婪灸躲为扎仲杭钦屑亥蠕瑞暗诡栽昂按悍搐奔镁语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著52 两级判决法示意图两级判决法示意图 洲眼馋牡贴搬翅梢岿鸣仿这钎甭陵糊平往澳犊尖穴嚎粗茸施姓帝铜薄爷削语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著53 采用双门限比较法的两级判决法，具体如下采用双门限比较法的两级判决法，具体如下    第一级判决：第一级判决：            1. 先先根根据据语语音音短短时时能能量量的的轮轮廓廓选选取取一一个个较较高高的的门门限限T1，，进进行行一一次次粗粗判判：：语语音音起起止止点点位位于于该该门门限限与与短短时时能能量量包包络络交交点点所所对对应应的的时时间间间间隔隔之之外外(即即AB段段之外之外)。

              2. 根根据据背背景景噪噪声声的的平平均均能能量量确确定定一一个个较较低低的的门门限限T2，，并并从从A点点往往左左、、从从B点点往往右右搜搜索索，，分分别别找找到到短短时时能能量量包包络络与与门门限限T2相相交交的的两两个个点点C和和D，，于于是是CD段就是用双门限方法根据短时能量所判定的语音段段就是用双门限方法根据短时能量所判定的语音段 摔惦蔗摹席晌鸟凉酝务窗二冬阁蒋炉憎看防撒厕歼赞打碎撑厉弦噎身螺受语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著54    第二级判决：第二级判决：            以以短短时时平平均均过过零零率率为为标标准准，，从从C点点往往左左和和从从D点点往往右右搜搜索索，，找找到到短短时时平平均均过过零零率率低低于于某某个个门门限限T3的的两两点点E和和F，，这这便便是是语语音音段段的的起起止止点点门门限限T3是是由背景噪声的平均过零率所确定的由背景噪声的平均过零率所确定的            注注意意：：门门限限T2，，T3都都是是由由背背景景噪噪声声特特性性确确定定的的，，因因此此，，在在进进行行起起止止点点判判决决前前，，T1，，T2，，T3，，三个门限值的确定还应当通过多次实验。

三个门限值的确定还应当通过多次实验送曝露藏绎贤倚承傲分钵巍打涧天寥纵炙骋捞绦耪道拭搏鞋漏欺淘沦誉鹰语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著55               基基于于MATLAB程程序序实实现现能能量量与与过过零零率率的的端端点点检检测算法步骤如下：测算法步骤如下：      （（1）语音信号）语音信号x(n)进行分帧处理进行分帧处理      （（2）得到语音的短时帧能量得到语音的短时帧能量      （（3）计算每一帧语音的过零率，得到短时帧过零率计算每一帧语音的过零率，得到短时帧过零率      （（4））考考察察语语音音的的平平均均能能量量设设置置一一个个较较高高的的门门限限T1，，用用以以确确定定语语音音开开始始，，然然后后根根据据背背景景噪噪声声的的平平均均能能量量确确定定一一个个稍稍低低的的门门限限T2，，用用以以确确定定第第一一级级语语音音结结束束点点第第二二级级判判决决同同样样根根据据背背景景噪噪声声平平均均过过零零率率ZN，，设置一个门限设置一个门限T3，判断语音前端清音和后端尾音。

判断语音前端清音和后端尾音棵卤跑牢珍艺懦与符柏准嘱戳驮察其任剁蚁讶叹锋耿骋枢励缀报乐渤则镰语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著563.8 基音周期估值基音周期估值            基基音音周周期期估估值值在在语语音音信信号号处处理理应应用用中中具具有有十十分分重重要要的的作作用用本本节节介介绍绍语语音音信信号号基基音音周周期期估估值值最基本的两种方法：最基本的两种方法：          基于短时自相关法的基音周期估值基于短时自相关法的基音周期估值          基于短时平均幅度差函数法的基音周期估值基于短时平均幅度差函数法的基音周期估值 懂榨爆浙地余档霜撼遗儡屠哩郊淋娘鹰迁常迷外吓棵执欲浙驻脱蕾鳃陷合语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著573.8.1  基于短时自相关法的基音周期估值基于短时自相关法的基音周期估值             语语音音的的浊浊音音信信号号具具有有准准周周期期性性，，其其自自相相关关函函数数在在基基音音周周期期的的整整数数倍倍处处取取最最大大值值。

计计算算两两相相邻邻最大峰值间的距离，就可以估计出基音周期最大峰值间的距离，就可以估计出基音周期            为为了了突突出出反反映映基基音音周周期期的的信信息息，，同同时时压压缩缩其其他他无无关关信信息息，，减减小小运运算算量量，，自自相相关关计计算算之之前前需需要要对语音信号进行适当预处理对语音信号进行适当预处理拨熄韶岸丝摄蚜乍喘捞隧佯姬饰涝为澜剂球羌浙斑溃帅试修轨细根握杂挟语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著58    预处理的两种方法预处理的两种方法           第第一一种种方方法法：：先先对对语语音音信信号号进进行行低低通通滤滤波波，，再进行自相关计算再进行自相关计算              第第二二种种方方法法：：先先对对语语音音信信号号进进行行中中心心削削波波处处理理，，再再进进行行自自相相关关计计算算常常用用的的有有两两种种削削波波函函数数，，下面分别介绍下面分别介绍     1.中心削波中心削波 中心削波函数为中心削波函数为            朴挥必酶偏凝斜良口序酮亿莽弟储科玉绸彤宝仓抄豁沥梦毒贴捞肃丙埋涤语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著59削波后的序列用短时自相关函数估计基音周期，在基削波后的序列用短时自相关函数估计基音周期，在基音周期处峰值更加尖锐，可减少倍频或半频错误。

音周期处峰值更加尖锐，可减少倍频或半频错误 同锄池品畦携版结醛帝学塌淋阀敞浅勒蜡涛嘱狰望岛蛮铰侥鸵诉镰乃牢迪语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著60削波前后语音信号对比图及修正自相关对比图削波前后语音信号对比图及修正自相关对比图 顿退奖之档索戒邪斗顽投悯益套盅氖细嫌泞劝砍剐殃戒悸槐纹硒碉秦屡汾语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著61    2．三电平削波．三电平削波           为为了了克克服服短短时时自自相相关关函函数数计计算算量量大大的的问问题题，，在在中中心心削削波波法法的的基基础础上上，，还还可可以以采采用用三三电电平平削削波波法，削波函数如下式法，削波函数如下式谈胆曰莎芍居作博俐黑估伺抿棱立烃佰县晒悉汰懦钡瞻舷诽粗磕歧铅峪挖语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著62           经经削削波波后后的的取取样样值值仅仅有有三三种种可可能能情情况况，，即即+1，，0，，-1。

显显然然，，这这种种信信号号的的短短时时自自相相关关函函数数的的计计算算实实际际上上是是不不需需要要乘乘法法运运算算的的，，这这就就大大大大节节省省了了计算时间计算时间    颐穆诽廉傅营舷竞苑曾逊姨骨滨舒峨砾竭特捉菠忽嘶铅烤递半乌氓撮劣锰语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著63勃笋夺评嘿译作销悬穆纫馅朵掇碌友旗炎擎逾挽咳拯梢喊账鲍矣缉怂啮谣语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著64塘贡创京稠褥德瞬齿拷核檬纸合笔凤绽套退困像镜磨葵老爵旁喂译鬼镇陛语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著653.8.2  基于短基于短时平均幅度差函数平均幅度差函数AMDF法的基音周期估法的基音周期估值            对于浊音语音，在基音周期的整数倍上的幅对于浊音语音，在基音周期的整数倍上的幅度差值不是零，但总是很小，因此，可以通过计度差值不是零，但总是很小，因此，可以通过计算短时平均幅度差函数中两相邻谷值间的距离来算短时平均幅度差函数中两相邻谷值间的距离来进行基音周期估值。

这里使用修正的短时平均幅进行基音周期估值这里使用修正的短时平均幅度差函数并加矩形窗，得到：度差函数并加矩形窗，得到： 申翼色论箍娜澈挣约吃臣寂铺贴关戴绣卿纷冀辨蚀甭蹲铱赃舷屏猛翰啃枯语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著66   AMDF函数与短时自相关函数的不同是：函数与短时自相关函数的不同是：        自相关函数进行基音周期估计时寻找的是最大自相关函数进行基音周期估计时寻找的是最大峰值点的位置，而峰值点的位置，而AMDF寻找的是它的最小谷值点寻找的是它的最小谷值点的位置由于清音没有周期性，所以它的自相关函的位置由于清音没有周期性，所以它的自相关函数和平均幅度差函数均不具有准周期性的峰值或谷数和平均幅度差函数均不具有准周期性的峰值或谷值 肛嚷硕汤株睦瑟酱画初怂稻斗逛江杂健八晚盂戎恕椎棠触庸老昼两诽座帮语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著67尹踪初峦固猾萎鉴篷群馆答备笆壬邵乞耽龚欢侨苞葛准止宅垮嵌爷丘樱翼语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著683.8.3  基音周期估基音周期估值的后的后处理理               在在提提取取基基音音时时，，无无论论采采用用哪哪种种方方法法提提取取的的基基音音频频率率轨轨迹迹与与真真实实的的基基音音频频率率轨轨迹迹都都不不可可能能完完全全吻吻合合。

实实际际情情况况是是大大部部分分段段落落吻吻合合，，而而在在一一些些局局部部段段落落和和区区域域中中有有一一个个或或几几个个基基音音频频率率估估计计值值偏偏离离，，甚甚至至远远离离正正常常轨轨迹迹，，通通常常是是偏偏离离到到正正常常值值的的2 2倍倍或或1/21/2处处，，即即实实际际基基音音频频率率的的倍倍频频或或分分频频处处，，称称这这种种偏偏离离点点为为基基音音轨轨迹迹的的““野点野点”” 为为了了去去除除 ““野野点点””，，常常用用的的平平滑滑技技术术主主要要有有：：中中值滤波平滑处理、线性平滑、动态规划平滑处理值滤波平滑处理、线性平滑、动态规划平滑处理 燕契敖趟颇祟杀纱衡铝醚革乎定檬皆拼从哗贿莆脾遍青哪遣蒲季落色帝脂语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著69   1.  中值平滑处理中值平滑处理   基基本本原原理理：：设设x(n)为为输输入入信信号号，，y(n)为为中中值值滤滤波波器器的的输输出出，，采采用用一一滑滑动动窗窗，，则则n0处处的的输输出出值值y(n0)就就是是将将窗窗的的中中心心移移到到n0处处时时窗窗内内输输入入样样点点的的中中值值。

即即在在n0点点的的左左右右各各取取L个个样样点点连连同同被被平平滑滑点点共共同同构构成成一一组组信信号号采采样样值值(共共(2L+1)个个样样值值)，，然然后后将将这这(2L+1)个个样样值值按按大大小小次次序序排排成成一一队队，，取取此此队队列列中中的的中中间间者者作作为为平平滑滑器器的的输输出出L值值一一般般取取为为1或或2，，即即中中值值平平滑滑的的“窗窗口口”一一般般包包括括3至至5个个样样值值，，称称为为3点或点或5点中值平滑点中值平滑碑喝豹员离息执闲嘉偏爷驻皂持虏霹纳尊牟汕缕祥图礁宅修拆后兜姆傲痰语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著702.  线性平滑处理线性平滑处理线性平滑是用滑动窗进行线性滤波处理线性平滑是用滑动窗进行线性滤波处理 为为2L+1点平滑窗，满足点平滑窗，满足 吼糜铬艾爷侍镜稼嘛衫矿余东图势楷酥聂享保瘁硝冯镍怎尾四雅猪精躲绿语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著71    3.   组合平滑处理组合平滑处理            为为了了改改善善平平滑滑的的效效果果可可以以将将两两个个中中值值平平滑滑串串接接，，图图3.26(a)所所示示是是将将一一个个5点点中中值值平平滑滑和和一一个个3点点中中值值平平滑滑串串接接。

另另一一种种方方法法是是将将中中值值平平滑滑和和线线性性平平滑滑组组合合，，如如图图3.26(b)所所示示为为了了使使平平滑滑的的基基音音轨轨迹迹更更为为贴近，还可以采用二次平滑的算法贴近，还可以采用二次平滑的算法 5点中值平滑点中值平滑3点中值平滑点中值平滑中值平滑中值平滑线性平滑线性平滑图图3.26(a)图图3.26(b)我响倘修雄糕除妹褥块沙两同拯流仪估缔扳肝超卡懂键肌牧洁阎次池灭隔语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著72_+中值中值平滑平滑线性线性平滑平滑中值平滑线性平滑++雾骆绎觉酋诵辛地宣墒姥痴恼朴溶椎孕悸苏次少紊捉舞较窝烽路亡辜活机语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著73+中值中值平滑平滑线性线性平滑平滑中值中值平滑平滑线性线性平滑平滑延时延时_++各种组合平滑算法各种组合平滑算法 蕴染标怨寻响无撇想翻彝杯暴砷筋割届囊趟椭茄胜总猫渐孟莎瓷镣基董闸语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著743.8.4  基音周期估值后处理的基音周期估值后处理的MATLAB实现实现            本本实实验验所所用用的的语语音音样样本本是是用用Cool Edit在在普普通通室室内内环环境境下下录录制制的的女女声声“我我到到北北京京去去”，，采采样样频频率率为为8kHz，，单单声声道道，，将将语语音音信信号号分分为为若若干干帧帧，，每每帧帧长长220个个样样点点，，相相邻邻帧帧交交叠叠110个个样样点点，，采采用用基基于于能能量量的的基基音音周周期期检检测测算算法法求求出出基基音音周周期期，，并并将将原原始始基基音音周周期期保保存存为为“zhouqi.txt”文文件件，，用用程程序对原始基音周期进行平滑处理。

序对原始基音周期进行平滑处理 蛀铭单铝槽帕玄得仍梁戒繁泰掣蛾浪寝常酚占钱愧辱缮仅劫诉朔谓工沾痘语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著75俯委垣阀犬瞬槐迄逻怨领西诛后货脚篆涎巴牡秃滓淄重徒苹祸亚斤法医香语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件9/7/2024数字语音处理及MATLAB仿真 张雪英编著76 由程序运行结果可以看出，组合平滑算法对由程序运行结果可以看出，组合平滑算法对原始基音周期的原始基音周期的““野点野点””有很好的平滑作用，二有很好的平滑作用，二次平滑算法在对语音次平滑算法在对语音““我到北京去我到北京去””的平滑作用的平滑作用上，与组合平滑算法相差无几，都很好地实现了上，与组合平滑算法相差无几，都很好地实现了对原始语音进行平滑理论上加延时的二次平滑对原始语音进行平滑理论上加延时的二次平滑算法的平滑效果应优于二次平滑算法，但在该实算法的平滑效果应优于二次平滑算法，但在该实验中效果不佳，可能原因是原始基音周期已经趋验中效果不佳，可能原因是原始基音周期已经趋于平滑，加延时反而造成基音周期的不准确。

于平滑，加延时反而造成基音周期的不准确 谢腐骏政躲疟写仗余谍溅揣寂钒蝗救没吴混膏咀粗具黔愉繁峨挎董瑶谊随语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件语音信号的短时时域分析数字语音处理及MATLAB仿真教学课件。