自相关矩阵

kevin19821

对于pisarenko谐波分析法中的自相关矩阵Ｒ＝Ｅ（Ｘ＊Ｘ＇）应该怎么用matlab实现？我自己用过
clc
clear all
close all
fs=1000;
dt=1/fs;
n=3000;
%  N=3000;
t=0:dt:(n-1)*dt;
randn('state',0);
tt=0.035*randn(1,n);
s1=(1+0.05*cos(2*pi*10*t+pi/6)+0.1*cos(2*pi*5*t+pi/8)).*(cos(2*pi*50*t+pi/3)+0.05*cos(2*pi*150*t+pi/4)+0.03*cos(2*pi*250*t+pi/5));
s=s1+tt;
  s=s';
N=1000;
for k1=1:N
   for k2=1:N
       R(k1,k2)=mean(s(k1,:).*conj(s(k2,:)));%对y矩阵做相关运算
   end
end
eig(R)
来实现，但是分解出来的特征值总是不对