语音实验报告2.docVIP

太原理工大学现代科技学院 数字语音处理 课程 实验报告 专业班级 通信14-4 学 号 2014101716 姓 名 杨蓉 指导教师 贾海蓉 太原理工大学现代科技学院实验报告 ……………………………………装………………………………………订………………………………………… ……………………………………装………………………………………订…………………………………………线……………………………………… 专业班级 通信14-4 学号 2014101716 姓名 杨蓉 成绩 实验二：语音修正的短时自相关的实现 一、实验目的 熟悉语音修正自相关的意义。 充分理解取取不同窗长时的语音的修正自相关的变化情况。 熟悉Matlab编程语言在语音信号处理中的作用。 能够实现对程序的重新编制。 二、实验原理 对于语音来说，采用短时分析方法，语音短时自相关函数为 但是，在计算短时自相关时，窗选语音段为有限长度N，而求和上限为N-1-k，因此当k增加时可用于计算的数据就越来越少了，从而导致k增加时自相关函数的幅度减小。为了解决这个问题，提出了语音修正的短时自相关。 修正的短时自相关函数，其定义如下： （1式） 若令，代入上式中得到 （2式） 定义： 则有： （3式） 式中，K为k的最大值，即0≤k≤K。 （4式） 由4式可知，要使为非零值，必须使，考虑到，可得，故3式可以写成： 三、实验要求 实验前自己用Cool Edit 音频编辑软件录制声音，并把它保存为.txt文件. 编程实现不同矩形窗长N＝320，160，70的短时修正自相关。 用Matlab画出短时修正自相关的图形。 写出实验报告，分析实验结果。 四、实验设备 计算机 Matlab软件 实验程序及运行结果 fid=fopen(zhuoyin.txt,rt);b=fscanf(fid,%f);?b1=b(1:320);?N=320;?A=[];?for?k=1:320;?????sum=0;?for?m=1:N-k+1;??????sum=sum+b1(m)*b1(m+k-1);?end?A(k)=sum;?end?for?k=1:320?????A1(k)=A(k)/A(1);?end?figure(1)?subplot(3,1,1)?plot(A1);?xlabel(yan?shi?k)?ylabel(R(k))?legend(N=320)?axis([0,320,-1,1]);?b2=b(1:320);??N=160;?B=[];?for?k=1:160;?sum=0;?for?m=1:N;??sum=sum+b2(m)*b2(m+k-1);?end??B(k)=sum;?end?for?k=1:160;?B1(k)=B(k)/B(1);?end?figure(1)?subplot(3,1,2)?plot(B1);?xlabel(yan?shi?k)?ylabel(R(k))?legend(N=160)axis([0,160,-1.5,1.5]);?b3=b(1:140);?N=70;?C=[];?for?k=1:70;?sum=0;?for?m=1:N;??sum=sum+b3(m)*b3(m+k-1);?end?C(k)=sum;?end?for?k=1:70C1(k)=C(k)/C(1);??end?figure(1)?subplot(3,1,3)?plot(C1);??xlabel(yan?shi?k)?ylabel(R(k))?legend(N=70)?axis([0,70,-1.5,1.5]); 六、思考题 1、在相同的实验环境下，用Matlab程序实现语音的短时自相关的图形，并与修正自相关进行比较，加深对修正自相关的理解。 答：通过采用短时自相关，可以看出：如果长基音周期用窄的窗，将得不到预期的基音 周期；但是如果是短的基音周期用长的窗，自相关函数将对多个基音周期作平均计算， 从而模糊语音的短时特性，这是胡希望的。采用短时的修正自相关函数可以采用较窄的 窗，同时避免了短时自相关函数随k值增加而衰减的不足。 2、改变不同的录制对象，比如清音，分析语音短时修正自相关的图形将发生什么