📄 xcorral.m
字号:
[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
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -