论坛 › 信号处理方法 › 信号自相关矩阵特征值

kevin19821 发表于 2008-4-11 16:26

信号自相关矩阵特征值

小弟在编写张贤达现代信号处理一书中得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 发表于 2008-4-21 19:52

你是自动化系哪个所的?我也选了这课

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

mylovekite 发表于 2008-5-1 08:59

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

查看完整版本: 信号自相关矩阵特征值