limesprit.m

功率谱估计程序

%Esprit算法。选用电力系统谐波信号作为输入信号。参考书目：现代信号处理(第二版），张贤达，清华大学出版社。141页
clear all;close all;
f1=50;ts=1/1000;n=1:256;N=256;fs=1000;
x=sin(2*pi*f1*n*ts)+0.3*sin(2*pi*350*n*ts)+0.3*sin(2*pi*360*n*ts);D=0;p=6;
y=x';
noise=randn(1,256);
noise=noise';
dB=10;
y=sigmerge(y,noise,dB);
x=y';
for m=0:N-1
    s=0;
    for n=0:N-1-m
        s=s+x(n+1)*x(n+m+1)/N;
    end
    r(m+1)=s;
end
r=toeplitz([r,0]);
r1=r(1:N,1:N);
r2=r(1:N,2:N+1);

e=eig(r1);
e=min(e);
C1=r1(1:p,1:p)-e*eye(p);
a=eye(p+1);
a=a(1:p,2:end);
C2=r2(1:p,1:p)-e*a;
l=eig(C1,C2);
f=angle(l)/2/pi;
stem(f,ones(1,length(f)))