gccphat.m

本次的算法是相位周期展开聚类

function y=gccphat(t1,t2);
%[sx,fs,nibits]=wavread('d:\aaa\ttt.wav');
%x=sx(:,1);
%x=x';
%fs=16000;
%x=randn(1,30);
%x1=delay(x,2,0);
%x2=delay(x,10,0);
%snr=10;
%t1=awgn(x1,snr);%加入SNR的正态分布白噪声随机序列
%t2=awgn(x2,snr);
fx1=fft(t1);
fx1=fx1(1:256);
fx2=fft(t2);
fx2=fx2(1:256);
Pxy=fx1.*conj(fx2);
Rxy=Pxy./abs(Pxy);
R=ifft(Rxy);
d1=real(R);
d2=max(abs(R));
a=find(abs(R)==d2); %找出最大点
Rl2=fix(length(R)/2);
figure(3);
time=(-Rl2:Rl2)/fs;
plot(time,R); 
grid;