短时平均能量的实现

1、语音信号处理实验报告实验目的：1、熟悉Matlab基本程序的运用；2、充分理解取不同窗长时的语音短时平均能量的变化情况;3、熟悉Matlab编程语言在语音信号处理中的作用；4、能够实现程序的重新编制。实验原理：定义n时刻某语音信号的短时平均能量En为：E 二 堂 x(m) w(n - m) 2 = x(m)w(n - m) 2nm=-3m=n-(N-1)式中N为窗长，可见短时平均能量为一帧样点值的加权平方和。特殊地，当窗函数为矩形窗时，有E = x 2(m)nm=n-(N-1)三、实验要求:1、实验前自己用Cool Edit音频编辑软件录制声音“我到北京 去”，并把它保存为.txt文件;2、编

2、程实现不同矩形窗长N=50、100、400、800的短时平均能量;3、用Matlab画出不同窗长的短时平均能量的图形;4、写出实验报告，分析实验结果。四、实验步骤:1、用Cool Edit读入语音“我到北京去”。设置采样率为8kHz，16位，单声道；2、将读入的语音wav文件保存为txt文件；3、把保存的文件speech.txt读入Matlab；4、对采集到的语音样点值进行分帧；5、对照取不同矩形窗长N的短时平均能量，画出图形。五、实验程序：分帧子程序为fra.m，如下所示：取不同窗长N时的短时能量对照程序为nengliang.m,如下所示:, Edi-t or C: DDcment s an

3、d Sett ingsAdKinisrat orI. - - I _ II II XF_ile Idi +_ Text Go Cell Tool s De Lug Ilesktop jVimlow Help蜜 ? X? 口心.|箫嘻s 驾，胸；危丙 ,曲 事njv :十冒掐| - i. 口 + + 1.1 *疫必适 id=f openC zq.q. -feKt1 , ? rt?):-x=fscarf (f id n %f?):fclose (fid);紧计算N=50!帧移二5Cl时的语音能重s=f r a (50j 50j x)-e2=s.、-ener gy=sum(s2j 2)-subp 1

4、 ot (2, 2 1)-plot(energy)|-wlabel C 顿数】-刃皿1虑时能蚩E)legend N=50,)-axisttO 1500, 0,2*1010)酣算N=ioa,顿移=id。时的语音能量-s=fra(lOOj LOOj-s2=s. 2IT-日riQtgy=sumEL N)16-subplot (2j 2, 2)-plot(energy)-slabel (5 帧数)-yL由bel短时能重I )-legend N=100,)-axisttO.VSO, 0,41010)vscriptLil 9 C ol 13 Editor - C: Docuents and Setting

5、sAdiiiistratorM. -.Fils Edi t Tewt Gc Cell T*als Debug 暴Mkt* Tfiikdow Hlp华X? D a 苗岭冒F* 何IB曲日 气Vf I - 1. 0 _ + I 1.1 _ x 麟必24豺计算职如0,帔移=如。时的语音能重425a=fr a (400,400)26s2=s. 227energysumtESj 2)2Ssubploi: (2j 2j 3)29plot(energy)30觅Lab ml帧数)31ylabel短时能重E )32legend11=如)33gwUj 190, 1勃LET 11)34务计算N=soa,帧移罚叩时的

6、语音能量35s=f ra (800j 800, x)Mi3632=5. ,37ener gy=nijjiL i_32j 2.)38subplot(2, 2S 4)39plot(energy)40slab el顿数)41ylabel短时能量E)42legendN=800J)43axis (0,95, 0, 3*lCiLl)VscriptLn. 9Col 13VE 六、实验结果及分析:)Figiire 1匚Fils Edi + Vi h-a1 IrLE&rt mlw Dskt r.pHlp e q |叫鼠鬼娘坚晟j直| 口005001000帧数Lu啊盟tel00200Lu蜩噩您蜜N=1001500400600帧