語音實驗(語音短時平均能量的實現(xiàn))

1、學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱專業(yè)班級通信0802學號 實驗成績學生姓名同組人姓名實驗日期課程名稱語音信號處理實驗題目實驗一語音短時平均能量的實現(xiàn)一、實驗?zāi)康?、熟悉matble基本程序的運用。2、充分理解取不同窗長時的語音短時平均能量的變化情況。3、熟悉matble編程語言在語音信號處理中的作用。4、能夠?qū)崿F(xiàn)程序的重新編制。二、實驗要求1、實驗前自己用cool edit音頻編輯軟件錄制聲音“我到北京去”，并把它保存為.txt文件。2、編程實現(xiàn)不同矩形窗長n=50、100、400、800

2、的短時平均能量。3、用matlab畫出不同窗長的短時平均能量的圖形。4、寫出實驗報告，分析實驗結(jié)果。三、 實驗步驟1、用cool edit讀入語音“我到北京去”。設(shè)置采樣率為8khz,16位，單聲道。2、將讀入的語音wav文件保存為txt文件。3、把保存的文件speech.txt讀入matlab。四、程序及運行結(jié)果fid=fopen(zqq.txt,rt);x=fscanf(fid,%f);fclose(fid);s=fra(50,50,x);s2=s.2;energy=sum(s2,2);subplot(2,2,1);plot(energy);xlabel(幀數(shù));ylabel(短時能量 e

3、a);legend(n=50);axis(0,1500,0,2*1010);s=fra(100,100,x);s2=s.2;energy=sum(s2,2);subplot(2,2,2);plot(energy);實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗一xlabel(幀數(shù));ylabel(短時能量 eb);legend(n=100);axis(0,750,0,4*1010);s=fra(400,400,x);s2=s.2;energy=sum(s2,2);subplot(2,2,3);plot(en

4、ergy);xlabel(幀數(shù));ylabel(短時能量 ec);legend(n=400);axis(0,190,0,1.5*1011);s=fra(800,800,x);s2=s.2;energy=sum(s2,2);subplot(2,2,4);plot(energy);xlabel(幀數(shù));ylabel(短時能量 ed);legend(n=800);axis(0,95,0,3*1011);實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗一五、思考題1、通過改動不同的窗長和窗移，觀看短時平均能量的變化。

5、答：窗長越小短時平均能量越小，細節(jié)越多，隨著窗長的增大短時平均能量越大，曲線越光滑，個別細節(jié)被平滑掉了。2、在相同的實驗環(huán)境下，把本文中的矩形窗改成hamming窗，來實現(xiàn)語音短時平均能量。實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗（二）思考題五、思考題1、在相同的實驗環(huán)境下，用matlab程序?qū)崿F(xiàn)語音的短時自相關(guān)的圖形，并與修正自相關(guān)進行比較，加深對修正自相關(guān)的理解。答：由于在計算短時自相關(guān)時，窗選語音段為有限長度n，而求和上限是n-1-k，因此當k增加時可用于計算的數(shù)據(jù)就越來越少了，從而導致k增加時

6、自相關(guān)函數(shù)的幅度減小了。采用修正的自相關(guān)函數(shù)可以避免這個問題。2、通過實驗的圖形，試分析它在語音信號處理中的應(yīng)用。答：1）可以作為區(qū)分清音濁音的特征參數(shù)； 2）在信噪比較高的情況下，短時能量還可以作為區(qū)分有聲和無聲的依據(jù)； 3）可以作為輔助的特征參數(shù)用于語音識別中。實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗（三）圖形實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱語音信號處理實驗題目 實驗二 語音修正的短時自相關(guān)的實現(xiàn)1、 實驗?zāi)?/p>

7、的1、熟悉語言修正自相關(guān)的意義。2、充分理解取不同窗長時的語言的修正自相關(guān)的變化情況。3、熟悉matlab編程語言在語言信號處理中的作用。4、 能夠?qū)崿F(xiàn)程序的重新編制。二、實驗要求1、實驗前自己用cool edit音頻編輯軟件錄制聲音，并把它保存為.txt文件。 2、編程實現(xiàn)不同矩形窗長n=320,160,70的短時修正自相關(guān)。3、用matlab畫出短時修正自相關(guān)的圖形。4、寫出實驗報告，分析實驗結(jié)果。3、 實驗步驟1、用cool edit讀入濁音語音，設(shè)置采樣率為8khz，16位，單聲道。2、將讀入的語音wav文件保存為txt文件。3、讀入matlab中，并且對照取不同矩形窗長n的短時修正自

8、相關(guān)函數(shù)，畫出圖形。 四、 程序及運行結(jié)果fid=fopen(voice.txt,rt);b=fscanf(fid,%f);b1=b(1:300);n=150;a=;for k=1:150; sum=0; for m=1:n; sum=sum+b1(m)*b1(m+k-1); end a(k)=sum;endfor k=1:150;endfor k=1:150;a1(k)=a(k)/a(1);endfigure(1);實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗二subplot(3,1,1);plot(a

9、1);xlabel(延時k);ylabel(r(k);legend(n=150);axis(0,320,-0.5,1);b2=b(1:320);n=160;b=;for k=1:160; sum=0; for m=1:n; sum=sum+b2(m)*b2(m+k-1); end b(k)=sum;endfor k=1:160;b1(k)=b(k)/b(1);endfigure(1);subplot(3,1,2);plot(b1);xlabel(延時k);ylabel(r(k);legend(n=80);axis(0,320,-0.5,1);b3=b(1:140);n=70;c=;for k=

10、1:70;sum=0;for m=1:n;sum=sum+b3(m)*b3(m+k-1);end實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗二c(k)=sum;endfor k=1:70;c1(k)=c(k)/c(1);endfigure(1);subplot(3,1,3);plot(c1);xlabel(延時k);ylabel(r(k);legend(n=70);axis(0,320,-0.5,1);實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班級通信0802學號 實驗成績學生

11、姓名同組人姓名實驗日期課程名稱語音信號處理實驗題目實驗三用修正的短時自相關(guān)檢測語言的基音周期1、 實驗?zāi)康?、 熟悉前一個實驗程序以及中心削波的意義。2、用matlab實現(xiàn)用修正的短時自相關(guān)檢測語音的基音周期。3、分析修正的短時自相關(guān)在基音周期檢測中的應(yīng)用。4、能夠?qū)Τ绦蜻M行重新編制。2、 實驗要求1、實驗前自己用cool edit音頻編輯軟件錄制濁音聲音，并把它保存為.txt文件。2、分別取長度n=160和n=n+k的矩形窗函數(shù)作用于語音信號上，首先對其進行中心削波，比較削波前和削波后語音信號波形，并使得削波后的信號進行乘積并求和，計算延遲為0k=th0 a(k)=a(k)-th0; els

12、eif a(k)=(-th0); a(k)=a(k)+th0; else a(k)=0; endendm=max(a);for i=1:l; a(i)=a(i)/m;endsubplot(2,1,2);plot(a,k);axis(0,388,-1,1);title(中心削波后語音波形);xlabel(樣點數(shù));ylabel(幅度);fclose(fid);n=50;a=;for k=1:50;sum=0;for m=1:n; sum=sum+ b(m)*b(m+k-1);enda(k)=sum;endfor k=1:50實驗室名稱指導老師簽名太原理工大學學生實驗報告學院名稱現(xiàn)代科技學院專業(yè)班

13、級學號 實驗成績學生姓名同組人姓名實驗日期課程名稱實驗題目實驗三 b(k)=a(k)/a(1);endfigure(2);subplot(2,1,1);plot(b,k);title(中心削波前修正自相關(guān));xlabel(延時k);ylabel(幅度);axis(0,50,-1,1);n=50;a=;for k=1:50;sum=0;for m=1:n; sum=sum+a(m)*a(m+k-1);enda(k)=sum;endfor k=1:50 c(k)=a(k)/a(1);endsubplot(2,1,2);plot(c,k);title(中心削波后修正自相關(guān));xlabel(延時k);ylabel(幅度);a