pisar1.m

谱估计 (1)修正协方差法 (2)多重信号分类（MUSIC）算法 (3)ESPRIT算法 (4)皮萨论科（Pisarenko）谐波分解法 并对各算法进行分析。 解压后即可在MATLAB上运

xt=xn'; 
sin_num=3;
N=length(xt);
rxx=xcorr(xt,'biased');
rxx=rxx(N:(2*sin_num)+N);
%Frequencies estimation
Rxx=toeplitz(rxx);
ev=eig(Rxx);
[S i]=min(ev);
[V D]=eig(Rxx);
a=V(:,i)
rts=roots(a)
w_est=[];
for i=1:6
   w_est(i)= angle(rts(i));
end
F=(w_est/(2*pi))';
figure
  stem(F,ones(1,length(F)));
  title('Pisarenko算法');