为什么我做scot广义互相关不对？求助

tony6257

我做广义互相关求出的时延是正确的，但我在原广义互相关上加权后，时延估计在0点处，找不到原因，请各位指教，谢谢。

1. clear all;
   
2. close all;
   
3. clc;
   
4. 
   
5. N=1024;  %长度
   
6. Fs=500;  %采样频率
   
7. n=0:N-1;
   
8. t=n/Fs;   %时间序列
   
9. a1=5;     %信号幅度
   
10. a2=5;
    
11. d=8;     %延迟点数
    
12. 
    
13. x1=a1*cos(2*pi*10*n/Fs);     %信号1
    
14. x1=x1+randn(size(x1));      %加噪声
    
15. snr=-5;
    
16. x1=awgn(x1,snr);
    
17. x2=a2*cos(2*pi*10*(n+d)/Fs); %信号2
    
18. x2=x2+randn(size(x2));
    
19. snr=-5;
    
20. x2=awgn(x2,snr);
    
21. 
    
22. figure;
    
23. subplot(211);
    
24. plot(t,x1,'r');
    
25. axis([-0.2 1.5 -6 6]);
    
26. hold on;
    
27. plot(t,x2,':');
    
28. axis([-0.2 1.5 -6 6]);
    
29. legend('x1信号', 'x2信号');
    
30. xlabel('时间/s');ylabel('x1(t) x2(t)');
    
31. title('原始信号');grid on;
    
32. hold off
    
33. 
    
34. %互相关函数
    
35. X1=fft(x1,2*N-1);
    
36. XX1=abs(X1).^2;
    
37. X2=fft(x2,2*N-1);
    
38. XX2=abs(X2).^2;
    
39. Sxy=X1.*conj(X2;
    
40. 
    
41. Cxy_modulus=abs(Sxy);
    
42. ratio=Sxy./Cxy_modulus;
    
43. Cxy=fftshift(real(ifft(ratio)));
    
44. 
    
45. subplot(212);
    
46. t1=(-(N-1):(N-1))/Fs;                        %注意
    
47. plot(t1,Cxy,'b');
    
48. [max,location]=max(Cxy);%求出最大值max,及最大值所在的位置（第几行）location;
    

复制代码