自相关法基音提取的matlab程序,自相关函数法基音周期提取(matlab版)

function nmax=find_maxn(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)

end

end

r(1:jiaocha)=0; %祛除影响

maxn=max(r); %找最大值

temp=find(r==maxn);%返回第一个最大值

nmax=temp(1);

function jiyinzhouqi(filename,shift)

%短时自相关分析

%filename语音文件*.wav

%zhouqi基音周期

shift=10;

[signal,fs]=wavread('f:/mywork/1.wav');

shift=round(fs*shift); %帧移

n1=fix(fs*0.97)+1; %分析起点970ms,帧长30ms

n2=fix(fs*1)+1;

ii=1;

for ii=1:(length(signal)-n1)/shift %分析次数

if n2

data=signal(n1:n2);

N=n2-n1+1;