怎么画出一下互功率谱函数啊

715739802

load sumsin;
x=sumsin;
load  noissin;
y=noissin;
Rxy=xcorr(x,y,'coeff');
plot(Rxy);
R=abs(fft(Rxy));
plot(R);
结果不对啊。
对于自功率谱函数我尝试过用pyulear直接画和先求出自相关函数然后将自相关函数傅里叶变换，两种方法结果完全不一样
怎么回事啊
例
load sumsin;
x=sumsin;
Rx=xcorr(x,'coeff');
plot(Rx);
R=abs(fft(Rx));
plot(R);
这是用pyulear直接画的
load sumsin;
x=sumsin;
R=pyulear(x,4,[]);
plot(R);