matlab 音调跟踪Pitch Tracking

% Pitch Tracking  
% 
[xfs]=wavread（‘c1_sentence1‘）;
N=600;
OL_rate=0.5;
N_shift=N*（1-OL_rate）;
framenumber=fix（length（x）/N_shift-N/N_shift）+1;
for i=1:framenumber
     k=N_shift*（i-1）+1;
     acf=xcorr（x（k:k+N-1））;
     [rmax1imax1]=max（acf）;  %  1st peak is at N
     [rmin1imin1]=min（acf）;
     acf（imin1:2*imax1-imin1）=0;
     [rmax2imax2]=max（acf）;
     pitch=abs（imax1-imax2）/fs;
     pm（i）=pitch;
     
end


clear all