清浊音检测.txt

欢迎光临程序员联合开发网(www.pudn.com)

clear; 
[x,fs,bits]=wavread('2.wav'); 
 
N=256; 
n=fix(length(x)/N); 
 
t=[]; 
for m=1:n 
     
xd=x((m-1)*N+1:m*N);  
 
[C,L]=wavedec(xd,4,'db4'); 
ca=appcoef(C,L,'db4'); 
[cd1,cd2,cd3,cd4]=detcoef(C,L,[1,2,3,4]); 
E1=sum(cd1.*cd1)/length(cd1); 
E2=sum(cd2.*cd2)/length(cd2); 
E3=sum(cd3.*cd3)/length(cd3); 
E4=sum(cd4.*cd4)/length(cd4); 
E=[1,E2/E1,E3/E1,E4/E1]; 
Em=max(E); 
 
 if (Em==1)&(E4/E1<0.9)  
     f=1; 
 else  
     f=0; 
  t=[t,f]; 
end;  
 
end