- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
华南理工大学
?
?
《语音信号处理》实验报告
实验名称:基音周期估计
?姓名:?
学号:
班级:11级电信6班?
日期:2014年3?月
????
1.实验目的?
本次试验的目的是通过matlab编程,验证课本中基音周期估计的方法,本实验采用的方法是自相关法。?
2.?实验原理?
1、基音周期?
基音是发浊音时声带震动所引起的周期性,而基音周期是指声带震动频率的倒数。基音周期是语音信号的重要的参数之一,它描述语音激励源的一个重要特征,基音周期信息在多个领域有着广泛的应用,如语音识别、说话人识别、语音分析与综合以及低码率语音编码,发音系统疾病诊断、听觉残障者的语音指导等。因为汉语是一种有调语言,基音的变化模式称为声调,它携带着非常重要的具有辨意作用的信息,有区别意义的功能,所以,基音的提取和估计对汉语更是一个十分重要的问题。?
???? 由于人的声道的易变性及其声道持征的因人而异,而基音周期的范围又很宽,而同—个人在不同情态下发音的基音周期也不同,加之基音周期还受到单词发音音调的影响,因而基音周期的精确检测实际上是一件比较困难的事情。基音提取的主要困难反映在:①声门激励信号并不是一个完全周期的序列,在语音的头、尾部并不具有声带振动那样的周期性,有些清音和浊音的过渡帧是很难准确地判断是周期性还是非周期性的。②声道共振峰有时会严重影响激励信号的谐波结构,所以,从语音信号中直接取出仅和声带振动有关的激励信号的信息并不容?易。③语音信号本身是准周期性的(即音调是有变化的),而且其波形的峰值点或过零点受共振峰的结构、噪声等的影响。④基音周期变化范围大,从老年男性的50Hz到儿童和女性的450Hz,接近三个倍频程,给基音检测带来了一定的困难。由于这些困难,所以迄今为止尚未找到一个完善的方法可以对于各类人群(包括男、女、儿童及不向语种)、各类应用领域和各种环境条件情况下都能获得满意的检测结果。?
尽管基音检测有许多困难,但因为它的重要性,基音的检测提取一直是一个研究的课题,为此提出了各种各样的基音检测算法,如自相关函数(ACF)法、峰值提取算法(PPA)、平均幅度差函数(AMDF)法、并行处理技术、倒谱法、SIFT、谱图法、小波法等等。?
2、自相关函数?
对于离散的语音信号x(n),它的自相关函数定义为:??
R(k)=Σx(n)x(n-k),?
如果信号x(n))具有周期性,那么它的自相关函数也具有周期性,而且周期与信号x(n)的周期性相同。自相关函数提供了一种获取周期信号周期的方法。在周期信号周期的整数倍上,它的自相关函数可以达到最大值,因此可以不考虑起始时间,而从自相关函数的第一个最大值的位置估计出信号的基音周期,这使自相关函数成为信号基音周期估计的一种工具。
?3、短时自相关函数?
语音信号是非平稳的信号,所以对信号的处理都使用短时自相关函数。短时自相关函数是在信号的第N个样本点附近用短时窗截取一段信号,做自相关计算所得的结果?
?? Rm(k)=Σx(n)x(n-k)?
式中,n表示窗函数是从第n点开始加入。?
?
3.?实验数据及平台?
本实验所采用的数据语音文件isolated?word.WAV,平台是MATLAB2010。?
4.?实验过程(步骤)?
1、读取wav文件;?
2、利用自相关法进行基音周期估计;
3、去除野点。??????????????
5.?实验结果及讨论?
1、读取wav文件?
function?pitch?
x=wavread(E:\yuuyin\isolated?word.WAV); %读取声音文件??
figure(1);?
stem(x,.);??????????????????? %显示声音信号的波形?
得到的波形如下:?
?
2、利用自相关法进行基音周期估计??
?n=160;?????????????????????????%取20ms的声音片段,即160个点?
?for?m=1:length(x)/n;???????????%对每一帧求短时自相关数??????
?for?k=1:n;?????????
? ?Rm(k)=0;?
????????? ?for?i=(k+1):n;?
????????????? ?Rm(k)=Rm(k)+x(i+(m-1)*n)*x(i-k+(m-1)*n);??????????
?end???????
end?
?????? p=Rm(10:n);???? %防止误判,去掉前边10个数值较大的点?????
?????? [Rmax,N(m)]=max(p);???????%读取第一个自相关函数的最点??
end????????????????????????????%补
您可能关注的文档
最近下载
- 零星维修工程施工技术方案(可编辑Word完整版).pdf
- [2018年最新整理]3电气安装强制性条文(变电站).doc
- 历代名画记翻译历代名画记片段翻译.doc VIP
- 年产5万吨有机肥建设项目可行性研究报告模板-立项拿地.doc
- 流程管理与流程成熟度培训(含附表).pptx VIP
- 全球及中国飞机辅助动力装置行业市场现状供需分析及市场深度研究发展前景及规划可行性分析研究报告(2024-2030).docx
- 血气分析及肺功能检查.pptx VIP
- 五邑大学2022-2023学年第2学期《高等数学(下)》期末试卷(A卷)附标准答案.pdf
- 《亚当理论》.pdf
- 食品检测对食品安全的重要意义研究.docx
文档评论(0)