语音实验报告1.docVIP

太原理工大学现代科技学院 数字语音处理 课程 实验报告 专业班级 通信14-4 学 号 2014101716 姓 名 杨蓉 指导教师 贾海蓉 太原理工大学现代科技学院实验报告 ……………………………………装………………………………………订…………………………………………线………………………………………实验名称 语音短时平均能量的实现 ……………………………………装………………………………………订…………………………………………线……………………………………… 专业班级 通信14-4 学号 2014101716 姓名 杨蓉 成绩 实验一 语音短时平均能量的实现 一、实验目的 熟悉Matlab基本程序的运用。 充分理解取不同窗长时的语音短时平均能量的变化情况。 熟悉Matlab编程语言在语音信号处理中的作用。 能够实现程序的重新编制。 二、实验原理 定义n时刻某语音信号的短时平均能量En为： 式中N为窗长，可见短时平均能量为一帧样点值的加权平方和。特殊地，当窗函数为矩形窗时，有 三、实验要求 实验前自己用Cool Edit 音频编辑软件录制声音“我到北京去”，并把它保存为.txt文件。 编程实现不同矩形窗长N＝50、100、400、800的短时平均能量。 用Matlab画出不同窗长的短时平均能量的图形。 写出实验报告，分析实验结果。 四、实验步骤 1、用Cool Edit读入语音“我到北京去”。设置采样率为8kHz,16位，单声道。 2、将读入的语音wav文件保存为txt文件。 3、把保存的文件speech.txt读入Matlab。 五、实验设备 计算机 Matlab软件 实验程序及运行结果 fid=fopen(C:\Users\a\Desktop\beijing.txt,rt); x=fscanf(fid,%f); fclose(fid); s=fra(100,100,x); s2=s.^2; energy=sum(s2,2); subplot(2,2,1) plot(energy) xlabel(zhenshu) ylabel(duanshinengliang) legend(N=50) axis([0,500,0,0.5* 10^10]) s=fra(100,100,x); s2=s.^2; energy=sum(s2,2); subplot(2,2,2) plot(energy) xlabel(zhenshu) ylabel(duanshinengliang) legend(N=100) axis([0,500,0,0.5* 10^10]) s=fra(400,400,x); s2=s.^2; energy=sum(s2,2); subplot(2,2,3) plot(energy) xlabel(zhenshu) ylabel(duanshinengliang) legend(N=400) axis([0,100,0,0.2* 10^11]) s=fra(800,800,x); s2=s.^2; energy=sum(s2,2); subplot(2,2,4) plot(energy) xlabel(zhenshu) ylabel(duanshinengliang) legend(N=800) axis([0,60,0,0.4* 10^11]) function f=fra(len,inc,x) fh=fix(((size(x,1)-len)/inc)+1); f=zeros(fh,len); i=1;n=1; while i=fh j=1; while j=len f(i,j)=x(n); j=j+1;n=n+1; end n=n-len+inc; i=i+1; end 思考题 通过改动不同的窗长和窗移，观看短时平均能量的变化。 答：通过实验结果，我们发现，当窗长N越大时，平滑作用非常显著，使得短时能量变化不大，所以反映不出能量的变化。当窗长N越小时，没有平滑作用，反应了能量的快变细节，而看不出包络的变化。所以，窗长的选择很重要，一般选择100-200. 2、在相同的实验环境下，把本文中的矩形窗改成hamming窗，来实现语音短时平均能量。 答：通过实验，可知汉明