%例程3-2  計算隨機信號的相關函數

x=[1 3 2 6 2 1 -2 0 1 5 3 2 -3 0 1 2 0 3 1];
y=[0 0 0 1 3 2 6 2 1 -2 0 1 5 3 2 -3 0 1 2 0 3 1];
k=length(y);
e=randn(1,k);
y=y+e;
xk=fft(x,2*k);
yk=fft(y,2*k);
rm=real(ifft(conj(xk).*yk));
rm=[rm(k+2:2*k) rm(1:k)];
m=(-k+1):(k-1);
stem(m,rm) 
xlabel('m');
ylabel('幅度');