基音周期提取 matlab

自己编写的用自相关法提取基音周期，最后得到的是平均基音周期

function nmax=find_maxn（r）
%寻找峰值最大的n值及基音周期
%r自相关序列
%maxn为峰值最大的n
zer=find（r==0）; %找第一个零点如果存在
jiaocha=0; %找第一近零点
ii=1;
while （jiaocha<=0）
if（r（ii）>0 && r（ii+1）<0 && （ii+1）jiaocha=ii;
end

ii=ii+1;
if ii==length（r） %没有找到符合要求的点
jiaocha=1;
end
end

if length（zer）>0 %检查是否存在零点
if zer（1）jiaocha=zer（1）;
end
end

r（1:jiaocha）=0; %祛除影响
maxn=max（r）; %找最大值
temp=find（r==maxn）;%返回第一个最大值
nmax=temp（1）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        582  2008-10-13 22:16  find_maxn.m

     文件       2177  2008-10-14 15:53  jiyinzhouqi.m

     文件       3800  2008-10-13 22:07  test2（228hz）.wav

     文件       4392  2008-10-13 22:08  test3（350hz）.wav

     文件       4236  2008-10-13 22:09  test4（308hz）.wav

----------- ---------  ---------- -----  ----

                15187                    5