% 求输入xn和yn的互相关函数和双谱
clc
N=100;
n=1:N;
xn=exp(j*pi*n);            
yn=xn;            
Rxy=xcorr(xn,yn,'biased'); % 求输入xn和yn的互相关
subplot(211)
plot(abs(Rxy))
title('xn和yn的互相关')
Sxy=fft(Rxy);              % 对xn和yn的互相关求傅氏变换，即为互谱
subplot(212)
n1=1:size(Rxy,2);
f=n1/(2*N-1);
plot(f,abs(Sxy))
title('xn和yn的互谱')

cpsd()这个函数什么意思呢，
是不是也是互功率谱？
这两个有什么不同。

[quote]
% 求输入xn和yn的互相关函数和双谱
clc
N=100;
n=1:N;
xn=exp(j*pi*n);            
yn=xn;            
Rxy=xcorr(xn,yn,'biased'); % 求输入xn和yn的互相关
subplot(211)
plot(abs(Rxy))
title('xn和yn的互相关')
Sxy=fft(Rxy);              % 对xn和yn的互相关求傅氏变换，即为互谱
subplot(212)
n1=1:size(Rxy,2);
f=n1/(2*N-1);
plot(f,abs(Sxy))
title('xn和yn的互谱')[/quote]
谢谢您
请问：用互相关的话，长度是2N（假设原信号都是N），fft过后也是2N。
假如我直接用cpsd（）函数，并且最后参数设为N（就是nfft），结果只有一半N/2个点，请问，这个有什么关系，联系呢。谢谢。