短時(shí)平均能量的實(shí)現(xiàn)

語音信號(hào)處理實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)?zāi)康模?、 熟悉Matlab基本程序的運(yùn)用；2、 充分理解取不同窗長時(shí)的語音短時(shí)平均能量的變化情況3、 熟悉Matlab編程語言在語音信號(hào)處理中的作用；4、 能夠?qū)崿F(xiàn)程序的重新編制。實(shí)驗(yàn)原理：定義n時(shí)刻某語音信號(hào)的短時(shí)平均能量En為：E=芳n[x(m)w(n一m)]2=工 [x(m)w(n一m)]2m=s m=n_(N_1)式中N為窗長，可見短時(shí)平均能量為一幀樣點(diǎn)值的加權(quán)平方和。特殊地，當(dāng)窗函數(shù)為矩形窗時(shí)，有E二乙 x2(m)nm=n-(N-1)三、 實(shí)驗(yàn)要求：1、 實(shí)驗(yàn)前自己用CoolEdit音頻編輯軟件錄制聲音“我到北京去”并把它保存為.txt文件；2、 編程實(shí)現(xiàn)不同矩形窗長N=50、100、400、800的短時(shí)平均能量;3、 用Matlab畫出不同窗長的短時(shí)平均能量的圖形；4、 寫出實(shí)驗(yàn)報(bào)告，分析實(shí)驗(yàn)結(jié)果。四、 實(shí)驗(yàn)步驟：1、用CoolEdit讀入語音“我到北京去”設(shè)置采樣率為8kHz,16位，單聲道;2、 將讀入的語音wav文件保存為txt文件；3、 把保存的文件speech.txt讀入Matlab；4、 對(duì)采集到的語音樣點(diǎn)值進(jìn)行分幀；5、 對(duì)照取不同矩形窗長N的短時(shí)平均能量，畫出圖形。五、實(shí)驗(yàn)程序：分幀子程序?yàn)閒ra.m，如下所示：

取不同窗長N時(shí)的短時(shí)能量對(duì)照程序?yàn)閚engliang.m,如下所示:

■Editor--C:\Docu>entsandSettings\Ad>iiiistratorM.-.■問XFilsEditTewtGcCellT*alsDebugIisEktapTfiikdowH^lp-1*X髓?■相IB”宿回?Vf]I-1.0_+I1.1_x疇逆24軸計(jì)算N=40a,幀移=4(JCl時(shí)的語音罪量4□25—a=fra(400,400^)26—s2=s.^2;27—energy^sumtESj2)—2S—subp1o-fc(2j2j3)29—plot(energy)30—Klabel幀歎)31—ylabel「短時(shí)能童E‘)—32—legendC!1=4口『)33—asisftO^190,0,1.5*1031]^34場計(jì)算N=soa,幀移二遡時(shí)的語音能量—35—s=fra(800j800,x)一36—32=5."2;37—已nergy=sijjTLi_32j2.)38—subp1ot⑵2,4)—39—plot(energy)40—slabel「幀數(shù)’)41—ylabel('短時(shí)能量E‘)42—legend「N=800J) 43—axis([0,95,山3+10'11])VscriptLn.9Col13■■■■六、實(shí)驗(yàn)結(jié)果及分析:FilsEditVih-a-IrLE&rtTools VfindowH^lpoeflfe聲頭心⑥原搖"凰l?oeflfe聲頭心⑥原搖"凰l?o1oJI—X50-NLuw攝tel?5o.N=100500 1000幀數(shù)1500500 1000幀數(shù)150000 200 400 600幀數(shù)ii?ii N=4OO--■oii?ii N=4OO--■o050HI50“x10113|~~n~H~1 1 N=30080Ba數(shù)40幀分析：從