语音短时平均能量

太原理工大学语音信号处理 课程 实验报告专业班级 学 号 姓 名 指导教师 太原理工大学现代科技学院实验报告装订线实验名称 语音短时平均能量 同组人 专业班级 学号 姓名 成绩 一、实验目的： 1.熟悉Matlab基本程序的运用。 2.充分理解取不同窗长时的语音短时平均能量的变化情况。 3.熟悉Matlab编程语言在语音信号处理中的作用。 4.能够实现程序的重新编制。 二、实验原理： 定义n时刻某语音信号的短时平均能量为： 式中N为窗长，可见短时平均能量为为一帧样点值的甲醛平方和。特殊地，当窗函数为矩形窗时，有 三、实验要求： 1.实验前自己用Cool Edit音频编辑软件录制声音“我到北京去”，并把它保存为.txt文件。 2.编程实现不同矩形窗长N=50、100、200、800的短时平均能量。 3.用Matlab画出不同窗长的短时平均能量的图形。 4.写出实验报告，分析实验结果。 四、实验难点： 对采集到的语音样点进行分帧。 装订线五、实验条件： 计算机 Matlab软件 ，2学时 六、实验步骤 1、用Cooledit读入语音“我到北京去”。设置采样率为8kHz，16位，单声道。 2、将读入的语音wav文件保存为“我到北京去.txt”文件。 装订线3、把保存的文件“我到北京去.txt”读入Matlab。 4、对采集到得语音样点值进行分帧。 5、对照取不同矩形窗N的短时平均能量，画出图形。取不同窗长N时的短时能量对照程序 七、实验源程序及结果图 x=wavread( .wav); %x=fscanf(fid,% f); %fclose(fid); s=fra(50,25,x) s2=s.2; energy=sum(s2,2) subplot(2,2,1) plot(energy); xlabel(帧数) ylabel(短时能量 E) legend(N=50) %axis(0,1500,0,10*105) s=fra(100,50,x) s2=s.2; energy=sum(s2,2) subplot(2,2,2) plot(energy); xlabel(帧数) ylabel(短时能量 E) legend(N=100) %axis(0,750,0,2*106) s=fra(400,200,x) s2=s.2; energy=sum(s2,2) subplot(2,2,3) plot(energy); xlabel(帧数) ylabel(短时能量 E) legend(N=400) %axis(0,190,0,7*106) s=fra(800,400,x) s2=s.2; energy=sum(s2,2) subplot(2,2,4) plot(energy); xlabel(帧数) ylabel(短时能量 E) legend(N=800) %axis(0,95,0,14*106) 装订线八、思考题 1、通过改动不同的窗长和窗移，观看短时平均能量的变化。 答：通过实验结果，我们发现，当窗长N越大时，平滑作用非常显著，使得短时能量变化不大，所以反映不出能量的变化。当窗长N越小，没有平滑作用，反映了能量的快变细节，而看不出包络的变化。所以窗长的选择很重要，一般选择100-200。 2、在相同的实验环境下，把本文中的矩形窗改成hamming窗，来实验语音短时平均能量。 答：通过实验，可知汉明窗的主瓣宽度较宽，但它的旁瓣衰减较大，具有更平滑的低通特性，能够在较高的程度上反映短时信号的频率特性。而矩形窗的主瓣宽度小，频谱泄漏严重。所以一般选用汉明窗好。 九、实验心得 本实验是通过MATLAB软件来实现的，目的在于观察在不同窗长时的语音短时平均能量的变化情况，从而更加熟悉MATLAB编程语言在语音信号处理中的作用，以及窗长对语音信号处理时的影响。 试验中，我知道窗长的选择很重要，一般选择100-200。 本实验是通过MATLAB软件来实现的，目的在于观察在不同窗长时的语音短时平均能量的变化情况，从而更加熟悉MATLAB编程语言，也懂得了窗长的重要性。 试验使我们会更加了解短时平均能量的主要用途。在做实验的过程中最重要的是程序的正确编写，尤其是MATLAB必须不