短時平均能量的實現

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幀