请教利用广义互相关算法求时延的问题

[sixsun]

各位前辈好，毕设要用到广义相关算法，根据自己的理解，写了一段代码，但是运行结果不对，由于个人能力有限，不知道问题出在哪，还烦请大家帮忙分析一下！谢谢！
function R = my_GCC(x,y)
NFFT = length(x)+length(y)-1;
xx = [x;zeros(NFFT-length(x),1)];
yy = [y;zeros(NFFT-length(y),1)];
xx = fft(xx,NFFT);
yy = fft(yy,NFFT);
Pxx = abs(xx).^2/(2*length(x));
Pyy = abs(yy).^2/(2*length(y));
xy = xcorr(x,y);

Pxy = fft(xy,NFFT);
r = Pxy./sqrt(Pxx.*Pyy);
R = ifft(r,NFFT);
以上是基于SCOT方法写的代码！
附件里是一篇关于广义相关的英文文献，具体用到的公式是17，18，19，请大家帮忙看看吧~谢谢了~