Code covered by the BSD License
function [f0] = pitch3(len,sr,xin) xin=hamming(len).*xin; cn1=rceps(xin); LF=floor(sr/500); HF=floor(sr/70); cn=cn1(LF:HF); [mx_cep ind]=max(cn); if mx_cep > 0.09 & ind >LF f0= sr/(LF+ind); else f0=0; end
Contact us at files@mathworks.com