本科实验报告
课程名称：
数字语音处理
实验项目：语音短时平均能量的实现
实验地点：
跨越机房
专业班级：学号：
学生姓名：
指导教师：
年月
f实验一语音短时平均能量的实现
一、实验目的
1熟悉Matlab基本程序的运用。2充分理解取不同窗长时的语音短时平均能量的变化情况。3熟悉Matlab编程语言在语音信号处理中的作用。4能够实现程序的重新编制。
二、实验原理
定义
时刻某语音信号的短时平均能量E
为：



E
xmw
m2
xmw
m2
m
m
N1
式中N为窗长，可见短时平均能量为为一帧样点值的甲醛平方和。特殊地，当窗函数为矩形窗时，有


E

x2m
m
N1
三、实验要求
1实验前自己用CoolEdit音频编辑软件录制声音“我到北京去”，并把它保存为txt文件。
2编程实现不同矩形窗长N50、100、200、800的短时平均能量。3用Matlab画出不同窗长的短时平均能量的图形。4写出实验报告，分析实验结果。
四、实验难点
对采集到的语音阳电只进行分帧。
五、实验条件
计算机，Matlab软件，2学时
六、实验步骤
1用Cooledit读入语音“我到北京去”。设置采样率为8kHz，16位，单声道。2将读入的语音wav文件保存为txt文件。3把保存的文件speechtxt读入Matlab。fidfope
speechtxtrt
fxfsca
ffidffclosefid4对采集到得语音样点值进行分帧。
分帧子程序为framfu
ctio
ffrale
i
cxfhfixsizex1le
i
c1fzerosfhle
i1
1whileifh
j1whilejle
fijx
jj1
1e
d
le
i
cii1e
d5对照取不同矩形窗N的短时平均能量，画出图形。取不同窗长N时的短时能量对照程序e
ergymsfra5025xs2s2e
ergysums22subplot221plote
ergyxlabel帧数ylabel短时能量Elege
dN50axis01500010105sfra10050xs2s2e
ergysums22subplot222plote
ergyxlabel帧数ylabel短时能量Elege
dN100axis075002106sfra400200xs2s2e
ergysums22subplot223plote
ergyxlabel帧数ylabel短时能量E
flege
dN400axis019007106sfra800400xs2s2e
ergysums22subplot224plote
ergyxlabel帧数ylabel短时能量Elege
dN800axis095014106
七、实验结果
f本科实验报告
课程名称：
数字语音处理
实验项目：语音修正的短时自相关的实现
实验地点：
跨越机房
专业班级：学号：
学生姓名：
指导教师：
年月
f实验二语音修正的短时自相关的实现
一、实验目的
1、熟悉语音修正自相关的意义。2、充分理解取不同窗长时的语音的修正自相关的变化情况。3、熟悉Matlab编程语音在语音r