pitch.m

来自「语音处理平台 可以分析语音能量 第一振峰频率等数据」· M 代码 · 共 24 行

function [f0] = pitch(len,sr,xin)

% Copyright (c) 1995 Philipos C. Loizou
%


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