语音实验报告1.docx

太原理工大学现代科技学院
数字语音处理课程实验报告
专业班级 通信14-4
学 号 2014101716
姓 名 杨蓉
指导教师贾海蓉
: 实验名称 同组人
j 专业班级通信14-4 学号 2014101716 姓名杨蓉 成绩太原理工大学现代科技学院
数字语音处理课程实验报告
专业班级 通信14-4
学 号 2014101716
姓 名 杨蓉
指导教师贾海蓉
: 实验名称 同组人
j 专业班级通信14-4 学号 2014101716 姓名杨蓉 成绩
[ 实验一语音短时平均能量的实现
:一、实验目的
i 1、熟悉Mat lab基本程序的运用。
i 2、充分理解取不同窗长时的语音短时平均能量的变化情况。
: 3、熟悉Matlab编程语言在语音信号处理中的作用。
: 4、能够实现程序的重新编制。
: 二、实验原理
: 定义n时刻某语音信号的短时平均能量En为：
: —» «
^, = E [x(w)w(n-m)]: = 2 [x(m)w(n - m)]:
• m «-se iw-n —(3T—1)
4 式中N为窗长，可见短时平均能量为一帧样点值的加权平方和。特殊地，当窗函数为矩形 订
: 窗时，有
:E..=% r(M)
: 三、实验要求
: 1、实验前自己用Cool Edit音频编辑软件录制声音“我到北京去”，
: 文件。
； 2、编程实现不同矩形窗长N=50、100、400、800的短时平均能量。
窜 3、用Mat lab画出不同窗长的短时平均能量的图形。
: 4、写出实验报告，分析实验结果。
: 四、实验步骤
: 1、用Cool Edit读入语音“我到北京去”。设置采样率为8kHz, 16位，单声道。
: 2、将读入的语音wav文件保存为txt文件。
: 3、把保存的文件speech, txt读入Mat lab □
: 五、实验设备
' 计算机 Mat lab软件
六、实验程序及运行结果
fid=fbpen('C :\U sers\a\Desktop\beij ing. txt', Tf);
x=fscanf(fid；%f);
fclose(fid);
s=fra(100,100,x);
s2=;
energy=sum(s2,2);
subplot(2,2,l)
plot(energy)
xlabel('zhenshu')
ylabel('duanshinengliang')
legend('N=50‘)
axis([0,500,0,* 10A10])
s=fra( 100,100,x);
s2=;
energy=sum(s2,2);
subplot(2,2,2)
plot(energy)
xlabel('zhenshu)
ylabel('duanshinengliang)
legend('N=100')
axis([0,500,0,* 10A10])
s=fra(400,400,x);
s2=;
energy=sum(s2,2);
subplot(2,2,3)
plot(energy)
xlabel('zhenshu')
ylabel('duanshinen