用Matlab編程實(shí)現(xiàn)語音信號的短時(shí)分析

1、用Matlab編程實(shí)現(xiàn)語音信號的短時(shí)分析課程設(shè)計(jì)目錄一、目的2二、原理及方法2三、實(shí)驗(yàn)步驟2四、實(shí)驗(yàn)結(jié)論與分析2五、心得體會(huì)4六、附錄4用Matlab編程實(shí)現(xiàn)語音信號的短時(shí)分析一、目的1.在理論學(xué)習(xí)的基礎(chǔ)上，進(jìn)一步地理解和掌握語音信號短時(shí)分析的意義，短時(shí)時(shí)域分析的基本方法。2.進(jìn)一步理解和掌握語音信號短時(shí)平均能量函數(shù)及短時(shí)平均過零數(shù)的計(jì)算方法和重要意義。二、原理及方法一定時(shí)寬的語音信號，其能量的大小隨時(shí)間有明顯的變化。其中清音段（以清音為主要成份的語音段），其能量比濁音段小得多。短時(shí)過零數(shù)也可用于語音信號分析中，發(fā)濁音時(shí)，其語音能量約集中于3kHz以下，而發(fā)清音時(shí)，多數(shù)能量出現(xiàn)在較高頻率上，可

2、認(rèn)為濁音時(shí)具有較低的平均過零數(shù)，而清音時(shí)具有較高的平均過零數(shù)，因而，對一短時(shí)語音段計(jì)算其短時(shí)平均能量及短時(shí)平均過零數(shù)，就可以較好地區(qū)分其中的清音段和濁音段，從而可判別句中清、濁音轉(zhuǎn)變時(shí)刻，聲母韻母的分界以及無聲與有聲的分界。這在語音識(shí)別中有重要意義。三、實(shí)驗(yàn)步驟(1)從網(wǎng)上找了一段錄音，采樣率為8KHZ，量化精度為16比特線性碼。(2)利用設(shè)計(jì)方案中給定的公式分別編程計(jì)算這段語音信號的短時(shí)能量、短時(shí)平均幅度、短時(shí)過零率，然后分別畫出它們的曲線；(3)然后畫出短時(shí)零能比曲線。(4)根據(jù)上述結(jié)果判斷找出其中的一幀濁音信號和一幀清音信號。判斷依據(jù)是，濁音：短時(shí)能量大、短時(shí)平均幅度大、短時(shí)過零率低；清

3、音：短時(shí)能量小、短時(shí)平均幅度小、短時(shí)過零率高。濁音，取13270-13510個(gè)點(diǎn)，清音，取12120-12360個(gè)點(diǎn)。分別計(jì)算他們的短時(shí)自相關(guān)函數(shù)和平均幅度差函數(shù)；(5)根據(jù)圖形，分析濁音段語音的基音周期。四、實(shí)驗(yàn)結(jié)論與分析（1）從圖中明顯可以看出，濁音信號的具有明顯的周期性，其自相關(guān)函數(shù)和平均幅度差函數(shù)也表現(xiàn)出周期性。清音信號稍微差一點(diǎn)，但不是很明顯，主要是由于清音信號的位置找的不是很好。（2）濁音：短時(shí)能量大、短時(shí)平均幅度大、短時(shí)過零率低；清音：短時(shí)能量小、短時(shí)平均幅度小、短時(shí)過零率高。（3）基音提取消除共振峰的影響。從畫出的圖形中間可以看出，沒有加濾波器的聲音信號處理后共振峰的影響很大。

4、但是加一個(gè)60500Hz的帶通濾波器，利用濾波后的信號進(jìn)行基因估計(jì)，這樣可除去大部分共振峰的影響，自相關(guān)函數(shù)和短時(shí)平均幅度差函數(shù)具有更尖銳地峰值，有利于判決地準(zhǔn)確性。（4）畫出的波形如下所示：圖1為沒有加窗nowindow.wav，無濾波，圖2為加窗后的語音信號為 havewindow.wav，有濾波。圖1沒有濾波的第一列的圖形從上到下分別表示：短時(shí)能量、短時(shí)平均幅度、短時(shí)過零率和零能比曲線，第二列從上到下分別為濁音的短時(shí)自相關(guān)函數(shù)、短時(shí)平均幅度差函數(shù)和輕清音的短時(shí)自相關(guān)函數(shù)、短時(shí)平均幅度差函數(shù)。圖1 圖2五、心得體會(huì)在本次的課程設(shè)計(jì)中，通過自己查找資料，修改程序，最后得出了基音頻率。這種以實(shí)

5、驗(yàn)方式學(xué)習(xí)的過程，加深了我對語音信號理論知識(shí)的理解，在設(shè)計(jì)中，也使我獲得了很多新的知識(shí)，熟練對matlab的使用，對MATLAB的應(yīng)用也更加熟練了。總的說來，這個(gè)實(shí)驗(yàn)在周老師的指導(dǎo)下完成，感謝周老師在這個(gè)課程教會(huì)了我們很多關(guān)于語音處理的知識(shí)。六、附錄8 / 8源程序：% 沒有加窗nowindow.wavN=240Y=WAVREAD('nowindow',1 18000);L=length(Y)%30秒，每秒8000個(gè)點(diǎn)，一共240000個(gè)點(diǎn)LL=length(Y)/N %一共 1000幀figure(1)set(1,'Position',10,35,350,65

6、0)%短時(shí)能量 Em=zeros(1,(LL-1)*240);for ii=1:(LL-1)*240, temp=Y(ii:ii+240); Em(ii)=sum(temp.*temp);endsubplot(4,1,1)jj=1:(LL-1)*240; plot(jj, Em,'b');grid% axis(12000,16000,0.15,0.4) %短時(shí)平均幅度Mn=sum(abs(Y)/NMn=zeros(1,(LL-1)*240);for ii=1:(LL-1)*240, temp=Y(ii:ii+240); Mn(ii)=sum(abs(temp)/N;endfig

7、ure(1) subplot(4,1,2)jj=1:(LL-1)*240; plot(jj, Mn,'b');grid% axis(12000,16000,0.02,0.04) %短時(shí)過零率 Zn=zeros(1,(LL-1)*240);for ii=2:(LL-1)*240, temp1=sign(Y(ii:ii+240); temp=sign(Y(ii-1:ii+240-1); Zn(ii)=sum(abs(temp1-temp);endfigure(1) subplot(4,1,3)jj=1:(LL-1)*240; plot(jj, Zn,'b');gri

8、d% axis(12000,16000,0,20)%零能比曲線figure(1)subplot(4,1,4)jj=1:(LL-1)*240; plot(jj, Zn./Em,'b');grid% axis(12000,16000,0,100)% % Wn=60/8000 ,500/8000% n=1:240;% B = FIR1(239,Wn,'bandpass');%濁音，取13270-13510個(gè)點(diǎn) %短時(shí)自相關(guān)函數(shù) temp=Y(13271:13510);% % wc1 = 0.0075; wc2 = 0.9625; % wn = wc1/pi,wc2/p

9、i;% h = fir1(239,wn,blackman(240);% h=reshape(h,240,1);% temp= temp.*abs(h);% % %temp2=reshape(temp,1,240); %temp=temp2.*B;Rn1=zeros(1,240); for nn=1:240, for ii=1:240-nn, Rn1(nn) =Rn1(nn)+ temp(ii)*temp(nn+ii); endend figure(2)set(2,'Position',400,35,350,650) subplot(4,1,1)jj=1:240;plot(jj,

10、 Rn1,'b');grid%axis(9,50,-0.01,0.02)%短時(shí)平均幅度差函數(shù)Yn1=zeros(1,240); for nn=1:240, for ii=1:240-nn, Yn1(nn) =Yn1(nn)+ abs(temp(ii)-temp(nn+ii); endend figure(2) subplot(4,1,2)jj=1:240;plot(jj, Yn1,'b');grid%清音，取12120-12360個(gè)點(diǎn) %短時(shí)自相關(guān)函數(shù) temp=Y(12121:12360); % %temp2=reshape(temp,1,240); %tem

11、p=temp2.*B;Rn2=zeros(1,240); for nn=1:240, for ii=1:240-nn, Rn2(nn) =Rn2(nn)+ temp(ii)*temp(nn+ii); endend figure(2) subplot(4,1,3)jj=1:240;plot(jj, Rn2,'b');grid%axis(9,50,-0.01,0.02)%短時(shí)平均幅度差函數(shù)Yn2=zeros(1,240); for nn=1:240, for ii=1:240-nn, Yn2(nn) =Yn2(nn)+ abs(temp(ii)-temp(nn+ii); enden

12、d figure(2) subplot(4,1,4)jj=1:240;plot(jj, Yn2,'b');grid%加窗后的語音信號為 havewindow.wavN=240Y=WAVREAD('havewindow',1 18000);L=length(Y)%30秒，每秒8000個(gè)點(diǎn)，一共240000個(gè)點(diǎn)LL=length(Y)/N %一共 1000幀figure(1)set(1,'Position',10,35,350,650)%短時(shí)能量 Em=zeros(1,(LL-1)*240);for ii=1:(LL-1)*240, temp=Y(i

13、i:ii+240); Em(ii)=sum(temp.*temp);endsubplot(4,1,1)jj=1:(LL-1)*240; plot(jj, Em,'b');grid% axis(0,(LL-1)*240,0,0.6)% axis(12000,16000,0,0.06) %短時(shí)平均幅度Mn=sum(abs(Y)/NMn=zeros(1,(LL-1)*240);for ii=1:(LL-1)*240, temp=Y(ii:ii+240); Mn(ii)=sum(abs(temp)/N;endfigure(1) subplot(4,1,2)jj=1:(LL-1)*240

14、; plot(jj, Mn,'b');grid% axis(12000,16000,0,0.015) %短時(shí)過零率 Zn=zeros(1,(LL-1)*240);for ii=2:(LL-1)*240, temp1=sign(Y(ii:ii+240); temp=sign(Y(ii-1:ii+240-1); Zn(ii)=sum(abs(temp1-temp);endfigure(1) subplot(4,1,3)jj=1:(LL-1)*240; plot(jj, Zn,'b');grid% axis(12000,16000,0,60)%零能比曲線figure(

15、1)subplot(4,1,4)jj=1:(LL-1)*240; plot(jj, Zn./Em,'b');grid%axis(12000,16000,0,3*107)% % Wn=60/8000 ,500/8000% n=1:240;% B = FIR1(239,Wn,'bandpass');%濁音，取13270-13510個(gè)點(diǎn) %短時(shí)自相關(guān)函數(shù) temp=Y(13271:13510); % %temp2=reshape(temp,1,240); %temp=temp2.*B;Rn1=zeros(1,240); for nn=1:240, for ii=1:2

16、40-nn, Rn1(nn) =Rn1(nn)+ temp(ii)*temp(nn+ii); endend figure(2)set(2,'Position',400,35,350,650) subplot(4,1,1)jj=1:240;plot(jj, Rn1,'b');grid%axis(9,50,-0.01,0.02)%短時(shí)平均幅度差函數(shù)Yn1=zeros(1,240); for nn=1:240, for ii=1:240-nn, Yn1(nn) =Yn1(nn)+ abs(temp(ii)-temp(nn+ii); endend figure(2) subplot(4,1,2)jj=1:240;plot(jj, Yn1,'b');grid%清音，取12120-12360個(gè)點(diǎn) %短時(shí)自相關(guān)函數(shù) temp=Y(12121:12360); % %temp2=reshape(temp,1,240); %temp=temp2.*B;Rn2=zeros(1,240); for nn=1:240, for ii=1:240-nn, Rn2(nn) =Rn2(nn)+ temp(ii)*temp(nn+ii);