www.pudn.com > mel_endpointdetection.rar > teager2.m

data1=double(data1); 
%data=filter([1 -0.9375], 1, data); 
data1=data1/max(data1); 
FrameLen=256; 
FrameInc=128; 
xe=enframe(data1, FrameLen, FrameInc); 
L=256; 
w1=window(@hamming,L);  
[m,n]=size(xe); 
xw=zeros(m,n); 
for i=1:m 
    xw(i,:)=xe(i,:).*w1'; 
    xf(i,:)=fft(xw(i,:),256); 
end 
xfe=abs(xf); 
w=2.5+([1:32]-1)*4; 
for i=1:32 
    for j=1:m 
        Eb(i,j)=w(i)*sum(xfe(j,(1+(i-1)*4+1):(1+(i-1)*4+3+1))')'; 
    end 
end 
figure,plot(sum(Eb));