function short_time_analysis()
%該函數(shù)的功能是將語音信號進行分幀,frame_num為所分幀的數(shù)目,xn為二維向量,存放個幀數(shù)據(jù)
%Number of samples: 17280
%Sampling frequency: 8000 Hz
%精度16位
%窗函數(shù)為海明窗Hamming Window
x=wavread('01.wav');           %讀取數(shù)據(jù)
%plot(x);
%x=[1,2,-1,1,1,-1,-1,0,1,1,1,1,3,1,1,1,1,1,6,-1,1,1,1,-1,1,1,1,1,-1,1,2,3];%測試向量
%x=x';
N=160;                          %設(shè)定幀長
M=160*1/4;                      %設(shè)定個幀交疊部分的長度
xn=enframe(x,N,M);              %分幀函數(shù)   注：要根據(jù)具體情況進行分幀 20ms
[En,Mn]=energy(xn);             %短時能量En及短時平均幅度函數(shù)Mn 
Zn=zero_crossing_rate(xn);      %求短時過零率
Zn2=zero_crossing_rate2(xn);    %設(shè)置門限的短時過零率，具有更強的抗噪性
Xn=fft(xn);                     %求各幀的短時DFT
N=400;                          %設(shè)定幀長
M=400*1/4;                      %設(shè)定個幀交疊部分的長度
xn=enframe(x,N,M);              %分幀函數(shù)   注：根據(jù)自相關(guān)的要求進行分幀  50ms
[Rn,Fn]=auto_correlation(xn);   %短時自相關(guān)函數(shù)Rn和短時平均幅度差函數(shù)AMDF Fn

subplot(2,2,1), plot(x);
subplot(2,2,2), plot(En);
subplot(2,2,3), plot(Mn);
subplot(2,2,4), plot(Zn);
%subplot(2,3,4), plot(Rn);
%subplot(2,3,5), plot(Fn)