信号自相关矩阵特征值

kevin19821

小弟在编写张贤达现代信号处理一书中得ESPRIT算法时,对信号形成得自相关矩阵进行特征值分解,为什么前19个都为0,间下面得运行结果(见张贤达书139页)
clc
clear all
clf
fs=1000;
n=2048;
dt=1/fs;
t=0:dt:(n-1)*dt;
s=cos(2*pi*50*t)+cos(2*pi*150*t)+0.1*randn(1,n);
plot(s)
N=20;
y=s(1:N)';
% %%%%%======求自相关矩阵
R=zeros(N,N);
RR=zeros(N,N);
for k=1:N
  for k1=1:N
    for k2=1:N
RR(k1,k2)=(y(k1,:).*conj(y(k2,:)));%相关运算
    end
  end
R=R+RR;
end
eig(R)
运行结果为ans =
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
   -0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
  386.3844

generalstaff

你经常来这个论坛吗?以后互相照应欧.

mylovekite

我认为只对一个信号的自相关矩阵进行特征分解得到的结果就是这样的，也就是说对一个向量的相关矩阵R=X*X‘/N得到的肯定是你所说的形式。