太原理工大学语音信号实验报告.doc

太原理工大学现代科技学院 语音信号处理 课程 实验报告 太原理工大学现代科技学院实验报告装订线实验名称 语音短时平均能量的实现 同组人 专业班级 学号 姓名 成绩 一、实验目的 1、熟悉matlab基本程序的运用。 2、充分理解取不同窗长时的语音短时平均能量的变化情况。 3、熟悉matlab编程语言在语音信号处理中的作用。 4、能够实现程序的重新编制。 二、实验要求 1、实验前自己用cool edit音频编辑软件录制声音“我到北京去”，并把它保存为.txt文件。 2、编程实现不同矩形窗长n=50、100、400、800的短时平均能量。 3、用matlab画出不同窗长的短时平均能量的图形。 4、写出实验报告，分析实验结果。 三、实验步骤 1、用cool edit读入语音“我到北京去”。设置采样率为8khz,16位，单声道。 2、将读入的语音wav文件保存为txt文件。 3、把保存的文件speech.txt读入matlab。 四、程序及运行成果 fid=fopen(zxq.txt,rt); 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,1000,0,4* 108) 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,500,0,8* 108) 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,150,0,3* 109) 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,75,0,6* 109) 五、思考题 装订线1、通过改动不同的窗长和窗移，观看短时平均能量的变化。 装订线答：窗长越小短时平均能量越小，细节越多，随着窗长的增大短时平均能量越大，曲线越光滑，个别细节被平滑掉了。 2、在相同的实验环境下，把本文中的矩形窗改成hamming窗，来实现语音短时平均能量。 装订线实验名称 语音修正的短时自相关的实现 同组人 专业班级 学号 姓名 成绩 一、实验目的 1、熟悉语言修正自相关的意义。 2、充分理解取不同窗长时的语言的修正自相关的变化情况。 3、熟悉matlab编程语言在语言信号处理中的作用。 4、 能够实现程序的重新编制。 二、实验要求 1、实验前自己用cool edit音频编辑软件录制声音，并把它保存为.txt文件。 2、编程实现不同矩形窗长n=320,160,70的短时修正自相关。 3、用matlab画出短时修正自相关的图形。 4、写出实验报告，分析实验结果。 三、实验步骤 1、用cool edit读入浊音语音，设置采样率为8khz，16位，单声道。 2、将读入的语音wav文件保存为txt文件。 3、读入matlab中，并且对照取不同矩形窗长n的短时修正自相关函数，画出图形。 四、程序及运行结果 fid=fopen(zhuoyin.txt,rt); b=fscanf(fid,%f); b1=b(1:320); n=160; a=; for k=1:160; sum=0; for m=1:n; sum=sum+b1(m)*b1(m+k-1); end a(k)=sum; end for k=1:160 a1(k)=a(k)/a(1); end figure(1) subplot(3,1,1) plot(a1); ylabel(r(k) legend(n=160) axis(0,320,-1,1); b2=b(1:200); n=100; 装订线b=; for k=1:100; sum=0; for m=1:n; sum=sum+b2(m)*b2(m+k-1); end b(k)=sum; end for k=1:100 b1(k)=b(k)/b(1); end figure(1) subplot(3,1,2) plot(b1); ylabel(r(k) legend(n=100) axis(0,320,-1,1); l=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); e