語音訊號短時域分析之短時平均能量（四）

標籤：語音訊號短時平均能量

       由於語音訊號的能量隨時間而變化，清音和濁音之間的能量差別相當顯著，因此對短時能量和短時平均幅度進行分析，可以描述語音的這種特徵變換情況。定義n時刻某語音訊號的短時平均能量E為：

                             

式中，N為窗長，可見短時能量為一幀樣點值的加權平方和。當窗函數為矩形窗時

                                                                                       

        短時平均能量用途：

       (1) 可以作為區分清音和濁音的特徵參數；

       (2) 在信噪比比較高的情況下，短時能量可以作為區分有聲和無聲的依據；

       (3) 可以作為輔助的特徵參數用於語音辨識中。

       其中短時能量matlab代碼如下： 

x = wavread('beijing.wav');%計算N=50，幀移=20時的語音能量s=fra(50,20,x);               s2=s.^2;                    %一幀內各樣點的能量energy=sum(s2,2);            %求一幀能量subplot(2,2,1);               %定義畫圖數量和布局 plot(energy);                %畫N=50時的語音能量圖xlabel('幀數');               %橫座標ylabel('短時能量 E');         %縱座標legend('N=50');              %曲線標識axis([0,1500,0,2*10^10]);      %定義橫縱座標範圍

取不同幀時能量如：

                                                                          

其中，fra()為分幀函數，matlab代碼如下：

function f=fra(len,inc,x)fh=fix(((size(x,1)-len)/inc)+1)f=zeros(fh,len);i=1;n=1;while i<=fh    j=1;    while j<=len        f(i,j)=x(n);        j=j+1;n=n+1;    end    n=n-len+inc;    i=i+1;end