语音实验语音短时平均能量的实现

1、学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验室名称指导老师签名太原理工大学学生实验报告学院名称专业班级通信0802学号 实验成绩学生姓名同组人姓名实验日期课程名称语音信号处理实验题目实验一语音短时平均能量的实现一、实验目的1、熟悉matble基本程序的运用。2、充分理解取不同窗长时的语音短时平均能量的变化情况。3、熟悉matble编程语言在语音信号处理中的作用。4、能够实现程序的重新编制。二、实验要求1、实验前自己用cool edit音频编辑软件录制声音“我到北京去”，并把它保存为.txt文件。2、编程实现不同矩形窗长n=50、100、400、800

2、的短时平均能量。3、用matlab画出不同窗长的短时平均能量的图形。4、写出实验报告，分析实验结果。三、 实验步骤1、用cool edit读入语音“我到北京去”。设置采样率为8khz,16位，单声道。2、将读入的语音wav文件保存为txt文件。3、把保存的文件speech.txt读入matlab。四、程序及运行结果fid=fopen('zqq.txt','rt');x=fscanf(fid,'%f');fclose(fid);s=fra(50,50,x);s2=s.2;energy=sum(s2,2);subplot(2,2,1);plot(en

3、ergy);xlabel('帧数');ylabel('短时能量 ea');legend('n=50');axis(0,1500,0,2*1010);s=fra(100,100,x);s2=s.2;energy=sum(s2,2);subplot(2,2,2);plot(energy);实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验一xlabel('帧数');ylabel('短时能量 eb');legend('n=100&

4、#39;);axis(0,750,0,4*1010);s=fra(400,400,x);s2=s.2;energy=sum(s2,2);subplot(2,2,3);plot(energy);xlabel('帧数');ylabel('短时能量 ec');legend('n=400');axis(0,190,0,1.5*1011);s=fra(800,800,x);s2=s.2;energy=sum(s2,2);subplot(2,2,4);plot(energy);xlabel('帧数');ylabel('短时能量 ed&

5、#39;);legend('n=800');axis(0,95,0,3*1011);实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验一五、思考题1、通过改动不同的窗长和窗移，观看短时平均能量的变化。答：窗长越小短时平均能量越小，细节越多，随着窗长的增大短时平均能量越大，曲线越光滑，个别细节被平滑掉了。2、在相同的实验环境下，把本文中的矩形窗改成hamming窗，来实现语音短时平均能量。实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验

6、日期课程名称实验题目实验（二）思考题五、思考题1、在相同的实验环境下，用matlab程序实现语音的短时自相关的图形，并与修正自相关进行比较，加深对修正自相关的理解。答：由于在计算短时自相关时，窗选语音段为有限长度n，而求和上限是n-1-k，因此当k增加时可用于计算的数据就越来越少了，从而导致k增加时自相关函数的幅度减小了。采用修正的自相关函数可以避免这个问题。2、通过实验的图形，试分析它在语音信号处理中的应用。答：1）可以作为区分清音浊音的特征参数； 2）在信噪比较高的情况下，短时能量还可以作为区分有声和无声的依据； 3）可以作为辅助的特征参数用于语音识别中。实验室名称指导老师签名太原理工大学

7、学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验（三）图形实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称语音信号处理实验题目 实验二 语音修正的短时自相关的实现1、 实验目的1、熟悉语言修正自相关的意义。2、充分理解取不同窗长时的语言的修正自相关的变化情况。3、熟悉matlab编程语言在语言信号处理中的作用。4、 能够实现程序的重新编制。二、实验要求1、实验前自己用cool edit音频编辑软件录制声音，并把它保存为.txt文件。 2、编程实现不同矩形窗长n=320,

8、160,70的短时修正自相关。3、用matlab画出短时修正自相关的图形。4、写出实验报告，分析实验结果。3、 实验步骤1、用cool edit读入浊音语音，设置采样率为8khz，16位，单声道。2、将读入的语音wav文件保存为txt文件。3、读入matlab中，并且对照取不同矩形窗长n的短时修正自相关函数，画出图形。 四、 程序及运行结果fid=fopen('voice.txt','rt');b=fscanf(fid,'%f');b1=b(1:300);n=150;a=;for k=1:150; sum=0; for m=1:n; sum=su

9、m+b1(m)*b1(m+k-1); end a(k)=sum;endfor k=1:150;endfor k=1:150;a1(k)=a(k)/a(1);endfigure(1);实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验二subplot(3,1,1);plot(a1);xlabel('延时k');ylabel('r(k)');legend('n=150');axis(0,320,-0.5,1);b2=b(1:320);n=160;b=;for k=1:

10、160; sum=0; for m=1:n; sum=sum+b2(m)*b2(m+k-1); end b(k)=sum;endfor k=1:160;b1(k)=b(k)/b(1);endfigure(1);subplot(3,1,2);plot(b1);xlabel('延时k');ylabel('r(k)');legend('n=80');axis(0,320,-0.5,1);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实验室名称指导老师签

11、名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验二c(k)=sum;endfor k=1:70;c1(k)=c(k)/c(1);endfigure(1);subplot(3,1,3);plot(c1);xlabel('延时k');ylabel('r(k)');legend('n=70');axis(0,320,-0.5,1);实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级通信0802学号 实验成绩学生姓名同组人姓名实验日期课程名称语音信号处理实验题目实验

12、三用修正的短时自相关检测语言的基音周期1、 实验目的1、 熟悉前一个实验程序以及中心削波的意义。2、用matlab实现用修正的短时自相关检测语音的基音周期。3、分析修正的短时自相关在基音周期检测中的应用。4、能够对程序进行重新编制。2、 实验要求1、实验前自己用cool edit音频编辑软件录制浊音声音，并把它保存为.txt文件。2、分别取长度n=160和n'=n+k的矩形窗函数作用于语音信号上，首先对其进行中心削波，比较削波前和削波后语音信号波形，并使得削波后的信号进行乘积并求和，计算延迟为0<k<160时的相关值，并且用得到的相关值来检测语音信号的基音周期，用matla

13、b画出图形。3、写出实验报告，分析实验结果。3、 实验步骤1、用cool edit读入浊音语音。设置采样率为8khz，16位，单声道。2、将读入的语音wav文件保存为txt文件。3、编程实现中心削波前后的语音波形。四、程序及运行结果fid=fopen('voice.txt','rt');a,count=fscanf(fid,'%f',1,inf);l=length(a);m=max(a);for i=1:l a(i)=a(i)/m;endm=max(a);n=min(a);ht=(m+n)/2;for i=1:l; a(i)=a(i)-ht;en

14、dfigure(1);subplot(2,1,1);plot(a,'k');axis(0,388,-1,1);title('中心削波前语音波形');实验室名称指导老师签名太原理工大学学生实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验三xlabel('样点数');ylabel('幅度');coeff=0.7;th0=max(a)*coeff;for k=1:l; if a(k)>=th0 a(k)=a(k)-th0; elseif a(k)<=(-th0); a(k)=a

15、(k)+th0; else a(k)=0; endendm=max(a);for i=1:l; a(i)=a(i)/m;endsubplot(2,1,2);plot(a,'k');axis(0,388,-1,1);title('中心削波后语音波形');xlabel('样点数');ylabel('幅度');fclose(fid);n=50;a=;for k=1:50;sum=0;for m=1:n; sum=sum+ b(m)*b(m+k-1);enda(k)=sum;endfor k=1:50实验室名称指导老师签名太原理工大学学生

16、实验报告学院名称现代科技学院专业班级学号 实验成绩学生姓名同组人姓名实验日期课程名称实验题目实验三 b(k)=a(k)/a(1);endfigure(2);subplot(2,1,1);plot(b,'k');title('中心削波前修正自相关');xlabel('延时k');ylabel('幅度');axis(0,50,-1,1);n=50;a=;for k=1:50;sum=0;for m=1:n; sum=sum+a(m)*a(m+k-1);enda(k)=sum;endfor k=1:50 c(k)=a(k)/a(1);endsubplot(2,1,2);plot(c,'k');title('中心削波后修正自相关&