pisar1.m

来自「谱估计 (1)修正协方差法 (2)多重信号分类（MUSIC）算法 (3)E」· M 代码 · 共 21 行

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算法');