所求的互相关函数图像对吗？

yangcui

以下是我作的，遇到困难了，请大家帮忙解决一下。求两个信号序列的互相关函数，4.txt和6.txt中第二列是实测信号序列，第一列为采样时间。我按照一本书所说的，先求出两者的互功率谱，然后进行傅立叶逆变换，取傅立叶变换的实部为互相关函数，结果如下图所示。我感觉和预想的结果不一样。下面是用的程序，是不是建立负频率段的互功率谱错了，还是其他的地方出错了，请高手帮忙。万分感激！
clc
clear
data1=load('4.txt');
data2=load('6.txt');
x=data1(:,2)';
y=data2(:,2)';
np=length(x);
fs=1000;
nfft=1024;
t=0:1/fs:(np-1)/fs;
Pxy=CSD(x,y,nfft,fs,[]);  %计算互功率谱
Pxy(nfft/2+1)=real(Pxy(nfft/2));  %建立负频率段的互功率谱
Pxy(nfft/2+2:nfft)=conj(Pxy(nfft/2:-1:2));
g=ifft(Pxy);    %傅立叶逆变换
g=real(g);    %取傅立叶变换的实部为互相关函数
r=g(1:np);%绘制互相关函数时程曲线图
plot(t,r);
grid;