xcorral.m

matlab中的相关算法已通过验证

[x1,fs1,bits1]=wavread('C:\Documents and Settings\hp\My Documents\MATLAB\ak2.wav');
[x2,fs2,bits2]=wavread('C:\Documents and Settings\hp\My Documents\MATLAB\ak1.wav');
cc=xcorr(x1,x2);
subplot(2,1,1);plot(cc);title('result of xcorr')
cmax1=max(abs(cc));
maxpoint1=find(abs(cc)==cmax1);
tdelay1=(maxpoint1-length(cc)/2)/fs1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
s1=fft(x1,2048);
s2=fft(x2,2048);
s=s1.*conj(s2);
r1=real(fftshift(ifft(s)));
subplot(2,1,2);plot(r1);title('result of ifft')
cmax2=max(abs(r1));
maxpoint2=find(abs(r1)==cmax2);
tdelay2=(maxpoint2-length(r1)/2)/fs2