語音短時平均能量

1、太原理工大學(xué)語音信號處理 課程 實驗報告專業(yè)班級 學(xué) 號 姓 名 指導(dǎo)教師 推薦精選裝訂線實驗名稱 語音短時平均能量 同組人 專業(yè)班級 學(xué)號 姓名 成績 一、實驗?zāi)康模?1.熟悉Matlab基本程序的運用。 2.充分理解取不同窗長時的語音短時平均能量的變化情況。 3.熟悉Matlab編程語言在語音信號處理中的作用。 4.能夠?qū)崿F(xiàn)程序的重新編制。 二、實驗原理： 定義n時刻某語音信號的短時平均能量為： 式中N為窗長，可見短時平均能量為為一幀樣點值的甲醛平方和。特殊地，當(dāng)窗函數(shù)為矩形窗時，有 三、實驗要求： 1.實驗前自己用Cool Edit音頻編輯軟件錄制聲音“我到北京去”，并把它保存為.txt

2、文件。 2.編程實現(xiàn)不同矩形窗長N=50、100、200、800的短時平均能量。 3.用Matlab畫出不同窗長的短時平均能量的圖形。 4.寫出實驗報告，分析實驗結(jié)果。 四、實驗難點： 對采集到的語音樣點進行分幀。 裝訂線五、實驗條件： 計算機 Matlab軟件 ，2學(xué)時 六、實驗步驟 1、用Cooledit讀入語音“我到北京去”。設(shè)置采樣率為8kHz，16位，單聲道。 2、將讀入的語音wav文件保存為“我到北京去.txt”文件。 推薦精選 裝訂線3、把保存的文件“我到北京去.txt”讀入Matlab。 4、對采集到得語音樣點值進行分幀。 5、對照取不同矩形窗N的短時平均能量，畫出圖形。取不同

3、窗長N時的短時能量對照程序 七、實驗源程序及結(jié)果圖 x=wavread(' .wav'); %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('幀數(shù)') ylabel('短時能量 E') legend('N=50') %axis(0,1500,0,10*105) s=fra(100,50,x) s2=s.2; energy=sum(s2,2)

4、subplot(2,2,2) plot(energy); xlabel('幀數(shù)') ylabel('短時能量 E') legend('N=100') %axis(0,750,0,2*106) s=fra(400,200,x) s2=s.2; energy=sum(s2,2) 推薦精選 subplot(2,2,3) plot(energy); xlabel('幀數(shù)') ylabel('短時能量 E') legend('N=400') %axis(0,190,0,7*106) s=fra(800,400

5、,x) s2=s.2; energy=sum(s2,2) subplot(2,2,4) plot(energy); xlabel('幀數(shù)') ylabel('短時能量 E') legend('N=800') %axis(0,95,0,14*106) 推薦精選 裝訂線八、思考題 1、通過改動不同的窗長和窗移，觀看短時平均能量的變化。 答：通過實驗結(jié)果，我們發(fā)現(xiàn)，當(dāng)窗長N越大時，平滑作用非常顯著，使得短時能量變化不大，所以反映不出能量的變化。當(dāng)窗長N越小，沒有平滑作用，反映了能量的快變細(xì)節(jié)，而看不出包絡(luò)的變化。所以窗長的選擇很重要，一般選擇100-2

6、00。 2、在相同的實驗環(huán)境下，把本文中的矩形窗改成hamming窗，來實驗語音短時平均能量。 答：通過實驗，可知漢明窗的主瓣寬度較寬，但它的旁瓣衰減較大，具有更平滑的低通特性，能夠在較高的程度上反映短時信號的頻率特性。而矩形窗的主瓣寬度小，頻譜泄漏嚴(yán)重。所以一般選用漢明窗好。 九、實驗心得 本實驗是通過MATLAB軟件來實現(xiàn)的，目的在于觀察在不同窗長時的語音短時平均能量的變化情況，從而更加熟悉MATLAB編程語言在語音信號處理中的作用，以及窗長對語音信號處理時的影響。 試驗中，我知道窗長的選擇很重要，一般選擇100-200。 本實驗是通過MATLAB軟件來實現(xiàn)的，目的在于觀察在不同窗長時的語音短時平均能量的變化情況，從而更加熟悉MATLAB編程語言，也懂得了窗長的重要性。 試驗使我們會更加了解短時平均能量的主要用途。在做實驗的過程中最重要的是程序的正確編寫，尤其是MATLAB必須不能有絲毫偏差，實驗全程都必須細(xì)心認(rèn)真，我們學(xué)習(xí)、生活