資源簡介
用acf算出的F0 存至矩陣中 可以直接plot(x)出來
代碼片段和文件信息
%?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
評論
共有 條評論