qingzhuo.m

来自「语音信号的清浊音的提取是研究语音信号的一个必经的过程」· M 代码 · 共 30 行

clear;
[x,fs,bits]=wavread('e:\p1.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